Re^2: Removing nucleotide frm sequence

[space_monk]

local $/ = "";
while (<>) {
    # should have a whole record
    # split up 4 line record
    my (@line) = split( "\n");
    # change if necessary to check we have a valid record
    if ($line[0] =~ /^@/) {
       $line[1] =~ s/([ACGT]{5})$//;
       $line[3] =~ s/(\w{5})$//;
       print join("\n", @line);
    }
}