This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to extract de novo transcript from RNA long read alignement ?

Hi,

I have long reads from transcripts sequencing ( one amplicon sequenced wtih PACBIO ) that I mapped to a gene locus using minimap2. I would like to extract transcript structure ( in GTF format ? ) with their abundance. For example, in the screenshot bellow, you can see an alignment showing 2 kind of transcripts. One, with an intronic retention. I would like to get the structure of those transcripts and the amount .

Pacbio alignement

rna pacbio
• 1,574 views • link
updated by Juke34 Juke34 • written by sacha sacha

1 answer

If your minimap2 output is a bam format you can use agat_convert_minimap2_bam2gff.pl from AGAT to convert the data into GFF and then you can extract the sequences using agat_sp_extract_sequences.pl. You may also convert the GFF into GTF using agat_convert_sp_gff2gtf.pl

• 0 views • link
written by Juke34 Juke34

From a GFF file, can I remove duplicate using AGAT ? And get the count of each item ?

• 0 views • link
written by sacha sacha

agat_convert_minimap2_bam2gff.pl will not remove the duplicates but the other scripts (with _sp_ in their name) will remove duplicates automatically (when parsing the file). If you need a close look at the removed duplicates you can run agat_convert_sp_gxf2gxf.pl that will generate a log file.

• 0 views • link
written by Juke34 Juke34

Log in to answer this question.

More posts like this