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

GGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCG  >  NC_000913/4182763‑4182880                                                          |                                                             ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                                       <  1:1603696/65‑1 (MQ=255)   ttCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcgtc                                                     <  1:571942/65‑1 (MQ=255)        gccaacTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTAtc                                                <  1:2280994/65‑1 (MQ=255)                  tcAGTTTATGGACCAGAACGACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCg                                      <  1:2032805/65‑1 (MQ=255)                  tcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCg                                      <  1:1253942/65‑1 (MQ=255)                  tcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCg                                      <  1:1913983/65‑1 (MQ=255)                    aGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGc                                    <  1:101706/65‑1 (MQ=255)                        tATGGACCAGAACAACCCGCTGTCTGAAATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAg                                <  1:872158/65‑1 (MQ=255)                            gACCAGAACAACCTGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                            <  1:1090408/65‑1 (MQ=255)                            gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                            <  1:2136844/65‑1 (MQ=255)                            gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                            <  1:686170/65‑1 (MQ=255)                            gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                            >  1:358291/1‑65 (MQ=255)                            gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                            <  1:3566310/65‑1 (MQ=255)                            gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                            <  1:3086623/65‑1 (MQ=255)                            gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                            >  1:2082553/1‑65 (MQ=255)                            gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                            <  1:2030155/65‑1 (MQ=255)                                   acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAcc                     >  1:3418447/1‑65 (MQ=255)                                     aacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCg                   >  1:1076687/1‑65 (MQ=255)                                        ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGa                >  1:3581424/1‑65 (MQ=255)                                        ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGa                >  1:1167465/1‑65 (MQ=255)                                         cGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGaa               <  1:268559/65‑1 (MQ=255)                                            tGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGt            <  1:1724284/65‑1 (MQ=255)                                                      tACGTACAAACGTCGTATCTCCGCACTGGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCg  >  1:440502/1‑65 (MQ=255)                                                          |                                                             GGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCG  >  NC_000913/4182763‑4182880

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