changeset 8:2dffdc276fd8

Uploaded mmseq analysis pipeline tool2

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/mmseq.xml	Mon Aug 15 14:10:38 2011 -0400
@@ -0,0 +1,39 @@
+<tool id="mmseq" name="MMSEQ">
+  <description>Estimation of transcript expression</description>
+  <command interpreter="bash">mmseq.sh $input1 $input2 $output1 $output2 $output3</command>
+  <inputs>
+    <param format="fasta" name="input1" type="data" label="Reference genome"/>
+    <param format="bam" name="input2" type="data" label="Sorted bam file"/>
+  </inputs>
+  <outputs>
+    <data format="tabular" name="output1"/>
+    <data format="tabular" name="output2"/>
+    <data format="tabular" name="output3"/>
+  </outputs>
+  <help>
+
+**About MMSEQ**
+
+MMSEQ_ is a novel statistical RNA-seq analysis method for estimating haplotype, isoform and gene specific expression. It deconvolves the mapping of reads to multiple transcripts (isoforms or haplotype-specific isoforms). It can take into account non-uniform read generation and works with paired-end reads. Please cite: Turro, E.; Su, S-Y.; Goncalves, A.; Coin, L.J.M.; Richardson, S. and A., Lewin(2011). Haplotype and isoform specific expression estimation using multi-mapping RNA-seq reads. Genome Biology. 12:R13.
+
+.. _MMSEQ: http://www.bgx.org.uk/software/mmseq.html
+
+--------
+
+**Input formats**
+
+MMSEQ accepts sorted BAM file as input. The SAM files obtained as a result of Bowtie alignment can be converted to BAM files and sorted using the SAMTools.
+
+--------
+
+**Outputs**
+
+MMSEQ generates three output files with expectation maximization (EM) and Gibbs sampling (GS) expression estimates with associated Monte Carlo standard errors (MCSE) tabulated in each one of them (Turro et al., 2011).
+
+The first file provides estimates at the transcript/haplo-isoform level. The second file provides aggregate estimates for sets of transcripts that have been amalgamated due to having identical sequences (and so indistinguishable expression levels). The third file aggregates transcript estimates into genes, thus providing gene level estimates. Homozygous transcripts are aggregated together, while heterozygous transcripts are aggregated separately to produce 'haplo-gene' level estimates.	
+
+Out of the three outputs, the amalgamated estimates are recommended for use as the individual transcript estimates exhibit high variability and anti-correlation, but the total expression of two identical transcripts can be well estimated. 
+
+
+  </help>
+</tool>