comparison filter_gd_snp.xml @ 22:95a05c1ef5d5

update to devshed revision aaece207bd01

21:d6b961721037

<tool id="gd_filter_gd_snp" name="Filter SNPs" version="1.0.0">
  <description>: Discard some SNPs based on coverage or quality</description>

  <command interpreter="python">
    modify_snp_table.py "$input" "$p1_input" "$output" "$lo_coverage" "$hi_coverage" "$low_ind_cov" "$lo_quality"
    #for $individual, $individual_col in zip($input.dataset.metadata.individual_names, $input.dataset.metadata.individual_columns)
        #set $arg = '%s:%s' % ($individual_col, $individual)
        "$arg"
    #end for
  </command>

  <inputs>
    <param name="input" type="data" format="gd_snp" label="SNP dataset" />
    <param name="p1_input" type="data" format="gd_indivs" label="Population individuals" />
    <param name="lo_coverage" type="integer" min="0" value="0" label="Lower bound on total coverage" />
    <param name="hi_coverage" type="integer" min="0" value="1000" label="Upper bound on total coverage" />
    <param name="low_ind_cov" type="integer" min="0" value="0" label="Lower bound on individual coverage" />
    <param name="lo_quality" type="integer" min="0" value="0" label="Lower bound on individual quality values" />
  </inputs>

  <outputs>

"population", by supplying a list that has been previously created using the
Specify Individuals tool.  SNPs are then discarded if their total coverage
for the population is too low or too high, or if their coverage or quality
score for any individual in the population is too low.

-----

**Example**

- input gd_snp::

22:95a05c1ef5d5

<tool id="gd_filter_gd_snp" name="Filter SNPs" version="1.1.0">
  <description>: Discard some SNPs based on coverage or quality</description>

  <command interpreter="python">
    filter_gd_snp.py "$input" "$p1_input" "$output" "$lo_coverage" "$hi_coverage" "$low_ind_cov" "$lo_quality"
    #for $individual, $individual_col in zip($input.dataset.metadata.individual_names, $input.dataset.metadata.individual_columns)
        #set $arg = '%s:%s' % ($individual_col, $individual)
        "$arg"
    #end for
  </command>

  <inputs>
    <param name="input" type="data" format="gd_snp" label="SNP dataset" />
    <param name="p1_input" type="data" format="gd_indivs" label="Population individuals" />
    <param name="lo_coverage" type="text" value="0" label="Lower bound on total coverage">
      <sanitizer>
        <valid initial="string.digits">
          <!-- &#37; is the percent (%) character -->
          <add value="&#37;" />
        </valid>
      </sanitizer>
    </param>
    <param name="hi_coverage" type="text" value="1000" label="Upper bound on total coverage">
      <sanitizer>
        <valid initial="string.digits">
          <!-- &#37; is the percent (%) character -->
          <add value="&#37;" />
        </valid>
      </sanitizer>
    </param>
    <param name="low_ind_cov" type="integer" min="0" value="0" label="Lower bound on individual coverage" />
    <param name="lo_quality" type="integer" min="0" value="0" label="Lower bound on individual quality values" />
  </inputs>

  <outputs>

"population", by supplying a list that has been previously created using the
Specify Individuals tool.  SNPs are then discarded if their total coverage
for the population is too low or too high, or if their coverage or quality
score for any individual in the population is too low.

The upper and lower bounds on total population coverage can be specified
either as read counts or as percentiles (e.g. "5%", with no decimal places).
For percentile bounds the SNPs are ranked by read count, so for example, a
lower bound of "10%" means that the least-covered 10% of the SNPs will be
discarded, while an upper bound of, say, "80%" will discard all SNPs above
the 80% mark, i.e. the top 20%.  The threshold for the lower bound on
individual coverage can only be specified as a plain read count.

-----

**Example**

- input gd_snp::