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
1 answer
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)
Log in to answer this question.
-
How to interpret that overrepresented and duplicate sequences influence %GC content?
written by tvibhaps •
-
MiSeq Fastqc - Drop in the first position
written by andrebolerbarros •
-
Read length for single -ended RNA seq data analysis
written by vinaya •
-
How should I handle the raw reads with failed per base sequence content in fastQC
written by yujinlong000703 •
-
Illumina HiSeq generates abnormally high duplication levels and per tile sequence quality patterns …
written by edikedikedikedik •
-
Pick highest quality fastq reads
written by Biomonika (Noolean)
-
Why does this plot of human GC-content have a peak at around 60%
written by willj •
-
FastQ quality check : what can we correct ?
written by Rad
-
Illumina Rnaseq: Adaptor Sequence Followed By Polya?
written by Biomonika (Noolean)
-
Removing Illumina Adapters From Rna-Seq Data
written by Agatha •
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!!