This mode is used to write annotations corresponding to consensus peaks found in each sequence of count archives. It is implemented by edu.cornell.med.icb.goby.modes.CountArchiveToPeakAnnotationsMode.java.
Mode Parameters
The following options are available in this mode
| Flag | Arguments | Required | Description |
|---|---|---|---|
n/a | input | yes | The basenames of the count archives to read as input. If multiple count archives are provided, consensus peaks are called across the union of the input samples. The mode then outputs counts for these consensus peaks. This results in the same peak definitions being used across all the input samples. |
(-o|--output) | output | yes | The tab delimited output filename. The output will contain called peaks and peak count information. |
(-r|--include-reference-names) | include-reference-names | no | When provided, only write counts for reference identifiers listed in this comma separated list. To process only counts for chromosome 19 and 1, if sequences are identified by 1 and 19, use: --include-reference-names 1,19 |
(-a|--alternative-count-archive) | alternative-count-archive | no | Specify an alternative count archive within the basename. The default is to use the count archive that ends in ".counts" (equivalent to -a counts). Default value: counts |
--threshold | threshold | yes | The peak detection threshold. A peak boundary is called whenever the count at a specific base crosses the threshold. The default value of one (1) indicates that a peak is called for any contiguous segment of a reference sequence where the count is non-zero. Default value: 1 |
Leave a Comment