Clumpify then Salmon? · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Clumpify then Salmon?

Hi,

I have a dataset that I would like to analyze using Salmon. The data is paired-end Illumina and has a very high duplication rate of 90-95% as determined by Markduplicates in Picard (small genome, oversequenced). I would like to remove duplicates using clumpify before Salmon. I know clumpify will do some sorting to reduce file sizes and I just want to make sure this will not interfere with Salmon, which I know does not want data sorted by coordinates. I'm assuming clumpify is not sorting this way, but just want to make sure I'm not missing anything. Thanks!

salmon rna-seq clumpify

• 942 views • link

updated 2.5 years ago by Ram

• written 2.5 years ago by JMB •

1 answer

No, clumpify is not going to sort based on co-ordinates. It only looks at the sequence. It is an alignment-free dedupe program. Re-ordering of the clumps is false by default. It is only useful for compression when turned on.

• 0 views • link

updated 2.5 years ago by Ram

• written 2.5 years ago by GenoMax

Log in to answer this question.

More posts like this

Why am I getting the exact same MarkDuplicates results when passing a query or coordinate sorted b…
written 1.2 years ago by curious

I am following a [standard protocol][1] for alignment, which mentions: *If a primary alignment is marked as duplicate, then all supplementary alignments for that same …
Differential expression analysis with very high duplication rate
written 2.5 years ago by JMB •

Hi, I am doing a DE analysis on an RNA-seq dataset and have a question about PCR duplicates. The organism is a bacteria with a …
Clarification on sequence duplication levels and Qualimap estimated duplication rate
written 3.0 years ago by keremozdel9 •

Hello, I'm working with illumina paired-end reads. I have merged the data of 2 different runs of the same sample, in order to increase coverage …
Randomizing reads for Salmon
written 3.0 years ago by majd.abdulghani •

Hi Biostars community, I want to run Salmon with aligned samples. I only have sorted BAMs, and the way Salmon works, it needs the reads …
Creating a loop to mark duplicates using Picard
written 3.4 years ago by brandnewatthis •

Hi there, I am trying to create a loop to run multiple sorted bam files through picard to mark duplicates. My bash script is not …
distinguishing PCR and optical duplicates
written 4.1 years ago by evoecogen •

Hello, Mammalian samples were sequenced using Illumina NovaSeq. FastQC reported that 28-35% of reads were duplicates. I wanted to know if the duplication happened at …
salmon light weight alignment mode
written 6.4 years ago by evelyn

Hello, I have sorted bam files for paired end `RNA` seq fastq files using `Hisat2`, `Picard` and I want to use `salmon` alignment mode: hisat2 …
how to set metrics file when using MarkDuplicates in Picard
written 7.8 years ago by Yingzi Zhang

Dear all, I tried to use Picard MarkDuplicates to mark duplicates. I found `METRICS_FILE` need to be set. However, I know nothing about the "metrics"about …
How to interpret the "Matrics" file of MarkDuplicates (Picard)
written 7.9 years ago by ravipatel4 •

Hi, I am using Picard to remove duplicates from BAM files. The manual says that it matters how the BAM files are sorted and you …
Picard vs samtools rmdup
written 9.3 years ago by xd_d

Hey all, I want to remove duplicates from my bam file. I use picard MarkDuplicates to remove the duplicates. (REMOVE_DUPLICATES=true) After I run picard to …