Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,036 T→C V931A (GTT→GCT)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,0360TC97.4% 128.6 / ‑4.0 38V931A (GTT→GCT) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (13/24);  total (13/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

ATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGT  >  NC_000913/4183888‑4184183
                                                                                                                                                    |                                                                                                                                                   
aTTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTc                                                                                                                                                   >  2:452121/1‑151 (MQ=255)
      gttggtAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCt                                                                                                                                             >  1:523273/1‑151 (MQ=255)
            aaggtaaCGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTc                                                                                                                                       <  1:136739/151‑1 (MQ=255)
             aggtaaCGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTTCTcg                                                                                                                                      <  1:596390/151‑1 (MQ=255)
                 aaCGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGAt                                                                                                                                  <  2:221141/151‑1 (MQ=255)
                    gCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGc                                                                                                                               <  1:802609/151‑1 (MQ=255)
                    gCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGc                                                                                                                               >  2:688841/1‑151 (MQ=255)
                       gAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTa                                                                                                                            <  1:111217/151‑1 (MQ=255)
                        aaaGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTa                                                                                                                            <  2:457792/150‑1 (MQ=255)
                          aGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGaa                                                                                                                         <  1:758870/151‑1 (MQ=255)
                             tGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGaaaa                                                                                                                       >  1:366603/1‑150 (MQ=255)
                               aaaCTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGaaaa                                                                                                                       >  2:533041/1‑148 (MQ=255)
                                      gCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGt                                                                                                             >  2:194405/1‑151 (MQ=255)
                                         gACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTgcg                                                                                                          <  1:720389/151‑1 (MQ=255)
                                         gACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTg                                                                                                            <  1:553352/149‑1 (MQ=255)
                                                           gctgcGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGaaa                                                                                          <  1:736931/149‑1 (MQ=255)
                                                              gcgtgcgATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCa                                                                                      <  1:500521/150‑1 (MQ=255)
                                                               cgtgcgATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCACGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGc                                                                                    >  2:508218/1‑151 (MQ=255)
                                                                gtgcgATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGc                                                                                    >  1:256670/1‑150 (MQ=255)
                                                                 tgcgATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTc                                                                                  >  2:391897/1‑151 (MQ=255)
                                                                         tCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAAc                                                                              <  2:691470/147‑1 (MQ=255)
                                                                         tCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACa                                                                             <  2:505808/148‑1 (MQ=255)
                                                                         tCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAg                                                                            <  2:479151/149‑1 (MQ=255)
                                                                         tCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGc                                                                          >  1:85350/1‑151 (MQ=255)
                                                                          cGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCg                                                                         >  1:182805/1‑151 (MQ=255)
                                                                              gagaAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCgaaga                                                                     <  1:297909/151‑1 (MQ=255)
                                                                               agaAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCgaagaa                                                                    <  1:194382/151‑1 (MQ=255)
                                                                                          gACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTc                                                           <  1:208072/149‑1 (MQ=255)
                                                                                                 aaGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACt                                                  <  2:366626/151‑1 (MQ=255)
                                                                                                  aGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTg                                                 <  1:214825/151‑1 (MQ=255)
                                                                                                         tctcTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGa                                             <  2:370945/148‑1 (MQ=255)
                                                                                                              gcgcGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCACACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTc                                        <  1:139012/148‑1 (MQ=255)
                                                                                                              gcgcGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCg                                       <  2:338632/149‑1 (MQ=255)
                                                                                                              gcgcGTACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCg                                       <  2:528538/149‑1 (MQ=255)
                                                                                                                   tACCAAACGGTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCggg                                >  2:363117/1‑151 (MQ=255)
                                                                                                                           ggTGTATCCGGTACGGTTATCGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACGGCAGATCCTCGAAGCGGGTCTGTTCa                        <  1:583026/151‑1 (MQ=255)
                                                                                                                                               cGACGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGtgctgtgctg    <  2:772954/151‑1 (MQ=255)
                                                                                                                                                 aCGCTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGt  >  1:714145/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTTCGGTGAGAAAGCCTCTGACGTTAAAGACTCTTCTCTGCGCGTACCAAACGGTGTATCCGGTACGGTTATCGACGTTCAGGTCTTTACTCGCGATGGCGTAGAAAAAGACAAACGTGCGCTGGAAATCGAAGAAATGCAGCTCAAACAGGCGAAGAAAGACCTGTCTGAAGAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGT  >  NC_000913/4183888‑4184183

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑