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

GTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCC  >  W3110S.gb/3824219‑3824347
                                                             |                                                                   
gTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCt                                                            <  1:1471723/71‑1 (MQ=255)
  gCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTa                                                           >  1:1521117/1‑70 (MQ=255)
     tCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAg                                                        >  1:698289/1‑70 (MQ=255)
         gCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCt                                                   <  1:612936/71‑1 (MQ=255)
                gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGa                                             >  1:1324105/1‑70 (MQ=255)
                gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGa                                             >  1:1578339/1‑70 (MQ=255)
                gtggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGa                                             >  1:272438/1‑70 (MQ=255)
                 tggCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaa                                           <  1:861168/71‑1 (MQ=255)
                     aCGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACgg                                       >  1:1014017/1‑71 (MQ=255)
                      cgcgTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGc                                      <  1:302917/71‑1 (MQ=255)
                          tACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                 <  1:746261/56‑1 (MQ=255)
                            ccGCGTCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTg                                <  1:313472/71‑1 (MQ=255)
                                 tCGATTACGGGTTCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCtgg                                                 <  1:1749089/49‑1 (MQ=38)
                                 tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTg                                <  1:707054/66‑1 (MQ=255)
                                 tCGATTACGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTg                                <  1:104309/66‑1 (MQ=255)
                                        cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGc                                  >  1:1529558/1‑57 (MQ=255)
                                              ccTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGa                >  1:1669843/1‑69 (MQ=255)
                                                  cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaaa                                          >  1:1073880/1‑39 (MQ=255)
                                                  cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGaaa                                           >  1:558933/1‑38 (MQ=255)
                                                  cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAggg           <  1:304977/70‑1 (MQ=255)
                                                      gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg       >  1:1257381/1‑70 (MQ=255)
                                                        gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCaac                       >  1:155966/1‑51 (MQ=255)
                                                        gcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAg    <  1:211996/71‑1 (MQ=255)
                                                          tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGcc  >  1:942051/1‑71 (MQ=255)
                                                             |                                                                   
GTGCTTCAGGCTGATGGTGGCACGCGTACCGCGTCGATTACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCC  >  W3110S.gb/3824219‑3824347

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: