Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,184,542	C→A	P1100Q (CCG→CAG)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,184,542	0	C	A	100.0%	39.7 / NA	17	P1100Q (CCG→CAG)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (6/11);  total (6/11)

TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTC  >  NC_000913/4184402‑4184675                                                                                                                                             |                                                                                                                                      ttttgcttttAAACGCCGTAGCCAGCCTGGTGCCAAGAGGGGAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCCAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGc                                                                                                                                     <  1:44019/135‑1 (MQ=255)          tAAATGCGGTATTGATCCTGGGGACATGAGGGCAGGGTGTCACTGTAACGAGGGGGTACTTTCTAAGATCGACTCGCGCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAc                                                                                                                            <  2:90389/143‑1 (MQ=255)                ccGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                      >  2:89479/1‑143 (MQ=255)                  gTATCCATCCTGGTGACAAGATTGCAGGTCGTCACGGTAACAAGCGTGTATTTTCTATGATCAACCCGATCGAAGATATGCTTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGt                                                                                                                    <  1:48895/143‑1 (MQ=255)                       cAGCCTGGTGACAAGATGGCAGGTCGTCCCGGTATCACGGCCGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                               <  2:175365/143‑1 (MQ=255)                              gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATAAAAACGGTACGCCAGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                        >  1:67456/1‑143 (MQ=255)                               tGACAAGCTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGCTATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                                       <  2:229466/143‑1 (MQ=255)                                 aCAAGATTGCAGGTCGTCACGGTATCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAATCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAg                                                                                                     <  2:62446/143‑1 (MQ=255)                                   aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATGAACCCGATCGAAGATATGACTTCCGGTGAAAACGGTACGCCGGGAGAAACGGTACTGAAGCAGCTGGGGGGACCGGCTAGTACGAACATCGGTCAGAt                                                                                                   >  2:254843/1‑143 (MQ=255)                                   aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGGATGAACATCGGTCAGAt                                                                                                   >  2:380712/1‑143 (MQ=255)                                                             gggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTg                                                                         <  1:199728/143‑1 (MQ=255)                                                                tgtAATTTCTAAGATAAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACAACGTACTAAACCAGCTGGGCGTAACGTATCGTATGAACATAGGTAAGATCCTAGAAACCCACCTAGGTATGGCTGCGa                                                                      >  1:109267/1‑143 (MQ=255)                                                                 gtAATTGCTAAGTTCAACCTGATCGATGATATGCTTTTCTATTTTTACGTTACGTCGTTAGACATCGTACGGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                     <  1:53478/143‑1 (MQ=255)                                                                   tctttCTAAGATCGTCTCTCTCGAAGATATTCCTTACGATGGAAACGTTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAGCTTCTGTCAGATCCTCGATACCCACCTGTGTATGGCTGCGAAAg                                                                   <  2:202800/141‑1 (MQ=255)                                                                          aaGATCAGCCCGATGGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGACAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCgg                                                            <  2:305249/143‑1 (MQ=255)                                                                          aaGATCACCCCGATCGAAGATATGTCTTATGATGATAACGGTACGCCGGTAGACACCGTACTCAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAATCCCACCTGGGTATGGCTGCGAAAGGTATCgg                                                            <  2:265006/143‑1 (MQ=255)                                                                                                                                    tACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGTTAGGTATCGGCGACTAGATAAACGCCATGCTGAAAAAGCAGAAAGAAGTCGCGAAACTGAGCGAATTc  >  1:41109/1‑143 (MQ=255)                                                                                                                                             |                                                                                                                                      TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTC  >  NC_000913/4184402‑4184675

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