Question: Assembly using cufflinks
0
gravatar for Bioinfonext
2.3 years ago by
Bioinfonext140
Korea
Bioinfonext140 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 • 790 views
ADD COMMENTlink written 2.3 years ago by Bioinfonext140
1

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 2.3 years ago • written 2.3 years ago by genomax65k

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 2.3 years ago by Bioinfonext140
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: 1394 users visited in the last hour