Beefy Boxes and Bandwidth Generously Provided by pair Networks
XP is just a number

Re^2: Removing nucleotide frm sequence

by space_monk (Chaplain)
on Jun 07, 2013 at 07:42 UTC ( #1037609=note: print w/replies, xml ) Need Help??

in reply to Re: Removing nucleotide frm sequence
in thread Removing nucleotide frm sequence

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); } }
If you spot any bugs in my solutions, it's because I've deliberately left them in as an exercise for the reader! :-)

Log In?

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://1037609]
and the web crawler heard nothing...

How do I use this? | Other CB clients
Other Users?
Others romping around the Monastery: (5)
As of 2016-10-23 02:32 GMT
Find Nodes?
    Voting Booth?
    How many different varieties (color, size, etc) of socks do you have in your sock drawer?

    Results (300 votes). Check out past polls.