remove duplicates in Nanopore cDNA PCR data
Hi to all. I have a ONT cDNA PCR reads mapping result in bam format. I use picard to remove PCR duplicate reads:
I found that there are too many duplicates. I don't know whether picard is effective for ONT data. Does anyone have any experience in using other tools?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
~80% duplicate pairs in RNAseq data
written by Irvin •
-
I can not see the intron region aligned to ref in IGV from ONT direct RNA sequencing
written by Zhen •
-
Nanopore cDNA sequencing - primary and complementary strand in one read
written by ovovo •
-
How to do preprocessing from fastqc
written by Lulu •
-
Increasing per base G content in QC sequencing files
written by robertsr •
-
Use limma for correlation analysis of RNA-seqdata and continous trait
written by Yibin •
-
How to use ClusterProfiler
written by ochiengwycliffe5 •
-
CollectHsMetrics (Picard) HS_LIBRARY_SIZE is empty
written by Tao •
-
Why my `Percent duplicates` in a shallow sequencing is too high?
written by 631079064 •
-
Picard vs samtools rmdup
written by xd_d
I don't have much experience with ONT cDNA data processing, but is this really that surprising? If you've done a good job in the wet lab steps, you should have full transcripts being sequenced at a reasonable depth. Most transcripts can be spanned by a single ONT read, so you're expecting a lot of duplicated reads. This is especially true for short transcripts and if the system doesn't have that much variation.
That said, I would check the distribution across the transcriptome. If you have coverage across most transcripts this could be real. I'm more alarmed you managed to get 100% mapping rate. No artefacts or unannotated data?
My 100% mapping rate is the bam file filtered by samtools, and then used for Picard.
hello,
I would like to suggest that you could try barcode specific adapter trimming by Porechop. this sometimes reduce duplication caused by technical artifacts.