This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to "filter_low_expr_transcripts"?

Hello I denovo assembled my samples and then transcript abundance estimatin with RSEM. Now, I would like to filter low expressed transcripts. Fot this I am feeding TPM matrix expression file as well as the Trinity.fasta file to "filter_low_expr_transcripts.pl". I use the following command as mentioned in the mannual:

tb44227@lido-gw02:~/nobackup/Radula_data/trinity_a2b2_seq_111016_1333/RSEM/RSEM_BWT2>~/nobackup/bioinfo_package/trinityrnaseq-2.2.0/util/filter_low_expr_transcripts.pl  --matrix trans_counts.TPM.not_cross_norm --transcripts Trinity.fasta --min_expr_any 5 --highest_iso_only --trinity_mode

and I got the following error message

Error, --min_pct_iso and --highest_iso_only are mutually exclusive parameters.  at /home/tb44227/nobackup/bioinfo_package/trinityrnaseq-2.2.0/util/filter_low_expr_transcripts.pl line 90.

I tried to google it filter_loq_expression_reads but could not find the way to fix it.

Any suggestion to fix this issue:

Thanks

filter transcripts low expression rna-seq
• 3,667 views • link
written by tcf.hcdg tcf.hcdg

As the error states, the two parameters are mutually exclusive. You can set the threshold for min_pct_iso for all isoforms, or take the highest one only.

• 0 views • link
written by st.ph.n st.ph.n

Thanks. I changed it and It works

• 0 views • link
written by tcf.hcdg tcf.hcdg

0 answers

No answers yet.

Log in to answer this question.

More posts like this