Guest User

Untitled

a guest
Feb 20th, 2018
100
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Perl 1.34 KB | None | 0 0
  1. #!/usr/bin/perl
  2.  
  3. use warnings;
  4.  
  5. print ("Please input the first DNA-sequence: ");
  6. chomp($dna1 = <STDIN>);
  7.  
  8. print ("Please input the second DNA-sequence: ");
  9. chomp($dna2 = <STDIN>);
  10.  
  11. $dna1len = length($dna1);
  12. $dna2len = length($dna2);
  13. $itr = 0;
  14. $num = 1;
  15. if ($dna1len < $dna2len){
  16.     for ($itr=0; $itr != $dna2len; $itr++){
  17.     $sub1 = substr $dna1, $itr, $num;
  18.     $sub2 = substr $dna2, $itr, $num;
  19.     push(@dna1, $sub1);
  20.     push(@dna2, $sub2);
  21.  
  22. }
  23. }
  24. elsif ($dna1len > $dna2len){
  25.     for ($itr=0; $itr != $dna1len; $itr++){
  26.     $sub1 = substr $dna1, $itr, $num;
  27.     $sub2 = substr $dna2, $itr, $num;
  28.     push(@dna1, $sub1);
  29.     push(@dna2, $sub2);
  30.    
  31.     }
  32. }
  33. else {
  34.     for ($itr=0; $itr != $dna1len; $itr++){
  35.     $sub1 = substr $dna1, $itr, $num;
  36.     $sub2 = substr $dna2, $itr, $num;
  37.     push(@dna1, $sub1);
  38.     push(@dna2, $sub2);
  39.    
  40. }
  41. }
  42. print ("DNA1 has $dna1len nucleotides, DNA2 has $dna2len nucleotides.\n");
  43. print ("@dna1 and @dna2\n");
  44. $itr1 = 0;
  45.  
  46. if ($dna2len < $dna1len){
  47.     $lim = $dna2len;
  48.     }
  49.     else {
  50.         $lim = $dna1len;
  51.         }
  52.  
  53. while ($itr1 != $lim){
  54.    for ($dna1[$itr1] eq $dna2[$itr1]){
  55.     push(@dnaarray, $dna1[$itr1]);
  56.     $itr1++
  57.     }
  58. }
  59.  
  60. #for ($itr1=0; $dna1[$itr1] eq $dna2[$itr1]; $itr1++){
  61. #    push(@dnaarray, $dna1[$itr1]);
  62. #    }
  63. print ("lol $itr1 @dnaarray\n");
Add Comment
Please, Sign In to add comment