Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,410 T→A 100% pseudogene (543/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,4100TA100.0% 116.2 / NA 40pseudogene (543/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (21/19);  total (21/19)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGC  >  W3110S.gb/3824345‑3824468
                                                                 |                                                          
gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                                                        <  1:2147498/70‑1 (MQ=255)
       tttCTGTCTGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCag                                                     <  1:2460481/66‑1 (MQ=255)
         tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag                                                   >  1:1864132/1‑66 (MQ=255)
                 taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                      >  1:2190348/2‑71 (MQ=255)
                   ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                      >  1:773132/1‑69 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                  >  1:1045375/4‑71 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                  >  1:39173/4‑71 (MQ=255)
                     gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                  >  1:2339685/4‑71 (MQ=255)
                            gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                               >  1:2485466/3‑67 (MQ=255)
                            gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                               >  1:1741404/3‑67 (MQ=255)
                            gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                               >  1:2608129/3‑67 (MQ=255)
                            gtgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgat                               >  1:1124634/3‑67 (MQ=255)
                             tgAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGAcc                           >  1:1789278/2‑70 (MQ=255)
                                      ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACggg                  >  1:263998/1‑70 (MQ=255)
                                      ttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGgc                 >  1:309334/1‑71 (MQ=255)
                                             cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat          >  1:94061/1‑71 (MQ=255)
                                             cTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat          >  1:2176741/1‑71 (MQ=255)
                                                gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                      >  1:2053877/1‑56 (MQ=255)
                                                gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                      >  1:1837798/1‑56 (MQ=255)
                                                gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                      >  1:2589547/1‑56 (MQ=255)
                                                gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat          <  1:1190493/68‑1 (MQ=255)
                                                gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa       <  1:883300/71‑1 (MQ=255)
                                                gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa       <  1:429570/71‑1 (MQ=255)
                                                gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa       <  1:1260863/71‑1 (MQ=255)
                                                gAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa       <  1:1272901/71‑1 (MQ=255)
                                                gAATACCTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa       <  1:40925/71‑1 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa       <  1:2019078/68‑1 (MQ=255)
                                                   tACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa       <  1:1297076/68‑1 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGa       >  1:221347/1‑66 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTg   >  1:758493/1‑70 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGc  <  1:1670348/71‑1 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGc  <  1:704145/71‑1 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGc  <  1:1839464/71‑1 (MQ=255)
                                                     cGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGc  <  1:1168640/71‑1 (MQ=255)
                                                       ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                  <  1:1335494/37‑1 (MQ=255)
                                                       ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                  <  1:1474261/37‑1 (MQ=255)
                                                       ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                  <  1:23010/37‑1 (MQ=255)
                                                       ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                  <  1:596747/37‑1 (MQ=255)
                                                       ttGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                  <  1:641169/37‑1 (MQ=255)
                                                       ttGAAGACACAGCCGCAGAGACCGACATGAACGTAGttat                               >  1:2438344/1‑40 (MQ=25)
                                                                 |                                                          
GCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGC  >  W3110S.gb/3824345‑3824468

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: