group|g=s
priority|p=f
use_score|u
- visibility|v=s
+ visibility|V=s
color|c=s
chunk_size|C=s
);
push @options, qw(
stream_in|I=s
stream_out|O=s
- verbose
+ verbose|v
help|?
);
$phastcons_file = Maasha::Config::genome_phastcons( $options->{ "genome" } );
$phastcons_index = Maasha::Config::genome_phastcons_index( $options->{ "genome" } );
- $index = Maasha::UCSC::phastcons_index_retrieve( $phastcons_index );
+ $index = Maasha::UCSC::fixedstep_index_retrieve( $phastcons_index );
$fh_phastcons = Maasha::Common::read_open( $phastcons_file );
while ( $record = get_record( $in ) )
{
if ( $record->{ "CHR" } and $record->{ "CHR_BEG" } and $record->{ "CHR_END" } )
{
- $scores = Maasha::UCSC::phastcons_index_lookup( $index, $fh_phastcons, $record->{ "CHR" }, $record->{ "CHR_BEG" }, $record->{ "CHR_END" }, $options->{ "flank" } );
+ $scores = Maasha::UCSC::fixedstep_index_lookup( $index, $fh_phastcons, $record->{ "CHR" },
+ $record->{ "CHR_BEG" },
+ $record->{ "CHR_END" },
+ $options->{ "flank" } );
push @{ $AoA }, [ @{ $scores } ];
}