Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,184,543	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,543	0	C	A	100.0%	36.7 / NA	16	P1100Q (CCG→CAG)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (6/10);  total (6/10)

TAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTC  >  NC_000913/4184412‑4184676                                                                                                                                    |                                                                                                                                      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)                                                                                                                                    |                                                                                                                                      TAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTC  >  NC_000913/4184412‑4184676

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: ‑