Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,499,394	C→A	A217D (GCC→GAC)	SAUSA300_RS07270 →	sulfite reductase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,499,394	0	C	A	100.0%	41.6 / NA	16	A217D (GCC→GAC)	SAUSA300_RS07270	sulfite reductase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (14/2);  total (14/2)

ATTCCAGTATTTGTATTAGTCTTAACTGGATTGCCATGGTCTGGTTTTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGCCCTTGGTCAAGGTGAATTAGCGATTAATCCACCTAAATCAGATTTGAACGAATTACCTTGGGCTACACGTAAAAATAAACAGCCAG  >  NC_007793/1499306‑1499480                                                                                         |                                                                                       aTTCCAGTATTTGTATTAGTCTTAACTGGATTGCCATGGTCTGGTTTTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGACCTTGGTCAAg                                                                             >  1:2887279/1‑100 (MQ=255)            tGTATTAGTCTTAATTGGATTGCCATGGTCTGGTTTTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGACCTTGGTCAAGGTGAATTAGCg                                                                  <  1:4989737/100‑1 (MQ=255)                       tAACTGGATTGCCATGGTCTGGTTTTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCAcc                                                       >  1:2141150/1‑100 (MQ=255)                       tAACTGGATTGCCATGGTCTGGTTTTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCAcc                                                       >  1:5954315/1‑100 (MQ=255)                        gaCTGGATTGCCATGGTCTGGTTTTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCACCt                                                      >  1:2950654/2‑100 (MQ=255)                         tcTGGATTGCCATGGTCTGGTTTTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCACCTa                                                     >  1:2186861/2‑100 (MQ=255)                                tGCCATGGTCTGGTTTTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCACCTAAATCAGa                                              >  1:403920/1‑100 (MQ=255)                                 gCCATGGTCTGGTTTTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCACCTAAATCAGAt                                             >  1:4384467/1‑100 (MQ=255)                                 gCCATGGTCTGGTTTTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCACCTAAATCAGAt                                             >  1:5754630/1‑100 (MQ=255)                                        acTGGTTTTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCACCTAAATCAGATTTGAACg                                      >  1:499915/2‑100 (MQ=255)                                         cTGGTTCTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCACCTAAATCAGATTTGAACGa                                     >  1:5082131/1‑100 (MQ=255)                                           ggTTTTATGGGCAGTAAAATTGCCGGTATGATGGCCCCAACCGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCACCTAAATCAGCTTTGAACGAAt                                   >  1:2705584/1‑100 (MQ=255)                                                                      atgatgGACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCACCTAAATCAGATTTGAACGAATTACCTTGGGCTACACGTAAAAATAAAc        >  1:1061944/1‑100 (MQ=255)                                                                         gctgACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCACCTAAATCAGATTTGAACGAATTACCTTGGGCTACACGTAAAAATAAAcagc     >  1:6053987/4‑100 (MQ=255)                                                                           ggACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCACCTAAATCAGATTTGAACGAATTACCTTGGGCTACACGTAAAAATAAAcagcca   >  1:5518174/1‑100 (MQ=255)                                                                            gACACAAACGGTGACCTTGGTCAAGGTGAATTAGCGATTAATCCACCTAAATCAGGTTTGAACGAATTACCTTGGGCTACACGTAAAAATAAAcagccag  <  1:3320784/100‑1 (MQ=255)                                                                                         |                                                                                       ATTCCAGTATTTGTATTAGTCTTAACTGGATTGCCATGGTCTGGTTTTATGGGCAGTAAAATTGCCGGTATGATGGACACAAACGGTGCCCTTGGTCAAGGTGAATTAGCGATTAATCCACCTAAATCAGATTTGAACGAATTACCTTGGGCTACACGTAAAAATAAACAGCCAG  >  NC_007793/1499306‑1499480

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