--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/PsiCLASS-1.0.2/FilterSplice.pl	Fri Mar 26 16:52:45 2021 +0000
@@ -0,0 +1,48 @@
+#!/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 ;
+