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.7% 127.8 / NA 44P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (25/18);  minor base G (0/1);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAT  >  NC_000913/4184409‑4184677
                                                                                                                                     |                                                                                                                                       
gTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                >  1:789766/1‑143 (MQ=255)
  tAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAc                                                                                                                              <  2:132109/143‑1 (MQ=255)
  tAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAc                                                                                                                              <  1:401003/143‑1 (MQ=255)
         cGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCg                                                                                                                       >  2:352530/1‑143 (MQ=255)
          gTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAAAAAGCTGGAAGTAACGTCTCGt                                                                                                                      >  2:178760/1‑143 (MQ=255)
           tATCCAGCTGTGTGACAAGATGGCAGGTCGGCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTa                                                                                                                     <  1:165913/143‑1 (MQ=255)
                     ggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCgg                                                                                                           >  1:230536/1‑143 (MQ=255)
                       tGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                                         >  2:665773/1‑143 (MQ=255)
                       tGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAGGAACATCGGTc                                                                                                         >  1:307088/1‑143 (MQ=255)
                           aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGCTATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                                                     >  1:461059/1‑143 (MQ=255)
                                         tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcc                                                                                       >  2:13470/1‑143 (MQ=255)
                                          cACGGTACCAAGGGTGGAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCt                                                                                      <  2:533265/143‑1 (MQ=255)
                                          cACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAGCt                                                                                      >  2:200981/1‑143 (MQ=255)
                                                aaCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTAt                                                                                <  2:852927/143‑1 (MQ=255)
                                                aaCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTAt                                                                                <  2:496306/143‑1 (MQ=255)
                                                   aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGc                                                                             >  2:804134/1‑143 (MQ=255)
                                                   aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGc                                                                             >  2:802084/1‑143 (MQ=255)
                                                            aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGCAAACGGTACACCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCTTATGAACATCGGTCACCAAAACGAAATCCACCTTGCAATGGCTGCGAAAgg                                                                    >  2:667541/1‑143 (MQ=255)
                                                              ttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAAGTa                                                                  >  1:45120/1‑143 (MQ=255)
                                                                cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTAGGGCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGGATGAACATCGGTCAGAGCCTCTAAACCCACGTGGGTGTGGCGGCGAGAGGTAtg                                                                >  2:232446/1‑142 (MQ=255)
                                                                  aaGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCgg                                                              <  1:480009/143‑1 (MQ=255)
                                                                    gATCAACCCGATCGAAGATATGCCTTACGATGAAAACGATCCGACAGCCCGCATCTTACTGAACCAGCCAGGCGTACCGTCTCGTCTGAACATCGGTCAGACACTCGAAACCCACCTGCGGATCGCTGCGACGGGTATCgcag                                                            >  1:385974/1‑140 (MQ=255)
                                                                          cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                      >  2:682384/1‑143 (MQ=255)
                                                                             gATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGATACTCACCCGGGTATGGCTGCGGCATGTCTCGGCGACAAGATCa                                                   >  2:274839/1‑143 (MQ=255)
                                                                              aTCGAAGATATGCCTTACGATGAAAACGGTACGCGTATAATCAGCGTACTGAACCAGCTGGGCGTACCGTCTAGTATGAACATCGGTCAGATCCTGGAAAGCAAACTGGGTATGGCTGCGAAAGGTAGCGGCGACAAGATCaa                                                  >  2:165393/1‑143 (MQ=255)
                                                                                    gATATGCCTTACGATGGAACCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCAt                                            <  1:402505/143‑1 (MQ=255)
                                                                                     atatGCCTTCCGATGACAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATg                                           <  1:675855/143‑1 (MQ=255)
                                                                                      tatGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGc                                          <  2:310200/143‑1 (MQ=255)
                                                                                           cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTCTGGCTGCGAAAGGTATCGGCGACAAGAGCCGCGCCCTGCTGaca                                     >  2:520034/1‑141 (MQ=255)
                                                                                           cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                     <  1:281930/143‑1 (MQ=255)
                                                                                            ttccGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAc                                    <  2:552232/140‑1 (MQ=255)
                                                                                              aCGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag                                  <  2:605790/143‑1 (MQ=255)
                                                                                                   gAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGGGAAAGGTCTCGGCGAGAAGATCACCGGCATGCTGGAAcagcgcgc                             >  2:242559/1‑139 (MQ=255)
                                                                                                          gTACGCCGGTAGACATCGTACTGAGCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTc                      <  1:594721/143‑1 (MQ=255)
                                                                                                            aCGCCAGTCGCCATCGTATCGAACCAGCTGGGCGTACCGTCGCGTATGTACATCGGTCAGATCCGCGAAACCTACCTGGGTAGGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc                    <  1:778460/143‑1 (MQ=255)
                                                                                                                cGGTAGACATCGTACTGAAACAGCTGGGCGTACCGTGTCTTATGACCATCGGTCAGTTCCTCGCACACCACCTGGGTATGCATGCGAAAGGAACCGGCGACACGGGCTACGCCATGCTGAAACAGCCGCaagaacgagcgaga                >  1:693543/1‑134 (MQ=255)
                                                                                                                cGGTAGACATCGTACTCAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaaa                <  2:376782/143‑1 (MQ=255)
                                                                                                                   tAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACGCACCTGGGTATGGCCGCGAACGGTATCGGCCACAAGATAACAGCCATGCTGAAACAGCAGCACGAAGTCGCGAAACTg             >  1:269234/1‑143 (MQ=255)
                                                                                                                   tAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAGCAACGCCATGCTGAAACAGCAGCAAGGAGTCGCGAAACTg             >  2:489262/1‑143 (MQ=255)
                                                                                                                   tAGACAGCGTACTGAACCAGCTGGGCGGACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAGCGGCGACAAGAGCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTg             <  2:795905/143‑1 (MQ=255)
                                                                                                                         tCGTACTGAACCGGCGCGGCGTACCGTCTCGTATGCACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa       <  1:352530/143‑1 (MQ=255)
                                                                                                                         tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa       >  2:682804/1‑143 (MQ=255)
                                                                                                                             aCTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCa   <  1:749033/143‑1 (MQ=255)
                                                                                                                              cTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAACAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGCAGTCag  >  1:111407/1‑142 (MQ=255)
                                                                                                                                     |                                                                                                                                       
GTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAT  >  NC_000913/4184409‑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: