no reads pseudoalign with stranded option in Kallisto
1
0
Entering edit mode
4.1 years ago

I am trying to use Kallisto to get transcript counts. I know my RNA Seq data is stranded. However, when I use either the -fr-stranded or -rf-stranded option in Kallisto, I see that 0 reads get pseudoaligned. However, if I dont choose the strandedness, I get ~ 75% of the reads aligning. These are the commands I have been using

kallisto quant --rf-stranded -i kallistoIdx.idx -o kallisto_aligned_rfStranded R1.gz R2.gz

kallisto quant --fr-stranded -i kallistoIdx.idx -o kallisto_aligned_frStranded R1.gz R2.gz

kallisto quant -i kallistoIdx.idx -o kallisto_aligned R1.gz R2.gz

Am I doing something wrong in the stranded command? I see this is how it is indicated in the kallisto manual. I would really appreciate any help.

I tried with both kallisto/0.43 and the latest version.

RNA-Seq Kallisto • 2.1k views
ADD COMMENT
0
Entering edit mode

Thanks! I realized that as I was editing the question. Works well now!

ADD REPLY
0
Entering edit mode

enter image description here

ADD REPLY
0
Entering edit mode
4.1 years ago

The second file should be R2. You gave it R1 twice.

ADD COMMENT

Login before adding your answer.

Traffic: 2072 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