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

CAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  W3110S.gb/3824468‑3824597
                                                                  |                                                                
cAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                                                 >  1:461268/1‑67 (MQ=255)
 aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                            <  1:2357366/71‑3 (MQ=255)
 aGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggga                                                            <  1:723407/71‑3 (MQ=255)
         gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                     <  1:821239/71‑1 (MQ=255)
         gCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaatc                                                     <  1:1714270/71‑1 (MQ=255)
            gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                                                             >  1:2595932/1‑58 (MQ=255)
            gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                                                             >  1:3509509/1‑58 (MQ=255)
            gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                                                             >  1:2550553/1‑58 (MQ=255)
            gAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAggggg                                                             >  1:1675742/1‑58 (MQ=255)
               ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:1984047/1‑70 (MQ=255)
               ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:3641532/1‑70 (MQ=255)
               ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:667357/1‑70 (MQ=255)
               ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGt                                               <  1:1503848/71‑1 (MQ=255)
               ggCGAGCCGTTCACCCATGAAGAGATACTCACCTTGTTGGCGCTGGACCGAgcg                                                               >  1:3489216/1‑52 (MQ=255)
                     ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGAc                                         >  1:755158/1‑71 (MQ=255)
                     ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGAc                                         >  1:2986264/1‑71 (MQ=255)
                      cGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACg                                        <  1:3033738/71‑1 (MQ=255)
                          cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     >  1:727954/1‑70 (MQ=255)
                               aTGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg                               >  1:2727104/1‑71 (MQ=255)
                                tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                                                  <  1:1281463/51‑1 (MQ=255)
                                tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAg                                     <  1:2281106/64‑1 (MQ=255)
                                tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcg                               >  1:3707253/1‑70 (MQ=255)
                                tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAGTGTAGc                                            <  1:523119/57‑1 (MQ=255)
                                 gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggcggc                             >  1:1166520/1‑71 (MQ=255)
                                  aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg                            >  1:1172906/1‑71 (MQ=255)
                                  aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGgcggcg                            >  1:3269975/1‑71 (MQ=255)
                                    gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:4076294/1‑49 (MQ=255)
                                    gagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTg                                                >  1:784009/1‑49 (MQ=255)
                                         aCTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGa                                          >  1:958847/1‑50 (MQ=255)
                                             accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                            >  1:819289/3‑44 (MQ=37)
                                             accTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt                  >  1:617599/3‑70 (MQ=255)
                                             accTTGCTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt                  >  1:1618231/3‑70 (MQ=255)
                                                ttgttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCa                                      <  1:1941479/47‑1 (MQ=255)
                                                   ttgGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGAttttt           <  1:549979/71‑1 (MQ=255)
                                                        gctGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc      <  1:2299076/70‑1 (MQ=255)
                                                            gCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACt  <  1:2322687/71‑1 (MQ=255)
                                                                  |                                                                
CAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGCGACT  >  W3110S.gb/3824468‑3824597

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: