Gene Expression Analysis
cuffdiff [options]* <transcripts.gtf> <sample1_replicate1.sam[,â€¦,sample1_replicateM.sam]> <sample2_replicate1.sam[,â€¦,sample2_replicateM.sam]> â€¦ [sampleN.sam_replicate1.sam[,â€¦,sample2_replicateM.sam]]
Prints the help message and exits
Sets the name of the directory in which Cuffdiff will write all of its output. The default is â€œ./â€.
Specify a label for each sample, which will be included in various output files produced by Cuffdiff.
Use this many threads to align reads. The default is 1.
Instructs Cuffdiff to analyze the provided samples as a time series, rather than testing for differences between all pairs of samples. Samples should be provided in increasing time order at the command line (e.g first time point SAM, second timepoint SAM, etc.)
With this option, Cufflinks counts all fragments, including those not compatible with any reference transcript, towards the number of mapped fragments used in the FPKM denominator. It is inactive by default.
With this option, Cufflinks counts only those fragments compatible with some reference transcript towards the number of mapped fragments used in the FPKM denominator. Using this mode is generally recommended in Cuffdiff to reduce certain types of bias caused by differential amounts of ribosomal reads which can create the impression of falsely differentially expressed genes. It is active by default.
Providing Cufflinks with the multifasta file your reads were mapped to via this option instructs it to run our bias detection and correction algorithm which can significantly improve accuracy of transcript abundance estimates. See How Cufflinks Workshow_it_works/index.html#) for more details.
Tells Cufflinks to do an initial estimation procedure to more accurately weight reads mapping to multiple locations in the genome. See How Cufflinks Works for more details.
The minimum number of alignments in a locus for needed to conduct significance testing on changes in that locus observed between samples. If no testing is performed, changes in the locus are deemed not signficant, and the locusâ€™ observed changes donâ€™t contribute to correction for multiple testing. The default is 10 fragment alignments.
Tells Cuffdiff to ignore all reads that could have come from transcripts in this GTF file. We recommend including any annotated rRNA, mitochondrial transcripts other abundant transcripts you wish to ignore in your analysis in this file. Due to variable efficiency of mRNA enrichment methods and rRNA depletion kits, masking these transcripts often improves the overall robustness of transcript abundance estimates.
The allowed false discovery rate. The default is 0.05.
This is the expected (mean) fragment length. The default is 200bp.
Note: Cuffdiff now learns the fragment length mean for each SAM file, so using this option is no longer recommended with paired-end reads.
The standard deviation for the distribution on fragment lengths. The default is 80bp.
Note: Cuffdiff now learns the fragment length standard deviation for each SAM file, so using this option is no longer recommended with paired-end reads.
Sets the number of iterations allowed during maximum likelihood estimation of abundances. Default: 5000
Print lots of status updates and other diagnostic information.
Suppress messages other than serious warnings and errors.
Turns off the automatic routine that contacts the Cufflinks server to check for a more recent version.
Use the Poisson fragment dispersion model instead of learning one in each condition.
Cuffdiff will output a file for each condition (called _counts.txt) containing the fragment counts, fragment count variances, and fitted variance model. For internal debugging only. This option will be removed in a future version of Cuffdiff.
Cuffdiff will round down to zero the abundance of alternative isoforms quantified at below the specified fraction of the major isoforms. This is done after MLE estimation but before MAP estimation to improve robustness of confidence interval generation and differential expression analysis. The default is 1e-5, and we recommend you not alter this parameter.
Sets the maximum number of fragments a locus may have before being skipped. Skipped loci are marked with status HIDATA. Default: 1000000
Cuffdiff will make this many draws from each transcriptâ€™s predicted negative binomial random numbder generator. Each draw is a number of fragments that will be probabilistically assigned to the transcripts in the transcriptome. Used to estimate the variance-covariance matrix on assigned fragment counts. Default: 100.
For each fragment drawn from a transcript, Cuffdiff will assign it this many times (probabilistically), thus estimating the assignment uncertainty for each transcript. Used to estimate the variance-covariance matrix on assigned fragment counts. Default: 50.
Cuffdiff wonâ€™t test genes for differential regulation unless the conditions in question have at least this many replicates. Default: 3.
Cuffdiff will not employ its â€œeffectiveâ€ length normalization to transcript FPKM.
Cuffdiff will not normalize fragment counts by transcript length at all. Use this option when fragment count is independent of the size of the features being quantified (e.g. for small RNA libraries, where no fragmentation takes place, or 3 prime end sequencing, where sampled RNA fragments are all essentially the same length). Experimental option, use with caution.
Share your experience or ask a question