Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,399,569	A→G	K415E (AAA→GAA)	SAUSA300_RS06895 →	aminoacyltransferase FemB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,399,569	0	A	G	100.0%	37.5 / NA	15	K415E (AAA→GAA)	SAUSA300_RS06895	aminoacyltransferase FemB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/3);  total (12/3)

TGGATTTAATGTACAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTAAAAAATTAAAGAAATAGATGTAAATGATGTTAGGGCTTATAGTTATTGATACTATAGGCTCTTTTTTATATG  >  NC_007793/1399472‑1399641                                                                                                  |                                                                         tGGATTTAATGTACAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTgaa                                                                        <  1:816153/100‑1 (MQ=255)  ggATTTAATGTACAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTgaaa                                                                       >  1:460391/1‑100 (MQ=255)       tAATGTACAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTGAAAAATTa                                                                  >  1:5614955/1‑100 (MQ=255)        aaTGTACAAGTCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTGAAAAATTaa                                                                 >  1:295549/1‑100 (MQ=255)          tGTACAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTGAAAAATTaagg                                                               >  1:2126003/1‑98 (MQ=255)          agTACAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTGAAAAATTAAAg                                                               >  1:4177701/2‑100 (MQ=255)           gTACAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTGAAAAATTAAGGa                                                              >  1:1654598/1‑100 (MQ=255)           gTACAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTGAAAAATTAAAGa                                                              >  1:3086723/1‑100 (MQ=255)           gTACAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTGAAAAATTAAAGa                                                              >  1:887880/1‑100 (MQ=255)           gTACAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTGAAAAATTAAAGa                                                              >  1:4677330/1‑100 (MQ=255)             aagcagcGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTGAAAAATTAAAGaaa                                                            >  1:2306022/7‑100 (MQ=255)              cAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTGAAAAATTAAAGAAAt                                                           <  1:1048582/100‑1 (MQ=255)              cAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTGAAAAATTAAAGAAAt                                                           <  1:2129176/100‑1 (MQ=255)                   agacgcaTTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACGTTGGATAAATTACGTGAAAAATTAAAGAAATAGATg                                                      >  1:3092417/7‑100 (MQ=255)                            aaTAGGGGATTTCTATAAACCAATTCATTATGTGAATTTTTGTTTGTTCACAACTTTGGTTAAATTACGTGAAAAATTAAAGAAATAGAtcggaagagca                                             >  1:781368/1‑90 (MQ=255)                                                       tAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTGAAAAATTAAAGAAATAGATGTAAATGATGTTAGGGCTTATAGTTATTGATACtatt                  <  1:1672979/100‑2 (MQ=255)                                                                       ctgttCACAACATTGGATAAATTACGTGAAAAATTAAAGAAATAGATGTAAATGATGTTAGGGCTTATAGTTATTGATACTATAGGCTCTTTTTTATATg  >  1:1214348/2‑100 (MQ=255)                                                                                                  |                                                                         TGGATTTAATGTACAAATCGAAGAATTAATAGGGGATTTCTATAAACCAATTCATAAAGTGAAATATTGGTTGTTCACAACATTGGATAAATTACGTAAAAAATTAAAGAAATAGATGTAAATGATGTTAGGGCTTATAGTTATTGATACTATAGGCTCTTTTTTATATG  >  NC_007793/1399472‑1399641

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