Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,037,513	T→C	intergenic (+178/‑6)	alaT → / → rrlA	tRNA‑Ala/23S ribosomal RNA of rrnA operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,037,513	0	T	C	100.0%	26.9 / NA	12	intergenic (+178/‑6)	alaT/rrlA	tRNA‑Ala/23S ribosomal RNA of rrnA operon
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/6);  total (6/6)

GTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGT  >  NC_000913/4037371‑4037543                                                                                                                                               |                               gtACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGCTGTGAGGt                        >  1:240493/1‑151 (MQ=25) gtACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGCTGTGAGGt                        <  2:658972/151‑1 (MQ=25)        cAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGCTGTGAGGTTAAGCg                  >  1:1029257/1‑150 (MQ=25)            ggTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGCTGTGAGGTTAAGCGACTa              >  1:245958/1‑150 (MQ=25)                cACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGCTGTGAGGTTAAGCGACTAAGCGt         <  2:220101/151‑1 (MQ=25)                 aCTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTATTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGCTGTGAGGTTAAGCGACTAAGCGTa        <  1:320459/151‑1 (MQ=25)                    gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGCTGTGAGGTTAAGCGACTAAGCGTacac     <  1:909834/151‑1 (MQ=21)                    gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGCTGTGAGGTTAAGCGACTAAGCGTacac     >  2:359114/1‑151 (MQ=21)                    gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGCTGTGAGGTTAAGCGACTAAGCGTacac     <  2:661363/151‑1 (MQ=21)                    gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGCTGTGAGGTTAAGCGACTAAGCGTac       <  2:117525/149‑1 (MQ=21)                     cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGCTGTGAGGTTAAGCGACTAAGCGTACACg    >  2:986118/1‑151 (MQ=18)                       aaGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGCTGTGAGGTTAAGCGACTAAGCGTACACGGt  >  2:914517/1‑151 (MQ=17)                                                                                                                                               |                               GTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGT  >  NC_000913/4037371‑4037543

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