Question: what is the input for cuffmerge???
3
gravatar for F
3.6 years ago by
F3.4k
Iran
F3.4k wrote:

hi friends,

i RNA-seq protocol i read...

tophat -p 8 -G genes.gtf -o C1_R1_thout genome C1_R1_1.fq C1_R1_2.fq

by which i have accepted_hits.bam

cufflinks -p 8 -o C1_R1_clout C1_R1_thout/accepted_hits.bam

then saying that :

Create a file called assemblies.txt that lists the assembly file for each sample. The file should contain the following lines:
./C1_R1_clout/transcripts.gtf
./C2_R2_clout/transcripts.gtf
./C1_R2_clout/transcripts.gtf
./C2_R1_clout/transcripts.gtf
./C1_R3_clout/transcripts.gtf
./C2_R3_clout/transcripts.gtf

then this command

cuffmerge -g genes.gtf -s genome.fa -p 8 assemblies.txt

what is the assemblies.txt please??? how i can create such a file by cufflinks output that as i saw is a file named QC_filtered.fastq

thanks for your simpathy

 

myposts cuffmerge cufflinks • 5.7k views
ADD COMMENTlink modified 3.6 years ago by geek_y9.4k • written 3.6 years ago by F3.4k
6
gravatar for geek_y
3.6 years ago by
geek_y9.4k
Barcelona/CRG/London/Imperial
geek_y9.4k wrote:

The output of cufflinks is transcripts.gtf . If you have many samples, cufflinks creates transcripts.gtf ( i.e assembles transcripts) file for each samples and these needs to be merged by cuffmerge,

ADD COMMENTlink modified 3.6 years ago • written 3.6 years ago by geek_y9.4k

then how i can create the assemblies.txt file from transcripts.gtf  files?

ADD REPLYlink written 3.6 years ago by F3.4k
3

Open a text editor. copy paste the path to transcripts.gft of each sample per line. Then save the file as assemblies.txt. Or you could also write a short command that does it for you.

ADD REPLYlink written 3.6 years ago by geek_y9.4k

thank you Goutham

ADD REPLYlink written 3.6 years ago by F3.4k

thank you for your response, would you please write a command line for cuffmerge and explain each item separately thanks

ADD REPLYlink written 14 months ago by hamid_fiuji0

Sorry, I don't have access to Linux or mac I then can't remember what I had done 2 years ago but the instruction of cuffmerge was somehow straightforward as I remember.

ADD REPLYlink written 14 months ago by F3.4k
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: 1399 users visited in the last hour