Question: Assembly using cufflinks
gravatar for Bioinfonext
3.8 years ago by
Bioinfonext250 wrote:

I am using this command to do assembly of mapped read on refernece using the hisat2:

/home/yog/software/cufflinks-2.2.1.Linux_x86_64/cufflinks out_mapped.bam -G gene.gtf -p 8 -o /data/SNU_work/Analysis/unmapped/unmapped/hisat2/assembly/

But it did not working. any suggestion?

rna-seq • 1.2k views
ADD COMMENTlink written 3.8 years ago by Bioinfonext250

Move all other options before the bam file. TopHat suite programs used to be sensitive to order of things on command line. Be more specific than it did not working in terms of errors you see (if any). -G is not a valid option for cufflinks (if I recall right, should be -g).

ADD REPLYlink modified 3.8 years ago • written 3.8 years ago by genomax92k

Does it require sorted sam only? or any sorted file format it accept whether it is sam or bam. if it is sam only than please share the samtool command that can extract only mapped read from the alinment sam file.

ADD REPLYlink written 3.8 years ago by Bioinfonext250
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1897 users visited in the last hour