Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,113,266	G→A	D33D (GAC→GAT)	SAUSA300_RS10710 ←	acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,113,266	0	G	A	100.0%	91.4 / NA	33	D33D (GAC→GAT)	SAUSA300_RS10710	acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (3/30);  total (3/30)

CCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCTTCCATTACGCCTATTTGGAATCCTGTTTCATAGTTCACTTTC  >  NC_007793/2113217‑2113362                                                  |                                                                                                 ctcgCTCTTTGTTGCGTTCTCTTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCt                                                <  1:3020436/97‑1 (MQ=255) ccAGCTCTTTGTTGCGTTCTCTTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCt                                                >  1:4482760/1‑100 (MQ=255) ccAGCTCTTTGTTGCGTTCTCTTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCt                                                <  1:2713815/100‑1 (MQ=255) ccAGCTCTTTGTTGCGTTCTCTTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCt                                                <  1:4000838/100‑1 (MQ=255) ccAGCTCTTTGTTGCGTTCTCTTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTGGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCt                                                <  1:688685/100‑1 (MQ=255) ccAGCACTTTGTTGCGTTCTCTTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCt                                                <  1:5690221/100‑1 (MQ=255)   aGCTCTTTGTTGCGTTCTCTTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAg                                              <  1:1566697/100‑1 (MQ=255)    gCTCTTTGTTGCGTTCTCTTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAgt                                             <  1:900139/100‑2 (MQ=255)     ctctTTGTTGCGTTCTCTTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAgag                                            <  1:4369498/100‑3 (MQ=255)     ctctTTGTTGCGTTCTCTTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGcg                                            <  1:5327704/100‑2 (MQ=255)     ctctTTGTTGCGTTCTCTTAACTTACCTATATCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGcg                                            <  1:2272290/100‑2 (MQ=255)          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:2123368/100‑1 (MQ=255)                                    tCCCCAATAAGCTCATCTCGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCGTCCATTACACCTAGTTGGAACCCTGTTTCAt             <  1:2679683/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)                                                  |                                                                                                 CCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTCACGTTGCTTCTTGTACTCATCACGTTGTTTTCTCATCTTCTTCAACCTAGCTTCCATTACGCCTATTTGGAATCCTGTTTCATAGTTCACTTTC  >  NC_007793/2113217‑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: ‑