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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5430CA97.7% 120.1 / NA 44P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (28/14);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCC  >  NC_000913/4184403‑4184680
                                                                                                                                            |                                                                                                                                         
tCTGGCGGTTAAACGCCGTATCCAGCCTGGAGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTCAGATCGACCCGATCGAAGGTATGCCTTACGATGCAAACGGTACGCCGGTGGACAGCGTACTGAACCAGc                                                                                                                                         >  1:320735/1‑143 (MQ=255)
  tGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTg                                                                                                                                       <  1:413636/143‑1 (MQ=255)
   ggCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGCACCAGCTgg                                                                                                                                      >  2:389919/1‑143 (MQ=255)
     cGGTTAAACGCCGTATCCAGCCTGGTGGCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                    <  1:298517/143‑1 (MQ=255)
     cGGTTAAACGCCGTATCCAGACTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTCGACATCGTACTAAACCAGCTGGGc                                                                                                                                    >  2:145159/1‑143 (MQ=255)
      ggTTAAACGCCGTATCCAGCATGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGCTATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACGAGCTgggag                                                                                                                                   >  2:502971/1‑141 (MQ=255)
             cGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtc                                                                                                                            >  1:65571/1‑143 (MQ=255)
                cGTATCCAGCCTGGTGACAAGATGGCAGGGCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCg                                                                                                                         <  2:87656/143‑1 (MQ=255)
                  taaccAGCCTGATGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTa                                                                                                                       <  1:453783/140‑1 (MQ=255)
                       aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTCCGCCGGTAGACATCGTCCTGACCCACCCGGGCGTACCGCCTCGTATGAAc                                                                                                                  >  2:118944/1‑143 (MQ=255)
                            ggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCgg                                                                                                             >  1:334546/1‑143 (MQ=255)
                            ggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCgg                                                                                                             >  1:533057/1‑143 (MQ=255)
                                aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGCCTCGTATGAACATCGGGCAg                                                                                                         >  2:247584/1‑143 (MQ=255)
                                  aaGATGGCAGGTCGTCACGGTAACAAGGGTGTCATTTCTAAGATCAACCCGATCGAAGATATGCCAAACGAAAAAAACAGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAGCGGTCAGAt                                                                                                       >  1:444880/1‑143 (MQ=255)
                                    gagggCAGGTCGTCACGGTAATAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATcc                                                                                                     <  2:454850/140‑1 (MQ=255)
                                          aGGTCGTCACGGTAACAAGTGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                               <  2:55548/143‑1 (MQ=255)
                                                tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACACCGGTCAGACCATCGAAACCGAcc                                                                                         >  2:316512/1‑143 (MQ=255)
                                                tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTACGAACATCGGCCACATCCGCGAAACCCAcc                                                                                         >  2:484889/1‑143 (MQ=255)
                                                               tgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCCGACCCTCGACCCCCACCGGGGTATGCATGCGa                                                                          >  2:611846/1‑143 (MQ=255)
                                                                    tttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCAAACTTGTAATGGACGCGGAAGGt                                                                     >  2:616553/1‑143 (MQ=255)
                                                                    gttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGt                                                                     <  2:539577/142‑1 (MQ=255)
                                                                      tCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCGCGTATGAACATCGGTCAGATCCTCGAGACCGACCTGGGTATGGCTGCGAAAGGTAt                                                                   >  1:281358/1‑143 (MQ=255)
                                                                           gATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGACTCGTATGAACATCGGTCAGATCATCGAAACCCACCTGGGCATGGCTGCGAAAGGTATCGGCg                                                              >  2:643792/1‑143 (MQ=255)
                                                                           gATCAACCCGATCGAAGACAGGCGGTACGATGAAATCGGTACGCCGGTAGACATCGGACTGAACCAGCTGGGCGTACCGGCGCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCg                                                              <  2:398060/143‑1 (MQ=255)
                                                                             tCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGCACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCACGAAACCCACGTGGGTATGGCTGCGAAGGGAATCGGCggc                                                            >  2:609037/1‑141 (MQ=255)
                                                                                 cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCCGCGTATGGCTGCGAAAGGTATCGGCGACACGa                                                        >  2:367540/1‑143 (MQ=255)
                                                                                  ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACAACGTACTGAACCAGCTGGGCGTACCGTCGCGCATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCAGCGAAAGGTCTCGGCGACAAGAt                                                       >  2:27929/1‑143 (MQ=255)
                                                                                     gtCGAAGATATGCCTTACGATGAGAGGGGTGCGCGGGTAGACAGCGTACTGAACCGGCTGGGCGTAGCGCAACACTTTATTATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCaa                                                    <  2:326293/142‑1 (MQ=255)
                                                                                         aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                                >  1:160866/1‑143 (MQ=255)
                                                                                                 ccTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaa                                        <  1:261712/143‑1 (MQ=255)
                                                                                                 ccTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACACACCCGGGTAGGGCTGCGAAAGGTATCGGCGACAAGCTCACCGACATGcaggc                                        >  2:574390/1‑139 (MQ=255)
                                                                                                  cTTACGATGAAACCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                       <  1:161059/143‑1 (MQ=255)
                                                                                                  cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTTAGCGTACCGTCTCGTATAAACATCAGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAATGTATCGGCGACAAGATCAACGCCATGCTGaaa                                       >  1:31185/1‑143 (MQ=255)
                                                                                                  cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTGTCGTATGAACATCGGTCAGATCCTCGAACCCCCCCTGGGCATGGCTGCGACAGGCCTCGGCGACAAGCCCACCGCCGTGCCGaac                                       >  2:193671/1‑142 (MQ=255)
                                                                                                  cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCAAAAAGATCAAAGACATGCTGaaa                                       >  1:528368/1‑143 (MQ=255)
                                                                                                   ttACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAc                                      >  2:384163/1‑143 (MQ=255)
                                                                                                     accaTGAAAACGGTATGCGGGTAGGCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTGAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag                                    <  2:363639/140‑1 (MQ=255)
                                                                                                        aTGAAAACGGTACGCCGGTAGACATCGGACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagcag                                 <  1:533522/143‑1 (MQ=255)
                                                                                                            aaaCGGGACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                             <  1:100830/143‑1 (MQ=255)
                                                                                                              ccGGGACGCCGGTAGACAGCGTACTGAACCAGCTGCGCGGACCGTCTCGTAGGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                           <  2:452620/142‑1 (MQ=255)
                                                                                                               cGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCCGGGTATGGCTGCGAAAGGTATCGGCGACTAGATCAACGCCATGCTGAAACAGCAGCaagaga                          >  2:473148/1‑141 (MQ=255)
                                                                                                                            gACAGCGGACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcg             <  1:261188/143‑1 (MQ=255)
                                                                                                                                      tGAACCAGCTGGGCGTACCGTCTCGTACGAACCTCCATGGTATGTGCGAAAGCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGTTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATc   >  1:285423/1‑143 (MQ=255)
                                                                                                                                       gAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATcc  >  2:264173/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                         
TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCC  >  NC_000913/4184403‑4184680

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: ‑