# HG changeset patch # User eugen # Date 1345038394 14400 # Node ID 8aec1cc3b71fdc3d1c51b0fb34bf59ad197f6627 # Parent 120a27b53be7ed7c4c8e21f929ed20dd1c019e8e Deleted selected files diff -r 120a27b53be7 -r 8aec1cc3b71f bismark.xml --- a/bismark.xml Wed Aug 15 08:53:14 2012 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,175 +0,0 @@ - - - - Bismark - - - - bismark_wrapper.sh - ##Reference genome - ref="${indices.fields.path}" - ##Output files (SAM output, Bismark summary) - mapped=$mapped - summary=$summary - ##Temp directory - tempdir=$mapped.files_path - #if str($singlePaired.sPaired) == "single": - library=single - mate1=$singlePaired.sInput1 - #if $singlePaired.sParams.sSettingsType == "full": - fullparam=true - qual=$singlePaired.sParams.qual - seedmms=$singlePaired.sParams.seedmms - seedlen=$singlePaired.sParams.seedlen - maqerr=$singlePaired.sParams.maqerr - directional=$singlePaired.sParams.non_directional - header=$singlePaired.sParams.sam_no_hd - #end if - #else: - library=paired - mate1=$singlePaired.pInput1 - mate2=$singlePaired.pInput2 - #if $singlePaired.pParams.pSettingsType == "full": - fullparam="true" - qual=$singlePaired.pParams.qual - seedmms=$singlePaired.pParams.seedmms - seedlen=$singlePaired.pParams.seedlen - maqerr=$singlePaired.pParams.maqerr - directional=$singlePaired.pParams.non_directional - header=$singlePaired.pParams.sam_no_hd - minins=$singlePaired.pParams.minins - maxins=$singlePaired.pParams.maxins - #end if - #end if - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -**What it does** - -Bismark is a program to map bisulfite treated sequencing reads to a genome of interest and perform methylation calls in a single step. The output can be easily imported into a genome viewer, such as SeqMonk, and enables a researcher to analyse the methylation levels of their samples straight away. It's main features are: - - - Bisulfite mapping and methylation calling in one single step - - - Supports single-end and paired-end read alignments - - - Supports ungapped and gapped alignments - - - Alignment seed length, number of mismatches etc. are adjustable - - - Output discriminates between cytosine methylation in CpG, CHG and CHH context - -.. _Bismark: http://www.bioinformatics.babraham.ac.uk/projects/bismark/ - -**Input formats** - -Bismark accepts files in Sanger FASTQ format. - -**Outputs** - -The output is in SAM format, and has the following columns:: - - Column Description - -------- -------------------------------------------------------- - 1 QNAME seq-ID - 2 FLAG this flag tries to take the strand a bisulfite read originated from into account (this is different from ordinary DNA alignment flags!) - 3 RNAME chromosome - 4 POS start position - 5 MAPQ always 255 - 6 CIGAR - 7 RNEXT - 8 PNEXT - 9 TLEN - 10 SEQ - 11 QUAL Phred33 scale - 12 NM-tag edit distance to the reference - 13 XX-tag base-by-base mismatches to the reference, not including indels - 14 XM-tag methylation call string - 15 XR-tag read conversion state for the alignment - 16 XG-tag genome conversion state for the alignment - - - -