|
|
Way to extract detail and its content above specific value problem asking
03-16-2010
>position_10 sample:68711 coords:5453-8666 number:3 type:complete len:344 MSINQYSSDFHYHSLMWQQQQQQQQHQNDVVEEKEALFEKPLTPSDVGKLNRLVIPKQHA ERYFPLAAAAADAVEKGLLLCFEDEEGKPWRFRYSYWNSSQSYVLTKGWSRYVKEKHLDA NRTS* >position_4 sample:68711 coords:553-866 number:4 type:partial len:483 MSGVVRSSPGSSQPPPPPPHHPPSSPVPVTSTPVIPPIRRHLAFASTKPPFHPSDDYHRF KITPSDVENDESDYWLLSNAEISMTDIWKTDSGIDWDYGIADVSTPPPGMGEIAPTAVDS TPR* >position_7 sample:68711 coords:453-86 number:2 type:partial len:214 KAAETLEVQKRRIYDITNVLEGIDLIEKPFKNRILWKGVDACPGDEDADVSVLQLQAEIE NLALEEQALDNQIRWLFVTEEDIKSLPGFQNQTLIAVKAPHGTTLEVPDPDEAADHPQRR TDSGIDWDYGIADVSTPPPGMGEIAPTAVDSTPR* >position_11 sample:68711 coords:53-86 number:1 type:complete len:558 MLGDFIIRLLVLILGYTYPAFECFKTVEKNKVDIEELRFWCQYWILLALISSFERVGDFF RAPRPLNKSLSALRSLEKQTSRGRKWPPPTPPPTPGRDSAGTFNGDDGVNIPDTIPGSPL TDARAKLRRSNSRTQPAA* . .
>position_10 sample:68711 coords:5453-8666 number:3 type:complete len:344 MSINQYSSDFHYHSLMWQQQQQQQQHQNDVVEEKEALFEKPLTPSDVGKLNRLVIPKQHA ERYFPLAAAAADAVEKGLLLCFEDEEGKPWRFRYSYWNSSQSYVLTKGWSRYVKEKHLDA NRTS* >position_11 sample:68711 coords:53-86 number:1 type:complete len:558 MLGDFIIRLLVLILGYTYPAFECFKTVEKNKVDIEELRFWCQYWILLALISSFERVGDFF RAPRPLNKSLSALRSLEKQTSRGRKWPPPTPPPTPGRDSAGTFNGDDGVNIPDTIPGSPL TDARAKLRRSNSRTQPAA* . .
03-16-2010
solution
use strict;
use warnings;
open FH , "testfile" or die "$0:$!" ;
while (<FH> )
{
my $line = $_ ;
# check for a type complete
if ( $line =~ />type complete/ )
{
$line =~ /([0-9]+)\s*$/ ;
# check for greater than 300
if ( $1 < 300)
{
next ;
}
print "$line" ;
while ( <FH> )
{
# leave the type partial
if ( $_ !~ />type partial/ )
{
print "$_" ;
}
else{
last ;
}
}
}
}
03-16-2010
|
|
|
|