search

ZW-xjtlu / exomePeak2

Peak calling and differential methylation for MeRIP-Seq
25 stars 5 forks source link

Problem about strand information #15

Open CTLife opened 2 years ago

CTLife commented 2 years ago

I split BAM into two files: plus (forward) strand and minus (reverse) strand.

However, ~50% peaks on minus strand for plus-strand BAM file, and ~50% peaks on plus strand for minus-strand BAM file.

Why? I am sure that the two bam files have no problem.

ZW-xjtlu commented 2 years ago

I split BAM into two files: plus (forward) strand and minus (reverse) strand.

However, ~50% peaks on minus strand for plus-strand BAM file, and ~50% peaks on plus strand for minus-strand BAM file.

Why? I am sure that the two bam files have no problem.

Please check if the parameter strandness in exomePeak2() function is set at either “1st_strand” or “2nd_strand”.