A program to pull out overlapping features from the database.
get_intersecting_features.pl [--options] <filename>
File options: -i --in <filename> input file -o --out <filename> optionally output file Database options: -d --db <database> database to search: name or sqlite -f --feature <text> db feature to search Modify search range: -b --begin --start <integer> adjust relative search start coordinate -e --end --stop <integer> adjust relative search stop coordinate -p --pos [5 | m | 3] relative position of search coordinate -x --extend <integer> extend search in both directions -r --ref [start | mid] measure distance from which coordinate General options: -z --gz compress output -v --version print version and exit -h --help show extended documentation
The command line flags and descriptions:
- --in <filename>
Specify an input file containing either a list of database features or genomic coordinates for which to collect data. The file should be a tab-delimited text file, one row per feature, with columns representing feature identifiers, attributes, coordinates, and/or data values. The first row should be column headers. Bed files are acceptable, as are text files generated by other BioToolBox scripts. Files may be gzipped compressed.
- --out <filename>
Optionally specify a new filename. A standard tim data text file is written. The default is to rewrite the input file.
- --db <database>
Specify the name of a Bio::DB::SeqFeature::Store annotation database from which gene or feature annotation may be derived. A database is required for generating new data files with features. This option may skipped when using coordinate information from an input file (e.g. BED file), or when using an existing input file with the database indicated in the metadata.
- --feature <text>
Specify the name of the target features to search for in the database that intersect with the list of reference features. The type may be a either a GFF "type" or a "type:method" string. If not specifed, then the database will be queried for potential GFF types and a list presented to the user to select one.
- --start <integer>
- --stop <integer>
- --begin <integer>
- --end <integer>
Optionally specify the relative start and stop positions from the 5' end (default) or the end specified by the "--pos" option with which to restrict the search region for target features. For example, specify "--start=-200 --stop=0" to restrict to the promoter region of genes. Both positions must be specified. Default is to take the entire region of the reference feature.
- --pos [ 5 | m | 3 ]
Indicate the relative position from which to make the adjustments to the search window. Both start and stop adjustments may be made from the respective 5 prime, 3 prime, or middle position as dictated by the feature's strand value.
- --extend <integer>
Optionally specify the number of bp to extend the reference feature's region on each side. Useful when you have small reference regions and you want to include a larger search region.
- --ref [start | mid]
Indicate the reference point from which to calculate the distance between the reference and target features. The same reference point is used for both features. Valid options include "start" (or 5' end for stranded features) and "mid" (for midpoint). Default is "start".
Specify whether the output file should (not) be compressed with gzip.
Print the version number.
Display the POD documentation
This program will take a list of reference features and identify target features which intersect them. The reference features may be either named features (name and type) or genomic regions (chromosome, start, stop). By default, the search region for each reference feature is the entire feature, but may be restricted or expanded in size with appropriate modifiers (--start, --stop, --extend). The target features are specifed as specific types.
Several attributes of the found features are appended to the original input file data. First, the number of target features are reported. If more than one are found, the feature with the most overlap with the reference feature is preferentially listed. The name, type, and strand of the selected target feature is reported. Finally, the distance from the reference feature to the target feature is reported. The reference points for measuring the distance is by default the start or 5' end of the features, or optionally the midpoints. Note that the distance measurement is relative to the coordinates after adjustment with the --start, --stop, and --extend options.
A standard BioToolBox data text file is written.
Timothy J. Parnell, PhD Howard Hughes Medical Institute Dept of Oncological Sciences Huntsman Cancer Institute University of Utah Salt Lake City, UT, 84112
This package is free software; you can redistribute it and/or modify it under the terms of the Artistic License 2.0.