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

ATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTT  >  W3110S.gb/3824459‑3824588
                                                                 |                                                                 
atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                              >  1:796824/1‑71 (MQ=255)
    aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                           <  1:1634410/70‑1 (MQ=255)
     aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGACGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                           <  1:1250714/69‑1 (MQ=255)
     aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                        <  1:1466172/71‑1 (MQ=255)
     aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                        <  1:1703102/71‑1 (MQ=255)
          aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                   <  1:703582/71‑3 (MQ=255)
                cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaa                                              <  1:342671/71‑1 (MQ=255)
                        ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                                <  1:67331/61‑1 (MQ=255)
                                gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                        >  1:840105/1‑44 (MQ=38)
                                gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                        >  1:1440321/1‑44 (MQ=38)
                                gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                        >  1:597342/1‑44 (MQ=38)
                                 ttCAGCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                          <  1:1630188/58‑1 (MQ=255)
                                 ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                                          <  1:55101/58‑1 (MQ=255)
                                    aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                           <  1:1725879/54‑1 (MQ=255)
                                    aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                           <  1:1790946/54‑1 (MQ=255)
                                    aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                           <  1:1639763/54‑1 (MQ=255)
                                                   cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                      <  1:219798/44‑1 (MQ=25)
                                                   cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                      <  1:1269989/44‑1 (MQ=25)
                                                   cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa             <  1:1676478/69‑1 (MQ=255)
                                                        cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAt      >  1:1575225/1‑71 (MQ=255)
                                                            ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt  <  1:310670/71‑1 (MQ=255)
                                                                 |                                                                 
ATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTT  >  W3110S.gb/3824459‑3824588

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: