Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007790	1,671	G→T	L31F (TTG→TTT)	SAUSA300_RS14710 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007790	1,671	0	G	T	100.0%	61.1 / NA	22	L31F (TTG→TTT)	SAUSA300_RS14710	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (21/1);  total (21/1)

AAAAACACAGTGTTAGATGAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTGGAAAATTTATATATATTAATTCCGTATGTTATTTTGATGGGAATAAATGCTATTTATTTGGTTGTTAAGGTTATGAATTTCAAAAAGAATAATT  >  NC_007790/1585‑1765                                                                                       |                                                                                               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:39520/1‑100 (MQ=255)       acaGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTaa                                                                             >  1:881286/1‑100 (MQ=255)       acaGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTaa                                                                             >  1:4252480/1‑100 (MQ=255)       acaGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTaa                                                                             >  1:4041095/1‑100 (MQ=255)                  taaaGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTAATTCTGTATGtt                                                                  >  1:2401180/3‑100 (MQ=255)                     aGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTAATTCTGTATGTTAtt                                                               >  1:353628/1‑100 (MQ=255)                                taGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTAATTCTGTATGTTATTTTGGTGGGAAt                                                    >  1:269815/1‑100 (MQ=255)                                taGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTAATTCTGTATGTTATTTTGATGGGAAt                                                    >  1:3629243/1‑100 (MQ=255)                                 aGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTAATTCTGTATGTTATTTTGATGGGAATa                                                   >  1:4562075/1‑100 (MQ=255)                                  gCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTAATTCTGTATGTTATTTTGATGGGagatc                                                  >  1:4562677/1‑96 (MQ=255)                                                                                  cATTTTGAAAATTTATATATATTAATTCTGTATGTTATTTTGATGGGAATAAATGCTATTTATTTGGTTGTTAAGGTTTTGAATTTCAAAAAGAATAAtt  >  1:3983873/1‑100 (MQ=255)                                                                                  cATTTTGAAAATTTATATATATTAATTCTGTATGTTATTTTGATGGGAATAAATGCTATTTATTTGGTTGTTAAGGTTTTGAATTTCAAAAAGAATAAtt  >  1:1536452/1‑100 (MQ=255)                                                                                  cATTTTGAAAATTTATATATATTAATTCTGTATGTTATTTTGATGGGAATAAATGCTATTTATTTGGTTGTTAAGGTTTTGAATTTCAAAAAGAATAAtt  >  1:3720321/1‑100 (MQ=255)                                                                                  cATTTTGAAAATTTATATATATTAATTCTGTATGTTATTTTGATGGGAATAAATGCTATTTATTTGGTTGTTAAGGTTTTGAATTTCAAAAAGAATAAtt  >  1:1887255/1‑100 (MQ=255)                                                                                       |                                                                                               AAAAACACAGTGTTAGATGAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTGGAAAATTTATATATATTAATTCCGTATGTTATTTTGATGGGAATAAATGCTATTTATTTGGTTGTTAAGGTTATGAATTTCAAAAAGAATAATT  >  NC_007790/1585‑1765

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: ‑