Skip to content

Commit

Permalink
added POD
Browse files Browse the repository at this point in the history
  • Loading branch information
daisieh committed Jul 16, 2013
1 parent f05fffe commit 8536a32
Showing 1 changed file with 33 additions and 2 deletions.
35 changes: 33 additions & 2 deletions parsing/trim_missing.pl
Expand Up @@ -11,10 +11,10 @@
my $runline = "running " . basename($0) . " " . join (" ", @ARGV) . "\n";

my $inputfile = "";
my $outname = "melded";
my $outname = "trimmed";
my ($help, $row_thresh, $col_thresh) = 0;

GetOptions ('files|input=s' => \$inputfile,
GetOptions ('fasta|input=s' => \$inputfile,
'outputfile=s' => \$outname,
'rowthreshold=f' => \$row_thresh,
'colthreshold=f' => \$col_thresh,
Expand Down Expand Up @@ -85,3 +85,34 @@
print FH ">$rowids[$i]\n$row\n";
}
close FH;


__END__
=head1 NAME
trim_missing
=head1 SYNOPSIS
trim_missing -input fastafile -output outputfile -row row_missing_frac -col col_missing_frac
=head1 OPTIONS
GetOptions ('fasta|input=s' => \$inputfile,
'outputfile=s' => \$outname,
'rowthreshold=f' => \$row_thresh,
'colthreshold=f' => \$col_thresh,
'help' => \$help) or pod2usage(-msg => "GetOptions failed.", -exitval => 2);
-fasta|input: fasta file of aligned sequences.
-outputfile: output file name.
-row: fraction of missing data allowed per row
-col: fraction of missing data allowed per column
=head1 DESCRIPTION
Takes a fasta file and removes rows and columns that contain more missing data than the specified thresholds.
=cut

0 comments on commit 8536a32

Please sign in to comment.