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%	111.5 / NA	31	H526Y (CAC→TAC)	rpoB	RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/24);  total (7/24)

TCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAA  >  NC_000913/4182761‑4182867                                                            |                                                tcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACaaa                                            >  1:3380917/1‑65 (MQ=255)   ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                          <  1:1978220/65‑1 (MQ=255)   ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                          <  1:371589/65‑1 (MQ=255)   ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                          <  1:3819354/65‑1 (MQ=255)   ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                          <  1:1368632/65‑1 (MQ=255)      tCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcgtcg                                       >  1:1954744/1‑65 (MQ=255)          gccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTAtc                                   <  1:2405338/65‑1 (MQ=255)          gccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTAtc                                   <  1:2920410/65‑1 (MQ=255)           ccagcTGTCTCAGTTTATGGACCAGAACAGCCCGCTGTCTGAGATTACGTACAAACGTCGTAtct                                  <  1:2852623/65‑1 (MQ=255)           ccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTAtct                                  <  1:2177076/65‑1 (MQ=255)           ccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTAtct                                  <  1:1280166/65‑1 (MQ=255)           ccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTAtct                                  <  1:1232633/65‑1 (MQ=255)               cTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGc                              <  1:1276814/65‑1 (MQ=255)                  tctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACt                           <  1:3599394/65‑1 (MQ=255)                    tcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCg                         <  1:1030636/65‑1 (MQ=255)                     cAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCgg                        <  1:938832/65‑1 (MQ=255)                      aGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGc                       <  1:3624307/65‑1 (MQ=255)                        tttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGccc                     >  1:1940805/1‑65 (MQ=255)                           caggACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAgg                  <  1:1813492/63‑1 (MQ=255)                              gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg               <  1:947228/65‑1 (MQ=255)                              gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg               <  1:3227848/65‑1 (MQ=255)                              gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg               <  1:3126097/65‑1 (MQ=255)                              gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg               <  1:3078390/65‑1 (MQ=255)                              gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg               <  1:2957012/65‑1 (MQ=255)                               aCCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGt              <  1:1649706/65‑1 (MQ=255)                                   gAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGa          <  1:443329/65‑1 (MQ=255)                                     acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTgtac        >  1:559585/1‑62 (MQ=255)                                     acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAct        >  1:3647920/1‑64 (MQ=255)                                      caacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAccc       <  1:1657861/65‑1 (MQ=255)                                          ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGa   >  1:2908573/1‑65 (MQ=255)                                          ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGa   >  1:414656/1‑65 (MQ=255)                                           cGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGaa  <  1:574/65‑1 (MQ=255)                                                            |                                                TCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAA  >  NC_000913/4182761‑4182867

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