Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,958,045	C→T	intergenic (‑806/‑77)	SAUSA300_RS09685 ← / → SAUSA300_RS09690	gamma‑hemolysin subunit A/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,958,045	0	C	T	100.0%	78.4 / NA	28	intergenic (‑806/‑77)	SAUSA300_RS09685/SAUSA300_RS09690	gamma‑hemolysin subunit A/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/13);  total (15/13)

TGTTGATGAGCTACATCTTTAAATCTTATTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGCGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATACAGGAACAGACATGGCTAAACTTGATTT  >  NC_007793/1957957‑1958138                                                                                         |                                                                                              tGTTGATGAGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCAtctc                                                                                    <  1:811230/100‑1 (MQ=255)     gatgaGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCt                                                                                <  1:1635120/100‑1 (MQ=255)       tgaGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTag                                                                              <  1:3757817/100‑2 (MQ=255)       tgaGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTAt                                                                              <  1:5188193/100‑1 (MQ=255)        gaGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTAta                                                                             <  1:3694810/100‑2 (MQ=255)        gaGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTAta                                                                             <  1:5671949/100‑2 (MQ=255)        gaGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATc                                                                             <  1:5347069/100‑1 (MQ=255)                   ttAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCttt                                                                  <  1:4032719/100‑1 (MQ=255)                                                ctgtggCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAAt                                     >  1:1003190/6‑100 (MQ=255)                                                cgcgggCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAAt                                     >  1:4375954/6‑100 (MQ=255)                                                 cagggCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATg                                    >  1:5154873/5‑100 (MQ=255)                                                  acggCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGCATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGc                                   >  1:3261891/4‑100 (MQ=255)                                                   ctGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGcc                                  >  1:2114212/2‑100 (MQ=255)                                                     gCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCta                                >  1:13761/1‑100 (MQ=255)                                                     gCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCta                                >  1:437921/1‑100 (MQ=255)                                                     gCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCta                                >  1:5009034/1‑100 (MQ=255)                                                     gCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCta                                >  1:2581018/1‑100 (MQ=255)                                                     gCGAGGAAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCta                                >  1:3249864/1‑100 (MQ=255)                                                       gAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCtata                              >  1:3011908/1‑100 (MQ=255)                                                        aGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATAc                             <  1:887515/100‑1 (MQ=255)                                                            agtaAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATACAGGa                         >  1:4593138/1‑100 (MQ=255)                                                                        ttatACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATACAGGAACAGACATGGCt             <  1:3748971/100‑1 (MQ=255)                                                                        ttatACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGGATTAAAAATTAATGCCTATACAGGAACAGACATGGCt             <  1:1241264/100‑1 (MQ=255)                                                                         tataCACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATACTGGAACAGatcggaaga            >  1:2775971/1‑92 (MQ=255)                                                                         tataCACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATACAGGAACAGACATGGCTa            <  1:3878150/100‑1 (MQ=255)                                                                          ataCACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATACAGGAACAGACATGGCTag           <  1:5783422/100‑2 (MQ=255)                                                                           taCACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATACAGGAACAGACATGGCTaaa          >  1:3894022/1‑100 (MQ=255)                                                                                   tGGTTGTGTATGCATCTCACCTATCATTTTTGCTTTTTTCAAATAAGAAAGCATTAAAAATTAATGCCTATACAGGAACAGACATGGCTAAACTTGAttt  >  1:151747/1‑100 (MQ=255)                                                                                         |                                                                                              TGTTGATGAGCTACATCTTTAAATCTTATTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGCGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATACAGGAACAGACATGGCTAAACTTGATTT  >  NC_007793/1957957‑1958138

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