Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	502,496	G→A	C193Y (TGT→TAT)	SAUSA300_RS02385 →	glutamate synthase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	502,496	0	G	A	100.0%	72.1 / NA	27	C193Y (TGT→TAT)	SAUSA300_RS02385	glutamate synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/4);  total (23/4)

AGCAGGATTAGCTGCTGCTGAAGAACTTAATCTACTAGGATATCAAGTAACTATTTATGAACGTGCTAGAGAATCAGGCGGTTTATTAATGTGTGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGCATTACTTTCATTAATGGTGTTG  >  NC_007793/502404‑502588                                                                                             |                                                                                             agcagGATTAGCTGCTGCTGAAGAACTTAATCTACTAGGATATCAAGTAACTATTTATGAACGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTaag                                                                                       <  1:4436142/100‑3 (MQ=255) agcagGATTAGCTGCTGCTGAAGAACTTAATCTACTAGGATATCAAGTAACTATTTATGAACGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTAtt                                                                                       >  1:1376681/1‑100 (MQ=255)       cttAGCTGCTGCTGAAGAACTTAATCTACTAGGATATCAAGTAACTATTTATGAACGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAat                                                                                 >  1:5173051/2‑100 (MQ=255)          aactgctgctgAAGAACTTAATCTACTAGGATATCAAGTAACTATTTATGAACGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGGATTCCGAATATg                                                                              >  1:1240571/3‑100 (MQ=255)           gctgctgctGAAGAACTTAATCTACTAGGATATCAAGTAACTATTTATGAACGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGa                                                                             >  1:873969/1‑100 (MQ=255)              gctgctgAAGAACTTAATCTACTAGGATATCAAGTAACTATTTATGAACGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAAc                                                                          >  1:2634629/1‑100 (MQ=255)                                  actaGGATATCAAGTAACTATTTATGAACGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGa                                                      >  1:4949031/1‑100 (MQ=255)                                  actaGGATATCAAGTAACTATTTATGAACGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGa                                                      >  1:3886260/1‑100 (MQ=255)                                            cAAGTAACTATTTATGAACGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTa                                            >  1:4561352/1‑100 (MQ=255)                                                         aTGAACGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGAATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGgaaga                               >  1:1714188/1‑100 (MQ=255)                                                             gcGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGaagaagat                           >  1:2784076/2‑98 (MQ=255)                                                             aCGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAGTATGAAACGTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCg                           >  1:3934081/1‑100 (MQ=255)                                                             aCGTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCg                           >  1:4756020/1‑100 (MQ=255)                                                               gTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCggg                         >  1:2070221/1‑100 (MQ=255)                                                               gTGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCgag                         >  1:4833768/1‑98 (MQ=255)                                                               atGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAACGgaagagatcgg                         >  1:1676847/2‑94 (MQ=255)                                                                tGCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGc                        >  1:933506/1‑100 (MQ=255)                                                                 gCTAGAGAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGCa                       >  1:4717383/1‑100 (MQ=255)                                                                    agagAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGCATTa                    >  1:3315899/1‑100 (MQ=255)                                                                    agagAATCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGCATTa                    >  1:2062821/1‑100 (MQ=255)                                                                        gaTCTGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGCATTACttt                <  1:4325310/99‑1 (MQ=255)                                                                        aaTCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGCATTACttt                >  1:1077812/1‑100 (MQ=255)                                                                         ctCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGCATTACTTTc               >  1:272517/2‑100 (MQ=255)                                                                          tCAGGCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGCagatcggaag              >  1:3543648/1‑91 (MQ=255)                                                                              gCGGTTTATTAATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGCATTACTTTCATTaa          >  1:73567/1‑100 (MQ=255)                                                                                      ttaATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGCATTACTTTCATTAATGGTGTTg  <  1:663578/100‑1 (MQ=255)                                                                                      ttaATGTATGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGCATTACTTTCATTAATGGTGTTg  <  1:1012234/100‑1 (MQ=255)                                                                                             |                                                                                             AGCAGGATTAGCTGCTGCTGAAGAACTTAATCTACTAGGATATCAAGTAACTATTTATGAACGTGCTAGAGAATCAGGCGGTTTATTAATGTGTGGTATTCCGAATATGAAACTTGATAAAGATGTGGTTCGACGTCGTATTAAGTTAATGGAAGAAGCGGGCATTACTTTCATTAATGGTGTTG  >  NC_007793/502404‑502588

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: ‑