sugar has asked for the wisdom of the Perl Monks concerning the following question:
#!/usr/bin/perl use strict; use warnings; if ( scalar @ARGV < 2 ) { die("Usage:\n\t$0 file1 file2\n\n"); } my %data; my $fn = $ARGV[0]; { my $sequence; open DF, $fn or die $!; while ( my $line = <DF> ) { chomp $line; next if $line =~ m/^\s*$/; if ( substr( $line, 0, 1 ) eq '>' ) { ( $sequence, undef ) = split /\s+/, $line; } else { @{ $data{$sequence} } = split /\s+/, $line; } } close DF; } $fn = $ARGV[1]; { my $sequence; my $side; my $len; open DF, $fn or die $!; while ( my $line = <DF> ) { chomp $line; next if $line =~ m/^\s*$/; if ( substr( $line, 0, 1 ) eq '>' ) { ( $sequence, $side, undef, $len ) = split /(?:_|=|\s)+/, $line; if ( defined @{ $data{$sequence} } ) { my @range = ( 0, $len - 1 ); if ( $side ne "left" ) { @range = ( -1 * $len, -1 ); } print $line, "\n"; print join( ' ', @{ $data{$sequence} } [ $range[0] .. $range[1] ] ), "\n"; } else { warn "No data read from " . $ARGV[0] . " for sequence $sequence\n"; } } else { print join( " ", split //, $line ), "\n"; } } close DF; }
|
|---|
| Replies are listed 'Best First'. | |
|---|---|
|
Re: When the input file is huge !!!
by tilly (Archbishop) on Jan 06, 2009 at 03:53 UTC | |
by Marshall (Canon) on Jan 06, 2009 at 05:28 UTC | |
by tilly (Archbishop) on Jan 06, 2009 at 19:54 UTC | |
by BrowserUk (Patriarch) on Jan 06, 2009 at 21:46 UTC | |
by tilly (Archbishop) on Jan 07, 2009 at 03:01 UTC | |
| |
|
Re: When the input file is huge !!!
by BrowserUk (Patriarch) on Jan 06, 2009 at 03:55 UTC | |
by sugar (Beadle) on Jan 06, 2009 at 04:18 UTC | |
by BrowserUk (Patriarch) on Jan 06, 2009 at 04:37 UTC | |
by sugar (Beadle) on Jan 06, 2009 at 09:32 UTC | |
by sugar (Beadle) on Jan 06, 2009 at 05:21 UTC | |
by BrowserUk (Patriarch) on Jan 06, 2009 at 06:03 UTC |