Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,579,338	A→G	N163N (AAT→AAC)	SAUSA300_RS07730 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,579,338	0	A	G	100.0%	82.4 / NA	33	N163N (AAT→AAC)	SAUSA300_RS07730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/21);  total (12/21)

CCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCATGATTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCTGGTCGATAGAATATAAGTTGAAAGCTATGTCTATCGCTATTACAATTGCCAAT  >  NC_007793/1579257‑1579428                                                                                   |                                                                                           cttccgatctcctTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCttt                                                                          <  1:5689309/92‑1 (MQ=38)  ccAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCcttt                                                                         >  1:4564413/1‑100 (MQ=38)     aaaCCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGt                                                                      <  1:319618/100‑1 (MQ=255)     aaaCCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGt                                                                      <  1:2132451/100‑1 (MQ=255)              attcTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCggata                                                             <  1:3772203/97‑5 (MQ=255)                   cGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAgt                                                        <  1:4781778/100‑2 (MQ=255)                    gATCTAATGTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAgtt                                                       <  1:6626718/100‑3 (MQ=255)                    gATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAgtt                                                       <  1:2160347/100‑3 (MQ=255)                    gATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAgtt                                                       <  1:1828343/100‑3 (MQ=255)                    gATCTAAAATGTCTGTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAgtt                                                       <  1:4732390/100‑3 (MQ=255)                         aataTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCAGATCg                                                  >  1:132888/4‑100 (MQ=255)                         aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAgatagatc                                                  >  1:3227921/4‑93 (MQ=255)                         aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTct                                                  <  1:6616270/97‑2 (MQ=255)                         aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCg                                                  >  1:722819/4‑100 (MQ=255)                         aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCg                                                  >  1:3111560/4‑100 (MQ=255)                         aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCg                                                  >  1:3280892/4‑100 (MQ=255)                         aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCg                                                  >  1:1652475/4‑100 (MQ=255)                         aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATATTGTCTTTCGTGGTTCCTTTTCATTTTCCTTTTGTATGCGTCAATGAGTTGGTCg                                                  >  1:6650409/4‑100 (MQ=255)                         aaGATGTCTTTACATGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCAGATCg                                                  >  1:5540304/1‑100 (MQ=255)                         aaGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAgtcagatc                                                  >  1:6141111/1‑93 (MQ=255)                         aaGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCtagatc                                                  >  1:2069662/1‑95 (MQ=255)                                cTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGATCGATAgtgt                                           <  1:2171842/100‑4 (MQ=255)                                  ttACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAatat                                         <  1:1606397/100‑1 (MQ=255)                                         tATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTgt                                  <  1:4394623/100‑2 (MQ=255)                                         tATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTgt                                  <  1:5644890/100‑2 (MQ=255)                                           tCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTgtag                                <  1:4759293/100‑4 (MQ=255)                                                ttCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTAt                           >  1:3056748/1‑100 (MQ=255)                                                  ccTGCGTCGGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGt                         <  1:963394/100‑1 (MQ=255)                                                  cagacGTGTGCTCTTCCGATC‑TGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGt                         <  1:3852389/96‑1 (MQ=255)                                                                  c‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTagctgggc         <  1:3300474/100‑8 (MQ=255)                                                                          tCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGcttgt  <  1:4233818/100‑5 (MQ=255)                                                                          ggcttCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGcttgt  <  1:3104356/97‑5 (MQ=255)                                                                          ccTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGAATATAAGTTGTAAGCTATGTCTATCACTATAACAATTGcttgt  <  1:3510153/99‑5 (MQ=255)                                                                                   |                                                                                           CCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCATGATTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCTGGTCGATAGAATATAAGTTGAAAGCTATGTCTATCGCTATTACAATTGCCAAT  >  NC_007793/1579257‑1579428

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