in reply to Splitting only on internal pattern, not at start or end of string

Hello,

Since you already know that your fasta string can only contain characters A,T,G,C other than N, the simplest way is to just use that bit of information. :-)

Your regex is fine, but the problem is that it can match 'N' anywhere in the string. Here's how we can use our little tidbit to advantage. Change your regex to: [ATGC]N+[ATGC].

Here's some sample code demonstrating it:

use strict; use warnings; while (my $line = <DATA>) { chomp $line; # The below regex tells perl to look for # any of A,T,G,C followed by a string of # one or more Ns, followed by A,T,G,C. my @info = split /[ATGC]N+[ATGC]/, $line; print join(", ", @info), "\n"; } __DATA__ NNNAAAATATGACAAAGGGGTTNNNNNNNNNNNNNNGATGTCTGGTCAATAGGAT CGCAGCCATTAACATCTCAACAAGCCAAAAATTCCTTCTCAGAAATTCGGNNN AAAATATGACAAAGGGGTTNNNNNNNNNNNNNNGATGTCTGGTCAATAGGAT

Update: My split solution has a problem in that it loses one of [ATGC] on either side of the internal pattern. Please use this solution by hdb or johngg's extractive matching. If you prefer to use lookaround assertions, here's one by oiskuu.

Replies are listed 'Best First'.
Re^2: Splitting only on internal pattern, not at start or end of string
by hdb (Monsignor) on Jan 16, 2014 at 10:10 UTC

    Did you realize that you lose one letter each side of the Ns from your sequence?

      Crap! What was I thinking? Yes, split is not the right solution for this situation.

      OP, apologies - please take johngg's solution. A global match is a better solution than mine.

      Update: added link to solution I was referring to

        If you want to use split, you need to apply lookbehind and lookahead assertions in the regex, to keep the letters A, G, C and T out of the match:

        my @info = split /(?<=[ATGC])N+(?=[ATGC])/, $line;
        (Waiting for AnomalousMonks expert answer...)
Re^2: Splitting only on internal pattern, not at start or end of string
by BiologySwede (Initiate) on Jan 16, 2014 at 09:51 UTC
    That is totally awesome, many thanks!

      Unfortunately, also totally wrong as it will consume an A, C, G or T adjacent to the Ns. Rather than split do a global match for one or more of A, C, G or T.

      $ perl -E ' $seq = q{NNACGTNNNACGTNACGTNN}; say for split m{[ACGT]N+[ACGT]}, $seq; say q{-} x 10; say for $seq =~ m{[ACGT]+}g;' NNACG CG CGTNN ---------- ACGT ACGT ACGT $

      I hope this is helpful.

      Update: Corrected wording, s/more than one of/one or more of/

      Cheers,

      JohnGG