Greetings, I'm working on a homework and I'm stuck with a part I just can't figure out. Today I asked another question and managed to solve that part, thanks to all who helped me on that occasion. Now, my problem is that I have printed the two genes I was supposed to find and now I'm supposed to print probes of length 6 with one mutation OF these 2 genes. These probes are supposed to be only explicit content of the 2 genes. A probe can't exist in the whole text and both of the genes. For example, a gene has CACTCT substring in it and that CACTCT should not be existent in the whole text and the other gene. Also CACTCT should have one letter different from the entire text. -mutations- my code is as below:
$infile="input.txt";
open(IN, "$infile") or die "can'topen \"$infile\"";
$genom="";
while($line=<IN>)
{
chomp $line;
$genom.=$line;
}
@allprobes=();
for($x=0;$x<=$alltextlength-6;$x++)
{
$all= substr $genom, $j, 6;
push(@allprobes,$all);
}
@probe=();
while($genom=~m/TATAAT.*?(ATG(...)*?(TAA|TAG|TGA))/g)
{
push (@gen, $1);
}
$cnt = 0;
$k=0;
$m=0;
foreach $gen (@gen)
{
$length=length($gen);
for($j=0; $j<$length-5; $j++)
{
$altilik= substr $gen, $j, 6;
push (@probe,$altilik);
while($allprobes[$m]=~m/$probe[$k]/g)
{
#push (@gen, $1);
$cnt++;
$k++;
}
$m++;
#print "@probe\n\n\n";
#print @allprobes;
}
}
print "$cnt\n";
Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
Read Where should I post X? if you're not absolutely sure you're posting in the right place.
Please read these before you post! —
Posts may use any of the Perl Monks Approved HTML tags:
- a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
| |
For: |
|
Use: |
| & | | & |
| < | | < |
| > | | > |
| [ | | [ |
| ] | | ] |
Link using PerlMonks shortcuts! What shortcuts can I use for linking?
See Writeup Formatting Tips and other pages linked from there for more info.