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%	55.8 / NA	21	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 (9/12);  total (9/12)

CTATGAAAATGTTGATGAGCTACATCTTTAAATCTTATTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGCGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATACAGGAACAGACATGGCTAAACTTGATT  >  NC_007793/1957948‑1958137                                                                                                  |                                                                                             ttgtGAAAATGTTGATGAGGTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTAGACACTTTTGGTTgtgt                                                                                            <  1:2002363/97‑1 (MQ=255) ctgtGAAAATGTTGATGAGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTgtgt                                                                                            <  1:4851761/97‑1 (MQ=255) cTATGAAAATGTTGATGAGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACATTTTTGGTTgtgt                                                                                            <  1:961675/100‑1 (MQ=255) cTATGAAAATGTTGATGAGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTgtgt                                                                                            <  1:809546/100‑1 (MQ=255) cTATGAAAATGTTGATGAGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTgtgt                                                                                            <  1:5608692/100‑1 (MQ=255)  tgtGAAAATGTTGATGAGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTa                                                                                           <  1:5428282/98‑1 (MQ=255)   aTGAAAATGTTGATGAGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTag                                                                                          <  1:80432/100‑2 (MQ=255)     gAAAATGTTGATGAGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGc                                                                                        <  1:1841803/100‑1 (MQ=255)     gAAAATGTTGATGAGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGc                                                                                        <  1:5425413/100‑1 (MQ=255)     gAAAATGTTGATGAGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGc                                                                                        <  1:5158480/100‑1 (MQ=255)     gAAAATGTTGATGAGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATACTGCGAGGTAGTAAATTGATATTAAACACTTTTGGTTGTGTATGc                                                                                        <  1:4704981/100‑1 (MQ=255)              gatgaGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCt                                                                               <  1:5091605/100‑1 (MQ=255)              gatgaGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCt                                                                               <  1:6452079/100‑1 (MQ=255)              gatgaGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCt                                                                               <  1:172521/100‑1 (MQ=255)                tgaGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTAt                                                                             <  1:6029960/100‑1 (MQ=255)                 gaGCTACATCTTTAAATCTTTTTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTagg                                                                            <  1:256389/100‑3 (MQ=255)                                                            ctGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAGTAAGAAAGTATTAAAAATTAATGcc                                 >  1:3431774/2‑100 (MQ=255)                                                             tGCGAGGTAGTAAATTGATATTATGCACTTTTGGTTGTGTATGCATCTCGCCTATCTTTTTTGCTTTTTACAAATAAGAAAGTATTAAAAATTAATGCCt                                <  1:2574706/100‑1 (MQ=255)                                                             cgCGAGGTAGTAAATTGATAGTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCt                                >  1:26060/2‑100 (MQ=255)                                                              gCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCta                               >  1:3852513/1‑100 (MQ=255)                                                              gCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCta                               >  1:4624729/1‑100 (MQ=255)                                                              gCGAGGTAGTAAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCta                               >  1:3017593/1‑100 (MQ=255)                                                              gCGAGGCAGTAAATTGATATTATACACTTTTGCTTGTGTATGCATCTCTCCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCta                               >  1:1243981/1‑100 (MQ=255)                                                                   gtagtaAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGCATTAAAAATTAATGCCTATACAg                          >  1:3563984/1‑100 (MQ=255)                                                                     agtaAATTGATATTATACACTTTTGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATACTGGa                        >  1:6355551/1‑100 (MQ=255)                                                                                           ttGGTTGTGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATACAGGAACAGACATGGCTAAACCTGAtt  >  1:1811657/1‑100 (MQ=255)                                                                                                  |                                                                                             CTATGAAAATGTTGATGAGCTACATCTTTAAATCTTATTATACATTGTATAAAATTATATTGCGAGGTAGTAAATTGATATTATACACTTTTGGTTGCGTATGCATCTCACCTATCTTTTTTGCTTTTTTCAAATAAGAAAGTATTAAAAATTAATGCCTATACAGGAACAGACATGGCTAAACTTGATT  >  NC_007793/1957948‑1958137

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