Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	3,248,064	+A	intergenic (+15/‑203)	h16_A3001 → / → hemK1	transcriptional regulator, LysR‑family/methylase of polypeptide chain release factors

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	3,248,063	1	.	A	100.0%	48.2 / NA	14	intergenic (+14/‑204)	h16_A3001/hemK1	transcriptional regulator, LysR‑family/methylase of polypeptide chain release factors
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (6/8);  total (6/8)

AAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTC‑AGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGG  >  AM260479/3247927‑3248188                                                                                                                                          |                                                                                                                              aaGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTc‑a                                                                                                                              <  2:586596/133‑1 (MQ=255)             gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTc‑a                                                                                                                              <  1:203785/121‑1 (MQ=255)             gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTc‑a                                                                                                                              >  2:203785/1‑121 (MQ=255)               gaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTAcc                                                                                                                <  2:409425/134‑1 (MQ=255)                 tCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTg                                                                                                              <  1:536085/134‑1 (MQ=255)                      cTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTc‑a                                                                                                                              >  2:6573/1‑112 (MQ=255)                      cTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTc‑a                                                                                                                              <  1:6573/112‑1 (MQ=255)                                        gCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcg                                                                                                            >  2:910637/1‑113 (MQ=255)                                        gCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcg                                                                                                            <  1:910637/113‑1 (MQ=255)                                                                     cgaccGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAcc                                                             >  2:884170/4‑134 (MQ=255)                                                                                               gccggccACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTc                                    >  2:453050/1‑134 (MQ=255)                                                                                                       ccTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAac                            <  1:453050/134‑1 (MQ=255)                                                                                                                       cGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTtc                                                                             >  2:130625/1‑70 (MQ=255)                                                                                                                       cGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTtc                                                                             <  1:130625/70‑1 (MQ=255)                                                                                                                           aCGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTg         <  1:884170/134‑1 (MQ=255)                                                                                                                             gCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGcc                    <  1:833474/121‑1 (MQ=255)                                                                                                                             gCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGcc                    >  2:833474/1‑121 (MQ=255)                                                                                                                               aCCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGa     >  1:733461/1‑134 (MQ=255)                                                                                                                                    gCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACgg  <  1:950432/132‑1 (MQ=255)                                                                                                                                          |                                                                                                                              AAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTC‑AGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGG  >  AM260479/3247927‑3248188

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