Mercurial > repos > lsong10 > psiclass
view PsiCLASS-1.0.2/FilterSplice.pl @ 0:903fc43d6227 draft default tip
Uploaded
author | lsong10 |
---|---|
date | Fri, 26 Mar 2021 16:52:45 +0000 |
parents | |
children |
line wrap: on
line source
#!/bin/perl use strict ; use warnings ; die "usage: a.pl a.raw_splice trusted.plisce > filtered.splice\n" if ( @ARGV == 0 ) ; my %trustedSplices ; open FP1, $ARGV[1] ; while ( <FP1> ) { chomp ; my @cols = split /\s+/, $_ ; my $key = $cols[0]." ".$cols[1]." ".$cols[2] ; $trustedSplices{ $key } = $cols[4] ; } close FP1 ; open FP1, $ARGV[0] ; while ( <FP1> ) { chomp ; my @cols = split /\s+/, $_ ; my $key = $cols[0]." ".$cols[1]." ".$cols[2] ; next if ( !defined $trustedSplices{ $key } ) ; my $trustedStrand = $trustedSplices{ $key } ; if ( $cols[3] <= 0 ) { print $cols[0], " ", $cols[1], " ", $cols[2], " 1 ", $trustedStrand, " 1 0 0 0\n" ; } else { if ( $cols[4] eq $trustedStrand ) { print $_, "\n" ; } else { $cols[4] = $trustedStrand ; print join( " ", @cols ), "\n" ; } } } close FP1 ;