low alignment rate of methyl-seq data
0
0
Entering edit mode
23 months ago
Sara ▴ 220

I have 8 samples and paired-end methyl-seq data for all of them. I have a serious problem with alignment of these files. when I aligned as paired-end (R1 and R2 together), the alignment rate is 4-5 % but if I align each file separately (R1 and R2 separately), the alignment rate would be at least 70 %. I have checked everything including the reference file to make sure there is no mistake. I also tried different aligners and different parameters but got similar results. for sure, the problem happened in the lab when processing the samples or during the sequencing. I have no idea in which step the problem happened. do you have any idea what could be the problem?

alignment • 372 views
ADD COMMENT
1
Entering edit mode

You would need some details. Which tools, which command lines. Organism etc.

ADD REPLY
0
Entering edit mode

@ATpoint: I used bismark and data is from Pork. bismark --quiet --bowtie2 --multicore 12 -n 1 /bismark_aln/genomes/ -1 bone.lane7.R1_val_1.fq.gz -2 bone.lane7.R2_val_2.fq.gz

ADD REPLY

Login before adding your answer.

Traffic: 747 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6