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

                • kinpict.com
                  Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
                  by kinpict.com
                  The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
                  Today, 06:37 PM
                • SEQadmin2
                  16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
                  by SEQadmin2





                  NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
                  10-30-2026, 06:46 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by SEQadmin2, 11-12-2026, 07:59 AM
                0 responses
                27 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 11-05-2026, 09:25 AM
                0 responses
                35 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 11-03-2026, 11:54 AM
                2 responses
                56 views
                0 reactions
                Last Post kinpict.com  
                Started by SEQadmin2, 10-28-2026, 09:57 AM
                0 responses
                175 views
                0 reactions
                Last Post SEQadmin2  
                Working...