diff scripts/isoform_exists.sh @ 0:05c27700e5ca

initial commit
author biomonika <biomonika@psu.edu>
date Thu, 04 Sep 2014 18:24:19 -0400
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/scripts/isoform_exists.sh	Thu Sep 04 18:24:19 2014 -0400
@@ -0,0 +1,6 @@
+#!/bin/bash
+#we will search for isoforms of this contig
+contig_to_search=$1;
+reference=$2;
+#including only putative isoforms
+part=`echo $contig_to_search | sed -n "s/^\(.*\_c[0-9]\+\)\_seq[0-9]\+/\1/p"`; egrep "$part\_" $reference | grep -v $contig_to_search;
\ No newline at end of file