Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,315,139 T→A 29.9% 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,1390TA29.9% 40.5 / 20.4 30D244V (GAT→GTT) infBfused protein chain initiation factor 2, IF2
Reads supporting (aligned to +/- strand):  ref base T (9/12);  new base A (4/5);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAGT  >  W3110S.gb/3315080‑3315206
                                                           |                                                                   
tcttcGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGtc                                                           <  1:1632008/70‑1 (MQ=255)
 cttcGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGtct                                                          <  1:1739585/70‑1 (MQ=255)
          cgcTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGc                                                <  1:567493/71‑1 (MQ=255)
            ctgctTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGc                                                >  1:535065/1‑69 (MQ=255)
                      gTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                    <  1:752577/71‑1 (MQ=255)
                       ttACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTgaga                                   <  1:2601218/71‑1 (MQ=255)
                        tACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCt                                                   >  1:1250529/1‑54 (MQ=255)
                         aCGACCACGGCCAAGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCa                                           >  1:709454/1‑61 (MQ=255)
                          cGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAag                                 <  1:1219682/70‑1 (MQ=255)
                                   ccacggccGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTg                       <  1:474280/71‑1 (MQ=255)
                                   ccacggccGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                    <  1:1857743/58‑1 (MQ=255)
                                     acggccGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGAt                     <  1:2481978/71‑1 (MQ=255)
                                       ggccGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATaa                   >  1:2590115/1‑71 (MQ=255)
                                       ggccGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATaa                   <  1:1790502/71‑1 (MQ=255)
                                            ccTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGc                                            <  1:2532166/41‑1 (MQ=255)
                                              ttCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAagt                                <  1:1308217/51‑1 (MQ=255)
                                              ttCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATa                    <  1:762512/63‑1 (MQ=255)
                                              ttCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGATgag                                    >  1:345891/1‑47 (MQ=255)
                                              ttCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTg             >  1:915651/1‑70 (MQ=255)
                                              ttCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTg             >  1:1658740/1‑70 (MQ=255)
                                               tCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTgg            >  1:1577870/1‑70 (MQ=255)
                                               tCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCAAGTTGAGAAGTAGTGACGTGATAATCGCTgg            >  1:673507/1‑70 (MQ=255)
                                                cGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGa           <  1:1759306/70‑1 (MQ=255)
                                                cGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGaa          >  1:1477233/1‑71 (MQ=255)
                                                 gACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAAt         <  1:2329973/71‑1 (MQ=255)
                                                    ttCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGaa                                  >  1:1950752/1‑43 (MQ=255)
                                                     tCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGaa                                  >  1:1923904/1‑42 (MQ=255)
                                                      cACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAAt         >  1:363118/1‑66 (MQ=255)
                                                        cGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAGt  <  1:421066/71‑1 (MQ=255)
                                                         gATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGAc                          <  1:1215224/46‑1 (MQ=255)
                                                           |                                                                   
TCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAGT  >  W3110S.gb/3315080‑3315206

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: