This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Cufflink on multiple bam files

Im using tophat protocol so now my step one is done I have set of 4 bamfles.Now according to the protocol the next step is to run the cufflinks using the generated bam files , how can I run cufflinks from multiple bam files ,can anyone help me how do I do that ?Im not sure if its going to reduce my execution time or not but it will help me to do other stuff at the same time I suppose .

rna-seq
• 4,004 views • link
updated by Devon Ryan Devon Ryan • written by 1769mkc 1769mkc

Those 6 are from the same sample or multiple samples?

• 0 views • link
written by WouterDeCoster WouterDeCoster

Sorry its 4 bam files , my mistake, so Sample 1 and its replicate which is control and Sample 2 and its replicate which is my test sample. so in total i have 4 bam files

• 0 views • link
written by 1769mkc 1769mkc

1 answer

Run cufflinks separately on each of the 4 BAM files and feed the results into cuffMerge.

• 0 views • link
written by Devon Ryan Devon Ryan

finally i did that ,but is there a way we can out that into a shell script and automate the whole thing ?

• 0 views • link
written by 1769mkc 1769mkc 
for bamfile in *.bam
do
cufflinks $bamfile etc etc etc #do stuff
done
• 0 views • link
written by WouterDeCoster WouterDeCoster

this looks really simple but i have to give output directory into it , how can i specify that in the script all together ? as I want them to go into different output folders

• 0 views • link
written by 1769mkc 1769mkc 
for bamfile in *.bam
do
cufflinks $bamfile -o ${bamfile%.*} etc etc etc #do stuff
done

Note that I don't know if you need the -o flag, just a guess. But ${filename%.*} removes the extension of variable $filename. Check it for yourself using echo first:

for bamfile in *.bam
do
echo cufflinks $bamfile -o ${bamfile%.*} etc etc etc #do stuff
done
• 0 views • link
written by WouterDeCoster WouterDeCoster

thank you very much i will give it a try and let know if I have issues with it.

• 0 views • link
written by 1769mkc 1769mkc

Log in to answer this question.

More posts like this