Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,824,524 T→G 100% pseudogene (657/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,5240TG97.8% 141.0 / ‑5.1 43pseudogene (657/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (26/16);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  W3110S.gb/3824457‑3824589
                                                                   |                                                                  
tcatTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTg                                                                 >  1:389631/1‑71 (MQ=255)
  atTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGc                                                               >  1:2463570/1‑71 (MQ=255)
   ttGAAGTGCAGGGGACGGCAGACGGCGCGCCGTTCACCCATGAAGAGCTACTCTCCTTGTTGGCGCTGGcc                                                              >  1:1612025/1‑71 (MQ=255)
   ttGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGcc                                                              >  1:962690/1‑71 (MQ=255)
     gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                             <  1:4561512/70‑1 (MQ=255)
     gAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGccc                                                             >  1:962571/1‑70 (MQ=255)
      aaGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCg                                                            <  1:1359733/70‑1 (MQ=255)
         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:235365/71‑1 (MQ=255)
                             gAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGc                                    <  1:3731694/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:2608057/71‑1 (MQ=255)
                                         cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                        >  1:2841980/1‑71 (MQ=255)
                                         cATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAggc                        >  1:2850328/1‑71 (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)
                                                                   |                                                                  
TCATTGAAGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTA  >  W3110S.gb/3824457‑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: