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.0% 122.7 / NA 52P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (27/21);  total (28/24)
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

CTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAT  >  NC_000913/4184404‑4184678
                                                                                                                                           |                                                                                                                                       
cTGGCGGTTGAACGCCGTCACGAACCGGGTCACGAGATGGCAGGTCGTGATGGTGAGAGGCGTGGAACTTTCAAGAACGACGCGAGCGAAGATATGCCTTACGATGAAAACGGTAGGGCGGTAGACATCGTACTGAACCAGCt                                                                                                                                      <  2:682466/143‑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)
                         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:233010/1‑143 (MQ=255)
                                               tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcc                                                                                       >  1:869785/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                                                                                                 >  1:229504/1‑103 (MQ=255)
                                                                             cAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                 <  2:229504/103‑1 (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)
                                                                                                      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)
                                                                                                                                   aCTGAACCAGCTGGGCGTGGCGACTCGGATGAACATCGGTCAGATCGTCTAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCa   <  1:830598/143‑1 (MQ=255)
                                                                                                                                    cTGAACCGGCCGGGCGTACCGTCTCGTATGAGCAGCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAt  <  2:85416/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                       
CTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAT  >  NC_000913/4184404‑4184678

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