Question: Zero reads aligning from bisulfite treated paired end illumina reads
gravatar for ryanzapotocny
4.7 years ago by
ryanzapotocny0 wrote:

I am attempting to use Bismark to align my bisulfite sequencing reads to the reference sequence. In my experiment I used a patch PCR method to enrich for a ~40kb region of the genome so the genome_prep performed in Bismark was done with this sequence. Trimming of the reads was done using cutadapt. I was having trouble removing the adapter sequence in skewer because of SNP's in the adapter so I switched over to cutadapt and chopped of 25 bases at the 5' end and quality trimmed the 3' end. When I then attempt to align these reads in Bismark using the --bowtie2 and --non_directional options I see that zero of around 14k reads are aligning. Even when reducing the quality threshold to its lowest in every category still zero reads align. Has anyone else had this issue before? I am  hoping that this points more towards an issue with the program and not my sequence data since I would expect even a part of a few reads to align at low quality thresholds. Is there any way to tell if my data is junk when it has been bisulfite treated?


sequencing snp alignment next-gen • 1.5k views
ADD COMMENTlink written 4.7 years ago by ryanzapotocny0

This is unlikely to be a bismark problem and likely has more to do with how you prepared the genome. Try aligning against the whole genome and see what happens. It may be that you simply aligned against the wrong region.

ADD REPLYlink written 4.7 years ago by Devon Ryan96k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1495 users visited in the last hour