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

CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTA  >  NC_000913/4182743‑4182894                                                                              |                                                                           cagcagTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtaca                                                                         <  1:829573/81‑1 (MQ=255)   gcagTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACaaa                                                                       <  2:829574/81‑1 (MQ=255)    cagTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAc                                                                      <  2:829575/81‑1 (MQ=255)         aaaGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcgtcgt                                                                 <  1:829576/81‑1 (MQ=255)         aaaGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcgtcg                                                                  <  2:829577/80‑1 (MQ=255)            gagTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTAtc                                                              <  1:829578/81‑1 (MQ=255)            gagTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTAtc                                                              >  2:829657/1‑81 (MQ=255)               ttcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTcc                                                           >  1:829644/1‑81 (MQ=255)                  ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCa                                                        >  1:829652/1‑81 (MQ=255)                     ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTc                                                     <  2:829579/81‑1 (MQ=255)                      gTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCg                                                    <  2:829580/81‑1 (MQ=255)                       ttCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCgg                                                   >  2:829640/1‑81 (MQ=255)                       ttCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCgg                                                   <  1:829581/81‑1 (MQ=255)                         ccagccagCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGc                                                  <  1:829582/80‑1 (MQ=255)                          cagccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGccc                                                <  2:829584/81‑1 (MQ=255)                          cagccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGccc                                                <  1:829585/81‑1 (MQ=255)                          cagccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGcc                                                 <  2:829583/80‑1 (MQ=255)                           agccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGcc                                                 >  2:829664/1‑79 (MQ=255)                            gccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAg                                              >  1:829651/1‑81 (MQ=255)                                gcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                          >  2:829655/1‑81 (MQ=255)                                 cTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGt                                         <  1:829587/81‑1 (MQ=255)                                 cTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGt                                         <  2:829586/81‑1 (MQ=255)                                    tctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCt                                       <  2:829589/80‑1 (MQ=255)                                    tctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCt                                       <  1:829590/80‑1 (MQ=255)                                    tctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTg                                      <  1:829588/81‑1 (MQ=255)                                     ctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGa                                     >  2:829648/1‑81 (MQ=255)                                        aGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAcc                                   <  1:829592/80‑1 (MQ=255)                                        aGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAcc                                   <  1:829591/80‑1 (MQ=255)                                        aGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAc                                    <  2:829593/79‑1 (MQ=255)                                          tttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                                >  2:829666/1‑81 (MQ=255)                                             aTGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGaa                             >  1:829647/1‑81 (MQ=255)                                             aTGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGaa                             >  1:829668/1‑81 (MQ=255)                                                    aGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAg                      <  2:829594/81‑1 (MQ=255)                                                    aGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAg                      >  1:829663/1‑81 (MQ=255)                                                       acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCt                   <  2:829595/81‑1 (MQ=255)                                                       acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCt                   <  2:829596/81‑1 (MQ=255)                                                        caacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCtt                  >  1:829674/1‑81 (MQ=255)                                                          acCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTc                 >  2:829660/1‑80 (MQ=255)                                                          acCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCg                <  1:829597/81‑1 (MQ=255)                                                            ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCg                >  1:829646/1‑79 (MQ=255)                                                                 gTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTc          >  2:829661/1‑80 (MQ=255)                                                                 gTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTc          >  1:829654/1‑80 (MQ=255)                                                                  tCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCg         <  2:829598/80‑1 (MQ=255)                                                                      agaTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCgag       <  2:829601/78‑1 (MQ=255)                                                                      agaTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCgag       <  2:829600/78‑1 (MQ=255)                                                                      agaTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGAc     <  2:829599/80‑1 (MQ=255)                                                                       gaTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACg    >  2:829667/1‑80 (MQ=255)                                                                        aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa  >  2:829650/1‑81 (MQ=255)                                                                              |                                                                           CAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTA  >  NC_000913/4182743‑4182894

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