Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,113,329	G→A	G12G (GGC→GGT)	SAUSA300_RS10710 ←	acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,113,329	0	G	A	89.5%	35.3 / ‑5.5	19	G12G (GGC→GGT)	SAUSA300_RS10710	acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (2/15);  total (2/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCTTCCATTACGCCTATTTGGAATCCTGTTTCATAGTTCACTTTC  >  NC_007793/2113232‑2113362                                                                                                  |                                  gttctctTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCTTCCATTACAcc                                 <  1:6299136/98‑1 (MQ=255) gttctctTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACAcc                                 <  1:2781762/98‑1 (MQ=255)      ctTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCCTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAgtt                            <  1:1847414/100‑1 (MQ=255)           cTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGaaa                       <  1:5213454/100‑2 (MQ=255)                cTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTgg                  <  1:3041088/100‑2 (MQ=255)                cTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTgc                  <  1:516544/100‑2 (MQ=255)                cTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTga                  <  1:4125618/100‑2 (MQ=255)                 tgtatCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTGtt                 <  1:1901454/98‑1 (MQ=255)                     tCCCCATTAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTGTTTCAt             <  1:714117/100‑1 (MQ=255)                     tCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTGTTTCAt             <  1:2679683/100‑1 (MQ=255)                     tCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTGTTTCAt             <  1:2123368/100‑1 (MQ=255)                      ccccCATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGAAACCCTGTTTCATa            >  1:4977265/1‑100 (MQ=255)                       cccAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTGTTTCATAg           <  1:2010765/100‑1 (MQ=255)                        ccAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTGTTTCATAGt          >  1:5072436/1‑100 (MQ=255)                         cAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTGTTTCATAGtt         <  1:6355599/100‑1 (MQ=255)                           aTAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACGCCTAGTTGGAGCCCTGTTTCATAGTTCa       <  1:3590348/100‑1 (MQ=255)                           aTAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTGTTTCATAGTTCa       <  1:77492/100‑1 (MQ=255)                             aaGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTGTTTCATAGTTCatg     <  1:1277790/100‑3 (MQ=255)                                cTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACGCCTAGTTGGAACCCTGTTTCATAGTTCattctg  <  1:6447816/100‑6 (MQ=255)                                                                                                  |                                  GCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCTTCCATTACGCCTATTTGGAATCCTGTTTCATAGTTCACTTTC  >  NC_007793/2113232‑2113362

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