Anonymous Monk has asked for the wisdom of the Perl Monks concerning the following question:
and this is where i've got to with the code...>gb|AE008689|:2073051-2073458, Atu4895 TTGGTCACATATTTTCGCTATTCGAAACCCAAATATCCCTGCGGAACACGTTTTATTGGAGCCGCCGTTT TGCAGTTGCTCGGTGCTGGATTCTTGGCGTTCGTCTTTTGTCTGCTGGATGGGCTCACTGCAAAACCAAC GATAATTCTCGGTCAGTTTGTATCATGCCTAGTGGGGAGCGTCGCCGGCTTTCATTTCGTGGCTTTTCGT CGCCCAGGCACGGACGGCCAACTTTACCTTATCGCGACGTCGCTCTTGGCATTTGGGACCCATTATTGGC TGGTGTCATATTCATTACCTGACCTTTTGCTAGCAAGATTGATTTCAGGATTTGGATCGGGTGTGGTTGT TGCAGGAACTTTCCGACGTCGCTTTCTGGAAAATCCGGTAATTCCCTGCGTTCGATAG >gb|AE008689|:c2074151-2073549, Atu4896 GTGTTTGGACCATATTTTCTGTGCAAACTAAACGATGACATAGGGCGATTTTTAGTGGCGGACAAATACA GACTTCCCGAAGAGTTTTTTACCACTCGGTTTCTCGTTAGACGCATCGTACCCACAGACGCTGAAGCTAT TTTCGAAGGGTGGAACACCGATCCCGAGGTGACGAAGTACCTGACGTGGAAACCCCACTCCGAGCTTGGC CAGACACAGCGGGCGATTGAAGAAAATTATAGTGCGTGGAATGCAGGTACATCGTTTCCAGCTGTCATCT GCCATCGCGAACGGCCACATGAACTAATCGGCCGTATTGATGCACGTCCGATGGGCCACAAGGTCTCTTA CGGGTGGCTTGTCCGAAGAACCTGGTGGGGCCGGGGTGTTGCAAGCGAGGTCGTTCAACTCGCTGTAGAA CACGCGTTATCGCATCCGCGCATCTTTCGCACCGAAGCATCCTGCGACGTTCTGAACACGGCGTCAGCAA GAGTGATGGAAAAAGTAGGGATGACAAAGGAGGCCGTGCTTCGACGGTACCTTTTTCACCCCAATTTTTC GAATATGCCGCGAGACGCCTTCCTGTATTCCAAGGTACGTTAA >gb|AE008689|:c2074749-2074345, Atu4897 ATGAAACATACCATCGCAGTTCTCGGCCTGATCACCTTCTCCAGCCCGGCCTTCGCAGCATCGTGCGAGA AAAACTTCACCGTCTCAGGCGTACCGATGGTCACGGCTGTCTCTTACAAATCCTTTCAGGAACTGCCGAA AGCCAAAGCACCAGCTGTCCTTCAAAAGCTCGCCCAGGCCGTCGCGGCAGAAGGTTTTTCAGGTATCCAG ATCAACAAGGCACTGTCGTCAATCGATGCCCATCAGGAAACCAGCGGAAGTGGCAGGATTCAGACGCTGC GGGTTGTCGCCCGCCAGAAAGGCGCCGCTGTCCGGATCGATGCTGTCTTCAATATTCAGGCAGGACAGAT CGCCGACAAAGACGTCATCCGCAAGGGCATCTGCGACATCATAAAAGGCGCGTAA
#! /usr/bin/perl -w use strict; my $num_of_params; $num_of_params = @ARGV; if ($num_of_params < 2) { die ("\nYou haven't entered enough parameters!\n\nTo run the progr +am type:\n\t\tperl gc3s_content.pl infile.ffn outfile.gc3s \n\n\t\tWh +ere window_size is the size of each window for which the program will + calculate results\n"); } open (INFILE, $ARGV[0]) or die "unable to open file $!"; open (OUTFILE,">$ARGV[1]") or die "unable to open file $!"; my @fasta = <INFILE>; my @gene; my $gene; my $flag; foreach my $line (@fasta) { if ($line =~ /^>/) { # print "$line\n"; $flag = 0; } # elsif ($line !~ /^>/) { elsif ($line =~ /^\w+/) { $gene = $line; $flag = 1; } if ($flag) { chomp ($gene); $gene .= $line; } } print "$gene\n";
|
|---|
| Replies are listed 'Best First'. | |
|---|---|
|
Re: parsing question
by BrowserUk (Patriarch) on Jan 15, 2004 at 12:45 UTC | |
by flounder99 (Friar) on Jan 15, 2004 at 12:55 UTC | |
|
Re: parsing question
by Roger (Parson) on Jan 15, 2004 at 12:49 UTC | |
|
Re: parsing question
by chimni (Pilgrim) on Jan 15, 2004 at 12:28 UTC | |
|
Re: parsing question
by flounder99 (Friar) on Jan 15, 2004 at 12:46 UTC | |
|
Re: parsing question
by alexg (Beadle) on Jan 15, 2004 at 14:44 UTC | |
|
Re: parsing question
by Hena (Friar) on Jan 15, 2004 at 12:50 UTC | |
|
Re: parsing question
by mr_mischief (Monsignor) on Jan 15, 2004 at 15:03 UTC |