PASA pipeline: just compare to the reference annotation

105 views
Skip to first unread message

Amina Echchiki

unread,
Jul 3, 2017, 4:49:17 AM7/3/17
to pasapipeline-users
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.

Best,
Amina

Brian Haas

unread,
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.

best,

~brian


--
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 pasapipeline-users+unsub...@googlegroups.com.
To post to this group, send email to pasapipeline-users@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/pasapipeline-users/020e65f0-7709-4da8-8a5a-e605925c576d%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.



--
--
Brian J. Haas
The Broad Institute
http://broadinstitute.org/~bhaas

 
Reply all
Reply to author
Forward
0 new messages