This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to extract the longest isoform from trinity outfils

Dear all,

I want to filtered the Trinity outfile for RNAseq assemble. I know this is not the best choice, and there are some scripts online. but they are all not good to use in my data, Could you give me some suggestion on this?

my data like this: (Start with > for name)

TRINITY_DN66086_c0_g1_i1 len=1158 path=[1258:0-49 1259:50-66 1268:67-201 1262:202-225 1267:226-656 1257:657-1157] [-1, 1258, 1259, 1268, 1262, 1267, 1257, -2] CTAGCCTACGGGGCTTTATTGGGTCACTTTTTCCAATCGAGTTTGCGGAAAGGCCTCTTT CGGCGGAGCCCTTACGCCTCGTGAAACGAGATATTCAGGGCTCCCCGAAAGGGCACGGCT

TRINITY_DN66086_c0_g1_i2 len=1207 path=[1253:0-49 1254:50-115 1268:116-250 1262:251-274 1267:275-705 1257:706-1206] [-1, 1253, 1254, 1268, 1262, 1267, 1257, -2] CTAGCCTACGGGGCTTTATTGGGTCACTTTTTCCAATCGAGTTTGCGGAAGAGGCTTTTC CACCATTTGGCCCAGTCTAGCCGGGGCCACATGGCGGAAAGGCCTCTTTCGGCGGAGCCC TTACGCCTCGTGAAACGAGATATTCAGGGCTCCCCGAAAGGGCACGGCTGCTCTTCCTTG G

How can I get the longest one for the same gene cluster?

thanks

assembly rna-seq
• 7,826 views • link
updated by Shahzad • Shahzad • written by wu.zhiqiang.1020 • wu.zhiqiang.1020

Dear All

I am using "get_longest_isoform_seq_per_trinity_gene.pl" script

Please help to solve this error. I also have installed the perl module Bio::Index::Fasta but the error is still there. 

Can't locate Fasta_reader.pm in @INC (you may need to install the Fasta_reader module) (@INC contains: /home/ubuntu/Desktop/isoform/../../PerlLib /home/ubuntu/perl5/lib/perl5 /home/ubuntu/miniconda3/envs/longiso/lib/site_perl/5.26.2/x86_64-linux-thread-multi /home/ubuntu/miniconda3/envs/longiso/lib/site_perl/5.26.2 /home/ubuntu/miniconda3/envs/longiso/lib/5.26.2/x86_64-linux-thread-multi /home/ubuntu/miniconda3/envs/longiso/lib/5.26.2 .) at get_longest_isoform_seq_per_trinity_gene.pl line 7.
BEGIN failed--compilation aborted at get_longest_isoform_seq_per_trinity_gene.pl line 7.
• 0 views • link
updated by GenoMax GenoMax • written by Shahzad • Shahzad

Looks like your $PERL5LIB variable is not correctly set. Something like: export PERL5LIB=/path_to/trinity/:${PERL5LIB} would be needed.

• 0 views • link
written by GenoMax GenoMax

I have trinity installed by conda and trinity is working correctly for assembly

• 0 views • link
written by Shahzad • Shahzad

1 answer 

$TRINITY_HOME/util/misc/get_longest_isoform_seq_per_trinity_gene.pl
• 0 views • link
written by h.mon h.mon

It need other pm file to load and the same time. I find this link. But I have to search other more pm files also thanks

• 0 views • link
written by wu.zhiqiang.1020 • wu.zhiqiang.1020

If you install Trinity correctly and set $TRINITY_HOME, the script should work. Also, when reporting an error, please explain it in more detail - for example, copy the error message.

• 0 views • link
written by h.mon h.mon

Log in to answer this question.

More posts like this

No similar posts found.