Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,182,820 C→T H526Y (CAC→TAC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,182,8200CT97.9% 146.8 / NA 47H526Y (CAC→TAC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/21);  minor base A (1/0);  total (26/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

CTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAT  >  NC_000913_3_bme_pgi/4182679‑4182959
                                                                                                                                              |                                                                                                                                            
ctgtctCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACACCCCGCTGTCTAAGATTACGTACAAAcg                                                                                                                                       >  2:183682/1‑149 (MQ=255)
     tctGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGACCAAACGTCGTa                                                                                                                                  >  1:7551/1‑149 (MQ=255)
         ggCGATCTGGATACGCTGATGCCACAGTTTCTGGTCACGGCCCAGCCGATTTCCG‑CAGCAGTGTAAGACATCTTCGGTGCAAGCCAGCGGGCTGAGTTTGTGTACCAGAACAACCCGCTGTCTGAGTTTACGTACAAACGTCGTAtatc                                                                                                                              >  2:343060/1‑149 (MQ=255)
            gATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGc                                                                                                                           >  1:612303/1‑149 (MQ=255)
            gATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGAATCTCCGc                                                                                                                           >  1:329435/1‑149 (MQ=255)
              tCTGGATACCCTGAGGCCACAGGAGATGATCAACGCCAAGCCGATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCAc                                                                                                                         >  1:552699/1‑149 (MQ=255)
                    tACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGcc                                                                                                                   >  2:89003/1‑149 (MQ=255)
                       ccTGATGCCACAGGATATGATCAACGCCAGGAC‑ATTTCCGCCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCCCACTCGGCC‑CAg                                                                                                               >  1:474976/1‑149 (MQ=255)
                        cTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGGGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAgg                                                                                                              >  2:343130/1‑149 (MQ=255)
                        cTGATGCCACAGGATATGATCAACGCCAAGCAAATTTCCT‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTTTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGACATTACGTACAAACGTCGTATCTACGCACTCGGCC‑CAgg                                                                                                              >  2:300856/1‑149 (MQ=255)
                             gCCACAGGATATGATCAACGCCAAGCCGATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTc                                                                                                         <  1:89003/149‑1 (MQ=255)
                               cacaGGATATGATCAACGCCAAGCCGATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTg                                                                                                       >  2:308527/1‑149 (MQ=255)
                               cacaGGATATGATCAACGCCAAGCCGAGTTCCG‑CAGCCGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTCTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGCCC‑CAGGCgcgctg                                                                                                       >  2:329716/1‑144 (MQ=255)
                                      tatGATCAACGCCAAGCCGATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTg                                                                                                >  1:130582/1‑149 (MQ=255)
                                                        gATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCTCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAg                                                                              <  1:506188/149‑1 (MQ=255)
                                                           ttCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTc                                                                           <  1:29233/149‑1 (MQ=255)
                                                           ttCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTc                                                                           <  2:130582/149‑1 (MQ=255)
                                                             ccG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATCACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCga                                                                         <  1:89306/149‑1 (MQ=255)
                                                               g‑cagcaGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGGTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCgaga                                                                       <  1:567691/149‑1 (MQ=255)
                                                                        gAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAAAGCGCAGGCCGCGAAGTTCGAGCCGTACAcc                                                               >  2:265214/1‑149 (MQ=255)
                                                                         aaaGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAccc                                                              >  1:192271/1‑149 (MQ=255)
                                                                             agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACt                                                          >  1:509517/1‑149 (MQ=255)
                                                                               ttcttgGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGCACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCa                                                        <  1:59280/149‑1 (MQ=255)
                                                                               ttcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGaa                                                                               >  1:27844/1‑126 (MQ=255)
                                                                               ttcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGaa                                                                               <  2:27844/126‑1 (MQ=255)
                                                                                 cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTCCACCCGACTCACt                                                      >  2:449342/1‑149 (MQ=255)
                                                                                  ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTa                                                     <  1:324830/149‑1 (MQ=255)
                                                                                         ccagccagCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgc                                              <  1:5763/149‑1 (MQ=255)
                                                                                          tagcgagcTGTCTGAGTTTATGGACCAGAACAACCCGCTTTCTGAGATTACGTATAAACGTCGTATCTCCGCACTCGGCC‑CTGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgcg                                             <  2:570029/148‑1 (MQ=255)
                                                                                           agccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGt                                            >  1:323147/1‑149 (MQ=255)
                                                                                                gcTGTCTCACTTTACGGCCCAGAACTACCCGTTCTCTGTGATTTCGTACCCCCGTCGTATCTCCGC‑CCCGGCCGGAGGCAGTCTTGGATGGCGACTTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTc                                       <  1:627020/149‑1 (MQ=255)
                                                                                                   gTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCaa                                    <  2:323147/149‑1 (MQ=255)
                                                                                                     ctcAGTTTATGGAACAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGCGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATc                                  <  1:183682/149‑1 (MQ=255)
                                                                                                                gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGCCTCACTACGGTCGCGTATGTCCACTCGAAAccccggg                       >  1:12239/1‑146 (MQ=255)
                                                                                                                 aCCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGaa                      >  1:559194/1‑149 (MQ=255)
                                                                                                                  ccAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTAGCGTATGTCCAATCGAACCTCCTTAAg                     >  2:67437/1‑149 (MQ=255)
                                                                                                                    aGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGt                   <  1:253010/149‑1 (MQ=255)
                                                                                                                     gAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTGCCGTATGTCCAATCGAAACCCCTGAAGGTc                  >  1:250295/1‑149 (MQ=255)
                                                                                                                     gAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGGAGGTc                  >  2:385455/1‑149 (MQ=255)
                                                                                                                      accaacCCGCTGTGTGGGATTACGTACAAACGTCGGATGCCCGCACTCGGCC‑AAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGAGTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTcc                 <  2:107314/147‑1 (MQ=255)
                                                                                                                       acaacCCGCTGGCTGAGATTACGTACACACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCg                <  2:393161/149‑1 (MQ=255)
                                                                                                                            ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGAc                                                           >  2:318534/1‑101 (MQ=255)
                                                                                                                            ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGAc                                                           <  1:318534/101‑1 (MQ=255)
                                                                                                                             cGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATc          >  2:350554/1‑149 (MQ=255)
                                                                                                                                tGTCTGAGATTACGTACAAACGTCGTATCCCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAATGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGt       <  1:343060/149‑1 (MQ=255)
                                                                                                                                     gagaTTACGTACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAt  <  1:295273/149‑1 (MQ=255)
                                                                                                                                     gagaTTACGTACAAACGTCGTAACTCCGCCCTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTTCCTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAt  <  2:164693/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
CTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCG‑CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCC‑CAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAT  >  NC_000913_3_bme_pgi/4182679‑4182959

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