Mercurial > repos > jdv > b2b_summarize_assembly
annotate summarize_assembly @ 1:5def63878840 draft default tip
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
| author | jdv | 
|---|---|
| date | Tue, 28 Sep 2021 06:16:07 +0000 | 
| parents | |
| children | 
| rev | line source | 
|---|---|
| 
1
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
1 #!/usr/bin/env perl | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
2 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
3 use strict; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
4 use warnings; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
5 use 5.012; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
6 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
7 use Getopt::Long; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
8 use Pod::Usage; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
9 use List::Util qw/sum max min/; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
10 use List::MoreUtils qw/uniq/; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
11 use BioX::Seq::Stream; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
12 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
13 my $PROGRAM = 'summarize_assembly'; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
14 my $VERSION = 0.002; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
15 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
16 my $fn_fasta; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
17 my @cutoffs; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
18 my $strip = 0; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
19 my $split = 0; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
20 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
21 # Collect command-line parameters | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
22 my $err_msg = 'Syntax error: please check your syntax'; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
23 GetOptions( | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
24 'fasta=s' => \$fn_fasta, | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
25 'cutoffs:i{,}' => \@cutoffs, | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
26 'strip_N' => \$strip, | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
27 'split_N' => \$split, | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
28 'help' => sub{ pod2usage( -verbose => 2 ) }, | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
29 'version' => sub{ print "This is $PROGRAM v$VERSION\n";exit; }, | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
30 ) or pod2usage( -msg => $err_msg, -verbose => 1 ); | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
31 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
32 # Set default cutoffs if necessary and sort | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
33 if (! scalar @cutoffs) { | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
34 warn "No cutoff supplied, defaulting to N50\n"; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
35 push @cutoffs, 50; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
36 } | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
37 @cutoffs = sort {$a <=> $b} uniq @cutoffs; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
38 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
39 # Only one of 'strip_N' or 'split_N' is valid | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
40 if ($strip && $split) { | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
41 warn "Only one of --strip_N or --split_N is valid, ignoring --strip_N\n"; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
42 $strip = 0; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
43 } | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
44 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
45 # Check for a few necessary conditions | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
46 die "Can't open FASTA file for reading" | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
47 if (defined $fn_fasta && ! -r $fn_fasta); | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
48 die "One or more cutoffs outside valid range (1-99)" | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
49 if (grep {$_ < 1 || $_ > 99} @cutoffs); | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
50 die "Cutoffs must be integer values" | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
51 if (grep {$_ ne int($_)} @cutoffs); | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
52 my @lens; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
53 my $N_sum = 0; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
54 my $GC_sum = 0; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
55 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
56 # Read in sequences and calculate descriptive stats | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
57 my $stream = BioX::Seq::Stream->new( $fn_fasta ); #STDIN if undefined | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
58 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
59 SEQ: | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
60 while (my $seq = $stream->next_seq) { | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
61 my @parts = ($seq); | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
62 @parts = split(/n+/i, $seq) if $split; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
63 for (@parts) { | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
64 my $Ns = ($_ =~ tr/Nn//); | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
65 $N_sum += $Ns; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
66 $GC_sum += ($_ =~ tr/GCgc//); | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
67 push @lens, length($_) - $Ns * $strip; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
68 } | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
69 } | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
70 @lens = sort {$b <=> $a} @lens; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
71 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
72 # Calculate basic stats | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
73 my $scaffold_count = scalar @lens; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
74 my $total_len = sum @lens; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
75 my $N_fraction = round( $N_sum/($total_len + $N_sum*$strip), 2 )*100; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
76 my $max_length = max @lens; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
77 my $min_length = min @lens; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
78 my $mean_length = round($total_len/$scaffold_count, 0); | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
79 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
80 # GC percentage calculated from non-ambiguous bases only | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
81 my $GC_fraction = round( $GC_sum/($total_len + ($strip - 1)*$N_sum), 2 )*100; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
82 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
83 # Calculate Nx (N50, N80, etc) values | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
84 # For example, N50 is the size of the smallest contig for which it and all | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
85 # larger contigs contain 50% of the total nucleotides in the database | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
86 my $cum_length = 0; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
87 my @fractions = map {$_/100} @cutoffs; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
88 my @Nx; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
89 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
90 LEN: | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
91 for (@lens) { | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
92 $cum_length += $_; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
93 if ($cum_length/$total_len >= $fractions[0]) { | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
94 push @Nx, $_; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
95 shift @fractions; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
96 last LEN if (@fractions < 1); | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
97 } | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
98 } | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
99 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
100 # Print summary | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
101 print '-' x 40 . "\n" | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
102 . "Summary\n" | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
103 . '-' x 40 . "\n" | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
104 . "number of scaffolds: $scaffold_count\n" | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
105 . "total length: $total_len\n" | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
106 . "average length: $mean_length\n" | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
107 . "G/C content: $GC_fraction\%\n" | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
108 . "ambiguous content: $N_fraction\%\n" | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
109 . "longest scaffold: $max_length\n"; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
110 for (0..$#Nx) { | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
111 my $label = sprintf "N%02d", $cutoffs[$_]; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
112 print "$label: $Nx[$_]\n"; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
113 } | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
114 print "shortest scaffold: $min_length\n"; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
115 print "NOTE: Ns were stripped for above calculations\n" if ($strip); | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
116 print "NOTE: Scaffolds were split on Ns for above calculations\n" if ($split); | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
117 print '-' x 40 . "\n"; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
118 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
119 exit; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
120 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
121 sub round { | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
122 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
123 my ($val,$places) = @_; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
124 return int($val*10**$places+0.5)/10**$places; | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
125 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
126 } | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
127 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
128 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
129 __END__ | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
130 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
131 =head1 NAME | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
132 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
133 summarize_assembly - print basic summary info for a file of assembly scaffolds | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
134 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
135 =head1 SYNOPSIS | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
136 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
137 summarize_assembly [--cutoffs I<cutoff_1> I<cutoff_2> .. I<cutoff_N> --strip_N --split_N ] --fasta I<input_file>] | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
138 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
139 =head1 DESCRIPTION | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
140 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
141 This program takes a FASTA file and optionally a list of cutoff values as | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
142 input and prints out summary information about the contigs/scaffolds contained | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
143 in the file. You can, of course, supply a FASTA file of any sort of nucleic | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
144 acid sequences, but the summary information makes most sense for contigs from | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
145 genomic sequencing assemblies. | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
146 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
147 =head1 OPTIONS | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
148 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
149 =over | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
150 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
151 =item B<--fasta> I<filename> | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
152 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
153 Specify contig/scaffold file from which to read input (default: STDIN) | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
154 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
155 =item B<--cutoffs> | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
156 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
157 Space-separated integer list of cutoffs to calculate (e.g. '--cutoffs 50 90' | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
158 will output N50 and N90 values) (default: 50) | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
159 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
160 =item B<--strip_N> | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
161 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
162 If specified, Ns will be stripped from scaffold sequences before statistics | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
163 are calculated (default: FALSE) | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
164 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
165 =item B<--split_N> | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
166 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
167 If specified, scaffold sequences will be split at regions of one or more Ns | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
168 before statistics are calculated (e.g. to get contig-level stats from a | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
169 scaffold file). Note that if this flag is specified, the value of '--strip_N' | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
170 will be ignored. (default: FALSE) | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
171 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
172 =item B<--help> | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
173 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
174 Display this usage page | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
175 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
176 =item B<--version> | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
177 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
178 Print version information | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
179 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
180 =back | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
181 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
182 =head1 CAVEATS AND BUGS | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
183 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
184 Please submit bug reports to the issue tracker in the distribution repository. | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
185 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
186 =head1 AUTHOR | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
187 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
188 Jeremy Volkening (jdv@base2bio.com) | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
189 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
190 =head1 LICENSE AND COPYRIGHT | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
191 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
192 Copyright 2014-16 Jeremy Volkening | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
193 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
194 This program is free software: you can redistribute it and/or modify | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
195 it under the terms of the GNU General Public License as published by | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
196 the Free Software Foundation, either version 3 of the License, or | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
197 (at your option) any later version. | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
198 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
199 This program is distributed in the hope that it will be useful, | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
200 but WITHOUT ANY WARRANTY; without even the implied warranty of | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
201 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
202 GNU General Public License for more details. | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
203 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
204 You should have received a copy of the GNU General Public License | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
205 along with this program. If not, see <http://www.gnu.org/licenses/>. | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
206 | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
207 =cut | 
| 
 
5def63878840
"planemo upload for repository https://github.com/jvolkening/galaxy-tools/tree/master/tools/b2b_utils commit 9bf8a0462bd44f170c0371b6cae67dd0c3b3da9f-dirty"
 
jdv 
parents:  
diff
changeset
 | 
208 | 
