Assembly using cufflinks
0
0
Entering edit mode
7.3 years ago
Bioinfonext ▴ 460

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.7k views
ADD COMMENT
1
Entering edit mode

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 REPLY
0
Entering edit mode

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 REPLY

Login before adding your answer.

Traffic: 3458 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6