# HG changeset patch # User wolma # Date 1408024669 14400 # Node ID b2bc4e91846254191a141e03f8cc7d3128af6551 Uploaded diff -r 000000000000 -r b2bc4e918462 sam_header.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/sam_header.xml Thu Aug 14 09:57:49 2014 -0400 @@ -0,0 +1,129 @@ + + Create a SAM format header from run metadata for sample annotation. + + mimodd + + + mimodd header + + --rg_id "$rg_id" + --rg_sm "$rg_sm" + + #if $str($rg_cn): + --rg_cn "$rg_cn" + #end if + #if $str($rg_ds): + --rg_ds "$rg_ds" + #end if + #if $str($anno) and $str($month) and $str($day): + --rg_dt "$anno-$month-$day" + #end if + #if $str($rg_lb): + --rg_lb "$rg_lb" + #end if + #if $str($rg_pl): + --rg_pl "$rg_pl" + #end if + #if $str($rg_ds): + --rg_pi "$rg_pi" + #end if + #if $str($rg_pu): + --rg_pu "$rg_pu" + #end if + + --outputfile $outputfile + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +.. class:: infomark + + **What it does** + +This tool takes the user-provided information about a next-generation sequencing run and constructs a valid header in the SAM file format from it. + +The result file can be used by the tools *Convert* and *Reheader* or in the *SNAP Read Alignment* step to add run metadata to sequenced reads files (or to overwrite pre-existing information). + +**Note:** + +**MiModD requires run metadata for every input file at the Alignment step !** + +**Tip:** + +While you can do Alignments from fastq file format by providing a custom header file directly to the *SNAP Read Alignment* tool, the **recommended approach** is to first convert all input files to and archive all datasets in SAM/BAM format with appropriate header information prior to any downstream analysis. Although a bit more time-consuming this practice protects against information loss and ensures that the input datasets will remain useful for others in the future. + + + +