Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	3,248,010	C→G	A258A (GCC→GCG)	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,010	0	C	G	100.0%	65.4 / NA	19	A258A (GCC→GCG)	h16_A3001	transcriptional regulator, LysR‑family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (9/10);  total (9/10)

TCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTC‑AGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACC  >  AM260479/3247895‑3248129                                                                                                                    |                                                                                                                         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)                                                                                                                    |                                                                                                                         TCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTC‑AGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACC  >  AM260479/3247895‑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: ‑