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,5430CA93.8% 86.1 / NA 32P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (16/14);  total (17/15)
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

GGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAA  >  NC_000913/4184406‑4184673
                                                                                                                                         |                                                                                                                                  
ggCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                               <  2:428158/143‑1 (MQ=255)
   ggTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTCGACATCGTACTGCACCAGCTGGGCg                                                                                                                            >  1:548952/1‑143 (MQ=255)
        aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGCACTGAACCAGCTGGGCGTACCg                                                                                                                       >  2:628811/1‑143 (MQ=255)
         aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGTTGCTGGAGCAACACCGCTAGCCAACGCACTGGACCAGCTGGGCGTACCGt                                                                                                                      >  1:487583/1‑143 (MQ=255)
                aTCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGCGCGTACCGTCTCGTAt                                                                                                               >  1:364221/1‑143 (MQ=255)
                           tGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                                    <  1:666606/143‑1 (MQ=255)
                            gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGACa                                                                                                   >  2:175306/1‑143 (MQ=255)
                                  aTGGCAGGTCGTCACGGTAACAAGGGCGTAATTTCTAAGATCAACACGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGCACCAGCTGGACGTACCGTCTCGTATGAACATCGGTACGATctt                                                                                             >  2:75398/1‑141 (MQ=255)
                                   tggtatgtcgtcACGGTAACTAGGTTGTAATGTCTAAGATCAACCCGATCGAAGATATGCCTTGCGTTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                            <  2:34872/137‑1 (MQ=255)
                                        ggTCGTCACGGTAACAAGGCTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAc                                                                                       <  2:346771/143‑1 (MQ=255)
                                          tcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAccc                                                                                     <  1:212201/143‑1 (MQ=255)
                                               aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGAATGAACATCGGACAGATCCTCGAAACCCACCTg                                                                                >  2:30601/1‑143 (MQ=255)
                                               aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGCCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTg                                                                                >  1:301431/1‑143 (MQ=255)
                                                cGGTAACAAGGCTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgg                                                                               <  1:79147/143‑1 (MQ=255)
                                                 ggTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCAggg                                                                              >  2:175917/1‑143 (MQ=255)
                                                 ggTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACACGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACATgag                                                                              >  2:310275/1‑141 (MQ=255)
                                                        aGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCt                                                                       >  1:297550/1‑143 (MQ=255)
                                                              tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTGTCGTATGAACATCGGGCAGATCCTCGAAACCCCCCTGGGAATGGCAGCGaca                                                                 >  2:457088/1‑141 (MQ=255)
                                                               aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACAGGACGCCGAAAAACAGCGTGCTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                <  1:419469/143‑1 (MQ=255)
                                                                aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAgg                                                               <  1:450985/143‑1 (MQ=255)
                                                                         aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                      <  1:595681/143‑1 (MQ=255)
                                                                            aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACac                                                   >  1:280979/1‑142 (MQ=255)
                                                                              cccGAACGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACAACGAACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                 <  1:55235/143‑1 (MQ=255)
                                                                               ccGATCGAAGATATGCCTTACGATGAAAGCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                                <  1:431399/143‑1 (MQ=255)
                                                                               ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                                <  1:541338/143‑1 (MQ=255)
                                                                                   tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGCCTCGTATGAACAGCGGTCAGATCACCGAAACCCTCCGGGGTAAGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                            >  2:231658/1‑143 (MQ=255)
                                                                                                          aaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACACGATCAACGCCATGCTGCAACAGCAGCaaga                     >  2:418739/1‑143 (MQ=255)
                                                                                                           acggagcGCCGGCAGACATCGTACTGAACCAGCTGGGCGGACCGTCTCGTATGGACATCGGTCAGATCCTCGAACCCCCCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                    <  1:377707/137‑1 (MQ=255)
                                                                                                                tcgcgcgTACACCTCGTACTGAACCAGCTCCGCGTACCGTCCCGTATGAACATCCGTCAGATCCTCGAATCCCACCTCGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc               <  1:196693/137‑1 (MQ=255)
                                                                                                                        atatgtagTACTAAGCCCGCGGGGCGTACCGTCTCGTATGAACAGCGGTCAGATCCTCGGAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgc       <  1:291140/136‑1 (MQ=255)
                                                                                                                             tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGCCAAGATCAACGCCATGCTGCACCAGCAGCAAGAAGTCGCGCGCCTGCGCGaa  >  1:13832/1‑143 (MQ=255)
                                                                                                                             tCGGACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa  >  1:388992/1‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                  
GGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAA  >  NC_000913/4184406‑4184673

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