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,5420CA94.3% 128.6 / NA 56P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (27/23);  total (29/27)
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

TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAT  >  NC_000913/4184402‑4184677
                                                                                                                                            |                                                                                                                                       
gcgggcgtgtagaCGCCGGATCTAGCCTGGTGACAAGATGGTAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGc                                                                                                                                       <  2:878603/131‑1 (MQ=255)
 cTGGCGGTTGAACGCCGTCACGAACCGGGTCACGAGATGGCAGGTCGTGATGGTGAGAGGCGTGGAACTTTCAAGAACGACGCGAGCGAAGATATGCCTTACGATGAAAACGGTAGGGCGGTAGACATCGTACTGAACCAGCt                                                                                                                                      <  2:682466/143‑1 (MQ=255)
     cggtgcagggCCGTAGCCAGCCTGGGGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCGAAGAGCATCCCGATCGAAGATATGCCTTACGATGATAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                  <  2:504299/134‑1 (MQ=255)
      ggTGAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAAGGTCTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTAAAGCAGCTggggg                                                                                                                                 >  2:142389/1‑141 (MQ=255)
        ttAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTa                                                                                                                               >  1:567998/1‑143 (MQ=255)
             cGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGCAGACATCGTACTGAACCAGCTGGGCGTACCgcc                                                                                                                          >  2:63752/1‑141 (MQ=255)
              gtggTATGCAGCCTGGTGACGAGATGGCAGGTTGTCAAGGTAACAAGGGTGTAATTTCTAAGATCCACCCGATGGAAGATATGCCTGACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtct                                                                                                                         <  2:575919/140‑1 (MQ=255)
              gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGGACCAGCTGGGCGTACCGtct                                                                                                                         >  1:173818/1‑143 (MQ=255)
               ccGTAGCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                        <  1:180230/143‑1 (MQ=255)
                  tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTa                                                                                                                     <  1:609189/143‑1 (MQ=255)
                    gccAGCCTTGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATg                                                                                                                   <  1:516044/142‑1 (MQ=255)
                      cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTCACAAGGGTGTAATTTCTCAGATCAACCCGATCGAAGATATGCCTTACGATGAGAACGGTACGCCGGTAGACAGCGCACTGAACCAGCTGCGCGTACAGCATCGTACGac                                                                                                                 >  2:52510/1‑142 (MQ=255)
                      cAGACTGGTGACAAGATGGCACGTCGTCACGTTAACAAGGGTGTAATTTCGAAGAGCAACCCGATCGAAGATATGCCGTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                                 <  2:689232/143‑1 (MQ=255)
                       aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAc                                                                                                                >  1:562326/1‑143 (MQ=255)
                       aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGCAc                                                                                                                >  1:561945/1‑143 (MQ=255)
                         ccTGGTGACAAGATGGCAGGTCGGCACGGTAACAAGGGTGAAATTTCTAAGATCAACCCGATCGAGGGCATGCCTTACGTTGAAAACGGTACGCCGGTAGACCTCGGACGGAACGCGCTGGGCGCACCggcccgcaggcgcag                                                                                                              >  2:889645/1‑129 (MQ=255)
                          cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATc                                                                                                             >  1:791858/1‑143 (MQ=255)
                           tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAGCGTCTCGTATGAACAGCg                                                                                                            >  1:591152/1‑143 (MQ=255)
                               gACAAGATGGCAGGTCGTCACGGTAACAAGTGTGTAATTTTTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCa                                                                                                        <  1:863774/143‑1 (MQ=255)
                                 cAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAGCATCGGTCAGa                                                                                                      >  2:194937/1‑143 (MQ=255)
                                  aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTAAGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAGCGTCGCGTATGAACGTCGGTCAGac                                                                                                     >  2:385706/1‑142 (MQ=255)
                                   ggATGGCAGGGCGTCGCGGTAACAAGGGTGTAATTTCTAAGATCCACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                                                    <  2:481878/142‑1 (MQ=255)
                                     ctGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCt                                                                                                  <  2:174520/142‑1 (MQ=255)
                                     aTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGCAAACGGTACGCCGGTAGACATCGTACTGAACCAGCGGGGCGTACCGCCTCGTATGAACATCGGTCAGATCCt                                                                                                  >  1:632247/1‑143 (MQ=255)
                                              cgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGTCATCGTACTGACCCAGCTGGGAGTAAAGACTTGTATGAACATCGGTCAGATCCTCGAAACCCa                                                                                         <  1:640644/143‑1 (MQ=255)
                                                tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcc                                                                                       >  1:869785/1‑143 (MQ=255)
                                                tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcc                                                                                       >  1:233010/1‑143 (MQ=255)
                                                  aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTg                                                                                     >  1:563510/1‑143 (MQ=255)
                                                  aCGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTCCCGTCTCGTATGAACATCGGTCAGAACGTCGAAACCCTccgg                                                                                     >  1:934220/1‑141 (MQ=255)
                                                     gTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACTGTCTCGTATGAACATCGGTCAGATCCTCGACACCCACCTGGGt                                                                                  >  1:524631/1‑143 (MQ=255)
                                                          aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGATAAGCCGGAAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCAGGGTATGGc                                                                             >  2:863874/1‑143 (MQ=255)
                                                             ggTTTAATTTCTAAGATTAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCTTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGc                                                                          <  2:243566/143‑1 (MQ=255)
                                                                 tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                      >  2:632433/1‑143 (MQ=255)
                                                                  aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                     <  2:288734/143‑1 (MQ=255)
                                                                  aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGGATGAACATCGGTCAGAACCTGGAAACCCACCTGGGTATGGCTGCGacgg                                                                     >  1:839828/1‑140 (MQ=255)
                                                                  aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGAAATCGTACTGCTCCATCTGGGCGTCCCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCCCCTGGGTATGGCTGCgcgag                                                                     >  1:131645/1‑139 (MQ=255)
                                                                   aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAgg                                                                    >  1:457071/1‑143 (MQ=255)
                                                                   aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGAGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGACACCCACCGGGGTCTGGCTCCGAAAgg                                                                    >  2:547642/1‑143 (MQ=255)
                                                                             tCAACCAGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTCAACCAGCTGGGCGTACCGTCTCGTGTGAACATCGATCAGATCCTCGCAAGCCACACGGGGATGGCGGGGACGGGGATCGGCGAc                                                          >  2:460357/1‑143 (MQ=255)
                                                                              cAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                 <  2:229504/103‑1 (MQ=255)
                                                                              cAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                 >  1:229504/1‑103 (MQ=255)
                                                                               aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGCACCACGGACAGAACCTCGGAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACaa                                                        >  2:229494/1‑143 (MQ=255)
                                                                               aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTCTGGCTGCGAAAGGTATCGGCGACaa                                                        >  2:605771/1‑143 (MQ=255)
                                                                                     accGAAGATATGCCTCACGATGAAAGCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCaa                                                  <  1:264155/141‑1 (MQ=255)
                                                                                          tggtctGCCTTGCGATGAATGCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCa                                             <  2:221329/138‑1 (MQ=255)
                                                                                                  ctgaCGATGAGAGCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                     <  2:706901/140‑1 (MQ=255)
                                                                                                  cgtACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                     <  2:570395/141‑1 (MQ=255)
                                                                                                       gatgaagcgggTACGCCGGTAGACATCGTACTGAACCGGCTGGGGGCACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                                <  2:658375/134‑1 (MQ=255)
                                                                                                       gATGAATGCGGCACGCCGGTGAACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGCTCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                                <  2:410857/143‑1 (MQ=255)
                                                                                                             aaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTCTGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCCTGCTGCAACCGCAGCacga                          >  1:723928/1‑143 (MQ=255)
                                                                                                                   aCGCCGGGAGCCATCGTACTGAACCGGCTTGGCGTACCGTCTCGTATGAACATCGGTCATATCTTTGTAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc                    <  1:758075/143‑1 (MQ=255)
                                                                                                                     gCCGGTCGACAACGTACTGAACCCGCTCGGCGTTCCGTCTCGTCTGAACATCGGTCAGATCCTCGAAACCCACCTGGCTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGa                  <  2:802607/143‑1 (MQ=255)
                                                                                                                           aGACATCGTACTGAACCAGCGGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgc            <  1:669393/143‑1 (MQ=255)
                                                                                                                              cATCGTACTGATCCAGCTTGGTGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGcgcg         <  1:185972/143‑1 (MQ=255)
                                                                                                                                ccggactgcgcccGCTGCGCGCCCCGTCGCGGATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa       <  2:87284/133‑1 (MQ=255)
                                                                                                                                    aCTGAACCAGCTGGGCGTGGCGACTCGGATGAACATCGGTCAGATCGTCTAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCa   <  1:830598/143‑1 (MQ=255)
                                                                                                                                     cTGAACCGGCCGGGCGTACCGTCTCGTATGAGCAGCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAt  <  2:85416/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                       
TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAT  >  NC_000913/4184402‑4184677

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: