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%	75.3 / NA	19	coding (770/813 nt)	h16_A3001	transcriptional regulator, LysR‑family
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (9/10);  total (9/10)

CGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTC‑AGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACC  >  AM260479/3247871‑3248129                                                                                                                                        |                                                                                                                             cgcTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACCtg                                                                                                                             <  1:748132/134‑2 (MQ=255)                         tCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTg                                                                                                    >  2:536085/1‑134 (MQ=255)                              cTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTgg                                                                                               >  2:950432/1‑134 (MQ=255)                                       aTCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAAcgt                                                                                      <  2:298641/134‑3 (MQ=255)                                                         aaGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTc‑a                                                                   <  2:586596/133‑1 (MQ=255)                                                                    ggAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCa                                                                                                        <  1:395520/87‑1 (MQ=255)                                                                    ggAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCa                                                                                                        >  2:395520/1‑87 (MQ=255)                                                                     gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCAtc                                                                                                                   >  1:988843/1‑75 (MQ=255)                                                                     gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCAtc                                                                                                                   <  2:988843/75‑1 (MQ=255)                                                                     gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTc‑a                                                                   >  2:203785/1‑121 (MQ=255)                                                                     gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTc‑a                                                                   <  1:203785/121‑1 (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                                                                   <  1:6573/112‑1 (MQ=255)                                                                              cTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTc‑a                                                                   >  2:6573/1‑112 (MQ=255)                                                                                                gCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCAtct                                                                                                                  <  1:937406/49‑1 (MQ=12)                                                                                                gCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCAtct                                                                                                                  >  2:937406/1‑49 (MQ=9)                                                                                                gCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcg                                                 <  1:910637/113‑1 (MQ=255)                                                                                                gCCGGCGTCGCCTATGATG‑CCGAG‑CGCGA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcg                                                 >  2:910637/1‑113 (MQ=255)                                                                                                                             cgaccGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAcc  >  2:884170/4‑134 (MQ=255)                                                                                                                                        |                                                                                                                             CGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTC‑AGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACC  >  AM260479/3247871‑3248129

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