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

AGACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGC  >  W3110S.gb/3824446‑3824567
                                                                    |                                                      
aGACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCt                                                      >  1:2200497/1‑71 (MQ=255)
  aCGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttg                                                    >  1:2093754/1‑71 (MQ=255)
           tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                           >  1:185524/1‑71 (MQ=255)
               tGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgc                                             >  1:2266295/1‑63 (MQ=255)
                     gCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggg                                >  1:1048568/1‑71 (MQ=255)
                      cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                               <  1:2834136/71‑2 (MQ=255)
                       aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                              <  1:465506/71‑3 (MQ=255)
                        ggggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgggggaa                             <  1:1731962/71‑4 (MQ=255)
                           gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGcgc                                             >  1:319411/1‑51 (MQ=255)
                              ggCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                   <  1:1930678/44‑1 (MQ=255)
                               gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                       <  1:2690513/71‑1 (MQ=255)
                                    aGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                       >  1:812206/1‑66 (MQ=255)
                                     ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                  >  1:3065803/1‑70 (MQ=255)
                                        gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                 <  1:1418018/51‑1 (MQ=255)
                                         aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg             <  1:2648002/71‑1 (MQ=255)
                                             gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGATGGGGAATCGAATCCATTGTAGCGACGc         <  1:3021076/71‑1 (MQ=255)
                                              ttCACCCATGAAGAGCTACTGACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg          >  1:3299930/1‑69 (MQ=255)
                                              ttCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                       <  1:1472110/40‑1 (MQ=38)
                                               tCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                               <  1:2030430/46‑2 (MQ=255)
                                                cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg       <  1:3361837/70‑1 (MQ=255)
                                                cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg       >  1:43609/1‑70 (MQ=255)
                                                    cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc  <  1:1495939/71‑1 (MQ=255)
                                                    cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                   >  1:580136/1‑37 (MQ=38)
                                                    cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGgaa                              <  1:1805752/42‑1 (MQ=38)
                                                                cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                 >  1:3276853/1‑44 (MQ=25)
                                                                 tCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc              <  1:981878/46‑1 (MQ=255)
                                                                    |                                                      
AGACGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGC  >  W3110S.gb/3824446‑3824567

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