This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Convert STAR sjdb to gtf

The RNA-seq aligner STAR can use annotation data to generate genomes. These data can be in gtf, gff3 or sjdb format. My question is, can I convert the sjdb file to gtf?

I would like to see what 'transcriptome' annotation was used to generate the genome.

Thanks

star rnaseq splice junctions gtf rna-seq
• 3,994 views • link
updated by Devon Ryan Devon Ryan • written by lkmklsmn lkmklsmn

1 answer

A GTF file has information not present in the sjdb files, namely gene_id and transcript_id columns. The simplest route to proceed might be to just take the first 3 columns of the sjdb file and compare that to the introns in the various possible annotations. Note that the sjdb coordinates are 1-based. I suspect that a simple bedtools intersect -f 0.99 ... | wc -l or something like that would work.

• 0 views • link
written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this