Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • pony2001mx
    Member
    • Aug 2013
    • 32

    #1

    perl script to generate six frame sequences

    Hi everyone
    I wrote a script to generate siz frame sequences (+1 +2 +3 -1 -2 -3 frames). However, it only works on frames +1 +2 +3 -1. For example the following is -3 frame sequence generation. Where is wrong? Note after the script i make an output example. Thank you very much!

    Code:
    #!/bin/perl
    use strict;
    use warnings;
    
    my $frame = -3;  #when changing to -1 output sequences are ok; but -2 or -3 does not work;
    open IN, $ARGV[0];
    
    while (<IN>){
    	chomp;
    	my $m = length ($_);
    	
    	if(/^>(\w+)/){
    	 print ">$1"."_"."frame"."$frame"."\n";
    	 }else{
    		my $comprevseq = reverse $_;
    		$comprevseq =~ tr/[A,T,C,G]/[T,A,G,C]/; 
    		my $frameseq = substr ($comprevseq, abs($frame)-1, $m-abs($frame)+1);
    		
    		print "$frameseq\n";
    	}
    }
    close IN;
    INPUT file
    Code:
    >seq1
    ATGC
    >seq2
    ATGC
    >seq3
    TTTT
    OUTPUT result
    Code:
    >seq1_frame-3
    CAT
    >seq2_frame-3
    CAT
    >seq3_frame-3
    AA
  • SNPsaurus
    Registered Vendor
    • May 2013
    • 525

    #2
    This gives the correct answer (and is pretty much your code without the reading in the lines):
    #!/bin/perl

    my $frame = -3;
    $s = "ATGC";

    my $m = length ($s);


    my $comprevseq = reverse $s;
    $comprevseq =~ tr/[A,T,C,G]/[T,A,G,C]/;
    my $frameseq = substr ($comprevseq, abs($frame)-1, $m-abs($frame)+1);

    print "$frameseq\n";

    exit(0);
    output:
    AT
    The fact that the last input in your example works points to the problem being the linefeeds. You do chomp the input. But if the input has windows-type carriage returns and linefeeds there will be an extra invisible character. You can:
    local $/ = "\r\n";

    so that chomp removes both characters.
    Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

    Comment

    • pony2001mx
      Member
      • Aug 2013
      • 32

      #3
      Thank you very much, SNPsaurus. You comment is really helpful!

      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, 11:13 AM
      0 responses
      19 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, Yesterday, 10:48 AM
      0 responses
      18 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-15-2026, 10:52 AM
      0 responses
      36 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-14-2026, 05:39 AM
      0 responses
      42 views
      0 reactions
      Last Post SEQadmin2  
      Working...