#!/usr/bin/perl -w use strict; use warnings; my $expected = <<'END_RESULT'; contig00000001 506 540 7 5.0 7 100 0 70 42 42 0 14 1.45 AAACCCT AAACCCTAAACCCTAAACCCTAAACCCTAAACCCT contig00000001 664 691 3 9.3 3 100 0 56 32 35 32 0 1.58 CAG CAGCAGCAGCAGCAGCAGCAGCAGCAGC contig00000002 128 188 3 20.3 3 80 6 70 34 29 31 4 1.79 GCA GCAGCAGCAGCAATAGCAGCAGCAGCAGATGCAACTGCAGCAGCAGCAGCAGCAGCAGCAG contig00000002 313 357 3 15.0 3 81 9 56 35 33 28 2 1.70 AGC AGCAGCAACAGCACAAGCAGCAGCAGCAGCAGCAGCTGCAGCAGC END_RESULT my @lines = grep { /\S/ && ! /^Parameters/ } ; my $result = ''; my $contig; foreach my $line (@lines) { # skip data untill we have a $contig or a line that can give # us $contig next unless $contig || $line =~ /^Sequence/; if ($line =~ /^Sequence: (contig\d+)/) { $contig = $1; next; } $result .= "$contig $line"; } if ($result eq $expected) { print "Good\n"; } else { print "Bad\n"; } __DATA__ Sequence: contig00000001 Parameters: 2 7 7 80 10 50 50 506 540 7 5.0 7 100 0 70 42 42 0 14 1.45 AAACCCT AAACCCTAAACCCTAAACCCTAAACCCTAAACCCT 664 691 3 9.3 3 100 0 56 32 35 32 0 1.58 CAG CAGCAGCAGCAGCAGCAGCAGCAGCAGC Sequence: contig00000002 Parameters: 2 7 7 80 10 50 50 128 188 3 20.3 3 80 6 70 34 29 31 4 1.79 GCA GCAGCAGCAGCAATAGCAGCAGCAGCAGATGCAACTGCAGCAGCAGCAGCAGCAGCAGCAG 313 357 3 15.0 3 81 9 56 35 33 28 2 1.70 AGC AGCAGCAACAGCACAAGCAGCAGCAGCAGCAGCAGCTGCAGCAGC