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

CTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAT  >  NC_000913_3_bme_pgi/4182679‑4182959
                                                                                                                                              |                                                                                                                                            
ctgtctCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACACCCCGCTGTCTAAGATTACGTACAAAcg                                                                                                                                       >  2:183682/1‑149 (MQ=255)
     tctGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGACCAAACGTCGTa                                                                                                                                  >  1:7551/1‑149 (MQ=255)
         ggCGATCTGGATACGCTGATGCCACAGTTTCTGGTCACGGCCCAGCCGATTTCCGCAGCAGTGTAAGACATCTTCGGTGCAAGCCAGCGGGCTGAGTTTGTGTACCAGAACAACCCGCTGTCTGAGTTTACGTACAAACGTCGTAtatc                                                                                                                              >  2:343060/1‑149 (MQ=255)
            gATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGc                                                                                                                           >  1:612303/1‑149 (MQ=255)
            gATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGAATCTCCGc                                                                                                                           >  1:329435/1‑149 (MQ=255)
              tCTGGATACCCTGAGGCCACAGGAGATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCAc                                                                                                                         >  1:552699/1‑149 (MQ=255)
                    tACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGcc                                                                                                                   >  2:89003/1‑149 (MQ=255)
                       ccTGATGCCACAGGATATGATCAACGCCAGGACATTTCCGCCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCCCACTCGGCCCAg                                                                                                               >  1:474976/1‑149 (MQ=255)
                        cTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGGGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAgg                                                                                                              >  2:343130/1‑149 (MQ=255)
                        cTGATGCCACAGGATATGATCAACGCCAAGCAAATTTCCTCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTTTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGACATTACGTACAAACGTCGTATCTACGCACTCGGCCCAgg                                                                                                              >  2:300856/1‑149 (MQ=255)
                             gCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTc                                                                                                         <  1:89003/149‑1 (MQ=255)
                               cacaGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTg                                                                                                       >  2:308527/1‑149 (MQ=255)
                               cacaGGATATGATCAACGCCAAGCCGAGTTCCGCAGCCGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTCTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGCCCCAGGCgcgctg                                                                                                       >  2:329716/1‑144 (MQ=255)
                                      tatGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTg                                                                                                >  1:130582/1‑149 (MQ=255)
                                                        gATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCTCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAg                                                                              <  1:506188/149‑1 (MQ=255)
                                                           ttCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTc                                                                           <  1:29233/149‑1 (MQ=255)
                                                           ttCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTc                                                                           <  2:130582/149‑1 (MQ=255)
                                                             ccGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATCACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCga                                                                         <  1:89306/149‑1 (MQ=255)
                                                               gcagcaGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGGTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCgaga                                                                       <  1:567691/149‑1 (MQ=255)
                                                                        gAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAAAGCGCAGGCCGCGAAGTTCGAGCCGTACAcc                                                               >  2:265214/1‑149 (MQ=255)
                                                                         aaaGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAccc                                                              >  1:192271/1‑149 (MQ=255)
                                                                             agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACt                                                          >  1:509517/1‑149 (MQ=255)
                                                                               ttcttgGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGCACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCa                                                        <  1:59280/149‑1 (MQ=255)
                                                                               ttcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGaa                                                                               >  1:27844/1‑126 (MQ=255)
                                                                               ttcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGaa                                                                               <  2:27844/126‑1 (MQ=255)
                                                                                 cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTCCACCCGACTCACt                                                      >  2:449342/1‑149 (MQ=255)
                                                                                  ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTa                                                     <  1:324830/149‑1 (MQ=255)
                                                                                         ccagccagCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgc                                              <  1:5763/149‑1 (MQ=255)
                                                                                          tagcgagcTGTCTGAGTTTATGGACCAGAACAACCCGCTTTCTGAGATTACGTATAAACGTCGTATCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgcg                                             <  2:570029/148‑1 (MQ=255)
                                                                                           agccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGt                                            >  1:323147/1‑149 (MQ=255)
                                                                                                gcTGTCTCACTTTACGGCCCAGAACTACCCGTTCTCTGTGATTTCGTACCCCCGTCGTATCTCCGCCCCGGCCGGAGGCAGTCTTGGATGGCGACTTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTc                                       <  1:627020/149‑1 (MQ=255)
                                                                                                   gTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCaa                                    <  2:323147/149‑1 (MQ=255)
                                                                                                     ctcAGTTTATGGAACAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGCGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATc                                  <  1:183682/149‑1 (MQ=255)
                                                                                                                gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGCCTCACTACGGTCGCGTATGTCCACTCGAAAccccggg                       >  1:12239/1‑146 (MQ=255)
                                                                                                                 aCCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGaa                      >  1:559194/1‑149 (MQ=255)
                                                                                                                  ccAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTAGCGTATGTCCAATCGAACCTCCTTAAg                     >  2:67437/1‑149 (MQ=255)
                                                                                                                    aGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGt                   <  1:253010/149‑1 (MQ=255)
                                                                                                                     gAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTGCCGTATGTCCAATCGAAACCCCTGAAGGTc                  >  1:250295/1‑149 (MQ=255)
                                                                                                                     gAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGGAGGTc                  >  2:385455/1‑149 (MQ=255)
                                                                                                                      accaacCCGCTGTGTGGGATTACGTACAAACGTCGGATGCCCGCACTCGGCCAAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGAGTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTcc                 <  2:107314/147‑1 (MQ=255)
                                                                                                                       acaacCCGCTGGCTGAGATTACGTACACACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCg                <  2:393161/149‑1 (MQ=255)
                                                                                                                            ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGAc                                                           >  2:318534/1‑101 (MQ=255)
                                                                                                                            ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGAc                                                           <  1:318534/101‑1 (MQ=255)
                                                                                                                             cGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATc          >  2:350554/1‑149 (MQ=255)
                                                                                                                                tGTCTGAGATTACGTACAAACGTCGTATCCCCGCACTCGGCCCAGGCGGTCTGACCCGTGAATGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGt       <  1:343060/149‑1 (MQ=255)
                                                                                                                                     gagaTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAt  <  1:295273/149‑1 (MQ=255)
                                                                                                                                     gagaTTACGTACAAACGTCGTAACTCCGCCCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTTCCTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAt  <  2:164693/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
CTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAT  >  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: