Mercurial > repos > lionelguy > spades

comparison tools/spades_2_5/spades.pl @ 7:95ddc2380130 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author lionelguy
date Thu, 28 Nov 2013 05:29:32 -0500
parents b5ce24f34dd7
children
comparison
equal deleted inserted replaced
6:1b1af74a54ae 7:95ddc2380130
74 open TAB, '>', $out or die $!; 74 open TAB, '>', $out or die $!;
75 print TAB "#name\tlength\tcoverage\n"; 75 print TAB "#name\tlength\tcoverage\n";
76 while (<FASTA>){ 76 while (<FASTA>){
77 next unless /^>/; 77 next unless /^>/;
78 chomp; 78 chomp;
79 my @a = split(/\s/, $_); 79 die "Not all elements found in $_\n" if (! m/^>NODE_(\d+)_length_(\d+)_cov_(\d+\.*\d*)_/);
80 my ($NODE, $n, $LENGTH, $l, $COV, $cov) = split(/_/, $a[0]); 80 my ($n, $l, $cov) = ($1, $2, $3);
81 die "Not all elements found in $_\n" unless ($n && $l && $cov);
82 print TAB "NODE_$n\t$l\t$cov\n"; 81 print TAB "NODE_$n\t$l\t$cov\n";
83 } 82 }
84 close TAB; 83 close TAB;
85 } 84 }