Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
Run Kallisto on Windows 10 Brigobestia Bioinformatics 1 01-21-2021 09:11 AM
counting with kallisto bfp7 RNA Sequencing 2 12-03-2015 01:51 PM
Piping to Picard in Windows ckingsley Bioinformatics 0 09-14-2015 10:45 AM
samtools piping problem malbuquerque Bioinformatics 9 10-29-2013 10:52 AM
how to maintain headers when piping into samtools Kennels Bioinformatics 2 02-26-2013 08:32 AM

Thread Tools
Old 04-09-2019, 12:19 PM   #1
Registered Vendor
Location: San Diego

Join Date: Feb 2012
Posts: 25
Default No bam from kallisto or Salmon after piping into samtools

Hi all,

I am trying to generate the pseudo bam files from both Kallisto v45 and Salmon v13 and can't get either to pipe into Samtools and generate a non-emtpy file. I think I am more currently focused on Salmon since I am getting good mapping results with it, but I've tried with Kallisto also just to see if I can generate the file.

In my script:

salmon quant -i $SalmonIndex -l A -1 $fq1 -2 $fq2 -p 20 -o `pwd` --validateMappings --writeMappings | samtools view -bS - > SalmonBamTest.bam

kallisto quant -i $KallistoIndex -o ./$out --threads=40 --bias --single-overhang $SPECIFICITY --pseudobam $fq1 $fq2 | samtools view -bS - > kallistoPseudobam.bam

Kallisto will run but I get an empty bam file.
Salmon stops after 'Index contained 137,271 targets' and goes to the command prompt, also an empty file I believe

If I remove the pipe and add writeMappings=SalmonBamTest.sam, I generate a sam file that can then be converted with the samtools view command. I realize I can do this in my script and delete the sam, but it's killing me that I can't get the direct pipe into samtools working.
Does anyone have any clue what I am doing wrong?

Thank you!
JeremyDay is offline   Reply With Quote

kallisto, pseudobam, salmon, samtools view

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 04:30 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO