Question: how to extract the longest isoform from trinity outfils
0
gravatar for wu.zhiqiang.1020
2.7 years ago by
United States
wu.zhiqiang.102020 wrote:

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

rna-seq assembly • 2.0k views
ADD COMMENTlink modified 4 months ago by Shahzad10 • written 2.7 years ago by wu.zhiqiang.102020

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.
ADD REPLYlink modified 4 months ago by genomax74k • written 4 months ago by Shahzad10

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

ADD REPLYlink written 4 months ago by genomax74k

Or use conda to install the trinity, as it takes full care of dependencies and paths.

https://anaconda.org/bioconda/trinity

ADD REPLYlink written 4 months ago by ATpoint25k

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

ADD REPLYlink written 4 months ago by Shahzad10
1
gravatar for h.mon
2.7 years ago by
h.mon28k
Brazil
h.mon28k wrote:
$TRINITY_HOME/util/misc/get_longest_isoform_seq_per_trinity_gene.pl
ADD COMMENTlink written 2.7 years ago by h.mon28k

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

ADD REPLYlink modified 2.7 years ago • written 2.7 years ago by wu.zhiqiang.102020

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.

ADD REPLYlink written 2.7 years ago by h.mon28k

thanks, I find it and make it go through. thanks again

ADD REPLYlink written 2.7 years ago by wu.zhiqiang.102020
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 873 users visited in the last hour