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,5420CA97.6% 116.4 / NA 42P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (18/23);  total (18/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 = 9.95e-01

GCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGC  >  NC_000913/4184406‑4184669
                                                                                                                                        |                                                                                                                               
gcagTTAAACGCTGTATCCAACATGCAGACAACATGGCAAGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCCTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTggg                                                                                                                           <  1:696081/140‑1 (MQ=255)
    ttAAACGCCGTATCCAGCCGGGTGACAAGATGGGAGGTCGTCACGGTAACAAGGGTGTAAGTTCTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTa                                                                                                                       <  2:8343/143‑1 (MQ=255)
         cGCCGTATCCAGCCTGGTGATAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtc                                                                                                                  <  2:312782/143‑1 (MQ=255)
                  cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACCTCGTACTGAACCAGCTAGGCGAACCGACTCGTATGaa                                                                                                         >  2:37838/1‑143 (MQ=255)
                  cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTTAGCCAGCTAGGCGTACCGTCTCGTATGaa                                                                                                         >  1:193474/1‑143 (MQ=255)
                  cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                         >  2:112909/1‑143 (MQ=255)
                       tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCg                                                                                                    <  1:5420/143‑1 (MQ=255)
                        ggTGACAAGAGGGAAGGTGGTGGAGTGAACGCGGGTGTAATTTCTAAGCTCCCCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCgg                                                                                                   <  1:452676/143‑1 (MQ=255)
                          tGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                                 >  2:451852/1‑143 (MQ=255)
                          tGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAAGGTCGAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACACCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGGATGAACATCGGTc                                                                                                 >  2:87784/1‑143 (MQ=255)
                          tGACAAGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGCTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCgcac                                                                                                 >  2:201269/1‑140 (MQ=255)
                              aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                                             >  2:92507/1‑143 (MQ=255)
                                  tGGCAGGTCGTCAGGGTAACAAGGGGGTAATTTCGAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                         <  1:78116/143‑1 (MQ=255)
                                   ggCAGGTCGTCACGGTAACTGGGCTGTAATTTCTAATATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACACCGTACTGAACCCGCTGGGCGTCCCGTCTCGTATGAACATCGGTCAGATCCTCg                                                                                        <  1:702118/143‑1 (MQ=255)
                                   ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCg                                                                                        <  1:325337/143‑1 (MQ=255)
                                      aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                     >  2:302314/1‑143 (MQ=255)
                                      aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                     <  1:340489/143‑1 (MQ=255)
                                        gtcgtcACGGTACCAAGGGTGTAATTTCTAAGACCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAcc                                                                                   <  2:408240/143‑1 (MQ=255)
                                        gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGCTCCTCGAGAcc                                                                                   >  2:724347/1‑143 (MQ=255)
                                               cggtttctaGGGTTTATTTTCTAATAGCAATCCGCTCGGTGATATGCCTTATTATGGAAACGGTACGCCTCGAGACATCGTACTGAACCAGTTGGGCGTACTGTCTCGTATGAACTTCGGTCAGATCCTCGAAACCCACCTgg                                                                            <  2:23668/135‑1 (MQ=255)
                                                 atAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGt                                                                          <  2:427584/142‑1 (MQ=255)
                                                  tAACAAGGGTGTAATTTCTAGGATCAACCCGATCGATGTACCATCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTa                                                                         <  1:720000/143‑1 (MQ=255)
                                                       aGGGGGGAATTTCTAAGATCACCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCATATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCt                                                                    <  1:446871/143‑1 (MQ=255)
                                                           tgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGa                                                                >  1:654965/1‑143 (MQ=255)
                                                             gaaTGTCTAAGATCAAACCGATCGAAGAGATGCCTTACGATGAAAACGGTACGCCGGTAGACAGCGTACTGAACCAGCGGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                              <  1:455996/142‑1 (MQ=255)
                                                               aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTCTGGCTCCGACAgg                                                            >  2:683031/1‑143 (MQ=255)
                                                               aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAAAGGTCGGATGCTCGAAACCCACCTGGGTATGGCTGCGAAAgg                                                            >  1:98739/1‑143 (MQ=255)
                                                                tttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGGCAGGt                                                           >  1:213986/1‑143 (MQ=255)
                                                                  tCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAt                                                         <  1:447727/143‑1 (MQ=255)
                                                                        aTCAACCCGATCGAAGATAGGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCGCGAAACCCACCCGGGTATGGCTGCGAAGGGTACCGGCGa                                                   >  2:129011/1‑143 (MQ=255)
                                                                            aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATACTCGAAACCCACCTAGAAATGGGAGTGGACGGTATCGGCGACAAg                                               >  1:193754/1‑143 (MQ=255)
                                                                             cacctacGAAGATCTGCCCTACGATGACGTCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACACCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                              <  2:649151/137‑1 (MQ=255)
                                                                                    gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTAGTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGGTCCCCGACACCCCCCCGGGTCTCGTTCCGCCAGATATCCGCGACCAGaacaacac                                       >  1:633715/1‑141 (MQ=255)
                                                                                    gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAACCCCACCTGGGTCTGGCTGCGCCAGGTACCCGCGCCCAGGTCACCGc                                       >  1:283714/1‑143 (MQ=255)
                                                                                            gCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGa                               <  2:98777/143‑1 (MQ=255)
                                                                                                       gaagcGGTACGCCGGTAGACATCGCACGGAAGCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaa                    <  2:653568/139‑1 (MQ=255)
                                                                                                          gcggggtgCCGGGAGAAACCGGACTGAACCAGCTGGGCGGACCGGCTCGGATGAACATCGGTCAGATCCTCGAAACCCACCGGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                 <  2:661873/136‑1 (MQ=255)
                                                                                                              tACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcg             <  1:602649/143‑1 (MQ=255)
                                                                                                               aCGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc            <  2:193474/143‑1 (MQ=255)
                                                                                                                   cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaaa        <  1:112909/143‑1 (MQ=255)
                                                                                                                   cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaaa        <  1:694973/143‑1 (MQ=255)
                                                                                                                         aCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcgc  >  1:241020/1‑143 (MQ=255)
                                                                                                                                        |                                                                                                                               
GCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGC  >  NC_000913/4184406‑4184669

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