salmon kallisto sleuth results in a genome browser ?
1
0
Entering edit mode
4.9 years ago
carlmannfr • 0

For differentially-expressed transcripts identified by salmon/kallisto-sleuth, is there was any way to retrieve the alignments so that they could be viewed in a genome brower, or is it necessary to independently map in parallel all the seq reads to the genome in order to get the alignments ?

RNA-Seq • 1.2k views
ADD COMMENT
3
Entering edit mode
4.9 years ago

There is a --pseudobam option to kallisto quant that may produce the output that you want:

Usage: kallisto quant [arguments] FASTQ-files
   ....
   --pseudobam               Output pseudoalignments in SAM format to stdout
ADD COMMENT
3
Entering edit mode

Just for future reference, this can also be done with salmon using the --writeMappings flag. However, it's probably worth noting that the pseudo-SAM files produced by these programs are mapping records of fragments directly to the transcriptome (and before probabilistic assignment). Thus, there will be quite a lot of multi-mapping, and the coordinates are with respect to individual transcripts, not the genome.

ADD REPLY

Login before adding your answer.

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