Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,183,814 T→A V857E (GTG→GAG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,183,8140TA96.4% 103.3 / NA 28V857E (GTG→GAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (16/11);  minor base C (1/0);  total (17/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTT  >  NC_000913/4183666‑4183961
                                                                                                                                                    |                                                                                                                                                   
aaTGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAgg                                                                                                                                                   <  1:345803/151‑1 (MQ=255)
 aTGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAggg                                                                                                                                                  >  2:87756/1‑151 (MQ=255)
         aaCTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTg                                                                                                                                          <  2:38398/151‑1 (MQ=255)
               gAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGctct                                                                                                                                    >  2:503102/1‑151 (MQ=255)
                     tCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAAc                                                                                                                              >  2:295438/1‑151 (MQ=255)
                              gTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAAt                                                                                                                     >  2:208090/1‑151 (MQ=255)
                                   cGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGt                                                                                                                <  1:268249/151‑1 (MQ=255)
                                    gAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGt                                                                                                                >  2:470022/1‑150 (MQ=255)
                                       cGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCg                                                                                                            >  1:273245/1‑151 (MQ=255)
                                               tCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACAtt                                                                                                    >  2:184747/1‑151 (MQ=255)
                                                            ttCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTg                                                                                        >  2:361383/1‑150 (MQ=255)
                                                             tCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCGAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTg                                                                                        <  1:12692/149‑1 (MQ=255)
                                                             tCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGAc                                                                                      >  1:389784/1‑151 (MQ=255)
                                                               accaccATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg                                                                                    >  1:17765/1‑151 (MQ=255)
                                                                   ccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTg                                                                                 <  1:488646/150‑1 (MQ=255)
                                                                                 gAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg                                                                  >  2:261792/1‑151 (MQ=255)
                                                                                      ggCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggta                                                               <  1:118647/149‑1 (MQ=255)
                                                                                                     tGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGt                                              >  2:13651/1‑151 (MQ=255)
                                                                                                       acacCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTg                                             >  1:471740/1‑150 (MQ=255)
                                                                                                             aGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCa                                      <  1:23661/151‑1 (MQ=255)
                                                                                                                tGGGGCCGGAAGAGATTACCACTGGCATCACGAACGCGGGAGAAGCTGGGCTCTCCAAACTGGGTGAATCCGGAATCGGTTACATTGGTGCGGAAGTGACCCGAGGCGACATGGTGGGGGGTAAGGTAACGCCGAAAGGTGATGCTCAGCt                                   >  1:330415/1‑151 (MQ=255)
                                                                                                                           gagaTCCCCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCagaag                        <  1:437142/151‑1 (MQ=255)
                                                                                                                                 aCCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAAc                  <  2:197311/151‑1 (MQ=255)
                                                                                                                                    gCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAActg                >  1:364794/1‑150 (MQ=255)
                                                                                                                                     cTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAActgct              >  1:363645/1‑151 (MQ=255)
                                                                                                                                      tGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAActgctg             >  1:514006/1‑151 (MQ=255)
                                                                                                                                       gACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACtgctgc            <  2:79469/151‑1 (MQ=255)
                                                                                                                                                 aCGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCtt  <  1:295354/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTT  >  NC_000913/4183666‑4183961

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