Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,113,320	T→G	E15D (GAA→GAC)	SAUSA300_RS10710 ←	acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,113,320	0	T	G	95.5%	50.6 / ‑5.8	22	E15D (GAA→GAC)	SAUSA300_RS10710	acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (2/19);  total (2/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCTTCCATTACGCCTATTTGGAATCCTGTTTCATAGTTCACTTTC  >  NC_007793/2113226‑2113362                                                                                               |                                           tgTTGCGTTCTCTTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCAt                                       <  1:5273051/99‑1 (MQ=255)    tgcgttctctTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTAc                                    <  1:1546978/95‑1 (MQ=255)      cgttctctTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTAcac                                  <  1:1541887/97‑3 (MQ=255)       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)                                                                                               |                                           CGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCTTCCATTACGCCTATTTGGAATCCTGTTTCATAGTTCACTTTC  >  NC_007793/2113226‑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: ‑