# HG changeset patch # User stef # Date 1419948573 18000 # Node ID fbe7fe47991f0976b3b0459a2be310b9ea989397 # Parent 831c62364275c753aac4c0af1c08e8161f1b54c0 Uploaded diff -r 831c62364275 -r fbe7fe47991f QDNAseq.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/QDNAseq.xml Tue Dec 30 09:09:33 2014 -0500 @@ -0,0 +1,297 @@ + + + + + samtools + + + + QDNASEQ_PATH + + + + Quantitative copy number abberation detection + + + + QDNAseq.R + $qdnaseq_cfg + \$QDNASEQ_PATH + + + + + + + + + + + + + + + + + + .*\S + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +## Desc: this file was sourced in QDNAseq R wrapper script +## as means to pass all galaxy params to R + +## ----- +## required params +## ----- +TRUE -> inGalaxy +"${binSizes}" -> binSizesString +"${experimentType}" -> experimentType +"${jobName}" -> outputName + +## ----- +## extra main params +## ----- +"${htmlFile}" -> outputHtml +"${htmlFile.id}" -> outputId +"${__new_file_path__}" -> newFilePath + +"${htmlFile.files_path}" -> outputPath +as.logical( "${doSegment}" ) -> doSegment +as.logical( "${debug}" ) -> debug + +## ----- +## own bin-annotations file options +## ----- +"${binannotations_source.binannotation_file}" -> binAnnotations + +## ----- +## advanced options +## ----- +as.double( "${advanced.undoSD}" ) -> undoSD +as.logical( "${advanced.blacklist}" ) -> filterBlacklistedBins +as.integer( "${advanced.mappability}" ) -> mappabilityCutoff +"${advanced.undo_splits}" -> undoSplits +as.logical( "${advanced.copynumbers_igv}" ) -> doOutputCopynumbersIgv + +## #for binSize in $binSizes}.split(",")# +## "${binSize}kbp_${igvCopyNumbers}" -> copyNumbersIgvDatasetFile +## #end for + +## ----- +## plot options +## ----- +as.integer( "${advanced.plot_width}" ) -> PLOT_WIDTH +as.integer( "${advanced.plot_height}" ) -> PLOT_HEIGHT +"${advanced.exclude_chrs}" -> excludeChrsString + +## ----- +## input BAMs init +## ----- +c() -> bamsPaths +c() -> bamsNames + +#for bam in $bams# +c( bamsPaths, "${bam}" ) -> bamsPaths +c( bamsNames, "${bam.name}" ) -> bamsNames +#end for + + + + + + + + + + + + + + + + + + + + + + +.. class:: infomark + +**Introduction** + +This tool is a wrapper for the R Bioconductor package QDNAseq_ + +.. _QDNAseq: http://www.bioconductor.org/packages/release/bioc/html/QDNAseq.html + +It determines the copy number state of human chromosomes 1 - 22 for (shallow coverage) whole genome sequencing data. + +For questions/remarks about the galaxy part of this tool, see contact form here_ + +.. _here: http://www.stefs.nl/wp/contact + +You can **test this tool** with built-in data by selecting the option "Run with test data" and press execute. + +----- + +.. class:: warningmark + +As there is no R 3.1.0 galaxy-package yet (a requirement for QDNAseq), the **dependencies** need to be installed by hand and available to the user under which galaxy runs: R (3.1.0) and bioconductor package QDNAseq (>= 1.0.5). In case the default R is not 3.1.0, also the wrapper xml must be updated to include the correct path during installation of this tool. + +.. class:: warningmark + +The input BAMs are expected to be **single end reads of 50bp length** mapped to GRCh37/hg19 genome build. Other experiment setups are currently not tested or supported within galaxy. See the documentation of QDNAseq at bioconductor on how to deal with different setups (or keep fingers crossed ;) ) + +.. class:: warningmark + +Requires **internet access** for downloading bin-annotations from bitbucket and to show some styling (css) of the final report + +----- + +**Citation** + +For the underlying QDNAseq R package please cite: +Scheinin I, Sie D, Bengtsson H, van de Wiel MA, Olshen AB, van Thuijl HF, van Essen HF, Eijk PP, Rustenburg F, Meijer GA, Reijneveld JC, Wesseling P, Pinkel D, Albertson DG and Ylstra B (2014). “DNA copy number analysis of fresh and formalin-fixed specimens by shallow whole-genome sequencing with identification and exclusion of problematic regions in the genome assembly.” Genome Research. doi:10.1101/gr.175141.114. + +See also the bioconductor package_ documentation. + +.. _package: http://www.bioconductor.org/packages/release/bioc/html/QDNAseq.html + +.. image:: LGG150_copynumber.png +.. image:: LGG150_copynumberSegmented.png + + + +