Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	347,608	C→G	Q153E (CAA→GAA)	SAUSA300_RS01570 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	347,608	0	C	G	100.0%	52.0 / NA	20	Q153E (CAA→GAA)	SAUSA300_RS01570	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (18/2);  total (18/2)

AATTAGGCCAAGAAAACTACTATAACTATAGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAACAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGAACAATAGAGGCGATAACATGATTTTCGAAGAAAAACTAAATGAAATGTACAA  >  NC_007793/347528‑347697                                                                                 |                                                                                          aaTTAGGCCAAGAAAACTACTATAACTATAGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTaa                                                                        >  1:3018067/1‑100 (MQ=255)  aTTAGGCCAAGAAAACTACTATAACTATAGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTaaa                                                                       >  1:510048/1‑100 (MQ=255)   ccaGGCCAAGAAAACTACTATAACTATAGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAag                                                                      >  1:3136249/3‑100 (MQ=255)    tAGGCCAAGAAAACTACTATAACTATAGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAaga                                                                     >  1:4090453/1‑100 (MQ=255)     aGGCCAAGAAAACTACTATAACTATAGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAagag                                                                    >  1:3512461/1‑100 (MQ=255)        ccAAGAAAACTACTATAACTATAGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGGTAAAGAAATCGAGCAATATATTAAAGAGCaa                                                                 >  1:1556640/1‑100 (MQ=255)          aaGAAAACTACTATAACTATAGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCaaga                                                               <  1:5501031/100‑1 (MQ=255)                   actataactataGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGa                                                      >  1:2881400/1‑100 (MQ=255)                    tgctaactataGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGaa                                                     >  1:569166/4‑100 (MQ=255)                    aaataactataGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGccaga                                                     >  1:2852781/3‑96 (MQ=255)                        aactataGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGAACAAt                                                 >  1:1931195/1‑100 (MQ=255)                        aactataGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGAACAAt                                                 >  1:156901/1‑100 (MQ=255)                        aactataGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGAACAAt                                                 >  1:1380226/1‑100 (MQ=255)                              aGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGAACAATAGAGGc                                           <  1:2333107/100‑1 (MQ=255)                                            tataCCAGAAATGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGAACAATAGAGGCGATAACATGAtttt                             >  1:1617078/1‑100 (MQ=255)                                                     aaTGGAATATGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGAACAATAGAGGCGATAACATGATTTTCGAAGaaaa                    >  1:4953906/1‑100 (MQ=255)                                                             tggtaaTGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGAACAATAGAGGCGATAGCATGATTTTCGAAGAAAAACTaaatg            >  1:1703925/5‑100 (MQ=255)                                                             atGAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGAACAATAGAGGCGATAACATGATTTTCGAAGAAAAACTaaatg            >  1:3977851/1‑100 (MQ=255)                                                               gAAATGGAAGAGGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGAACAATAGAGGCGATAACATGATTTTCGAAGAAAAACTaaatgaa          >  1:2716338/1‑100 (MQ=255)                                                                       agagGTTAAAGAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGAACAATAGAGGCGATAACATGATTTTCGAAGAAAAACTAAATGAAATGTACaa  >  1:916188/1‑100 (MQ=255)                                                                                 |                                                                                          AATTAGGCCAAGAAAACTACTATAACTATAGAAAATTCGGAATTATACCAGAAATGGAATATGAAATGGAAGAGGTTAAACAAATCGAGCAATATATTAAAGAGCAAGAAGAAGCTGAACAATAGAGGCGATAACATGATTTTCGAAGAAAAACTAAATGAAATGTACAA  >  NC_007793/347528‑347697

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