Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,447,843 A→G 76.9% V1163V (GTT→GTC rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,447,8430AG76.9% 49.8 / 11.6 26V1163V (GTT→GTCrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (2/18);  total (5/21)
Fisher's exact test for biased strand distribution p-value = 6.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TACGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGAC  >  W3110S.gb/3447775‑3447886
                                                                    |                                           
tACGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGAc                                           <  1:144186/71‑1 (MQ=255)
tACGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGAc                                           <  1:1846552/71‑1 (MQ=255)
            ccGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGAt                               <  1:1308540/71‑1 (MQ=255)
            ccGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGAt                               <  1:529216/71‑1 (MQ=255)
            ccGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGAt                               <  1:1534695/71‑1 (MQ=255)
            ccGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGAt                               <  1:167243/71‑1 (MQ=255)
                    cGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGAt                               <  1:1885094/63‑1 (MQ=255)
                      gggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAgg                     <  1:1181442/71‑1 (MQ=255)
                      gggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAgg                     <  1:1535859/71‑1 (MQ=255)
                      gggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAgg                     <  1:179680/71‑1 (MQ=255)
                      gggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAgg                     <  1:386632/71‑1 (MQ=255)
                            tAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGcc               <  1:527325/71‑1 (MQ=255)
                              aCCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCgg             >  1:108548/1‑71 (MQ=255)
                              aCCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCgg             >  1:1073362/1‑71 (MQ=255)
                              aCCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAg                      <  1:779510/62‑1 (MQ=255)
                              aCCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAg                      <  1:484053/62‑1 (MQ=255)
                              aCCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAg                      <  1:1328736/62‑1 (MQ=255)
                               ccAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGAt                               <  1:50575/52‑1 (MQ=255)
                                cAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGct           >  1:1826474/1‑71 (MQ=255)
                                cAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGct           >  1:437241/1‑71 (MQ=255)
                                cAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGct           >  1:1187178/1‑71 (MQ=255)
                                      acgTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCg     <  1:1846456/71‑1 (MQ=255)
                                      acgTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCg     <  1:546277/71‑1 (MQ=255)
                                      acgTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCg     <  1:709820/71‑1 (MQ=255)
                                         tttACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgac  <  1:1667839/71‑1 (MQ=255)
                                         tttACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgac  <  1:1426793/71‑1 (MQ=255)
                                                                    |                                           
TACGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGAC  >  W3110S.gb/3447775‑3447886

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