Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,543 C→A P1100Q (CCG→CAG)  rpoB → RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5430CA97.7% 129.4 / NA 44P1100Q (CCG→CAG) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (24/19);  total (25/19)
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.96e-01

TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184405‑4184681
                                                                                                                                          |                                                                                                                                          
tGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTg                                                                                                                                        <  1:439944/143‑1 (MQ=255)
     gTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACTCACTGGGCGt                                                                                                                                   >  2:327210/1‑143 (MQ=255)
     gTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                   >  1:567541/1‑143 (MQ=255)
                tagccAGCCTGGTGACACGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGCCCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTa                                                                                                                        <  2:472281/140‑1 (MQ=255)
                tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGCGCCGTCTCGTa                                                                                                                        >  1:544313/1‑143 (MQ=255)
                     aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAc                                                                                                                   >  1:890155/1‑143 (MQ=255)
                             gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCACCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGGACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCa                                                                                                           <  2:496931/143‑1 (MQ=255)
                             gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCa                                                                                                           >  1:514348/1‑143 (MQ=255)
                                aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCATATGAAAAACGGTCAGAt                                                                                                        >  1:56966/1‑143 (MQ=255)
                                        aGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                <  2:3824/143‑1 (MQ=255)
                                           agagcACGGTAACACGGGTGGAATGTAGAGGATCAACCCGAGCGAAGTGAGGCC‑GAGGATGAAAACGGTACGCCGGTAGCCACCGTACTGGACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                                                            <  1:453970/139‑1 (MQ=255)
                                                 cGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgg                                                                                       <  1:442233/143‑1 (MQ=255)
                                                       cAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATgg                                                                                 <  1:155224/143‑1 (MQ=255)
                                                          gggTGTAATTTCTAAGATCAACCCGATCCAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTAGTATGAGCATCGGTCAGATCCTCGAAACCCCCCTGGGTATGGCTg                                                                              >  1:132086/1‑143 (MQ=255)
                                                           ggTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGTTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGc                                                                             <  2:464264/143‑1 (MQ=255)
                                                           ggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGAGTACCGTCTCGTATGAACATCGGACAGAACCTCGAAACCCACCCAGGTATGGCCGc                                                                             >  2:302132/1‑143 (MQ=255)
                                                                agtttCTAATATCACCCCGATCGAAGATATGTCTTACGATGACAACGGTACGCCGGTAGACATCGTACTGATCCAGCTGGGCGTACCGTGTCGTATGAACATCGGTCAGATCCTCGAATCCCACCTGTGTATGGCTGCGAAAg                                                                        <  2:301427/141‑1 (MQ=255)
                                                                 aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAgg                                                                       >  1:348913/1‑143 (MQ=255)
                                                                   ttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCCAAACCCACCTGGCTATGGCTGCGAGAGCTa                                                                     >  2:707502/1‑143 (MQ=255)
                                                                      tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATACTCGAAACCCACCCGGGTATGGCGGCGgaaacaatcg                                                                  >  2:256559/1‑136 (MQ=255)
                                                                      tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCATACTGAACCAGCTGGGCGTACCGTCTCGCATCAACATCGGTCAGATCCTCGAAACCCACCCGGGTATCACTGCCAAAGGTATCg                                                                  >  2:240894/1‑143 (MQ=255)
                                                                          aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                              <  1:732945/143‑1 (MQ=255)
                                                                               cccGCTCGAAGATATGCCTTACGACGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                         <  2:940789/143‑1 (MQ=255)
                                                                               cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGACCATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGCAACGGCGACAAGa                                                         >  2:243731/1‑143 (MQ=255)
                                                                                ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                                        >  1:807594/1‑143 (MQ=255)
                                                                                ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                                        >  1:253140/1‑143 (MQ=255)
                                                                                      gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGc                                                  >  2:507188/1‑143 (MQ=255)
                                                                                      gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGc                                                  <  2:514348/143‑1 (MQ=255)
                                                                                       aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                                 >  1:287277/1‑143 (MQ=255)
                                                                                       aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCGCGAGACCCACCTGGGTATGGCTGCGAAAGGTAACAGCGACAACACCAGCGcc                                                 >  2:654944/1‑143 (MQ=255)
                                                                                              gCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCGGCGAAAGGTATCGGCGACAAGATCAACGCCAGGCTGa                                          >  2:526739/1‑143 (MQ=255)
                                                                                               ccTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGACAGGTATCGGCGACAAGATCAACGCCATGCTGaa                                         >  2:410942/1‑143 (MQ=255)
                                                                                                   aCGATGAAAACGGTACGCCGGTAGACATAGTACTGAACCAGCTGGGCGTACCGTCGCGTATGAACATCGGACAGATCCTCGAAAAACACCTGGGAATGGACACGAAAGGTATCGCCAACAAGTACAACGCTATACTGAAAcag                                     >  2:526094/1‑143 (MQ=255)
                                                                                                    cGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCTGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATGAACGCCGTGCAGAAAcagc                                    >  2:631198/1‑143 (MQ=255)
                                                                                                            aCGGTACGCCGGTCGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                            <  2:465082/143‑1 (MQ=255)
                                                                                                                 aCGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc                       <  2:567541/143‑1 (MQ=255)
                                                                                                                 aCGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc                       <  1:345751/143‑1 (MQ=255)
                                                                                                                 aCGCCGGCAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc                       <  1:886489/143‑1 (MQ=255)
                                                                                                                    ccGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaa                    <  1:196299/143‑1 (MQ=255)
                                                                                                                               cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAACGGTATCCGCGACAAGATCAACGCCCTGATGACACCGCAGCAAGAAGACGCGACACTGCGCGacg         >  2:325987/1‑141 (MQ=255)
                                                                                                                               cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGCAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAt         >  1:578974/1‑143 (MQ=255)
                                                                                                                                gTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCGGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAtt        <  1:555368/143‑1 (MQ=255)
                                                                                                                                     gAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGCACTCAccc   >  2:454618/1‑143 (MQ=255)
                                                                                                                                      aaCCAGCTGGGCGGACCGTCTAGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCa  <  2:56966/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184405‑4184681

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