I wonder if there is an existing software studies duplicates from NGS bam file. I understand that Picard and samtools can both mark/remove duplicates, but what if I want to count the frequency of each duplicated fragment?
I could do some string comparison in Python and categorize the duplicates (or its aligned location) into dictionaries, but I am not sure this is optimal. Thank you in advance for reading.
1 answer
You can use clumpify for this if you have original read data available (or you could convert the bam back to fastq using reformat.sh from BBMap suite). Use the option addcount to get information about counts added to fastq headers.
By default clumpify allows two errors/substitutions in the reads when doing "clumping". You can use dupesubs=0 to allow only perfect matches. addcount appends a copies=N message to the fastq header. BTW: I am a user of BBMap suite, not author. That would be @Brian Bushnell. If you just run clumpify.sh on command line you will see the extensive in-line help.
Log in to answer this question.
-
BAI files after removing duplicates and multimappers
written by Maka •
-
Mark Duplicates
written by priya.bmg •
-
Issue with Duplicates reads in Exom Data
written by Ranan Jyoti Sarma
-
10x genomics linked reads-- preprocessing
written by nadiabeg.comsats •
-
Split the bam file by chromosome and speed up picard Markduplicates
written by lghust2011
-
Count Unique Alignments in BAM file
written by colin.kern
-
Question about Details in Picard MarkDuplicates
written by DVA
-
Will Picard MarkDuplicates also un-mark duplicates?
written by Chris Miller
-
Remove duplicates in featurecounts
written by lkmklsmn
-
Samtools Rmdup And Picard Mark Duplicates
written by Kssr
just use
samtools view -f 1024 -c in.bam "chr1:234-567"?Sorry I guess I was not very clear in the post. I want to count the frequency of the duplicated reads, not regions - I do not know what regions to look at yet, unless I scan through the bam file first. Thank you so much for the reply.
so what about sorting the reads by names and counting the number of duplicate for the same name ? (I think this data is provided by picard)Thank you very much for the reply. I'm not sure what do you mean by "name" - I thought every read has a unique name at the very start (starting with @)?
ah yes, sorry, I was wrong.
Maybe dupRadar?