Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	3,248,092	G→C	intergenic (+43/‑175)	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,092	0	G	C	100.0%	89.0 / NA	24	intergenic (+43/‑175)	h16_A3001/hemK1	transcriptional regulator, LysR‑family/methylase of polypeptide chain release factors
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (11/13);  total (11/13)

CCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTC‑AGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCC  >  AM260479/3247991‑3248203                                                                                                       |                                                                                                                gcgcgA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCa                                                                               <  2:869438/133‑1 (MQ=255)  cgcgA‑CCGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCaa                                                                              >  2:939984/1‑134 (MQ=255)     cgaccGCAACC‑TGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAAcc                                                                            <  1:633108/131‑1 (MQ=255)                       gCATCTCGGCCGGCCACCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAg                                                           >  2:188538/1‑134 (MQ=255)                                      aCCTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGaa                                            >  2:309443/1‑134 (MQ=255)                                       ccTGTTTGGCAGCAA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAac                                           <  1:309443/134‑1 (MQ=255)                                               gcagcaA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGcc                                   <  1:939984/134‑1 (MQ=255)                                               gcagcaA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGcc                                   <  2:867099/134‑1 (MQ=255)                                                cagcaA‑CGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCt                                  <  1:936139/134‑1 (MQ=255)                                                     aacGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGcc                                   <  1:931033/127‑1 (MQ=255)                                                     aacGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGcc                                   >  2:931033/3‑129 (MQ=255)                                                       cGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCa                                                                               <  1:496021/83‑1 (MQ=255)                                                       cGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCa                                                                               >  2:496021/1‑83 (MQ=255)                                                                       gtcaagcagGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgc            >  2:579075/5‑134 (MQ=255)                                                                       gtcaagcagGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgc            >  2:397885/5‑134 (MQ=255)                                                                             cagGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCctc      <  1:924503/134‑1 (MQ=255)                                                                             cagGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCctc      <  1:397885/134‑1 (MQ=255)                                                                               gggTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCt                                                                                 >  2:431049/1‑57 (MQ=255)                                                                               gggTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCt                                                                                 <  1:431049/57‑1 (MQ=255)                                                                                 gTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTccc  >  1:172751/1‑134 (MQ=255)                                                                                                gCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgc            >  2:392306/1‑109 (MQ=255)                                                                                                gCTTTGTCTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCccgc            <  1:392306/109‑1 (MQ=255)                                                                                                gCTTTGTCTACACCTTCATCGAGCTGTTCTCACCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCctc      >  1:974326/1‑115 (MQ=255)                                                                                                gCTTTGTCTACACCTTCATCGAGCTGTTCTCACCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCctc      <  2:974326/115‑1 (MQ=255)                                                                                                       |                                                                                                                CCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTC‑AGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCGAGCTGTTCTCGCCGACGCTCAACCGCAAGCTGGTTGAGCAGGCCATGTCCGGCGAACACGAAGCCTACGAACTCTGAGGACGGCCCGCCGCCTCTCCC  >  AM260479/3247991‑3248203

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