+
+
+
+
+
+
+
+
+
+
+ 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 ); | +