How to determine number of contigs used in Trinity assembly?
2
0
Entering edit mode
7.8 years ago

Hello,

How do I go about determining how many contigs were used in the de novo trinity assembly i assembled?

Thanks for the help!

Nikelle

trinity DE assembly rna-seq • 4.8k views
ADD COMMENT
1
Entering edit mode

Run the following command using your fasta assembly file

  grep -c "^>" trinity_fasta_assembly_file
ADD REPLY
1
Entering edit mode
7.8 years ago
st.ph.n ★ 2.7k

There is also a perl script called TrinityStats.pl in the 'util' directory that will give you stats for genes, transcripts, reconstruction size and N metrics. Use as 'perl /path/to/your/install/util/TrinityStats.pl your_assembly.fasta'

ADD COMMENT
0
Entering edit mode
7.8 years ago
Benn 8.3k

There is a perl script in the Trinity util folder called "TrinityStats.pl"

You can use it as follow:

perl /trinityrnaseq-2.2.0/util/TrinityStats.pl  /folder2file/Trinity.fasta

Edit: I see @st.ph.n just commented the same solution! Nevermind...

ADD COMMENT

Login before adding your answer.

Traffic: 2632 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6