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.G90.9% 27.4 / ‑2.7 11pseudogene (666/687 nt)rphECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (3/7);  total (4/7)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGC  >  W3110S.gb/3824473‑3824572
                                                             |                                        
gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                <  1:897506/71‑1 (MQ=255)
  cGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcgaa                              >  1:1556389/1‑71 (MQ=255)
      aGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCa                          <  1:704075/71‑1 (MQ=255)
                     cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAg                                         >  1:1061868/1‑41 (MQ=38)
                       cccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                <  1:103367/48‑1 (MQ=255)
                       cccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAg         >  1:1211786/1‑71 (MQ=255)
                        ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAgg        <  1:1530249/71‑1 (MQ=255)
                           tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                <  1:1425942/44‑1 (MQ=25)
                           tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGGgaatcg                                <  1:908625/44‑1 (MQ=25)
                              agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGgcgc  >  1:1337977/1‑71 (MQ=255)
                                   tACTCACCTTGTTGGCGCTGGCCCGAGGGGGAATCGAATCCAt                         <  1:1414918/43‑1 (MQ=25)
                                            tgttgGCGCTGGCCCGAGGGGAAATCGAATCCATTGTa                     >  1:19124/1‑38 (MQ=37)
                                                             |                                        
GACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGAGGGGAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGC  >  W3110S.gb/3824473‑3824572

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