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% 66.4 / NA 21pseudogene (647/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (4/17);  total (4/17)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCA  >  W3110S.gb/3824446‑3824577
                                                                    |                                                                
aGACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCt                                                                <  1:454578/71‑1 (MQ=255)
aGACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCt                                                                <  1:1775263/71‑1 (MQ=255)
             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:597342/1‑44 (MQ=38)
                                             gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                             >  1:1440321/1‑44 (MQ=38)
                                              ttCAGCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                               <  1:1630188/58‑1 (MQ=255)
                                              ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                               <  1:55101/58‑1 (MQ=255)
                                                 aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                <  1:1790946/54‑1 (MQ=255)
                                                 aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATcc                                <  1:1725879/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)
                                                                    |                                                                
AGACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCA  >  W3110S.gb/3824446‑3824577

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: