Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,542 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,5420CA97.7% 121.5 / NA 45P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (28/15);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

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

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