Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,499,927	A→G	I395V (ATA→GTA)	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,927	0	A	G	100.0%	49.4 / NA	20	I395V (ATA→GTA)	SAUSA300_RS07270	sulfite reductase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/13);  total (7/13)

AAATCATTAATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTA  >  NC_007793/1499838‑1500019                                                                                          |                                                                                             aaaTCATTAATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAGTAAAAGTAcc                                                                                    <  1:2965547/100‑1 (MQ=255)   aTCATTAATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAGTAAAAGTACcac                                                                                  <  1:6422168/100‑1 (MQ=255)          atttatttGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAGTAAAAGTACCACATCGCGt                                                                           <  1:47149/100‑1 (MQ=255)                       tatCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGGTAAAGCGTACAAAAAATACTGCAGTAAAAGTACCACATCGCGTAAAAAAACCAgac                                                              <  1:4753722/100‑3 (MQ=255)                       tatCGCTTAATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAGTAAAAGTACCACATCGCGTAAAAAAACCAGcc                                                              <  1:3618368/100‑2 (MQ=255)                        atCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAGTAAAAGTACCACATCGCGTAAAAAAACCAGCag                                                             <  1:4713806/100‑2 (MQ=255)                          cGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAGAATACTGCAGTAAAAGTACCACATCGCGTAAAAAAACCAGCATCt                                                           <  1:6201208/100‑1 (MQ=255)                          cGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAGTAAAAGTACCACAGCGCGTAAAAAAACCAGCATCt                                                           <  1:3370740/100‑1 (MQ=255)                            cTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAGTAAAAGTACCACATCGCGTAAAAAAACCAGCATCtat                                                         <  1:3401759/100‑1 (MQ=255)                                 gtagtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAGTAAAAGAACCACATCGCGTAAAAAAACCAGCATCTATATCAc                                                    >  1:3737292/2‑100 (MQ=255)                                   agtagCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAGTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTc                                                  >  1:4406356/1‑100 (MQ=255)                                            tGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAGTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGt                                         >  1:5022881/1‑100 (MQ=255)                                                      ccgatctcTTAGATAAAGCGTACAAAAAATACTGCAGTAAAAGTACCACATCGCGT‑AAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTAtt                              <  1:2808940/96‑1 (MQ=255)                                                        tGTCTCTTGGATAAAGCGTACAAAAAATACTGCAGTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTa                             >  1:2278041/1‑100 (MQ=255)                                                           ctctTGGATAAAGCGTACAAAAAATACTGCAGTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGa                          >  1:5364134/1‑100 (MQ=255)                                                               tGGATAAAGCGTACAAAAAATACTGCAGTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGattat                      >  1:59925/1‑100 (MQ=255)                                                                        cGTACAAAAAATACCGCAGTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCat             <  1:4991288/100‑1 (MQ=255)                                                                         gTACAAAAAATACTGCAGTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCatg            <  1:1414596/100‑2 (MQ=255)                                                                          tACAAAAAATACTGCAGTAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAAATGTATTTGGATTATTAATGCCatag           <  1:14287/100‑3 (MQ=255)                                                                                   tACTGCAGTAAAAGAACCACATCGCGTAAAAAAACCAGCATCTATACCACTCATAATATGCTTAATTGTACTAGGATTATTAATGCCATTATTTGGATTa  >  1:6042677/1‑100 (MQ=255)                                                                                          |                                                                                             AAATCATTAATTTATTTGTATGTATCGCTTTATTAGTAGCCATTGGCATGGGGTTTGTCTCTTGGATAAAGCGTACAAAAAATACTGCAATAAAAGTACCACATCGCGTAAAAAAACCAGCATCTATATCACTCATAATATGTTTAATTGTATTAGGATTATTAATGCCATTATTTGGATTA  >  NC_007793/1499838‑1500019

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