Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	49,106	A→G	intergenic (‑215/+19)	SAUSA300_RS00205 ← / ← SAUSA300_RS00210	recombinase RecA/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	49,106	0	A	G	100.0%	37.4 / NA	16	intergenic (‑215/+19)	SAUSA300_RS00205/SAUSA300_RS00210	recombinase RecA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (4/12);  total (4/12)

CTCGATAAGTTCATCCATGTTTTTCACTCCGTTATTTGTTCTATCTTCAATACGTCGATTATTCAGTTTGATGCTTCACAGTTGTATGATAAAGACAATTAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCTTCATTAGCGGATTTAAAAAGTTCTTTCGCAGCTTTGT  >  NC_007793/49027‑49198                                                                                |                                                                                             ctcGATAAGTTCATCCATGTTTTTCACTCCGTTATTTGTTCTATCTTCAATACGTCGATTATTCAGTTTGATGCTTCACGGTTGTATGATAAAGACAAtc                                                                          <  1:2605991/100‑2 (MQ=255)   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)                                                                                |                                                                                             CTCGATAAGTTCATCCATGTTTTTCACTCCGTTATTTGTTCTATCTTCAATACGTCGATTATTCAGTTTGATGCTTCACAGTTGTATGATAAAGACAATTAGAAATCTTCGTGAACTCCTGAAGGGCCTATCCCTTCATTAGCGGATTTAAAAAGTTCTTTCGCAGCTTTGT  >  NC_007793/49027‑49198

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