# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> DESCRIPTION <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
-# List all tracks availible in a local installation of the Biopieces Genome Browser.
+# Intersects two tracks in the Biopieces Genome Browser.
# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> DESCRIPTION <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
-# Write Biopiece records to the KISS Genome browser.
+# Write Biopiece records to the Biopieces Genome browser.
# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
$new_record->{ 'S_ID' } = $s_id;
$new_record->{ 'STRAND' } = $strand;
+ $new_record->{ 'S_BEG' } = $beg;
+ $new_record->{ 'S_END' } = $end;
$new_record->{ 'CONTIG_ID' } = $contig;
- $new_record->{ 'CONTIG_BEG' } = $beg;
- $new_record->{ 'CONTIG_END' } = $end;
$new_record->{ 'CONTIG_LEN' } = $end - $beg + 1;
$new_record->{ 'CONTIG_MIN' } = $min;
$new_record->{ 'CONTIG_MAX' } = $max;
Maasha::Common::run( "formatdb", "$arg -l /dev/null" );
}
-Maasha::Filesys::dir_create( $options->{ 'directory' } );
+Maasha::Filesys::dir_create_if_not_exists( $options->{ 'directory' } );
unlink $file_tmp;
# Returns a hashref
+ if ( not defined $record->{ 'S_ID' } and
+ not defined $record->{ 'S_BEG' } and
+ not defined $record->{ 'S_END' } )
+ {
+ return undef;
+ }
+
$record->{ 'SCORE' } ||= $record->{ 'E_VAL' } || ".";
$record->{ 'HITS' } ||= ".";
$record->{ 'BLOCK_COUNT' } ||= ".";