File: ensure_min_asm.pl

package info (click to toggle)
trinityrnaseq 2.15.2%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 468,004 kB
  • sloc: perl: 49,905; cpp: 17,993; java: 12,489; python: 3,282; sh: 1,989; ansic: 985; makefile: 717; xml: 62
file content (16 lines) | stat: -rwxr-xr-x 369 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
#!/usr/bin/env perl

use strict;
use warnings;

my $min_val = $ARGV[0] or die "Error, specify min full length";
my $trin_fa_file = $ARGV[1] or die "Error, need trin fa file";

my $asm_text = `cat $trin_fa_file.pslx.maps.summary`;
my @vals = split(/\t/, $asm_text);
if ($vals[1] < $min_val) {
    die "Error, $vals[1] < $min_val transcripts reconstructed";
}

exit(0);