SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   RNA Sequencing (http://seqanswers.com/forums/forumdisplay.php?f=26)
-   -   Cuffdiff not generating the fpkm values (http://seqanswers.com/forums/showthread.php?t=73045)

mavadhani 12-13-2016 09:35 AM

Cuffdiff not generating the fpkm values
 
Hello All,

I am trying to generate a python pipeline to analyze RNASeq data. I am able to build an index of the genome, align the fastq files to the genome and get a sam output and sort the sam file. My script for running cuffdiff works fine on my personal mac. I get the fpkm values of the query files. However, when I run the exact same thing on my windows computer at work, I do get the output files, but the genes.fpkm_tracking file appears to be incomplete. There are no fpkm values instead the fields show 0. I am confused. Any help in this regard is greatly appreciated.

Here is my script:

# Step 6: Compare the sorted sam files using cuffdiff
print("Comparing files")
process = subprocess.Popen([cuffdiff, '-o/cuff_output', gffFile2, sorted_output1, sorted_output2], stdout=subprocess.PIPE)
process.stdout.read()
process.wait()


Thanks,
-M


All times are GMT -8. The time now is 05:56 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.