Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	49,126	T→C	*596* (TAA→TGA)	SAUSA300_RS00210 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	49,126	0	T	C	100.0%	51.4 / NA	20	*596* (TAA→TGA)	SAUSA300_RS00210	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (4/16);  total (4/16)

CGATAAGTTCATCCATGTTTTTCACTCCGTTATTTGTTCTATCTTCAATACGTCGATTATTCAGTTTGATGCTTCACAGTTGTATGATAAAGACAATTAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCTTCATTAGCGGATTTAAAAAGTTCTTTCGCAGCTTTGTTATCATTTGACGGTGTCCA  >  NC_007793/49029‑49217                                                                                                  |                                                                                            cGATAAGTTCATCCATGTTTTTCACTCCGTTATTTGTTCTATCTTCAGTACGTCGATTATTCAGTTTGATGCTTCACGGTTGTATGATAAAGACAATCAg                                                                                           >  1:270185/1‑100 (MQ=255) cGATAAGTTCATCCATGTTTTTCACTCCGTTATTTGTTCTATCTTCAATACGTCGATTATTCAGTTTGATGCTTCACGGTTGTATGATAAAGACAATCAg                                                                                           >  1:3251203/1‑100 (MQ=255)   aTAAGTTCATCCATGTTTTTCACTCCGTTATTTGTTCTATCTTCAATACGTCGATTATTCAGTTTGATGCTTCACGGTTGTATGATAAAGACAATCAGaa                                                                                         <  1:3221550/100‑1 (MQ=255)               aTGTTTTTCACTCCGTTATTTGTTCTATCTTCAATACGTCGATTATTCAGTTTGATGCTTCACGGTTGTATGATAAAGACAATCAGAAATCTTCGTGaaa                                                                             <  1:617144/100‑2 (MQ=255)                           ccGTTATTTGTTCTATCTTCAATACGTCGATTATTCAGTTTGATGCTTCACGGTTGTATGATAAAGACAATCAGAAATCTTCGTGAACTCCTGAAGGGcc                                                                 <  1:1822939/100‑1 (MQ=255)                            cGTTATTTGTTCTATCTTCAATACGTCGATTATTCAGTTTGATGCTTCACGGTTGTATGATAAAGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCt                                                                <  1:2029531/100‑1 (MQ=255)                             cgatcttGTTCTATCTTCAATACGTCGATTATTCAGTTTGGTGCTTCACGGTTGTATGATAAAGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCtt                                                               <  1:1444628/95‑2 (MQ=255)                                    gTTCTATCTTCAATACGTCGATTATTCAGTTTGATGCTTCACGGTTGTATGATAAAGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCTATccccta                                                        <  1:4268782/100‑4 (MQ=255)                                       cTATCTTCAATACGTCGATTATTCAGTTTGATGCTTCACGGTTGTATGATAAAGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCCTCAtt                                                     <  1:1825114/100‑1 (MQ=255)                                            ccgatctagTCGATTATTCAGTCTGATGCTTCACGGTTGTATGATAAAGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCCTCATTAGCgg                                                <  1:3597846/92‑1 (MQ=255)                                                 tACGTCGATTATTCAGTTTGATGCTTCACGGTTGTATGATAAAGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCCTCATTAGCGGAttaa                                           <  1:3644077/100‑3 (MQ=255)                                                  aCGTCGATTATTCAGTTTGATGCTTCACGGTTGTATGATAAAGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCCTCATTAGCGGATTTaa                                          >  1:5568432/1‑100 (MQ=255)                                                             tCAGTTTGATGCTTCACGGTTGTATGATAAAGACAATCAGAGATCTTCGTGAACTCCTGAAGGGCCTATCCCCTCATTAGCGGATTTAAAAAGTTCTTTc                               >  1:4424057/1‑100 (MQ=255)                                                               aGTTTGATGCTTCACGGTTGTATGATAAAGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCCTCATTAGCGGATTTAAAAAGTTCTTTCGc                             <  1:4013360/100‑1 (MQ=255)                                                                       tcTTCACGGTTGTATGATAAAGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCCTCATTAGCGGATTTAAAAAGTTCTTTCGCAGCTTTGt                     <  1:4790723/99‑1 (MQ=255)                                                                                         aaaGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCCTCATTAGCGGATTTAAAAAGTTCTTTCGCAGCTTTGTTATCATTTGCCGGTGTcc   <  1:435462/100‑1 (MQ=255)                                                                                         aaaGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCCTCATTAGCGGATTTAAAAAGTTCTTTCGCAGCTTTGTTATCATTTGCCGGTGTcc   <  1:1640532/100‑1 (MQ=255)                                                                                          aaGACAATCGGAAATCTTCGTGAACTCCTGAAGGACCTATCCCCTCATTAGCGGATTTAAAAAGTTCTTTCGCAGCTTTGTTATCATTTGCCGGTGTCCa  <  1:1097196/100‑1 (MQ=255)                                                                                          aaGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCCTCATTAGCGGATTTAAAAAGTTCTTTCGCAGCTTTGTTATCATTTGCCGGTGTCCa  <  1:4025709/100‑1 (MQ=255)                                                                                          aaGACAATCAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCCTCATTAGCGGATTTAAAAAGTTCTTTCGCAGCTTTGTTATCATTTGCCGGTGTCCa  <  1:806172/100‑1 (MQ=255)                                                                                                  |                                                                                            CGATAAGTTCATCCATGTTTTTCACTCCGTTATTTGTTCTATCTTCAATACGTCGATTATTCAGTTTGATGCTTCACAGTTGTATGATAAAGACAATTAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCTTCATTAGCGGATTTAAAAAGTTCTTTCGCAGCTTTGTTATCATTTGACGGTGTCCA  >  NC_007793/49029‑49217

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