Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,579,335	A→G	H164H (CAT→CAC)	SAUSA300_RS07730 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,579,335	0	A	G	100.0%	85.8 / NA	34	H164H (CAT→CAC)	SAUSA300_RS07730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/22);  total (12/22)

AAGTATTGCGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCATGATTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCTGGTCGATAGAATATAAGTTGAAAGCTATGTCTATCGCTATTACAATTGCCAAT  >  NC_007793/1579240‑1579428                                                                                                 |                                                                                              aaGTATTGTGTCATTTCCCCAAAACCTCCTCGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTggt                                                                                           <  1:6414682/100‑3 (MQ=255)                  cttccgatctcctTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCttt                                                                          <  1:5689309/92‑1 (MQ=38)                   ccAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCcttt                                                                         >  1:4564413/1‑100 (MQ=38)                      aaaCCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGt                                                                      <  1:2132451/100‑1 (MQ=255)                      aaaCCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGt                                                                      <  1:319618/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:3280892/4‑100 (MQ=255)                                          aaaaTGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCg                                                  >  1:3111560/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)                                                                                                 |                                                                                              AAGTATTGCGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCC‑ATCTTGTCTTTCATGATTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCTGGTCGATAGAATATAAGTTGAAAGCTATGTCTATCGCTATTACAATTGCCAAT  >  NC_007793/1579240‑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: ‑