Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	806,587	T→C	Y207H (TAT→CAT)	SAUSA300_RS03910 →	glycerate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	806,587	0	T	C	100.0%	92.8 / NA	38	Y207H (TAT→CAT)	SAUSA300_RS03910	glycerate kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/4);  total (34/4)

ACACCCTAAAATGGAAACAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGTATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGTTATTTAAAAGTGAATTGAAAATTGCAATTGGTCCAG  >  NC_007793/806490‑806683                                                                                                  |                                                                                                 gcacCCTAAAATGGAAACAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCAt                                                                                                >  1:3809824/2‑100 (MQ=255) acccccTAAAATGGAAACAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCAt                                                                                                >  1:2938090/4‑100 (MQ=255) acacCCTAAAATGGAAACAGCAAGAATTCAAGTATTGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCAt                                                                                                >  1:3619035/1‑100 (MQ=255)   gcccTAAAATGGAAACAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCa                                                                                              >  1:4855963/2‑100 (MQ=255)            tGGAAACAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCa                                                                                     >  1:4263158/1‑100 (MQ=255)               aaaCAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATaa                                                                                  >  1:5991339/1‑100 (MQ=255)                taCAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAAt                                                                                 >  1:1982207/2‑100 (MQ=255)                ctcAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAAt                                                                                 >  1:5867065/3‑100 (MQ=255)                 ccAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATc                                                                                >  1:4569654/2‑100 (MQ=255)                 aCAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATc                                                                                >  1:5156463/1‑100 (MQ=255)                 aCAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATc                                                                                >  1:4340699/1‑100 (MQ=255)                 aCAGCAAGAATTCAAGTAATGTCGGA‑TTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAGTCATAATCa                                                                               >  1:867073/1‑100 (MQ=255)                  cAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCa                                                                               >  1:5189145/1‑100 (MQ=255)                    cctaGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCagg                                                                             >  1:4105221/4‑98 (MQ=255)                    acAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAag                                                                             >  1:2479151/2‑100 (MQ=255)                     cAAGAATTCAAGTAATGTCGGATTTCTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTAAATCATAATCAagc                                                                            >  1:4984168/1‑100 (MQ=255)                     ataGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAagc                                                                            >  1:5849890/3‑100 (MQ=255)                        gAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAagcagc                                                                         <  1:3050578/100‑1 (MQ=255)                        gAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAagcaac                                                                         <  1:1058734/100‑3 (MQ=255)                         gaTTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAgcagca                                                                        >  1:3044870/2‑100 (MQ=255)                         caTTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAGgcagca                                                                        >  1:5358896/2‑100 (MQ=255)                         aaTTCAAGTAATGTCGGATTTGTCAAGTCGATTATATGGGAAGGAAAGTGAAATCATGCAAGCTTATGATGCGCATCAGTTGAATCATAATCAagcagct                                                                        <  1:1690184/100‑2 (MQ=255)                          aTTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGcagcag                                                                       <  1:4936931/100‑1 (MQ=255)                          aTTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGcagcag                                                                       >  1:594968/1‑100 (MQ=255)                            tCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGC‑GCAGatc                                                                    >  1:5697852/1‑98 (MQ=255)                             cAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAagcagccagat                                                                    >  1:621514/1‑95 (MQ=255)                             cAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGCAGCAGatc                                                                    >  1:5758126/1‑98 (MQ=255)                              caGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGCAGCAGatcg                                                                   >  1:345299/2‑97 (MQ=255)                                  aaTGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACCTATGATGCGCATCTGTTGAATCATAATCAAGCAGCAGAAATCgct                                                               >  1:187621/1‑98 (MQ=255)                                               tCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGt                                                  >  1:564018/1‑100 (MQ=255)                                                         tataTGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGa                                        >  1:963544/1‑100 (MQ=255)                                                            atGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGtt                                     >  1:673536/1‑100 (MQ=255)                                                             tagTAAGCAAAGGGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGTTa                                    >  1:4774238/3‑100 (MQ=255)                                                              ggTAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGTTAt                                   >  1:1338215/1‑100 (MQ=255)                                                                tAAGCAAAGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGTTAttt                                 >  1:5825803/1‑100 (MQ=255)                                                                       aGTGAAATCATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGTTATTTAAAAGTg                          >  1:5288081/1‑100 (MQ=255)                                                                               cATGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGTTATTTAAAAGTGAATTGaaa                  >  1:5452817/1‑100 (MQ=255)                                                                                aTGCAAACTTATGATGCGCATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGTTATTTAAAAGTGAATTGaaaa                 >  1:5882538/1‑100 (MQ=255)                                                                                             atgCGCATCCGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGTTATTTAAAAGTGAATTGAAAATTGCAATTGGTcc    <  1:6604503/100‑1 (MQ=255)                                                                                               gcgcATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGTTATTTAAAAGTGAATTGAAAATTGCAATTGGTCCAg  >  1:4793061/1‑100 (MQ=255)                                                                                               gcgcATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGTTATTTAAAAGTGAATTGAAAATTGCAATCGGTCCAg  >  1:5281315/1‑100 (MQ=255)                                                                                               gcgcATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGTTATTTAAAAGTGAATTGAAAATTGCAATCGGCCCAg  >  1:2240785/1‑100 (MQ=255)                                                                                                  |                                                                                                 ACACCCTAAAATGGAAACAGCAAGAATTCAAGTAATGTCGGATTTTTCAAGTCGATTATATGGTAAGCAAAGTGAAATCATGCAAACTTATGATGCGTATCAGTTGAATCATAATCAAGCAGCAGAAATCGATAATTTAATTTGGTATTTTAGTGAGTTATTTAAAAGTGAATTGAAAATTGCAATTGGTCCAG  >  NC_007793/806490‑806683

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