Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • codecatcher
    Junior Member
    • Jan 2015
    • 2

    Help with a simple script in PERL

    sounds like homework
    Last edited by codecatcher; 01-23-2015, 03:12 PM.
  • Richard Finney
    Senior Member
    • Feb 2009
    • 701

    #2
    What does your first draft look like?

    Comment

    • Phoenix_ICE
      Member
      • Jan 2011
      • 14

      #3
      It sounds easy. But I really have no idea about the score! How to calculate the score?
      Which one is the most important, and which one is the second, of your 3 limitations
      Last edited by Phoenix_ICE; 01-22-2015, 12:54 PM.

      Comment

      • Bukowski
        Senior Member
        • Jan 2010
        • 388

        #4
        Perl is not an acronym and therefore never capitalised. Show some working, or try your luck on stackoverflow.com.

        Comment

        • rdeborja
          Member
          • Aug 2008
          • 42

          #5
          A quick starting point

          Here's a starting point for you. I copied the 23mer you provided 3 times with a slight modification to the first occurrence.

          === script.pl ======
          Code:
          #!/usr/bin/env perl
          use Data::Dumper;
          
          my $sequence = 'ACGATCTTTGCCCCGACGTGATCGAGGTTTTTTTTTTTTTTCAGAGACCGAGATACGATCCCCCGACGTGATCGAGGACGATCCCCCGACGTGATCGAGG';
          my @sequences = $sequence =~ /.{21}GG/g;
          print Dumper(\@sequences);
          
          __END__
          === output =======
          $ perl script.pl
          $VAR1 = [
          'TCTTTGCCCCGACGTGATCGAGG',
          'ACGATCCCCCGACGTGATCGAGG',
          'ACGATCCCCCGACGTGATCGAGG'
          ];
          Last edited by rdeborja; 01-23-2015, 05:19 AM.

          Comment

          • GenoMax
            Senior Member
            • Feb 2008
            • 7142

            #6
            @rdeborj: Please enclose your code in [CODE_] code here [/_CODE]. Remove the underscores when you use the tags.

            Comment

            • codecatcher
              Junior Member
              • Jan 2015
              • 2

              #7
              Thanks rdeborja and GenoMax for the comments. I will try it out and see if I can get it working with that.
              Phoenix ICE I assume they are all equal parameters that need to be parsed if possible

              Comment

              • sklages
                Senior Member
                • May 2008
                • 628

                #8
                The above
                Code:
                /.{21}GG/g
                won't catch overlapping 23mers. Just keep that in mind. E.g.:
                Code:
                $sequence = 'ACGATCTTTGCCCCGACGTGATCGAGGTTTGGTTGGTTTTTTTTTCAGAGACCGAGATACGATCCCCCGACGTGATCGAGGACGATCCCCCGACGTGATCGAGG';

                Comment

                Latest Articles

                Collapse

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by SEQadmin2, Today, 10:09 AM
                0 responses
                9 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, Yesterday, 08:59 AM
                0 responses
                16 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 06-02-2026, 12:03 PM
                0 responses
                24 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 06-02-2026, 11:40 AM
                0 responses
                21 views
                0 reactions
                Last Post SEQadmin2  
                Working...