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

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

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: