Getting libType without running quantification with Salmon
1
0
Entering edit mode
6 months ago
compuTE ▴ 90

I want to use Salmon just to check the strandness of my samples. Is there a way to get the guess of libtype (--libType A) without getting all the output from the quantification? Or skipping running the quantification all-together?

strandness salmon rnaseq python strand • 356 views
ADD COMMENT
3
Entering edit mode
ADD COMMENT
0
Entering edit mode

Whoops, I guess I missed reading this. Thanks.

ADD REPLY
3
Entering edit mode

Also, it's worth noting that to estimate the strandedness of the library, it's not necessary for salmon to look at all of the reads if you're not attempting to perform quantification. So, you could do something like passing only the first million reads to salmon along with the --skipQuant option, and then you can see what library type it guessed. Since salmon can read from named pipes, you can accomplish this with something like -1 <(gunzip -c fastq_1.fq.gz | head -n 4000000) -2 <(gunzip -c fastq_2.fq.gz | head -n 4000000) to only pass the first million records of the fastq files to salmon.

ADD REPLY

Login before adding your answer.

Traffic: 2883 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6