This is a test version of Biostars. For the public version, visit https://www.biostars.org.
PASA pipeline: just compare to the reference annotation

Hello,

I would like to ask if anyone has experience in running a subset of the PASA pipeline, in particular for the reconciliation of some experimental 'transcripts' with the reference annotation.

In more detail, I am working with RNA-seq data from D. melanogaster. I have reconstructed the 'transcipts' using Trinity. I have aligned these 'transcripts' to the reference genome using GMAP. Now I would like to match these 'experimental transcripts' with the reference annotation, to see how they compare.

I was wondering if I can run just a subset of the whole PASA pipeline, so basically skipping the step of alignment to the reference and using the bam file that I generated externally.

Many thanks in advance.

cross-posted: https://groups.google.com/forum/#!topic/pasapipeline-users/cyG9e8-HjG8

pasa rna-seq transcripts isoforms trinity
• 3,356 views • link
written by dovah • dovah

1 answer

Answer from the authors:

PASA ends up using gff3 format instead of sam format for uploading the alignments. If somebody has the gmap gff3 output, then can be done to just upload that directly using the custom alignment importer, but not with sam. They suggest just having PASA rerun GMAP as part of its regular routine.

• 0 views • link
written by dovah • dovah

Log in to answer this question.

More posts like this