'=' symbol in the BAM/CRAM SEQ field

OpenGene / gencore

Generate duplex/single consensus reads to reduce sequencing noises and remove duplications

MIT License

116 stars 31 forks source link

Closed ShannonDaddy closed 1 year ago

ShannonDaddy commented 2 years ago

Hi，after I use gencore 0.17.2 to process my bam file, I found '=' symbol in a consensus read sequence as follows: c0770a42d8da6271c3391ce7e453fcc

my command is like this:

please help me out, thnaks!

robinjugas commented 2 years ago

Hi, I have the same issue.

jammywwh commented 1 year ago

same issue someone else can help me?

biolxy commented 1 year ago

I have the same issue. When sequences in Reads1 and Reads2 are inconsistent, the inconsistent site(base) is set to '='