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%	85.3 / NA	31	E9K (GAA→AAA)	SAUSA300_RS14710	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/1);  total (30/1)

AAATAAAGGTATTTATAATTTGAGTAAGTCTTGTGATAAAGTTTTTTAAAAAAGGAGTGTTTATTTTGAATAAAAACACAGTGTTAGATGAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTGGAAAATTTATATATATTAATTCCGTAT  >  NC_007790/1514‑1698                                                                                          |                                                                                                aaaTAAAGGTATTTATAATTTGAGTAAGTCTTGTGATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGt                                                                                       <  1:950307/100‑1 (MQ=255)        ggTATTTATAATTTGAGTAAGTCTTGTGATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCa                                                                                >  1:2607006/1‑100 (MQ=255)                        gTAAGTCTTGTGATAAAGTTTTTTAAAAAAGGAGTCTTTATATTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGttt                                                                >  1:6066737/1‑100 (MQ=255)                          aaGTCTTGTGATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTc                                                              >  1:1228287/1‑100 (MQ=255)                              cTTGTGATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCtttt                                                          >  1:595129/1‑100 (MQ=255)                               ctgtgATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCttttt                                                         >  1:3979953/2‑100 (MQ=255)                                tgtgATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACGCCGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTa                                                        >  1:6102128/1‑100 (MQ=255)                                tgtgATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTa                                                        >  1:2182066/1‑100 (MQ=255)                                tgtgATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTa                                                        >  1:588003/1‑100 (MQ=255)                                tgtgATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTa                                                        >  1:4560922/1‑100 (MQ=255)                                 gtgATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTAc                                                       >  1:4162163/1‑100 (MQ=255)                                  tgATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACt                                                      >  1:734205/1‑100 (MQ=255)                                   gATAAAGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACtt                                                     >  1:1226096/1‑100 (MQ=255)                                      aaaGTTTTTTAAAAAAGGAGTCTTTATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGa                                                  >  1:3919042/1‑100 (MQ=255)                                                            gttATTTTGAATAAAGACACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCAttttg                            >  1:5791129/2‑98 (MQ=255)                                                                        atagacacAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTtatata                >  1:2050571/5‑100 (MQ=255)                                                                        aaagaCACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTtatata                >  1:4220468/1‑100 (MQ=255)                                                                        aaagaCACAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTtatata                >  1:4101120/1‑100 (MQ=255)                                                                         cagacacAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTtatatat               >  1:631161/4‑100 (MQ=255)                                                                          tgacacAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTtatatata              >  1:2695337/3‑100 (MQ=255)                                                                          agacacAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTtatatata              >  1:1495101/3‑100 (MQ=255)                                                                           gacacAGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTatatatat             >  1:5390205/2‑100 (MQ=255)                                                                             cacaGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTa           >  1:5853594/1‑100 (MQ=255)                                                                              acaGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGGTTAATTTTTATATCTCATTTTGAAAATTTATATATATTaa          >  1:5270955/1‑100 (MQ=255)                                                                              acaGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTaa          >  1:5508835/1‑100 (MQ=255)                                                                               caGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGCAAATTTATATATATTAAt         >  1:1200192/1‑100 (MQ=255)                                                                               caGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTAAt         >  1:5974105/1‑100 (MQ=255)                                                                               caGTGTTAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTAAt         >  1:1925805/1‑100 (MQ=255)                                                                                   cctAGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTAATTctg     >  1:5069649/3‑98 (MQ=255)                                                                                     tcgATAAAGGTTACGTTATAGCAACTGTTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTAATTCTGTa   >  1:2947514/3‑100 (MQ=255)                                                                                      aGATAAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTTGAAAATTTATATATATTAATTCTGTAt  >  1:238327/1‑100 (MQ=255)                                                                                          |                                                                                                AAATAAAGGTATTTATAATTTGAGTAAGTCTTGTGATAAAGTTTTTTAAAAAAGGAGTGTTTATTTTGAATAAAAACACAGTGTTAGATGAAGGTTACGTTATAGCAACTATTTTAAATGTTTTCTTTTTACTTGGATTAATTTTTATATCTCATTTGGAAAATTTATATATATTAATTCCGTAT  >  NC_007790/1514‑1698

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