Bismark reads failing to align
1
0
Entering edit mode
7.9 years ago
Anna S ▴ 500

Hello,

I'm using Bismark for the first time, and 99.99% of the reads are failing to align for my paired dataset of seq length 51 bp. I have tried using the default -l, as well as -l of 40, the default -n as well as -n of 1, the original fastq files as well as quality trimmed and adaptor cut ones, etc., but none of these variations have made much of a difference. The FastQC output shows that the per base and per sequence quality scores are very good. Do you have any suggestions on how to proceed?

Thank you!!

Anna

bismark fastqc • 2.3k views
ADD COMMENT
0
Entering edit mode

Please delete this question as it turns out to have been my misunderstanding of the datasets (2 lanes of the same end files vs paired-end files). Bismark is working great now!!

ADD REPLY
1
Entering edit mode
7.9 years ago

Without more detail is hard to tell where the problem is... Assuming the reference genome is the correct one, I would try to align only read 1 and see if the mapping goes up. Also make sure the paired reads have been trimmed in such way that the two fastq files in a pair remain in sync (i.e. by using trim_galore)

ADD COMMENT

Login before adding your answer.

Traffic: 1246 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