Unconfigured Ad

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

    #1

    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

                • SEQadmin2
                  Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
                  by SEQadmin2



                  Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

                  Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
                  ...
                  10-14-2026, 09:55 AM
                • SEQadmin2
                  New Genomics Technologies Take Aim at Long-Standing Limits
                  by SEQadmin2


                  Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

                  We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
                  ...
                  09-28-2026, 10:25 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by SEQadmin2, Yesterday, 10:52 AM
                0 responses
                10 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 10-14-2026, 05:39 AM
                0 responses
                19 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 10-13-2026, 09:56 AM
                0 responses
                17 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 10-09-2026, 11:14 AM
                0 responses
                27 views
                0 reactions
                Last Post SEQadmin2  
                Working...