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

GACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCA  >  W3110S.gb/3824447‑3824577
                                                                   |                                                                 
gACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCtt                                                                <  1:1864773/71‑1 (MQ=255)
  cGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                              <  1:2286272/71‑1 (MQ=255)
  cGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                              <  1:2283234/71‑1 (MQ=255)
     gcgcATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTgg                                                           <  1:1425829/71‑1 (MQ=255)
     gcgcATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTgg                                                           <  1:2382108/71‑1 (MQ=255)
        catcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgc                                                        <  1:736327/71‑3 (MQ=255)
          tcatTGAAGCGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                      <  1:1652209/71‑1 (MQ=255)
            atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTgg                                                           >  1:1663011/1‑64 (MQ=255)
              tGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcg                                                         <  1:1950651/64‑2 (MQ=255)
                     cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                                         <  1:2112123/71‑2 (MQ=255)
                              gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGct                                                         <  1:1945420/48‑1 (MQ=255)
                              gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                 >  1:1300405/1‑71 (MQ=255)
                               cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                              >  1:2429069/1‑57 (MQ=255)
                                   aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGTGAATCGAATCCATTg                            >  1:2041341/1‑71 (MQ=255)
                                    ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                           <  1:347950/71‑1 (MQ=255)
                                    ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                           <  1:2533899/71‑1 (MQ=255)
                                               cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                 <  1:192656/70‑1 (MQ=255)
                                                    aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg           >  1:1658272/1‑71 (MQ=255)
                                                     tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                              <  1:748659/35‑1 (MQ=37)
                                                     tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                              <  1:2168904/35‑1 (MQ=37)
                                                     tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                              <  1:2054062/35‑1 (MQ=37)
                                                                tCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa  <  1:1784528/68‑1 (MQ=255)
                                                                   |                                                                 
GACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCA  >  W3110S.gb/3824447‑3824577

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: