Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,447,843 A→G 77.5% 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,8430AG77.5% 80.6 / 19.4 40V1163V (GTT→GTCrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base G (12/19);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGACGACGTGCTTCGAAC  >  W3110S.gb/3447779‑3447900
                                                                |                                                         
gATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGa                                                       <  1:539765/69‑1 (MQ=255)
    gCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGa                                                       <  1:2548082/65‑1 (MQ=255)
    gCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGa                                                       >  1:1516592/1‑65 (MQ=255)
          gTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGa                                              <  1:3213289/68‑1 (MQ=255)
          gTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGa                                              <  1:2391102/68‑1 (MQ=255)
          gTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGAt                                             >  1:15359/1‑69 (MQ=255)
          gTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGAt                                             >  1:431037/1‑69 (MQ=255)
           tCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGa                                              <  1:2469705/67‑1 (MQ=255)
           tCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGa                                              <  1:2204722/67‑1 (MQ=255)
               ccGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGAc                                                         <  1:2813932/52‑1 (MQ=255)
               ccGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGAc                                                         <  1:157098/52‑1 (MQ=255)
               ccGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGAc                                                         <  1:1643833/52‑1 (MQ=255)
               ccGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCa                                         >  1:3558308/1‑68 (MQ=255)
                 ggggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGa                                                       >  1:857475/1‑52 (MQ=255)
                 ggggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGcc                                      >  1:3123578/1‑69 (MQ=255)
                   ggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAg                                    <  1:2098085/69‑1 (MQ=255)
                    gTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                    <  1:1053854/68‑1 (MQ=255)
                       aTAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCt                                                >  1:1869084/1‑53 (MQ=255)
                         aaCCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCa                                         >  1:1101653/1‑58 (MQ=255)
                          aCCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGcc                             <  1:315811/69‑1 (MQ=255)
                             aGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGc                          >  1:1162230/1‑69 (MQ=255)
                               acgacgTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGcag                                     >  1:210791/1‑54 (MQ=255)
                               acgacgTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                    >  1:225930/1‑57 (MQ=255)
                               acgacgTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                    >  1:854992/1‑57 (MQ=255)
                                   cgTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCttt                     >  1:3048024/1‑68 (MQ=255)
                                    gTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCg                   <  1:2048273/69‑1 (MQ=255)
                                    gTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCg                   <  1:1630963/69‑1 (MQ=255)
                                       tACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgac                <  1:2373179/69‑1 (MQ=255)
                                       tACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgac                <  1:3709876/69‑1 (MQ=255)
                                       tACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgac                >  1:3165641/1‑69 (MQ=255)
                                       tACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGa                                  >  1:1968473/1‑51 (MQ=255)
                                          cTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgacga              <  1:1001172/68‑1 (MQ=255)
                                          cTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgacga              <  1:2132979/68‑1 (MQ=255)
                                          cTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgacga              <  1:1411696/68‑1 (MQ=255)
                                           tttGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGacgacg            <  1:984406/69‑1 (MQ=255)
                                                 ttCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                    <  1:3163766/39‑1 (MQ=255)
                                                 ttCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                    <  1:2324266/39‑1 (MQ=255)
                                                 ttCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                    <  1:830171/39‑1 (MQ=255)
                                                  tCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                    >  1:362768/1‑38 (MQ=255)
                                                        ccGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGACGACGTGCTTCGAAc  <  1:1183900/66‑1 (MQ=255)
                                                                |                                                         
GATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGACGACGTGCTTCGAAC  >  W3110S.gb/3447779‑3447900

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: