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

GGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCG  >  W3110S.gb/3824261‑3824387
                                                                 |                                                             
ggTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCg                                                            <  1:2041035/69‑1 (MQ=255)
          tGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                    <  1:2332350/67‑1 (MQ=255)
              gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                              <  1:3050358/69‑1 (MQ=255)
              gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                              <  1:2917845/69‑1 (MQ=255)
              gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt                                              <  1:2837151/69‑1 (MQ=255)
                        gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                    >  1:3394623/1‑69 (MQ=255)
                        gcgcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt                                    >  1:1713390/1‑69 (MQ=255)
                         cgcTACAGAAGCTGATGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg                                                    >  1:3472158/1‑52 (MQ=39)
                           cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt                                 <  1:3180067/69‑1 (MQ=255)
                                gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATTAAAgg                                                    <  1:2267067/45‑1 (MQ=38)
                                gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                        >  1:1026692/1‑41 (MQ=255)
                                gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                                                        >  1:1465748/1‑41 (MQ=255)
                                   gCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAt                         <  1:2070653/69‑1 (MQ=255)
                                     tggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTAtt                        >  1:2257668/1‑68 (MQ=255)
                                         ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                     <  1:2385950/67‑2 (MQ=255)
                                         ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta                     <  1:1430028/67‑2 (MQ=255)
                                         ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                   >  1:2711419/1‑69 (MQ=255)
                                                          aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc   >  1:2665610/1‑68 (MQ=255)
                                                          aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg  >  1:3133106/1‑69 (MQ=255)
                                                          aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCg  >  1:1076050/1‑69 (MQ=255)
                                                           aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggt                           >  1:331044/1‑42 (MQ=255)
                                                           aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggt                           >  1:2935870/1‑42 (MQ=255)
                                                                 |                                                             
GGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCG  >  W3110S.gb/3824261‑3824387

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: