Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	3,248,007	Δ1 bp	coding (771/813 nt)	h16_A3001 →	transcriptional regulator, LysR‑family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	3,248,006	0	T	.	100.0%	68.1 / NA	16	coding (770/813 nt)	h16_A3001	transcriptional regulator, LysR‑family
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (6/10);  total (6/10)

AGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTC‑AGCAGGGTGCCTACCTGCGCAGCTT  >  AM260479/3247886‑3248088                                                                                                                         |                                                                                    aGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTcggcc                                                                    <  1:27629/134‑1 (MQ=255) aGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTcggcc                                                                    <  1:94936/134‑1 (MQ=255)                   aGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGgcag                                                   >  1:2694/1‑133 (MQ=255)                           gACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTcggccgg                                                                  <  1:337057/110‑1 (MQ=255)                           gACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTcggccgg                                                                  >  2:337057/1‑110 (MQ=255)                                        tCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTc                            <  1:1093797/134‑1 (MQ=255)                                                      gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTcggcc                                                                    >  1:612166/1‑81 (MQ=255)                                                      gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTcggcc                                                                    >  1:792981/1‑81 (MQ=255)                                                      gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTcggcc                                                                    <  2:612166/81‑1 (MQ=255)                                                      gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTcggcc                                                                    <  2:792981/81‑1 (MQ=255)                                                             ggCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCt                                   <  1:692705/106‑1 (MQ=255)                                                             ggCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCt                                   >  2:692705/1‑106 (MQ=255)                                                             ggCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcgc       <  2:2694/134‑1 (MQ=255)                                                             ggCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcg        >  2:1125456/1‑133 (MQ=255)                                                                  ggCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCtt  <  1:310863/134‑1 (MQ=255)                                                                     gTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcgc       <  2:253968/126‑1 (MQ=255)                                                                                                                         |                                                                                    AGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTC‑AGCAGGGTGCCTACCTGCGCAGCTT  >  AM260479/3247886‑3248088

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