small FASTa file from genome guided trinity assembly?

0

Entering edit mode

6.0 years ago

bsly • 0

Hello,

I have run Genome Guided Trinity with RSEM data from RNAseq of my organism and a closely related reference. Despite getting about 13,000 identified genes with RSEM, my FASTa file from ggTrinity only contains about 30 genes of short length. Any ideas of why this might be are greatly appreciated!

Thanks, Belinda

rna-seq assembly genome guided trinity • 947 views

ADD COMMENT • link 6.0 years ago by bsly • 0

0

Entering edit mode

Any code/ commands you can paste from your Trinity run?

Trinity uses fastq and outputs fasta. You then use your fastq reads and count the read alignments with RSEM. So I'm not sure why you are running RSEM data in Trinity...?? Do you mean you used a particular Trinity pipeline rather than the assembly Trinity.pl script?

Do a: grep -c '>' X with X being the name of your fasta file. Report back.

ADD REPLY • link 6.0 years ago by Biogeek ▴ 470

Login before adding your answer.