Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,362 A→T 100% pseudogene (495/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,824,3620AT100.0% 189.5 / NA 67pseudogene (495/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (34/33);  total (34/33)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGA  >  W3110S.gb/3824294‑3824424
                                                                    |                                                              
aaGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                                                              >  1:1163782/1‑71 (MQ=255)
 aGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTTTCGGTAtt                                                             >  1:647909/1‑71 (MQ=255)
     ggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                          >  1:1047394/1‑69 (MQ=255)
      gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                          <  1:1029863/69‑2 (MQ=255)
         gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg                                                      <  1:3073902/70‑1 (MQ=255)
           aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                          >  1:907218/1‑63 (MQ=255)
           aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                          >  1:1910478/1‑63 (MQ=255)
           aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                   <  1:3928350/71‑3 (MQ=255)
           aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                   <  1:2024184/71‑3 (MQ=255)
            aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                   >  1:274404/1‑68 (MQ=255)
             aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg                                                   >  1:4196372/1‑67 (MQ=255)
               ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                >  1:1330687/1‑70 (MQ=255)
                  aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                                                          <  1:3124732/57‑2 (MQ=255)
                  aaGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaat                                            <  1:4306097/71‑4 (MQ=255)
                       gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                        >  1:3694425/1‑70 (MQ=255)
                       gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                        >  1:232767/1‑70 (MQ=255)
                        aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                      >  1:2329814/1‑71 (MQ=255)
                        aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                      <  1:2884486/71‑1 (MQ=255)
                        aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa                                      <  1:1813908/71‑1 (MQ=255)
                            ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                  <  1:202790/71‑1 (MQ=255)
                             cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                 <  1:1846119/71‑1 (MQ=255)
                             cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                 <  1:3809293/71‑1 (MQ=255)
                             cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                 <  1:2392574/71‑1 (MQ=255)
                             cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTg                                  >  1:1456726/1‑70 (MQ=255)
                               acccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                >  1:3077161/3‑54 (MQ=255)
                                cccGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGaa                               >  1:2823998/2‑70 (MQ=255)
                                  cGATGACAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc                                               >  1:1887/1‑52 (MQ=255)
                                  cGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                        >  1:2842749/1‑59 (MQ=255)
                                   gATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                           <  1:4138096/71‑1 (MQ=255)
                                    aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                             <  1:3406528/52‑3 (MQ=255)
                                    aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                             <  1:2363009/52‑3 (MQ=255)
                                      gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                             >  1:3574559/1‑66 (MQ=255)
                                      gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                             >  1:4186287/1‑66 (MQ=255)
                                      gAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGtt                         >  1:195445/1‑70 (MQ=255)
                                       aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                              <  1:3848014/48‑2 (MQ=255)
                                       aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                        <  1:2759955/70‑1 (MQ=255)
                                       aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                        >  1:4378864/1‑70 (MQ=255)
                                       aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                        <  1:3877002/70‑1 (MQ=255)
                                       aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                        <  1:2387139/70‑1 (MQ=255)
                                       aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                        <  1:1080610/70‑1 (MQ=255)
                                        aaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                      <  1:3541650/71‑1 (MQ=255)
                                             aTGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActc                 <  1:4418413/71‑1 (MQ=255)
                                               ggTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGAct                  >  1:1583316/1‑68 (MQ=255)
                                                gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc              >  1:1650945/1‑71 (MQ=255)
                                                gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc              >  1:3423536/1‑71 (MQ=255)
                                                 tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc              >  1:1458063/1‑70 (MQ=255)
                                                  aGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg            >  1:2958070/1‑71 (MQ=255)
                                                   gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                 <  1:4310905/49‑1 (MQ=25)
                                                   gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg            >  1:3188825/1‑70 (MQ=255)
                                                   gCCGCAGTTTCTGTCGGTATTGTAAACCGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg            <  1:3261125/70‑1 (MQ=255)
                                                     cGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg            >  1:2442978/1‑68 (MQ=255)
                                                      gCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa          <  1:3875356/69‑1 (MQ=255)
                                                       cAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa          >  1:2427079/1‑68 (MQ=255)
                                                        aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg                                         >  1:49513/1‑36 (MQ=16)
                                                        aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                        <  1:487192/37‑1 (MQ=16)
                                                        aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc                                        <  1:1146740/37‑1 (MQ=16)
                                                        aGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAgaga      >  1:1288725/1‑71 (MQ=255)
                                                          tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                 <  1:1899496/42‑1 (MQ=19)
                                                          tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                             <  1:3957182/46‑1 (MQ=25)
                                                          tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATa                             <  1:695985/46‑1 (MQ=25)
                                                          tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa          <  1:660060/65‑1 (MQ=255)
                                                          tttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCa          <  1:1046803/65‑1 (MQ=255)
                                                           ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActcag               >  1:452878/1‑57 (MQ=255)
                                                           ttCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActcag               >  1:474601/1‑57 (MQ=255)
                                                            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGt                          >  1:225996/1‑47 (MQ=25)
                                                            tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCaga        >  1:269608/1‑65 (MQ=255)
                                                             ctgtctgtATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa  <  1:1316784/70‑1 (MQ=255)
                                                             cTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGa  <  1:3379863/70‑1 (MQ=255)
                                                                    |                                                              
AAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGA  >  W3110S.gb/3824294‑3824424

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: