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,5420CA98.3% 169.6 / NA 59P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (36/22);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCC  >  NC_000913/4184402‑4184679
                                                                                                                                            |                                                                                                                                         
tCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGc                                                                                                                                         >  1:687575/1‑143 (MQ=255)
    gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTggg                                                                                                                                     >  2:635721/1‑143 (MQ=255)
     tgcttAAAAGCCGTATCCAGCCTAGTGACAAGATGGCAGGTTGTCACGTTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                    <  1:30981/140‑1 (MQ=255)
      ggTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCg                                                                                                                                   >  1:636995/1‑143 (MQ=255)
        ttAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTa                                                                                                                                 >  2:516883/1‑143 (MQ=255)
               ccGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                          >  2:615212/1‑143 (MQ=255)
                    tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATAGAAGATATGCCTTACGCTGAAAACGGTACGCCGGTCGTCACCGTCCTGAACCAGCTGGGCGTGCCGGCTCCTaga                                                                                                                     >  2:306503/1‑141 (MQ=255)
                      cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                                   <  1:635700/143‑1 (MQ=255)
                      cAGCCGGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                                   <  1:541800/143‑1 (MQ=255)
                        gCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                                                                 >  2:593388/1‑143 (MQ=255)
                        gCCTGGTGAAAGGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGGCTGACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                                                                 <  1:386635/143‑1 (MQ=255)
                           tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCg                                                                                                              >  1:824060/1‑143 (MQ=255)
                           tGGCGACAAGATCGCAGGTCGTCCCGGTAACAAGGGTGTAATTTCTACGATCAACCCGATCGAAGATATGCCTTCCGATCAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCg                                                                                                              <  1:498750/143‑1 (MQ=255)
                            ggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCgg                                                                                                             >  1:13874/1‑143 (MQ=255)
                               gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCa                                                                                                          >  1:326591/1‑143 (MQ=255)
                                 cacgaggccAGGTCGTCACGCTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGa                                                                                                        <  1:511895/135‑1 (MQ=255)
                                  aaGATGGCAGGTCGTCACGGTAAGAAGGGTGTAATTTCTAAGGTCAACCCGATCGAAGATATGCCTTACGATGCAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                                                       >  1:841923/1‑143 (MQ=255)
                                  aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAGCGGTACGCCGGTAGACATCGTACTGCACCAGCTGGGCGTACCGTCCCGTATGAACATCGGTCAGAt                                                                                                       >  2:128154/1‑143 (MQ=255)
                                  aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                                                       >  1:629684/1‑143 (MQ=255)
                                  aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                                                       >  2:174836/1‑143 (MQ=255)
                                     agagCAGGTCGTCACGTTCACAAGGGTGTAATTTCTAAGCTCAACCCGATCGAAGATATGCCTTACGATGCAAACGGTACGCCTGTAGACAGCGTACTGAACCAGCTGGGAGTACCGTCTCGTATGAACATCGGTCAGATCCt                                                                                                    <  2:473267/140‑1 (MQ=255)
                                     aTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACGTCGGTCAGATCCt                                                                                                    >  1:737259/1‑143 (MQ=255)
                                        gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGa                                                                                                 >  1:633467/1‑143 (MQ=255)
                                             tcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGATCAGATCATCGAAAccc                                                                                            >  2:235339/1‑143 (MQ=255)
                                              cgtcTCGGTAACATGGGTGTAATTTCTAAGATCATCCCGCTCGAACATATGCCTTACGCTGAACCCGGTCCGCCGGTAGACATCGTACTGAACCCGCCGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                                                           <  2:559906/143‑1 (MQ=255)
                                                tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGCACCGTCACGTATGAACATCGGTCAGATCTGCGAAACCCAca                                                                                         >  1:754002/1‑142 (MQ=255)
                                                   cGGTAACCAGGGTGGCTTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAATATCGGTCAGATCCTCGAAACCCACCTgg                                                                                      <  2:251750/143‑1 (MQ=255)
                                                    cgcaaacaGGGTGTCGTTGCTAAGATCAACCCGATCGAAGAGATGNCTTACGATGGAACCGGTACGCCGGTAGCCATCGTACTGAACCAGCTGGGCGTACCGTCGCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                                                                     <  2:481223/136‑1 (MQ=255)
                                                        aCAAGGGTGTAATGTCTAAGACCAACCCGAGCGAAGATAGGCCTTACGATGAAAGCGGTCCGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGCCTCGGATGAACAGCGGGCAGATCCTCGAAACCCACCTGGGTATg                                                                                 <  2:628978/143‑1 (MQ=255)
                                                         cAAGGGTTTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACTATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATgg                                                                                <  2:745036/143‑1 (MQ=255)
                                                          aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGc                                                                               >  1:311197/1‑143 (MQ=255)
                                                             gttgtAATTTCTAAGATCAACCCGATCGAAGAGATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGc                                                                            <  1:294225/141‑1 (MQ=255)
                                                                 tAATTTCTTAGATCAACCCGATCGACGAGATGCCTTACGATGAAAACGGTACGCCGGTAGACATCCTACTGAACCATCTGGGCGTCCCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                        <  2:807769/143‑1 (MQ=255)
                                                                 tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCCGGGGATGGCTGCGagc                                                                        >  2:752169/1‑141 (MQ=255)
                                                                  agtttGTAAGCTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACGAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                       <  2:36834/141‑1 (MQ=255)
                                                                      tCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTCCTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCGGCGTACGGCTGCGAAAGGTAt                                                                   >  2:673357/1‑143 (MQ=255)
                                                                       cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATc                                                                  >  2:897628/1‑143 (MQ=255)
                                                                        tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                 >  1:761927/1‑143 (MQ=255)
                                                                            aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGAACCGGCGa                                                             >  2:714232/1‑143 (MQ=255)
                                                                             tCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                                                            >  1:775423/1‑143 (MQ=255)
                                                                                   cGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATc                                                      <  1:174836/143‑1 (MQ=255)
                                                                                         aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTGTCGGCGACAAGGTCCACGcc                                                >  1:33688/1‑143 (MQ=255)
                                                                                            gtatGCCCTACGGAGAAAGCGGTACGCCGGTAGACATCGTACTGAACCAGCCGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATg                                             <  2:875563/142‑1 (MQ=255)
                                                                                             tatGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACACGCTCAACGCCATGc                                            >  1:488392/1‑143 (MQ=255)
                                                                                                   ttACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCa                                               <  1:874149/134‑1 (MQ=255)
                                                                                                   ttACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCa                                               >  2:874149/1‑134 (MQ=255)
                                                                                                              aCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                           <  2:774438/143‑1 (MQ=255)
                                                                                                                  tACGCCGGCAGCCATCGGACTGAACCAGCTGGGCGTACCGTCTCGTATGAACACCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcg                       <  2:87226/143‑1 (MQ=255)
                                                                                                                   gcGCCGGTAGACCTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATTCTGAAACAGCAGCAAGAAGTcgc                      <  2:11149/142‑1 (MQ=255)
                                                                                                                       cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAACCAGCAGCAAGAAGTCGCGaaa                  >  1:932/1‑143 (MQ=255)
                                                                                                                       cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCCGGGTATGGCTGCGAAAGGTATCGGCGACCAGAGCAACGCCATGCCGAAACAGCAGCCAGGAGTCGCGaaa                  >  2:653438/1‑143 (MQ=255)
                                                                                                                          tAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTg               <  1:615212/143‑1 (MQ=255)
                                                                                                                             aCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcgc            >  1:55689/1‑143 (MQ=255)
                                                                                                                                 cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGCAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTCAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAt        >  1:271727/1‑143 (MQ=255)
                                                                                                                                 cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCCGAGTATGGATGCGAAAGGTATCGGCGACAAGATCAAAGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaaa        >  2:309419/1‑142 (MQ=255)
                                                                                                                                   tACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAGCTGCACGactcc      >  2:714332/1‑139 (MQ=255)
                                                                                                                                    aCTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCa     <  1:226537/143‑1 (MQ=255)
                                                                                                                                      tGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATc   >  1:710131/1‑143 (MQ=255)
                                                                                                                                       gAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATcc  >  1:619998/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: