Mercurial > repos > biomonika > linkyx

annotate scripts/isoform_exists.sh @ 0:1955f03f092e

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
initial commit
author biomonika <biomonika@psu.edu>
date Tue, 02 Sep 2014 14:16:04 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
1 #!/bin/bash
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
2 #we will search for isoforms of this contig
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
3 contig_to_search=$1;
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
4 reference=$2;
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
5 #including only putative isoforms
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
6 part=`echo $contig_to_search | sed -n "s/^\(.*\_c[0-9]\+\)\_seq[0-9]\+/\1/p"`; egrep "$part\_" $reference | grep -v $contig_to_search;