Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,533:1 +G 100% pseudogene (666/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,5331.G100.0% 59.7 / NA 18pseudogene (666/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (9/9);  total (9/9)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

AGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGG  >  W3110S.gb/3824464‑3824592
                                                                      |                                                           
aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                            <  1:1392140/71‑1 (MQ=255)
 gTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAgg                                                           >  1:1559632/1‑71 (MQ=255)
    cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                                                        <  1:1161823/71‑2 (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)
                       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)
                                                  ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGa           <  1:108036/70‑1 (MQ=255)
                                                     tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt         >  1:392764/1‑70 (MQ=255)
                                                             ctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAgg  >  1:1091769/1‑69 (MQ=255)
                                                                      |                                                           
AGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGG  >  W3110S.gb/3824464‑3824592

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: