http://www.perlmonks.org?node_id=960989

anisha3 has asked for the wisdom of the Perl Monks concerning the following question:

This node falls below the community's threshold of quality. You may see it by logging in.
  • Comment on programme to get the sequence with cut ends of 250 bp on both sides and also i need to get the cut pieces please suggest the rigght way

Replies are listed 'Best First'.
Re: perl programme to sort out snps from human genome ....plz suggest me the right way . iam confused how to write
by marto (Cardinal) on Mar 22, 2012 at 13:33 UTC

    "this is the one which i wrote by myself if its wrong plz letme knw"

    I'm sure if you ran your script you'd see it was wrong, in the sense that it doesn't work. Here is an example which does what I think you are trying to do:

    #!/usr/bin/perl

use strict;
use warnings;

my $word1 = 'Hcr1fa.fa'; # you missed the ;

if ( $word1 eq "chr1.subst.fa" ){ # I added quotes
    print "$word1\n"; 
} elsif ( $word1 eq "Hcr1fa.fa" ) { 
    print "$word1\n" 
}

    [download]

    Now I'm sure you know what you're trying to achieve with your code, I don't. My code prints the value of $word1 if the value is "chr1.subst.fa" or "Hcr1fa.fa". You'll see from the code and the comments that I added quotes around chr1.subst.fa, are you expecting this string or are you trying to match something else?

    When you previewed your post, you must have seen that something was wrong with the formatting. Your closing code tag should be </code> rather than <\code>.

    Here are some links I think you should read and understand:

[reply]
[d/l]
[select]

Back to Seekers of Perl Wisdom