## RNA-seq / Assemble reads into transcripts using Cufflinks

### Description

This tool assembles reads into transcripts using version Cufflinks 2.1.1.

### Parameters

- Chromosome names in my BAM file look like (chr1, 1) [1]
- Enable multi-mapped read correction (yes,no) [no]
- Correct for sequence-specific bias (yes, no) [no]
- Genome used for bias correction (Human (hg19), Mouse (mm10), Mouse (mm9), Rat (rn4)) [Human (hg19)]

### Details

Given aligned RNA-Seq reads, Cufflinks assembles the alignments into a parsimonious set of transcripts. It then estimates the relative abundances of these transcripts based on how many reads support each one, taking into account biases in library preparation protocols. It is recommended to use the TopHat aligner to map your reads to the reference genome. You can merge the resulting GTFs from several samples to one using the Cuffmerge tool, and use it in differential expression analysis using Cuffdiff.

Cufflinks can detect sequence-specific bias and correct for it in abundance estimation.

By default, Cufflinks will uniformly divide each multi-mapped read to all of the positions it maps to. If multi-mapped read correction is enabled, Cufflinks will re-estimate the transcript abundances dividing each multi-mapped read probabilistically based on the initial abundance estimation, the inferred fragment length and fragment bias, if bias correction is enabled.

### Output

The analysis output consists of the following:

- transcripts.gtf: This GTF file contains Cufflinks' assembled isoforms. The first 7 columns are standard GTF, and the last column contains attributes, some of which are also standardized ("gene_id", and "transcript_id"). There is one GTF record per row, and each record represents either a transcript or an exon within a transcript.
- genes.fpkm_tracking.tsv: This file contains the estimated gene-level expression values in the generic FPKM Tracking Format. Note, however that as there is only one sample, the "q" format is not used.
- isoforms.fpkm_tracking.tsv: This file contains the estimated isoform-level expression values in the generic FPKM Tracking Format. Note, however that as there is only one sample, the "q" format is not used.
- skipped.gtf: This file contains the skipped loci.

### References

This tool uses the Cufflinks package for statistical analysis. Please cite the following article:

Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010 May;28(5):511-5.