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%	67.3 / NA	27	H164H (CAT→CAC)	SAUSA300_RS07730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/12);  total (15/12)

AGTATTGCGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCATGATTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCTGGTCGATAGAATATAAGTTGAAAGCTATGTCTATCGCTATTACAATTGCCAATTG  >  NC_007793/1579241‑1579430                                                                                               |                                                                                                ccatgtcgctCACTTCCCCAAAACCTCCTTGACTCGATATAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGtt                                                                                            <  1:3721656/91‑1 (MQ=21)                     aaaCCTCCTTGACTCGATATAAGATGTCTTTACGCGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGt                                                                        <  1:1283397/100‑1 (MQ=39)                          gatctGACTCGATATAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGcc                                                                   <  1:3255413/96‑2 (MQ=255)                           cctTGACTCGATATAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGt                                                                  >  1:3324940/1‑100 (MQ=255)                           cctTGACTCGATATAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCCTTTTGTATGCGt                                                                  >  1:2884837/1‑100 (MQ=255)                               gACTCGATATAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAAt                                                              >  1:4762511/1‑100 (MQ=255)                                aCTCGATATAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATg                                                             <  1:5085023/100‑1 (MQ=255)                                    gATATAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAgtt                                                         <  1:72134/100‑3 (MQ=255)                                     atatAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAgttt                                                        <  1:1001331/97‑4 (MQ=255)                                        tAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGtg                                                     <  1:4975236/100‑2 (MQ=255)                                         aaGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCg                                                    >  1:744690/1‑100 (MQ=255)                                                  ttACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTatat                                           >  1:35024/1‑100 (MQ=255)                                                  ttACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTatat                                           >  1:4279777/1‑100 (MQ=255)                                                  ttACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCCTTTTGTATGCGTCAATGAGTTGGTCGATAGTatat                                           >  1:3470705/1‑100 (MQ=255)                                                  cgtgtggTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCCTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTatat                                           >  1:5298016/7‑100 (MQ=255)                                                   tACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGtatata                                          >  1:2173707/1‑100 (MQ=255)                                                         tATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTgt                                    <  1:1033537/100‑2 (MQ=255)                                                         tATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTgt                                    >  1:412143/1‑99 (MQ=255)                                                              ttttCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCt                               >  1:915041/1‑100 (MQ=255)                                                               tttCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTa                              >  1:1148100/1‑100 (MQ=255)                                                                ttCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTAt                             <  1:3841395/100‑1 (MQ=255)                                                                   cTGCGTCTGCTGTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTAAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTc                          <  1:5210808/100‑1 (MQ=255)                                                                         ctgctgTTCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGc                    >  1:1299253/1‑100 (MQ=255)                                                                                tCCATCTTGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATcaca             <  1:4577422/100‑1 (MQ=255)                                                                                       tGTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCa      >  1:2009892/1‑100 (MQ=255)                                                                                       ggTCTTTCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCa      >  1:3832926/2‑100 (MQ=255)                                                                                           tttCGTGGTTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGTTGGTCGATAGTATATAAGTTGTAAGCTATGTCTATCGCTATCACAATTGCCAAttt  <  1:656603/100‑2 (MQ=255)                                                                                               |                                                                                                AGTATTGCGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCATGATTCCTTTTCATTTTCTTTTTGTATGCGTCAATGAGCTGGTCGATAGAATATAAGTTGAAAGCTATGTCTATCGCTATTACAATTGCCAATTG  >  NC_007793/1579241‑1579430

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