Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	789,781	T→A	I489F (ATC→TTC)	SAUSA300_RS03825 ←	peptide ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	789,781	0	T	A	100.0%	16.1 / NA	8	I489F (ATC→TTC)	SAUSA300_RS03825	peptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (3/5);  total (3/5)

TACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATTGTTGTAACAATAATTGCAACAACACCTAGGAAAATAAAGTAATTTGTTTGACCTAGTGGTT  >  NC_007793/789712‑789855                                                                      |                                                                           tACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAATGCTAATACAATTGTTGTAACAATAATTGc                                              >  1:426139/1‑100 (MQ=255) tACTGCTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTGTGATGATTAAAGGTGAGAATGCTAATACAATTGTTGTAACAATAATTGc                                              <  1:3747091/100‑1 (MQ=255)           tGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAATGCTAATACAATTGTTGTAACAATAATTGCAACAACAccg                                    <  1:276020/100‑2 (MQ=255)           ccgatctggCCAGCAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAATGCTAATACAATTGTTGTAACAATAATTGCAACAACACCt                                    <  1:6336317/93‑1 (MQ=255)              aTTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAATGCTAATACAATTGTTGTAACAATAATTGCAACAACACCTAgg                                 <  1:6291646/100‑1 (MQ=255)                tgggCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAATGCTAATACAATTGTTGTAACAATAATTGCAACAACACCTAGGaa                               >  1:2567493/3‑100 (MQ=255)                                aCGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGAATGCTAATACAATTGTTGTAACAATAATTGCAACAACACCTAGGAAAATAAAGTAAtttgtt               <  1:927272/100‑1 (MQ=255)                                             aTCGCTTTGATGATTAAAGGTGAGAATGCTAATACAATTGTTGTAACAATAATTGCAACAACACCTAGGAAAATAAAGTAATTTGTTTGACCTAGTGGtt  >  1:58436/1‑100 (MQ=255)                                                                      |                                                                           TACTGTTTTTTGAATTAGGCCAACAATATTAACGTATACCTTTCATCGCTTTGATGATTAAAGGTGAGATTGCTAATACAATTGTTGTAACAATAATTGCAACAACACCTAGGAAAATAAAGTAATTTGTTTGACCTAGTGGTT  >  NC_007793/789712‑789855

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