Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,447,843 A→G 76.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,8430AG76.5% 99.4 / 28.0 51V1163V (GTT→GTCrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (4/8);  new base G (17/22);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TACGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGACGACGTGCTTCGAACAGGTCCGCAA  >  W3110S.gb/3447775‑3447910
                                                                    |                                                                   
tACGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGAc                                                                   <  1:1018742/71‑1 (MQ=255)
  cGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGa                                                                 <  1:1308855/71‑1 (MQ=255)
  cGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGa                                                                 <  1:2682508/71‑1 (MQ=255)
  cGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGa                                                                 <  1:2548153/71‑1 (MQ=255)
  cGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGa                                                                 <  1:3005727/71‑1 (MQ=255)
  cGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGa                                                                 <  1:1738025/71‑1 (MQ=255)
  cGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGa                                                                 <  1:1343014/71‑1 (MQ=255)
  cGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGa                                                                 <  1:1548160/71‑1 (MQ=255)
  cGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGa                                                                 <  1:1576260/71‑1 (MQ=255)
   ggATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGa                                                                 <  1:2210568/70‑1 (MQ=255)
         cTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCt                                                          <  1:740823/71‑1 (MQ=255)
         cTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCt                                                          <  1:2097926/71‑1 (MQ=255)
         cTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCt                                                          <  1:3043493/71‑1 (MQ=255)
               tCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTc                                                    <  1:2374225/71‑1 (MQ=255)
                     ggggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTAGCGAAGGAGACGATACCGCTGATTTCAGCCAg                                              >  1:700654/1‑71 (MQ=255)
                     ggggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                              <  1:2287215/71‑1 (MQ=255)
                     ggggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                              <  1:350197/71‑1 (MQ=255)
                     ggggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                              <  1:341890/71‑1 (MQ=255)
                     ggggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAg                                              <  1:2209004/71‑1 (MQ=255)
                      gggTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                              <  1:2193820/70‑1 (MQ=255)
                           aTAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATAcag                                                            >  1:197969/1‑49 (MQ=255)
                           aTAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                              >  1:187588/1‑65 (MQ=255)
                           aTAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                              >  1:2513635/1‑65 (MQ=255)
                           aTAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAg                                              >  1:3452886/1‑65 (MQ=255)
                           aTAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGc                                        >  1:2829320/1‑71 (MQ=255)
                              aCCAGACGACGTTTACCTTTGGTTTCTTTACCTAAGGAGAc                                                                   <  1:1203226/41‑1 (MQ=38)
                              aCCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGAc                                                                   <  1:667493/41‑1 (MQ=255)
                              aCCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGa                                            >  1:3596306/1‑64 (MQ=255)
                              aCCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGa                                            >  1:1844061/1‑64 (MQ=255)
                              aCCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGa                                            >  1:1557854/1‑64 (MQ=255)
                                   acgacgTTTACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCtt                                <  1:1668842/71‑1 (MQ=255)
                                         tttACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCt                                                          <  1:70011/39‑1 (MQ=255)
                                           tACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgacga                        >  1:2234517/1‑71 (MQ=255)
                                           tACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgacga                        >  1:1047697/1‑71 (MQ=255)
                                           tACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgacga                        >  1:1525708/1‑71 (MQ=255)
                                           tACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgacga                        >  1:3396618/1‑71 (MQ=255)
                                           tACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgacga                        >  1:3289231/1‑71 (MQ=255)
                                           tACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgacga                        >  1:1821218/1‑71 (MQ=255)
                                           tACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgacga                        >  1:2887963/1‑71 (MQ=255)
                                           tACCTTTGGTTTCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGgacga                        >  1:2749240/1‑71 (MQ=255)
                                           tACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGAt                                                       >  1:2773008/1‑40 (MQ=255)
                                           tACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGAt                                                       >  1:188373/1‑40 (MQ=255)
                                           tACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGAt                                                       >  1:223973/1‑40 (MQ=255)
                                                      tCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGc                                    <  1:2552807/48‑1 (MQ=255)
                                                      tCTTTACCGAAGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGc                                    <  1:572836/48‑1 (MQ=255)
                                                      tCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGACGACGTGCTTCGaa             >  1:1386404/1‑71 (MQ=255)
                                                               aaGGAGACGATACTGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGACGACGTGCTTCGAACAGGTCCGc    <  1:242906/71‑1 (MQ=255)
                                                               aaGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGACGACGTGCTTCGAACAGGTCCGc    <  1:3411231/71‑1 (MQ=255)
                                                               aaGGAGACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGACGACGTGCTTCGAACAGGTCCGc    <  1:96774/71‑1 (MQ=255)
                                                                 ggAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGACGACGTGCTTCGAACAGGTCCGCaa  <  1:3321145/71‑1 (MQ=255)
                                                                 ggAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGACGACGTGCTTCGAACAGGTCCGCaa  <  1:2449621/71‑1 (MQ=255)
                                                                    |                                                                   
TACGGATCGCTACCGTCTACCGGGGTGATAACCAGACGACGTTTACCTTTGGTTTCTTTACCGAAGGAAACGATACCGCTGATTTCAGCCAGGATTGCCGGCTCTTTCGGACGACGTGCTTCGAACAGGTCCGCAA  >  W3110S.gb/3447775‑3447910

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