Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007790	1,603	G→A	E9K (GAA→AAA)	SAUSA300_RS14710 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007790	1,603	0	G	A	100.0%	47.0 / NA	17	E9K (GAA→AAA)	SAUSA300_RS14710	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/1);  total (16/1)

AGGTATTTATAATTTGAGTAAGTCTTGTGATAAAGTTTTTTAAAAAAGGAGTGTTTATTTTGAATAAAAACACAGTGTTAGATGAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTGGAAAATTTATATATATTAA  >  NC_007790/1520‑1690                                                                                    |                                                                                        aGGTATTTATAATTTGAGTAAGTCTTGTGATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGc                                                                         >  1:337650/1‑100 (MQ=255)                  agaaGTCTTGTGATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGttt                                                        >  1:591011/3‑100 (MQ=255)                          tgtgATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTa                                                >  1:1928985/1‑100 (MQ=255)                          tgtgATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTa                                                >  1:3991439/1‑100 (MQ=255)                                aaaGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGa                                          >  1:2982867/1‑100 (MQ=255)                                                                  taagacacAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTtatata        >  1:2606373/5‑100 (MQ=255)                                                                   ctgacacAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTtatatat       >  1:2855425/4‑100 (MQ=255)                                                                   aagacacAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTtatatat       >  1:3525581/4‑100 (MQ=255)                                                                    cgacacAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTtatatata      >  1:2221741/3‑100 (MQ=255)                                                                    agacacAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTtatatata      >  1:896156/3‑100 (MQ=255)                                                                     gacacAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTatatatat     >  1:399492/2‑100 (MQ=255)                                                                      acacAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATAtt    >  1:977365/1‑100 (MQ=255)                                                                       cacaGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTa   <  1:705702/100‑1 (MQ=255)                                                                        acaGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTaa  >  1:4041095/1‑100 (MQ=255)                                                                        acaGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTaa  >  1:4252480/1‑100 (MQ=255)                                                                        acaGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTaa  >  1:39520/1‑100 (MQ=255)                                                                        acaGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTaa  >  1:881286/1‑100 (MQ=255)                                                                                    |                                                                                        AGGTATTTATAATTTGAGTAAGTCTTGTGATAAAGTTTTTTAAAAAAGGAGTGTTTATTTTGAATAAAAACACAGTGTTAGATGAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTGGAAAATTTATATATATTAA  >  NC_007790/1520‑1690

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