Question: RSEM error for for Transcript quantification
0
gravatar for tcf.hcdg
2.8 years ago by
tcf.hcdg60
European Union
tcf.hcdg60 wrote:

I have paired end Illumina highseq 2000 data from non model organism. I did denovo assembly with Trinity. Trinity assembly generated around 1.5 million transcripts. Then I align the read to the transcriptome with the help of bowtie2 and bam file has been generated. Now I would like to do transcript abundance estimation via alignment based RSEM method.

I follow the instrusction from the Trinity transcript quantification mannual. I tried to prepare the reference for alignemt with the following script and it gives me the error

 "Mapping Info is not correctt, cannot find TRINITY_DN134574_c1_g1_i12's gene_id!"

the script I used for this is:

/home/trinityrnaseq-2.2.0/util/align_and_estimate_abundance.pl --transcripts ~/nobackup/trinity_a2b2_seq_111016_1333/Trinity.fasta  --est_method RSEM --aln_method ~/nobackup/trinity_a2b2_seq_111016_1333/bowtie2.nameSorted.bam  --trinity_mode --thread_count 30 --prep_reference --output_dir ~/nobackup/trinity_a2b2_seq_111016_1333/RSEM/rsem_outdir

Could you please guide me what I did wrong and how can I fix it.

Thanks

ADD COMMENTlink modified 2.8 years ago • written 2.8 years ago by tcf.hcdg60
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: 1056 users visited in the last hour