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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,182,820	0	C	T	100.0%	116.1 / NA	32	H526Y (CAC→TAC)	rpoB	RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/18);  total (14/18)

TCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTT  >  NC_000913/4182758‑4182878                                                               |                                                           tcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                          >  1:1999410/1‑65 (MQ=255)  cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtaca                                                         <  1:3355352/65‑1 (MQ=255)    tcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACaaa                                                       >  1:1960049/1‑65 (MQ=255)      ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGTTGTCTGAGATTACGTACAAAcg                                                     <  1:1160071/65‑1 (MQ=255)      ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                                     <  1:3153670/65‑1 (MQ=255)      ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                                     <  1:587923/65‑1 (MQ=255)        ttCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcgtc                                                   <  1:150335/65‑1 (MQ=255)                agcTGTCTCATTTTATGGACCAGNACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTcc                                           <  1:2922160/65‑1 (MQ=255)                 gcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCg                                          >  1:2885503/1‑65 (MQ=255)                       tcAGTTTATGTACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCg                                    <  1:2934876/65‑1 (MQ=255)                         ggTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGc                                  <  1:2615375/64‑1 (MQ=255)                         aGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGc                                  <  1:3610044/65‑1 (MQ=255)                         aGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGc                                  <  1:2670449/65‑1 (MQ=255)                         aGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGc                                  >  1:785224/1‑65 (MQ=255)                          gTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGcc                                 <  1:3405367/65‑1 (MQ=255)                           tttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGccc                                >  1:3584229/1‑65 (MQ=255)                             tATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAg                              <  1:2894873/65‑1 (MQ=255)                                 gACCAGGACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                          <  1:188684/65‑1 (MQ=255)                                 gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCcggcgg                          <  1:2187340/65‑1 (MQ=255)                                 gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                          <  1:1104042/65‑1 (MQ=255)                                 gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                          >  1:1906010/1‑65 (MQ=255)                                  aCCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGt                         <  1:1094591/65‑1 (MQ=255)                                     aGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTg                      <  1:2654035/65‑1 (MQ=255)                                        acaacCNGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAcc                   >  1:2206488/1‑65 (MQ=255)                                        acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAcc                   >  1:3113488/1‑65 (MQ=255)                                        acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAcc                   >  1:3542376/1‑65 (MQ=255)                                        acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAcc                   >  1:1519052/1‑65 (MQ=255)                                        acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAcc                   >  1:603141/1‑65 (MQ=255)                                             ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGa              >  1:29480/1‑65 (MQ=255)                                             ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGa              >  1:582409/1‑65 (MQ=255)                                             ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGCCCCAGGCGGTCTCACCCGCGa              >  1:2670699/1‑65 (MQ=255)                                               gCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGTTCTGACCCGTGAAc            <  1:2819236/65‑1 (MQ=255)                                                 tGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGt          <  1:2701466/65‑1 (MQ=255)                                                         aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCtt  >  1:1829942/1‑65 (MQ=255)                                                               |                                                           TCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTT  >  NC_000913/4182758‑4182878

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