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,5430CA96.3% 148.7 / NA 54P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (25/27);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCC  >  NC_000913/4184405‑4184680
                                                                                                                                          |                                                                                                                                          
tGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACACGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTCCGCCGGCAGACATCGTATTCAGCCAGCTg                                                                                                                                        >  1:558700/1‑143 (MQ=255)
   cGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                     <  1:128370/143‑1 (MQ=255)
   cGGTGAAACGCCGTCGGCACCCTGGTGACGAGATGGCAGGGCGTCGCGGTAACAAGGGTGTACTTTCTCAGATCAGCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                     <  2:449709/143‑1 (MQ=255)
      ttAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTa                                                                                                                                  >  1:291080/1‑143 (MQ=255)
       taggcgcgGTATCCCGCCTCGTGCCAAGATGGCAGGTCGTGTCGGTAACAAGGGTGTAATCTCTAAGATCAACCCGATCGAAGACATGCCTTACGATGAAAACGGTACGCCGGTCGACATCGTACTGAACCAGCTGGGCGTAc                                                                                                                                 <  2:204949/139‑1 (MQ=255)
        gaaCGCCGTGTCCAGGCTGGTGACAAGATTGTAGGTCGTCACGGTAACAAGGGTGTAATTTCTAGGATCAGCGCGAGCGATGATATGCCTTAGGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAcc                                                                                                                                <  2:291359/142‑1 (MQ=255)
                  tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTAAGGAACCAGCTGGGCGTACCGTCTCGTagt                                                                                                                      >  1:222471/1‑141 (MQ=255)
                      gCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                                                                 >  2:761970/1‑143 (MQ=255)
                           gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                            <  2:240901/143‑1 (MQ=255)
                                 ataTGGCAGGTCGTCAGGGTAACAAGGGTGTAATTTCTAAGGTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                                                      <  2:134769/141‑1 (MQ=255)
                                 aGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATCAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCCGGACGTACCGGCTCGTAGAACGATCCCTCAGATc                                                                                                      >  2:248623/1‑143 (MQ=255)
                                 aGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCCTACGATGAAAACGGTACGCCGGTATACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                                                      <  1:612875/143‑1 (MQ=255)
                                  gATGGCAGGTCGTCACGGTAACAAGTGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTAGTATGAACATCGGTCAGCTcc                                                                                                     >  1:278244/1‑143 (MQ=255)
                                  gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTCACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCAGCTGGGCGTACCGTCTCGTCTGAACATCGGTCCGATcc                                                                                                     >  1:456892/1‑143 (MQ=255)
                                    ggctgggtggtCACGGTAAGAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAATCGGTACGCCGTTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                   <  2:302476/138‑1 (MQ=255)
                                      gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATAGTACTGAACCAGCTGGGCGTCCCGACTCGTATGAACATCGGGCAGATCCTCGa                                                                                                 >  2:103889/1‑143 (MQ=255)
                                       cAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTCCTGAACCAGCTGGGCGTACCGTCCCGTATGAACATCGGTCAGATCCTCGag                                                                                                >  2:775576/1‑142 (MQ=255)
                                        aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                               <  2:434561/143‑1 (MQ=255)
                                          gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAcc                                                                                             >  1:23213/1‑143 (MQ=255)
                                          gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGAACGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGGCAGATCCTCGAAAcc                                                                                             >  1:385871/1‑143 (MQ=255)
                                            cgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGGCAGATCCCCGAAACCCa                                                                                           >  2:744399/1‑143 (MQ=255)
                                             gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGGTGGGCGTACCGTCGCGTATGAAGATCGGTCAGGTCCTCGAAACCCag                                                                                          >  1:315378/1‑142 (MQ=255)
                                               cACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATAGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACcctcct                                                                                        >  1:657104/1‑143 (MQ=255)
                                                aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGCCACCCACCTg                                                                                       >  2:317905/1‑143 (MQ=255)
                                                aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGCCGTACCGCCTCGTATGAACATCGATCAGATCCTCGACACCCACCTg                                                                                       >  2:498566/1‑143 (MQ=255)
                                                 cGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATACTCTAAACCCACCCgg                                                                                      >  1:620570/1‑143 (MQ=255)
                                                 cGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAGGAACATCGGTCGGATCCTCGAAACCCACCTgg                                                                                      <  2:409441/143‑1 (MQ=255)
                                                       cAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGCGTagtg                                                                                >  1:296585/1‑140 (MQ=255)
                                                               tAATTTCTACGATCAACCCGATCGAAGATATGCCTTACGGTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                        <  2:180639/143‑1 (MQ=255)
                                                               tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                        <  1:569304/143‑1 (MQ=255)
                                                               tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                        >  1:153083/1‑143 (MQ=255)
                                                               gaatagcTAAGATAAACTGGATAGAAGAGAGGCGTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTTGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                        <  1:409562/137‑1 (MQ=255)
                                                                aaTTTCTAAACTCACCCCGATAGAAGAAATGCTTTAAGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                       <  1:488631/143‑1 (MQ=255)
                                                                     cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATc                                                                  >  1:756952/1‑143 (MQ=255)
                                                                      tAAGAACAACCCGATCGACGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                 <  1:123840/143‑1 (MQ=255)
                                                                          gtCAACGCGATCGAAGATATGCCTTACGATGACAACGGTACGCCGGTAGACACCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                             <  1:161367/142‑1 (MQ=255)
                                                                           gcAACCCGATCGAAGATATGCCTTACGATGAAATCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGCCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                                                            <  1:744399/142‑1 (MQ=255)
                                                                              gcccGATCGGAGACATGCGTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCGCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAg                                                         <  2:566073/142‑1 (MQ=255)
                                                                              aCCCGATCGACGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAg                                                         <  2:100247/143‑1 (MQ=255)
                                                                               cccGATCGAAGATATGCCTTACGATGAAGACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGCCTCGTATGAACATAGGTCAGATCCTCGAAACCAACATGGGAATGGCTGCGAAAGGTGCCGGCGCCaggc                                                        >  2:157438/1‑140 (MQ=255)
                                                                                 cGATCGAAGATATGCCTTACGAAGAAAACGGTACGCCGGTAGACATAGTACTGAACCAGATGGGAGTACCGCCTCGAACGAACATAGGTCAGAGCCTCGACACCCAACGGGGGAAAGCTGGGAACGGTATCGGAGACACGagc                                                      >  2:683465/1‑141 (MQ=255)
                                                                                    tCGAAGATATGCCTTACGATTAGATCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                   <  2:683936/143‑1 (MQ=255)
                                                                                      gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAGCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGc                                                 <  2:137283/143‑1 (MQ=255)
                                                                                          atatGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAGCATCGGTCAGATCCGCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATg                                             <  1:307095/143‑1 (MQ=255)
                                                                                                cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCCACAAGATCAACGCCATGCTGaaa                                       >  2:456500/1‑143 (MQ=255)
                                                                                                 ttACGATGAAAACGGTACGCCGGTAGACATCGGACTGAACCAGCTGGGCGTGCCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAc                                      >  2:412044/1‑143 (MQ=255)
                                                                                                         aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaa                              <  1:469514/143‑1 (MQ=255)
                                                                                                          caatggTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                             <  2:472347/139‑1 (MQ=255)
                                                                                                          aagcGGTACGCCGGTAGACAGCGTACTGGACCAGCTGGGCGTACCGTCTCGTATGAACCTCGGTCAGATCCTCGAAACCCCCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                             <  2:518238/140‑1 (MQ=255)
                                                                                                                    ccGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATGCTCGAAACCCAGCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCTCGaa                   >  1:384580/1‑143 (MQ=255)
                                                                                                                          gACATCGTACTGAATCAGCTGGGCGTACCGTCTCGGATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcg             <  2:291080/143‑1 (MQ=255)
                                                                                                                               ggTACTGAACGCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCGCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAt        <  2:595648/142‑1 (MQ=255)
                                                                                                                               cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCCTGCGGGAACAGCAGCACGCAGGCGCGAAACTGCGCGaag        >  1:231584/1‑142 (MQ=255)
                                                                                                                                     gAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATcc  <  2:456892/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCC  >  NC_000913/4184405‑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: