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,5420CA96.5% 158.4 / NA 57P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (25/30);  total (26/31)
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

TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCC  >  NC_000913/4184404‑4184679
                                                                                                                                          |                                                                                                                                          
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)
             gcgtttgtaggcTGGTGGCAAGATGGCGGGTCGTCACGGTAACCAGGGTGTAAGTTCTAAGATCAACCCGATCGAAGATCTGCCTTGCGATGAAAACGGTACGCCGGTAGACCTCGTACTGATCCAGCTGGGCGTACCGtctc                                                                                                                           <  2:754058/132‑1 (MQ=255)
                  tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTAAGGAACCAGCTGGGCGTACCGTCTCGTagt                                                                                                                      >  1:222471/1‑141 (MQ=255)
                      gCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑a                                                                                                                 >  2:761970/1‑143 (MQ=255)
                           gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGt                                                                                                            <  2:240901/143‑1 (MQ=255)
                                 ataTGGCAGGTCGTCAGGGTAACAAGGGTGTAATTTCTAAGGTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATc                                                                                                      <  2:134769/141‑1 (MQ=255)
                                 aGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATCAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCCGGACGTACCGGCTCGTA‑GAACGATCCCTCAGATc                                                                                                      >  2:248623/1‑143 (MQ=255)
                                 aGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCCTACGATGAAAACGGTACGCCGGTATACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATc                                                                                                      <  1:612875/143‑1 (MQ=255)
                                  gATGGCAGGTCGTCACGGTAACAAGTGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTAGTATGAAC‑ATCGGTCAGCTcc                                                                                                     >  1:278244/1‑143 (MQ=255)
                                  gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTCACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCAGCTGGGCGTACCGTCTCGTCTGAAC‑ATCGGTCCGATcc                                                                                                     >  1:456892/1‑143 (MQ=255)
                                    ggctgggtggtCACGGTAAGAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAATCGGTACGCCGTTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTc                                                                                                   <  2:302476/138‑1 (MQ=255)
                                      gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATAGTACTGAACCAGCTGGGCGTCCCGACTCGTATGAAC‑ATCGGGCAGATCCTCGa                                                                                                 >  2:103889/1‑143 (MQ=255)
                                       cAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTCCTGAACCAGCTGGGCGTACCGTCCCGTATGAAC‑ATCGGTCAGATCCTCGag                                                                                                >  2:775576/1‑142 (MQ=255)
                                        agggcttccctggaaCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAGCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGaaa                                                                                               <  2:435052/130‑1 (MQ=255)
                                        aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGaaa                                                                                               <  2:434561/143‑1 (MQ=255)
                                          gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAAcc                                                                                             >  1:23213/1‑143 (MQ=255)
                                          gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGAACGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGGCAGATCCTCGAAAcc                                                                                             >  1:385871/1‑143 (MQ=255)
                                            cgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGGCAGATCCCCGAAACCCa                                                                                           >  2:744399/1‑143 (MQ=255)
                                             gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGGTGGGCGTACCGTCGCGTATGAAG‑ATCGGTCAGGTCCTCGAAACCCag                                                                                          >  1:315378/1‑142 (MQ=255)
                                               cACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATAGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACcctcct                                                                                        >  1:657104/1‑143 (MQ=255)
                                                aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGCCACCCACCTg                                                                                       >  2:317905/1‑143 (MQ=255)
                                                aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGCCGTACCGCCTCGTATGAAC‑ATCGATCAGATCCTCGACACCCACCTg                                                                                       >  2:498566/1‑143 (MQ=255)
                                                 cGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATACTCTAAACCCACCCgg                                                                                      >  1:620570/1‑143 (MQ=255)
                                                 cGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAGGAAC‑ATCGGTCGGATCCTCGAAACCCACCTgg                                                                                      <  2:409441/143‑1 (MQ=255)
                                                       cAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGCGTagtg                                                                                >  1:296585/1‑140 (MQ=255)
                                                               tAATTTCTACGATCAACCCGATCGAAGATATGCCTTACGGTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                        <  2:180639/143‑1 (MQ=255)
                                                               tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                        <  1:569304/143‑1 (MQ=255)
                                                               tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                        >  1:153083/1‑143 (MQ=255)
                                                               gaatagcTAAGATAAACTGGATAGAAGAGAGGCGTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTTGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                        <  1:409562/137‑1 (MQ=255)
                                                                aaTTTCTAAACTCACCCCGATAGAAGAAATGCTTTAAGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                       <  1:488631/143‑1 (MQ=255)
                                                                     cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATc                                                                  >  1:756952/1‑143 (MQ=255)
                                                                      tAAGAACAACCCGATCGACGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                 <  1:123840/143‑1 (MQ=255)
                                                                          gtCAACGCGATCGAAGATATGCCTTACGATGACAACGGTACGCCGGTAGACACCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                             <  1:161367/142‑1 (MQ=255)
                                                                           gcAACCCGATCGAAGATATGCCTTACGATGAAATCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGCCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                                                            <  1:744399/142‑1 (MQ=255)
                                                                              gcccGATCGGAGACATGCGTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCGCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAg                                                         <  2:566073/142‑1 (MQ=255)
                                                                              aCCCGATCGACGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAg                                                         <  2:100247/143‑1 (MQ=255)
                                                                               cccGATCGAAGATATGCCTTACGATGAAGACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGCCTCGTATGAAC‑ATAGGTCAGATCCTCGAAACCAACATGGGAATGGCTGCGAAAGGTGCCGGCGCCaggc                                                        >  2:157438/1‑140 (MQ=255)
                                                                                 cGATCGAAGATATGCCTTACGAAGAAAACGGTACGCCGGTAGACATAGTACTGAACCAGATGGGAGTACCGCCTCGAACGAAC‑ATAGGTCAGAGCCTCGACACCCAACGGGGGAAAGCTGGGAACGGTATCGGAGACACGagc                                                      >  2:683465/1‑141 (MQ=255)
                                                                                    tCGAAGATATGCCTTACGATTAGATCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                   <  2:683936/143‑1 (MQ=255)
                                                                                      gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGAGCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGc                                                 <  2:137283/143‑1 (MQ=255)
                                                                                          atatGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAGC‑ATCGGTCAGATCCGCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATg                                             <  1:307095/143‑1 (MQ=255)
                                                                                                cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCCACAAGATCAACGCCATGCTGaaa                                       >  2:456500/1‑143 (MQ=255)
                                                                                                 ttACGATGAAAACGGTACGCCGGTAGACATCGGACTGAACCAGCTGGGCGTGCCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAc                                      >  2:412044/1‑143 (MQ=255)
                                                                                                         aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaa                              <  1:469514/143‑1 (MQ=255)
                                                                                                          caatggTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                             <  2:472347/139‑1 (MQ=255)
                                                                                                          aagcGGTACGCCGGTAGACAGCGTACTGGACCAGCTGGGCGTACCGTCTCGTATGAAC‑CTCGGTCAGATCCTCGAAACCCCCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                             <  2:518238/140‑1 (MQ=255)
                                                                                                                    ccGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATGCTCGAAACCCAGCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCTCGaa                   >  1:384580/1‑143 (MQ=255)
                                                                                                                          gACATCGTACTGAATCAGCTGGGCGTACCGTCTCGGATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcg             <  2:291080/143‑1 (MQ=255)
                                                                                                                             agtgTACTGAACCAGCTGGGCGTCCCGTCGCGTATGAAC‑ATCGGTCAGATCCTCCAAACCCACCTGGGAATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAATTctcgcagtttcgcga          <  1:98064/140‑12 (MQ=255)
                                                                                                                               ggTACTGAACGCGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCGCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAt        <  2:595648/142‑1 (MQ=255)
                                                                                                                               cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCCTGCGGGAACAGCAGCACGCAGGCGCGAAACTGCGCGaag        >  1:231584/1‑142 (MQ=255)
                                                                                                                                     gAACCAGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATcc  <  2:456892/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAAC‑ATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCC  >  NC_000913/4184404‑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: ‑