# HG changeset patch # User modencode-dcc # Date 1358542526 18000 # Node ID f5e6e4a22981dcff4572ec85e3ae9d07d6df4c41 # Parent 855ad707594f3a31309b9cf92850e04e6f959b6f Deleted selected files diff -r 855ad707594f -r f5e6e4a22981 bamedit.xml --- a/bamedit.xml Fri Jan 18 15:25:48 2013 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,123 +0,0 @@ - - - - SCRIPT_PATH - samtools - - Merging, splitting, filtering, and QC of BAM files - bamedit.py $options_file \$SCRIPT_PATH - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - input_type['input_type_selector'] == 'merge' - - - input_type['input_type_selector'] == 'filter' - - - input_type['input_type_selector'] == 'split' - - - input_type['input_type_selector'] == 'split' - - - input_type['input_type_selector'] == 'pileup' - - - - - <% -import simplejson -%> -#if str($input_type.input_type_selector) == 'merge': -#if str($input_type.input3) != 'None': -#set $__options = { 'action':str("merge"), 'bamout':str($bamoutm), 'input1':str($input1), 'input2':str($input2), 'input3':str($input3) } -#else: -#set $__options = { 'action':str("merge"), 'bamout':str($bamoutm), 'input1':str($input1), 'input2':str($input2) } -#end if -#end if -#if str($input_type.input_type_selector) == 'split': -#set $__options = { 'action':str("split"), 'bamout':str($bamouts1), 'bamout2':str($bamouts2), 'input1':str($input1) } -#end if -#if str($input_type.input_type_selector) == 'filter': -#set $__options = { 'action':str("filter"), 'bamout':str($bamoutf), 'input1':str($input1), 'quality':int($input_type.quality) } -#end if -#if str($input_type.input_type_selector) == 'pileup': -#set $__options = { 'action':str("pileup"), 'bamout':str($bamoutp), 'input1':str($input1), 'input2':str($input2), 'bamname':str($input1.name), 'refname':str($input2.name) } -#end if -${ simplejson.dumps( __options )} - - - - - - - - - -**What it does** - -Simple interface to manipulate BAM files through the use of SAMtools. - -View the original SAMtools documentation: http://samtools.sourceforge.net/ - ------- - -**Usage** - -**Merge BAM files**: Merges 2 or 3 input BAM files into a single BAM file. - -**Split BAM file**: Randomly split input BAM file into 2 BAM files. - -**Filter BAM file**: Removes aligned reads from input BAM files with mapping quality lower than entered Minimum Mapping Quality (MMQ) value. - -**Calculate coverage QC**: Calculates QC metrics such as percentage of genome covered and average coverage of bases covered for the input BAM and reference files. - - ------- - -**Citation** - -Li H.*, Handsaker B.*, Wysoker A., Fennell T., Ruan J., Homer N., Marth G., Abecasis G., Durbin R. and 1000 Genome Project Data Processing Subgroup (2009) The Sequence alignment/map (SAM) format and SAMtools. Bioinformatics, 25, 2078-9. [PMID: 19505943] - -Written by Ziru Zhou ( ziruzhou@gmail.com ). Please send your comments/questions to modENCODE DCC at help@modencode.org. - - - \ No newline at end of file