-
Notifications
You must be signed in to change notification settings - Fork 99
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Low unique alignment % #676
Comments
Note: One possibility I will investigate is perhaps these two samples had a higher than expected spike-in genome %, accounting for the unmapped reads. |
(As a general comment, if you run the Now for the problem at hand, I agree that all QC profiles you shared look very similar, and they also look good. Some standard trimming should get rid of the unwanted adapter, so it is not obvious why the samples would behave very differently. I have compiled a few FAQs regarding low mapping efficiency here: https://felixkrueger.github.io/Bismark/faq/low_mapping/ Maybe they can set you on the right path? |
Thank you. I had another idea from your FAQ website--these are NEB em-seq libraries. And I forgot to set the Max insert size to 1000. So perhaps those two samples have higher insert sizes and lost more reads due to that. I see also now there is an nf-core pipeline for bismark with an em-seq preset. So I will just switch to that. |
good point. Here are some trimming recommendations for EM-seq (https://felixkrueger.github.io/Bismark/bismark/library_types/#em-seq-neb), and there is preset for the nf-core/methylseq workflow, too (be sure to use the |
Thanks. How do I use the dev version exactly? |
on the command line it is |
Hi, It looks like the dev version is still broken, with at least two major bugs: nf-core/methylseq#406 Any suggestions? Thanks! |
Hi, Seems to be working now, I had to add this to the config: process.stageInMode = 'copy' |
Hi,
2 out of 23 samples, all prepared in the same batch, have low unique alignment % (38% and 45%) relative to other samples ~(75% average).
See attached bismark alignment reports for the 2 problematic samples (M3 and M5) relative to a good sample (M6).
M3.pdf
M5.pdf
M6.pdf
Representative FastQC of the samples does not indicate any abnormality in terms of excessive adapter or overrepresented sequences. I don't see any other reason why these two samples would have low unique alignment %. Do you have a suggestion of how to troubleshoot this?
M3.fastqc.pdf
M4.fastqc.pdf
M5.fastqc.pdf
Thanks
The text was updated successfully, but these errors were encountered: