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

CGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAA  >  W3110S.gb/3824449‑3824578
                                                                 |                                                                 
cGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                              <  1:346888/71‑1 (MQ=255)
 gggCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgtt                                                             >  1:1268968/1‑71 (MQ=255)
 gggCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgt                                                              <  1:2079678/70‑1 (MQ=255)
  ggCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCttgttg                                                            <  1:1282302/71‑1 (MQ=255)
     gcATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTCTTGGcg                                                         >  1:411836/1‑70 (MQ=255)
          atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTgg                                                     >  1:1881762/1‑70 (MQ=255)
              aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                 <  1:713783/70‑1 (MQ=255)
               aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                              <  1:1392140/71‑1 (MQ=255)
                gTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                             >  1:1559632/1‑71 (MQ=255)
                   cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                                          <  1:1161823/71‑2 (MQ=255)
                       ggACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGcc                                                   <  1:227999/59‑1 (MQ=255)
                         aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                      <  1:565452/70‑1 (MQ=255)
                                aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                <  1:2009166/53‑1 (MQ=255)
                                aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaat                                   >  1:1867536/1‑66 (MQ=255)
                                  ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                            <  1:167934/71‑1 (MQ=255)
                                   gCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                             <  1:637893/52‑1 (MQ=255)
                                    cGAGCCGTTCACCCATGAAGAGCTACTCACCttgtt                                                             >  1:1233423/1‑36 (MQ=255)
                                    cGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                    >  1:1432209/1‑45 (MQ=255)
                                      aGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCg                        <  1:1833455/71‑1 (MQ=255)
                                              aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                       <  1:435612/64‑1 (MQ=255)
                                                ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg              >  1:1733010/1‑71 (MQ=255)
                                                ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                         <  1:76499/43‑3 (MQ=38)
                                                   tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcgg           >  1:2001661/1‑71 (MQ=255)
                                                    gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                <  1:1508491/65‑1 (MQ=255)
                                                     aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc             >  1:102837/1‑67 (MQ=255)
                                                     aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg         >  1:623443/1‑71 (MQ=255)
                                                       gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                             <  1:1270355/49‑1 (MQ=255)
                                                        agCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAtt                              >  1:1177191/1‑47 (MQ=255)
                                                          cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg              >  1:238985/1‑61 (MQ=255)
                                                             cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaa  <  1:585966/70‑1 (MQ=255)
                                                             cTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCa   <  1:834322/69‑1 (MQ=255)
                                                                 |                                                                 
CGGGCGCATCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAA  >  W3110S.gb/3824449‑3824578

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