-
-
-
-
-
-
-
-
-
-
- Permalink
-
-
-
-
-
-
-
-
-
-
-
-
- Join GitHub today
-GitHub is home to over 36 million developers working together to host and review code, manage projects, and build software together.
- Sign up -
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- Find file
-
-
- Copy path
-
-
-
- - fasta-tool/fasta-tool.pl -
- -
-
- Find file
-
-
- Copy path
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- | - | #!/usr/bin/perl -w | -
| - | - | -
| - | # usage : perl toolExample.pl <FASTA file> <output file> | -
| - | - | -
| - | open (IN, "<$ARGV[0]"); | -
| - | open (OUT, ">$ARGV[1]"); | -
| - | while (<IN>) { | -
| - | chop; | -
| - | if (m/^>/) { | -
| - | s/^>//; | -
| - | if ($. > 1) { | -
| - | print OUT sprintf("%.3f", $gc/$length) . "\n"; | -
| - | } | -
| - | $gc = 0; | -
| - | $length = 0; | -
| - | } else { | -
| - | ++$gc while m/[gc]/ig; | -
| - | $length += length $_; | -
| - | } | -
| - | } | -
| - | print OUT sprintf("%.3f", $gc/$length) . "\n"; | -
| - | close( IN ); | -
| - | close( OUT ); | -