Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 541,440 A→G 10.0% Q1272R (CAG→CGG)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947541,4400AG10.0% 108.2 / 3.3 41Q1272R (CAG→CGG) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (20/16);  new base G (2/2);  total (22/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.15e-01

TGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAA  >  NC_002947/541299‑541587
                                                                                                                                             |                                                                                                                                                    
tGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTgg                                                                                                                                              <  2:420324/150‑1 (MQ=255)
  gTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGggt                                                                                                                                            >  1:735598/1‑150 (MQ=255)
   tGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGggtg                                                                                                                                           >  1:582034/1‑150 (MQ=255)
       gTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCcagca                                                                                                                                                      >  2:253868/1‑135 (MQ=255)
       gTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCcagca                                                                                                                                                      <  1:253868/135‑1 (MQ=255)
            aTGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGgcgc                                                                                                                                  >  2:427927/1‑150 (MQ=255)
                  gTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGGCGGTTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑Cg                                                                                                                           >  1:571964/1‑150 (MQ=255)
                    aCCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑Cggt                                                                                                                         <  2:546693/150‑1 (MQ=255)
                     ccGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑Cggtg                                                                                                                        <  1:36999/150‑1 (MQ=255)
                                  cGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCgg                                                                                                           <  1:281699/150‑1 (MQ=255)
                                       gTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGgaga                                                                                                      >  2:58101/1‑150 (MQ=255)
                                         ccTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATg                                                                                                    >  2:405886/1‑150 (MQ=255)
                                           tgtgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGCTGCTGggtaa                                                                                                                                          >  2:697392/1‑109 (MQ=255)
                                           tgtgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGCTGCTGggtaa                                                                                                                                          <  1:697392/111‑3 (MQ=255)
                                                      gTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCAAGCGAGTACGACAGGTTCCTACGCTCGGATGACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGgcgc                                                                                       >  1:467238/1‑150 (MQ=255)
                                                           atgtacatgCTCAAGCTGAACCACTTGGTGGCAACCAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTCGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                                  <  2:435306/150‑1 (MQ=255)
                                                           atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGACGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                                  <  1:408486/150‑1 (MQ=255)
                                                            tgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAg                                                                                 >  2:132283/1‑150 (MQ=255)
                                                                    cTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGgcg                                                                                                                                   >  1:511753/1‑93 (MQ=255)
                                                                    cTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGgcg                                                                                                                                   <  2:511753/93‑1 (MQ=255)
                                                                          cTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCa                                                                   <  1:132283/150‑1 (MQ=255)
                                                                              aCCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCCGCCATTCCG‑GATTCGATGGGTTTGAGGTGGTCAGCGTTTCGGGCAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATaca                                                               <  2:25908/150‑1 (MQ=255)
                                                                               ccACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑Cgg                                                                                                                          <  2:670626/90‑1 (MQ=255)
                                                                               ccACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑Cgg                                                                                                                          >  1:670626/1‑90 (MQ=255)
                                                                                       tggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAg                                                      >  1:162100/1‑150 (MQ=255)
                                                                                         gacgacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGaa                                                    >  1:506667/1‑150 (MQ=255)
                                                                                                 gATGCACGCGCGTTCCAATGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTg                                                                                                                                               >  1:389956/1‑52 (MQ=255)
                                                                                                 gATGCACGCGCGTTCCAATGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTg                                                                                                                                               <  2:389956/52‑1 (MQ=255)
                                                                                                             ttCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGATATGGAAGTGTGGGCGCTGGAAGCa                                                                               >  1:167108/1‑103 (MQ=255)
                                                                                                             ttCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCa                                                                               <  2:167108/103‑1 (MQ=255)
                                                                                                                cACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacga                             >  1:504658/1‑150 (MQ=255)
                                                                                                                cACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacga                             >  1:223681/1‑150 (MQ=255)
                                                                                                                 aCTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacgac                            <  2:735598/150‑1 (MQ=255)
                                                                                                                      ttCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGaa                       <  2:467238/150‑1 (MQ=255)
                                                                                                                       tCCTACAGCCTGGTTACGCAGCATCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTAAAc                      <  2:621035/150‑1 (MQ=255)
                                                                                                                           aCAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGcc                  >  2:517015/1‑150 (MQ=255)
                                                                                                                             aGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCG‑TTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAATTCGGACGACGTGAACGGCCGTa               <  1:492152/150‑1 (MQ=255)
                                                                                                                                cTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTAcc             >  2:113796/1‑150 (MQ=255)
                                                                                                                                   gTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAg          <  1:574729/150‑1 (MQ=255)
                                                                                                                                   atctccgagcagcCGCTTGGTTGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAg          <  2:232114/143‑1 (MQ=255)
                                                                                                                                        ccAGCGGCTGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGaa                                      <  2:223628/117‑1 (MQ=255)
                                                                                                                                        ccAGCGGCTGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGaa                                      >  1:223628/1‑117 (MQ=255)
                                                                                                                                        ccAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGaagat     >  1:232114/1‑146 (MQ=255)
                                                                                                                                        ccAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTa     >  2:180006/1‑150 (MQ=255)
                                                                                                                                         cagcagCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTAc    <  2:223681/150‑1 (MQ=255)
                                                                                                                                           gcagcCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACCGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACaa  <  1:427927/150‑1 (MQ=255)
                                                                                                                                             agcCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGaa                                      >  2:227238/1‑112 (MQ=255)
                                                                                                                                             agcCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGaa                                      <  1:227238/112‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
TGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAA  >  NC_002947/541299‑541587

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