Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,543 C→A P1100Q (CCG→CAG)  rpoB → RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5430CA96.9% 90.6 / NA 32P1100Q (CCG→CAG) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (18/13);  minor base G (0/1);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTC  >  NC_000913/4184403‑4184676
                                                                                                                                            |                                                                                                                                     
tCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTCCCATGAAAACGGTACGCCGGTAGACAACGTACTGAACCAGc                                                                                                                                     >  2:165625/1‑143 (MQ=255)
   ggCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCCgg                                                                                                                                  >  2:388009/1‑143 (MQ=255)
      gggtAAACGCCGTAGCCAGCCTGGGGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCg                                                                                                                               <  2:275968/140‑1 (MQ=255)
             cgatgTATCCAGGCTGGTGAGAAGCGGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGCTAGACATCGTAATGAACCAGCTGGGCGTACCGtc                                                                                                                        <  2:741481/139‑1 (MQ=255)
             cGCCGTATCCAGCCTGGCGACGAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTAAGCCGGGAGACAGCGTACTGAACCAGCTGGGGGTACGGtg                                                                                                                        >  2:786620/1‑142 (MQ=255)
                  tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGAGTAATTTCTAAGATCAACCCGCTCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACCCAACCAGCTCTGCGTACCGTCTCGtg                                                                                                                   >  1:52670/1‑142 (MQ=255)
                   aTCCAGCCTGGTGACAAGATGGCAGGTCGTCACGTTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                                  <  1:288496/143‑1 (MQ=255)
                    tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATg                                                                                                                 <  1:701250/143‑1 (MQ=255)
                     ccAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAGCGTCTCGTATgg                                                                                                                >  1:779130/1‑142 (MQ=255)
                                    gttGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGCACAACCCGATCGAAGATATTCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAGGCACATCGGTCAGATcc                                                                                                 >  1:436580/3‑143 (MQ=255)
                                    gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTACGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGACCATCGGTCAGATcc                                                                                                 >  1:498706/1‑143 (MQ=255)
                                      tGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGCATGTTATTAACATCGGTCAGATCCTc                                                                                               <  2:699538/143‑1 (MQ=255)
                                          aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                           >  2:203377/1‑143 (MQ=255)
                                                    ggTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                                                                 >  1:343316/1‑143 (MQ=255)
                                                               tgtAAGTTCTAAGATCAAGCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGa                                                                      <  2:305184/143‑1 (MQ=255)
                                                                      tCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAt                                                               <  1:396212/143‑1 (MQ=255)
                                                                         acgATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCgg                                                            <  2:322206/141‑1 (MQ=255)
                                                                            aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCCCGGAACCCACCCGGGGAGGGCTGCGCACGGTATGGGCGa                                                         >  2:89709/1‑143 (MQ=255)
                                                                                 cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAGAGGTATCGGCGACAAGa                                                    >  1:446145/1‑143 (MQ=255)
                                                                                                  cTTACGATGAAAACGGTCCGCCGATAAAAAACACACTGTTCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAACCCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                   <  1:751755/143‑1 (MQ=255)
                                                                                                     acgcggAAAGCGGTGCGCCGGGAGGCAGCGTGCTGAGCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCGGAGCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag                                <  2:798681/138‑1 (MQ=255)
                                                                                                      cgacgaACACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGGATGAACATCGGTCAGATCCGCGAACCCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc                               <  2:230087/143‑1 (MQ=255)
                                                                                                       gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAAAAGATCAACGCCAGGCAGAAAcagca                              >  1:536445/1‑143 (MQ=255)
                                                                                                          gAAAACGGTACGCCGGTAGACATCGGACTGAACCATCGGGGCGGACGGTCTCGTATGAACAGCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGGTCAACGCCATGCTGAAACAgcagca                           >  1:799817/1‑143 (MQ=255)
                                                                                                           aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTATCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTCTGGCTGCGAACGGTCTCGGCGACAAGCTCAAAGCCATGCTGAAACAgcagcac                          >  1:557474/1‑142 (MQ=255)
                                                                                                             aaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaaga                        >  1:500210/1‑143 (MQ=255)
                                                                                                                  tACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcg                   <  2:376812/143‑1 (MQ=255)
                                                                                                                             aCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcgc        <  1:308761/143‑1 (MQ=255)
                                                                                                                                 cGTACTGAGGCGGCTGGCCGTACCGTCTCGTATGACCATCGGTCAGTTCCTCTAAACCCACCTTGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAt    <  2:44526/143‑1 (MQ=255)
                                                                                                                                  gTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAACTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAACAAGAAGTAGCGAAACTGCCCGaaat   >  2:572005/1‑141 (MQ=255)
                                                                                                                                   tACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACCGCGCGAATTc  >  2:233858/1‑143 (MQ=255)
                                                                                                                                   tACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCGGAGACAGCAGCCACAAGTCGCGAAACTGCgagacgtg  >  1:210389/1‑139 (MQ=255)
                                                                                                                                            |                                                                                                                                     
TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTC  >  NC_000913/4184403‑4184676

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑