Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,182,877 T→G F545V (TTC→GTC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,182,8770TG97.6% 146.9 / NA 42F545V (TTC→GTC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (21/20);  minor base A (0/1);  total (21/21)
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

AAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAA  >  NC_000913/4182730‑4183022
                                                                                                                                                   |                                                                                                                                                 
aaGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCg                                                                                                                                                <  2:581744/151‑1 (MQ=255)
     gATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGtt                                                                                                                                           >  1:64754/1‑151 (MQ=255)
               gcagTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTac                                                                                                                                 >  2:873781/1‑151 (MQ=255)
                     aaaGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCg                                                                                                                            >  2:438375/1‑150 (MQ=255)
                     aaaGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGa                                                                                                                           <  1:86552/151‑1 (MQ=255)
                      aaGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGAc                                                                                                                          >  2:36213/1‑151 (MQ=255)
                        gagTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCCCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTc                                                                                                                        <  1:781820/151‑1 (MQ=255)
                        gagTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTACGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGCCTGACCCGTGAACGCGCAGGCGTAGAAGTTCGAGACGGACACCCGACTc                                                                                                                        >  2:778503/1‑151 (MQ=255)
                            tcttcGATTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACt                                                                                                                     >  1:782303/1‑150 (MQ=255)
                             cttcGGTTCCAGCCTGCTGTCTCAGTTTATGGACCAGAACACCCCGCTGTCTGAGATTACGCACAACCGGCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCATCGAAGTGCGAGACGTACACCCGACTCACTa                                                                                                                    <  2:134313/150‑1 (MQ=255)
                                 ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTc                                                                                                               <  1:621688/151‑1 (MQ=255)
                                 ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTc                                                                                                               >  2:700099/1‑151 (MQ=255)
                                  gTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcg                                                                                                              >  1:291872/1‑151 (MQ=255)
                                        gccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATg                                                                                                        >  2:710755/1‑151 (MQ=255)
                                              tGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAAt                                                                                                  >  1:741516/1‑151 (MQ=255)
                                                tctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCg                                                                                                >  2:647593/1‑151 (MQ=255)
                                                   cAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaa                                                                                              <  1:338060/150‑1 (MQ=255)
                                                                aGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTcc                                                                                <  2:71163/151‑1 (MQ=255)
                                                                  aacaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTTCAATCGAAACCCCTGAAGGTCCg                                                                               >  1:258713/1‑150 (MQ=255)
                                                                                 gagaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAt                                                                 <  1:153457/149‑1 (MQ=255)
                                                                                          cacaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGACGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTg                                                       <  2:457953/150‑1 (MQ=255)
                                                                                            caAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTcc                                                    >  2:767991/1‑151 (MQ=255)
                                                                                                cgtcgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtgt                                                <  1:673991/151‑1 (MQ=255)
                                                                                                   cgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTAc                                              <  1:750514/150‑1 (MQ=255)
                                                                                                    gtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGc                                            >  2:725836/1‑151 (MQ=255)
                                                                                                      aTCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                          >  2:610983/1‑151 (MQ=255)
                                                                                                          ccGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                          <  2:37188/147‑1 (MQ=255)
                                                                                                             cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTaa                                   <  1:621982/151‑1 (MQ=255)
                                                                                                                  ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                              >  2:280047/1‑151 (MQ=255)
                                                                                                                    cccAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                            >  2:202039/1‑151 (MQ=255)
                                                                                                                     ccAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                           <  2:620672/151‑1 (MQ=255)
                                                                                                                         gCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt                       <  2:879608/151‑1 (MQ=255)
                                                                                                                           ggTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTcc                     <  1:923205/151‑1 (MQ=255)
                                                                                                                            gTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                    >  1:380711/1‑151 (MQ=255)
                                                                                                                                gACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTga                 <  2:581454/150‑1 (MQ=255)
                                                                                                                                  cccGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGAc              <  1:562136/151‑1 (MQ=255)
                                                                                                                                    cGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTc            >  1:248376/1‑151 (MQ=255)
                                                                                                                                       gAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCg          >  2:887796/1‑150 (MQ=255)
                                                                                                                                             gCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGt    <  1:108869/150‑1 (MQ=255)
                                                                                                                                             gCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa   <  2:689817/151‑1 (MQ=255)
                                                                                                                                             gCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa   <  1:273299/151‑1 (MQ=255)
                                                                                                                                               aGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaa  >  2:678830/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAA  >  NC_000913/4182730‑4183022

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