PASA pipeline: just compare to the reference annotation

Skip to first unread message

Amina Echchiki

Jul 3, 2017, 4:49:17 AM7/3/17
to pasapipeline-users

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.


Brian Haas

Jul 3, 2017, 7:51:49 AM7/3/17
to Amina Echchiki, pasapipeline-users
Hi Amina,

PASA ends up using gff3 format instead of sam format for uploading the alignments.  If you had the gmap gff3 output, you could just upload that directly using the custom alignment importer, but not with sam.

I'd suggest just having PASA rerun GMAP as part of its regular routine.



You received this message because you are subscribed to the Google Groups "pasapipeline-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to
To post to this group, send email to
To view this discussion on the web visit
For more options, visit

Brian J. Haas
The Broad Institute

Reply all
Reply to author
0 new messages