Skip to content

Commit 7b45201

Browse files
author
Vyacheslav Brover
committed
#3: suffix
1 parent 5fe6307 commit 7b45201

File tree

1 file changed

+8
-20
lines changed

1 file changed

+8
-20
lines changed

phylogeny/database/genomeHash_make.sh

Lines changed: 8 additions & 20 deletions
Original file line numberDiff line numberDiff line change
@@ -5,12 +5,13 @@ if [ $# -ne 3 ]; then
55
echo "Create GenomeHash.{CDS,PRT,HMM}"
66
echo "#1: distance tree directory"
77
echo "#2: #1/../genome is large (0/1)"
8-
echo "#3: Create Genomehash.CDS (0/1)"
8+
echo "#3: Suffix: HMM, PRT, CDS"
99
exit 1
1010
fi
1111
INC=$1
1212
LARGE=$2
13-
CDS=$3
13+
SUF=$3
14+
1415

1516

1617
TMP=$( mktemp )
@@ -19,31 +20,18 @@ TMP=$( mktemp )
1920
GENOME=$INC/../genome
2021
$THIS/../../check_file.sh $GENOME 0
2122

22-
23-
# $TMP
23+
# $TMP, H
2424
if [ $LARGE -eq 1 ]; then
2525
$THIS/../../trav $GENOME "ls %d/%f" > $TMP
26+
H="%h/"
2627
else
2728
ls $GENOME > $TMP
29+
H=""
2830
fi
2931
wc -l $TMP
3032

31-
function run
32-
{
33-
local SUF=$1
34-
H=""
35-
if [ $LARGE -eq 1 ]; then
36-
H="%h/"
37-
fi
38-
section "GenomeHash.$SUF"
39-
$THIS/../../trav $TMP "cat $GENOME/$H%f/%f.hash-$SUF | sed 's/^/%f\t/1'" > GenomeHash.$SUF
40-
}
41-
42-
if [ $CDS -eq 1 ]; then
43-
run "CDS"
44-
fi
45-
run "PRT"
46-
run "HMM"
33+
$THIS/../../trav $TMP "cat $GENOME/$H%f/%f.hash-$SUF | sed 's/^/%f\t/1'" > GenomeHash.$SUF
34+
4735

4836

4937
rm $TMP*

0 commit comments

Comments
 (0)