Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 541,443 C→T 8.4% P1273L (CCG→CTG)  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,4430CT 8.4% 126.7 / 2.7 48P1273L (CCG→CTG) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (24/19);  new base T (2/2);  total (27/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 = 9.06e-01

GTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAAC  >  NC_002947/541301‑541591
                                                                                                                                              |                                                                                                                                                    
gTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGggt                                                                                                                                               >  1:678899/1‑150 (MQ=255)
      ttCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAg                                                                                                                                         >  1:280905/1‑150 (MQ=255)
      ttCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCAAGCAGCCGCTGGGTGGTAAg                                                                                                                                         <  2:572874/150‑1 (MQ=255)
         gATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGctgctg                                                                                                                                                  >  2:452713/1‑138 (MQ=255)
         gATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGctgctg                                                                                                                                                  <  1:452713/138‑1 (MQ=255)
                              ttCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGCTGCTGggtaa                                                                                                                                             <  1:128826/122‑3 (MQ=255)
                              ttCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGCTGCTGggtaa                                                                                                                                             >  2:128826/1‑120 (MQ=255)
                                cGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgg                                                                                                               <  1:10233/150‑1 (MQ=255)
                                      tCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGAAAAAATGCACGCGCGTTCGACTGGTTCCTACAGCCGGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGAt                                                                                                         >  1:297119/1‑150 (MQ=255)
                                              ctgTTGGTTACATGTACATGCTCAAGCTGAACCACTTGAAGTACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAgtgt                                                                                                 <  1:407229/148‑1 (MQ=255)
                                                 ttGGTTACATGTACATGATCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCATGGTTACCCAGCAGCAGACGGGTAATCAGTAGCATGTCGTTGGACAGCGGTTCGTGGATATGGAAGtgtgag                                                                                              >  2:407229/1‑148 (MQ=255)
                                                    gTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTggg                                                                                              <  2:743196/147‑1 (MQ=255)
                                                    gTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTggg                                                                                              >  1:743196/1‑147 (MQ=255)
                                                    gTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGgcgc                                                                                           >  1:471550/1‑150 (MQ=255)
                                                             acatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCAt                                                                                  >  2:759200/1‑150 (MQ=255)
                                                                  cTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGgc                                                                             >  1:577245/1‑150 (MQ=255)
                                                                   tCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGgcg                                                                            >  2:447854/1‑150 (MQ=255)
                                                                      aGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgg                                                                                                               >  1:455821/1‑112 (MQ=255)
                                                                      aGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgg                                                                                                               <  2:455821/112‑1 (MQ=255)
                                                                           aaCCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgg                                                                                                               >  1:466388/1‑107 (MQ=255)
                                                                           aaCCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgg                                                                                                               <  2:466388/107‑1 (MQ=255)
                                                                            aCCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                                              >  2:513892/1‑107 (MQ=255)
                                                                            aCCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                                              <  1:513892/107‑1 (MQ=255)
                                                                                 ttGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTg                                                              <  2:297119/150‑1 (MQ=255)
                                                                                     tggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATg                                                                                                        >  1:780340/1‑104 (MQ=255)
                                                                                     tggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATg                                                                                                        <  2:780340/104‑1 (MQ=255)
                                                                                          gacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCAGAGGCGGTCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATg                                                     >  1:251418/1‑150 (MQ=255)
                                                                                            cAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCt                                                   <  2:577245/150‑1 (MQ=255)
                                                                                             aaGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCa                                                                       <  1:177926/129‑1 (MQ=255)
                                                                                             aaGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCa                                                                       >  2:177926/1‑129 (MQ=255)
                                                                                                      gcgcgTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTgg                                                                                               <  1:335777/96‑1 (MQ=255)
                                                                                                      gcgcgTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTgg                                                                                               >  2:335777/1‑96 (MQ=255)
                                                                                                                cTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAgtgt                                                                                                 <  1:650987/84‑1 (MQ=255)
                                                                                                                cTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAgtgt                                                                                                 >  2:650987/1‑84 (MQ=255)
                                                                                                                   gTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGTATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGa                            <  2:471550/150‑1 (MQ=255)
                                                                                                                         aCAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTc                                       >  2:137669/1‑133 (MQ=255)
                                                                                                                         aCAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTc                                       <  1:137669/133‑1 (MQ=255)
                                                                                                                          cAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGTAGATGGAAGTGTGGTCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCg                     >  1:297504/1‑150 (MQ=255)
                                                                                                                             ccTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgac                                   <  2:279197/133‑1 (MQ=255)
                                                                                                                             ccTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgac                                   >  1:279197/1‑133 (MQ=255)
                                                                                                                             ccTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATAAAGATGAACATACACCCTGCAAGACATGTTCACTGTGAAGTCGGACGACGTGAACGGCCGTAc                  >  1:607173/1‑150 (MQ=255)
                                                                                                                               tGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATgg                                                                                                       <  2:422850/63‑1 (MQ=255)
                                                                                                                               tGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATgg                                                                                                       >  1:422850/1‑63 (MQ=255)
                                                                                                                               tGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCa                >  1:296090/1‑150 (MQ=255)
                                                                                                                                ggTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCaa               >  2:339714/1‑150 (MQ=255)
                                                                                                                                    aCCCAGCGGCTGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATg                                                     <  1:116683/108‑1 (MQ=255)
                                                                                                                                    aCCCAGCGGCTGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATg                                                     >  2:116683/1‑108 (MQ=255)
                                                                                                                                      ccAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTa         <  2:296090/150‑1 (MQ=255)
                                                                                                                                      ccAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTa         >  1:78269/1‑150 (MQ=255)
                                                                                                                                         gcagcCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACaa      <  2:78269/150‑1 (MQ=255)
                                                                                                                                             ccGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGAAGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAAc  >  1:382407/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAAC  >  NC_002947/541301‑541591

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