This is a test version of Biostars. For the public version, visit https://www.biostars.org.
find overlapping transcripts

Hi,

I used stringtie and scallop to assemble transcripts from a few RNAseq samples. Both tools generated annotated and novel transcripts in GTF format. I want to see the overlaps in annotated and novel transcripts between stringtie and scallop. Any recommendation for some handy tools? p.s. have looked at gffcompare, but it seems to work only on overlaps between a reference annotation and your sample gtf.

Thanks!

transcript rnaseq
• 1,511 views • link
updated by Alex Reynolds Alex Reynolds • written by robinycfang • robinycfang

2 answers

bedtools and bedops are the tools most often used to compare interval datatypes:

bedtools: a powerful toolset for genome arithmetic

BEDOPS: the fast, highly scalable and easily-parallelizable genome analysis toolkit

• 0 views • link
written by Istvan Albert Istvan Albert 
$ bedmap --echo --echo-map <(gtf2bed < annotated.gtf) <(gtf2bed < novel.gtf) > answer.bed

• https://bedops.readthedocs.io/en/latest/content/reference/statistics/bedmap.html • https://bedops.readthedocs.io/en/latest/content/reference/file-management/conversion/gtf2bed.html

• 0 views • link
written by Alex Reynolds Alex Reynolds

Log in to answer this question.

More posts like this