This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Mapping of De novo assembled transcripts to draft genome

I want to map de novo assembled transcript which is in full length (start codon to stop codon) to the draft genome and want to extract partially mapped or unmapped transcript.

I also want to discard chimeric transcripts.

Please help me to improve the gmap command line:

/home/software/gmap-2016-11-07/src/gmap -D /home/database -d Genome Trinity.fasta --chimera-margin=0 --min-trimmed-coverage=90 --min-identity=95 --failed-input /home/yog/failed.fasta

Will --chimera-margin=0 discard chimeric transcripts?

--min-trimmed-coverage=90 help to put all transcripts in failed.fasta which not able to cover 90% coding sequence of a gene?

--failed-input is these are the transcripts which do not completely map on genome according to above criteria?

**The genome which is I am using here is draft genome ...do not contain all the gene sequences and some of the gene sequences partially present in the draft genome.

Thanks

rna-seq

0 answers

No answers yet.

Log in to answer this question.