Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	374,582	A→G	intergenic (+95/+147)	SAUSA300_RS01705 → / ← SAUSA300_RS01710	lipase/alpha/beta hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	374,582	0	A	G	100.0%	65.6 / NA	25	intergenic (+95/+147)	SAUSA300_RS01705/SAUSA300_RS01710	lipase/alpha/beta hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/3);  total (22/3)

TCTGAATTTAATATGCTATGTAAATCGTGCTGTTATCACGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTAAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTGTATGAATAATATGGTTGATCG  >  NC_007793/374496‑374678                                                                                       |                                                                                                 tctGAATTTAATATGCTATGTAAATCGTGCTGTTATCACGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAg                                                                                     >  1:672968/1‑100 (MQ=255) tctGAATTTAATATGCTATGTAAATCGTGCTGTTATCACGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAg                                                                                     >  1:111873/1‑100 (MQ=255) tctGAATTTAATATGCTATGTAAATCGTGCTGTTATCACGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAg                                                                                     >  1:140620/1‑100 (MQ=255) tctGAATTTAATATGCTATGTAAATCGTGCTGTTATCACGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAg                                                                                     >  1:5683706/1‑100 (MQ=255) tctGAATTTAATATGCTATGTAAATCGTGCTGTTATCACGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAg                                                                                     >  1:2776824/1‑100 (MQ=255) tctGAATTTAATATGCTATGTAAATCGTGCTGTATTCACGGCACATCAGCTATAAGTAGCATCCCAGTGTTGAATTTAAAAATAGTGTTGTGAAATaatg                                                                                     >  1:3009743/1‑98 (MQ=255)         tAATATGCTATGTAAATCGTGCTGTTATCACGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGtc                                                                             <  1:4136870/100‑1 (MQ=255)          aaTATGCTATGTAAATCGTGCTGTTATCACGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGccagatc                                                                            >  1:6197173/1‑95 (MQ=255)               gCTATGTAAATCGTGCTGTTATCACGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGcccagatcggaa                                                                       >  1:5457043/1‑90 (MQ=255)                cTATGTAAATCGTGCTGTTATCACGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAg                                                                      >  1:4815938/1‑100 (MQ=255)                                     cACGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCg                                                 >  1:961446/1‑100 (MQ=255)                                      aCGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGt                                                >  1:4012996/1‑100 (MQ=255)                                        ggCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTAt                                              >  1:5083624/1‑100 (MQ=255)                                              acAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGttt                                        >  1:5026009/2‑100 (MQ=255)                                               cAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTa                                       >  1:5634436/1‑100 (MQ=255)                                                         tAGCATCACAGTGTTGAGTTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTtaat                             >  1:555778/1‑100 (MQ=255)                                                               tccaGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTg                       >  1:5355651/3‑100 (MQ=255)                                                               cacaGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTg                       >  1:1746436/1‑100 (MQ=255)                                                               cacaGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTg                       >  1:1149994/1‑100 (MQ=255)                                                               cacaGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTg                       >  1:744925/1‑100 (MQ=255)                                                               cacaGTGTTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTg                       >  1:1000006/1‑100 (MQ=255)                                                                    tgtTGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTGTATGa                  <  1:4820910/100‑1 (MQ=255)                                                                     cttGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTGTATGaa                 >  1:3607425/2‑100 (MQ=255)                                                                      ttGAATTTAAAAATAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTGTATGaat                >  1:5492752/1‑100 (MQ=255)                                                                                aaaTAGTGAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTGTATGAATAATATGGTTg      <  1:2902479/100‑1 (MQ=255)                                                                                    agtgaagtgaAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTGTATGAATAATATGGCAGATCg  >  1:1633875/1‑100 (MQ=255)                                                                                    agtgaagtgaAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTGTATGAATAATATGGAAGATCg  >  1:1070344/1‑100 (MQ=255)                                                                                       |                                                                                                 TCTGAATTTAATATGCTATGTAAATCGTGCTGTTATCACGGCACATCAGATATAAGTAGCATCACAGTGTTGAATTTAAAAATAGTAAAGTGAAATAAAGCGCCTGTCTCATTAGCGAAAACTAAAGGGACAGGCGTATCTGTTTATGAGCTTAATAAATTGTATGAATAATATGGTTGATCG  >  NC_007793/374496‑374678

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