Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Help with a simple script in PERL

    sounds like homework
    Last edited by codecatcher; 01-23-2015, 03:12 PM.

  • #2
    What does your first draft look like?

    Comment


    • #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


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

        Comment


        • #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


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

            Comment


            • #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


              • #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

                • seqadmin
                  Best Practices for Single-Cell Sequencing Analysis
                  by seqadmin



                  While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
                  Yesterday, 07:15 AM
                • seqadmin
                  Latest Developments in Precision Medicine
                  by seqadmin



                  Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

                  Somatic Genomics
                  “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
                  05-24-2024, 01:16 PM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, Today, 06:58 AM
                0 responses
                13 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, Yesterday, 08:18 AM
                0 responses
                19 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, Yesterday, 08:04 AM
                0 responses
                18 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-03-2024, 06:55 AM
                0 responses
                13 views
                0 likes
                Last Post seqadmin  
                Working...
                X