p3-feature-gap

Compute Gap Between Feature Pairs

p3-feature-gap [options] <featurePairFile

This script reads in a list of feature pairs and computes the gap between the features. If the features are on different contigs or belong to different genomes, the gap will be listed as 2 billion. This behavior can be overridden with a command-line option.

Parameters

There are no positional parameters.

The standard input can be specified using Input Options. The input column can be specified using Column Options. The following additional command-line options can be specified.

  • inf

Value to return for an infinite gap (different contigs or genomes). The default is 2000000000.

  • col2

Name or index of the column containing the ID of the second feature. The default is -1, which is the second-to-last column.

Example

p3-echo -t f1.patric_id -t f2.patric_id “fig|1302.21.peg.966” “fig|1302.21.peg.1019” | p3-feature-gap

f1.patric_id f2.patric_id gap fig|1302.21.peg.966 fig|1302.21.peg.1019 55253