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% 113.5 / NA 34pseudogene (666/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (20/14);  total (20/14)
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.

TGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  W3110S.gb/3824466‑3824589
                                                                    |                                                        
tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGa                                                           <  1:666498/68‑1 (MQ=255)
tGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggg                                                       >  1:169494/1‑71 (MQ=255)
        aCGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcga                                                <  1:609740/71‑1 (MQ=255)
          ggCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaat                                              >  1:1935389/1‑71 (MQ=255)
            cAGAAGGCGAGCCGTTCACCAATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                             >  1:1808500/1‑70 (MQ=255)
               aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGATGGCCCGAGGGGGAATCGAATCCAt                                          >  1:1088362/1‑70 (MQ=255)
                    gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                    <  1:3731694/71‑1 (MQ=255)
                    gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                    <  1:235365/71‑1 (MQ=255)
                            cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                        >  1:1138796/1‑59 (MQ=255)
                            cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                             >  1:1660170/1‑70 (MQ=255)
                             aCCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                           <  1:777572/71‑1 (MQ=255)
                                cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                        >  1:2850328/1‑71 (MQ=255)
                                cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                        >  1:2841980/1‑71 (MQ=255)
                                cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                        <  1:2608057/71‑1 (MQ=255)
                                  tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                             <  1:86225/48‑1 (MQ=255)
                                   gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg                       >  1:539004/1‑69 (MQ=255)
                                    aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                    >  1:3362324/1‑55 (MQ=255)
                                    aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc                     <  1:1986870/70‑1 (MQ=255)
                                    aaGAGCTACTCACCTTGTTGGCGCTGCCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc                     >  1:659279/1‑70 (MQ=255)
                                     agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                             <  1:3714994/61‑1 (MQ=255)
                                       tgCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                             <  1:1028504/58‑1 (MQ=255)
                                       agCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                              >  1:1297176/1‑58 (MQ=255)
                                       agCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                              >  1:1291733/1‑58 (MQ=255)
                                       agCTACTAACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGaa                           >  1:3995946/1‑61 (MQ=255)
                                         cTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGa                            >  1:4638611/1‑58 (MQ=255)
                                               accTTGTTGGCTCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTg         <  1:3058868/69‑1 (MQ=255)
                                               accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                  >  1:409568/3‑46 (MQ=255)
                                               accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt          >  1:3672084/3‑70 (MQ=255)
                                               accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTg         <  1:1037568/69‑1 (MQ=255)
                                                ccTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                       <  1:3866153/39‑1 (MQ=37)
                                                 cTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                       <  1:611246/39‑1 (MQ=37)
                                                   tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttt     >  1:3410898/1‑71 (MQ=255)
                                                   tgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAtt      >  1:1621828/1‑70 (MQ=255)
                                                    gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                  >  1:4026852/1‑41 (MQ=37)
                                                    gttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa            >  1:1304892/1‑63 (MQ=255)
                                                       ggCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa  <  1:1136973/70‑1 (MQ=255)
                                                                    |                                                        
TGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  W3110S.gb/3824466‑3824589

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: