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
                  Essential Discoveries and Tools in Epitranscriptomics
                  by seqadmin


                  The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...
                  Yesterday, 07:01 AM
                • seqadmin
                  Current Approaches to Protein Sequencing
                  by seqadmin


                  Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
                  04-04-2024, 04:25 PM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 04-11-2024, 12:08 PM
                0 responses
                55 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 04-10-2024, 10:19 PM
                0 responses
                51 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 04-10-2024, 09:21 AM
                0 responses
                45 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 04-04-2024, 09:00 AM
                0 responses
                55 views
                0 likes
                Last Post seqadmin  
                Working...
                X