Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,315,139 T→A 100% D244V (GAT→GTT)  infB ← fused protein chain initiation factor 2, IF2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,315,1390TA80.2% 48.1 / 6.8 25D244V (GAT→GTT) infBfused protein chain initiation factor 2, IF2
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base A (9/11);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTGCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAGTCGG  >  W3110S.gb/3315073‑3315209
                                                                  |                                                                      
ttgCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCTCGGCCGCCTTCGACTTCACGAACGc                                                                     <  1:532834/70‑1 (MQ=255)
ttgCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCt                                                                    <  1:1539533/71‑1 (MQ=255)
    cTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCg                                                                <  1:1922642/71‑1 (MQ=255)
        cttcGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGtct                                                             <  1:674747/70‑1 (MQ=255)
           cGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCtt                                                                   >  1:1199733/1‑61 (MQ=255)
           cGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTg                                                         >  1:395626/1‑71 (MQ=255)
           cGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTg                                                         >  1:2240543/1‑71 (MQ=255)
                gcgcTGCTTTGGCGTTACGTCCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTgg                                                    <  1:219485/71‑1 (MQ=255)
                gcgcTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTgg                                                    <  1:163849/71‑1 (MQ=255)
                 cgcTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTgg                                                    <  1:991933/70‑1 (MQ=255)
                  gctgctTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCg                                                  <  1:1108618/71‑1 (MQ=255)
                     gctTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAg                                                >  1:2271052/1‑70 (MQ=255)
                     gctTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAg                                                >  1:819247/1‑70 (MQ=255)
                                aCGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAag                                    >  1:19774/1‑71 (MQ=255)
                                  gACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAagta                                  <  1:2658023/71‑1 (MQ=255)
                                     cacggccacggcCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTg                               <  1:2386078/71‑1 (MQ=255)
                                        ggccacggccGCCTTCGACTTCACGAACGCTTTCGTCTTCTGc                                                        <  1:261321/43‑1 (MQ=255)
                                            acggccGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTg                                         >  1:2279642/1‑54 (MQ=255)
                                                 cgccTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCg                   <  1:2000355/71‑1 (MQ=255)
                                                 cgccTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCg                   <  1:577164/71‑1 (MQ=255)
                                                  gccTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCg                   <  1:2026938/70‑1 (MQ=255)
                                                   ccTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCg                                                  <  1:166837/38‑1 (MQ=255)
                                                             cACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCtt         >  1:2320614/1‑69 (MQ=255)
                                                             cACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCtt         >  1:1215595/1‑69 (MQ=255)
                                                               cGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAGt     <  1:754041/71‑1 (MQ=255)
                                                                  tCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAGTCgg  >  1:2159363/1‑71 (MQ=255)
                                                                  |                                                                      
TTGCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAGTCGG  >  W3110S.gb/3315073‑3315209

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: