crivaldi · sy-li · Nov 27, 2018 · Nov 27, 2018
diff --git a/Q1.sh b/Q1.sh
@@ -0,0 +1,18 @@
+# This script is for Question 1 of EX11
+# Usage: bash Q1.sh
+
+cat gene_sequences/sporecoat01.fasta > gene_sequences/sp.fasta
+for sp in gene_sequences/sporecoat0[234].fasta
+do
+	echo \n >> gene_sequences/sp.fasta
+	cat $sp >> gene_sequences/sp.fasta
+done 
+./bin/muscle -in gene_sequences/sp.fasta -out gene_sequences/sp.align
+
+cat gene_sequences/transporter01.fasta > gene_sequences/tr.fasta
+for tr in gene_sequences/transporter0[234].fasta
+do
+  	echo \n >> gene_sequences/tr.fasta
+        cat $tr >> gene_sequences/tr.fasta
+done
+./bin/muscle -in gene_sequences/tr.fasta -out gene_sequences/tr.align
diff --git a/Q2.sh b/Q2.sh
@@ -0,0 +1,19 @@
+# This script is for Question 2 of EX11
+# Usage: bash Q2.sh
+
+# build a HMM profile for the transporter gene
+./bin/hmmbuild proteomes/tr.hmm gene_sequences/tr.align
+
+# search the hits to HMM profile
+for bp in proteomes/*.fasta
+do
+	./bin/hmmsearch --tblout ${bp%.*}.hits proteomes/tr.hmm $bp
+done
+
+# count for number of hits
+echo "Number of hits to the HMM profile in each proteome:" > proteomes/hit_counts.txt
+for hits in proteomes/*.hits
+do
+	echo ${hits##*/} >> proteomes/hit_counts.txt
+	cat $hits | grep -v "#" | wc -l >> proteomes/hit_counts.txt
+done
diff --git a/bin/alimask b/bin/alimask
diff --git a/bin/hmmalign b/bin/hmmalign
diff --git a/bin/hmmbuild b/bin/hmmbuild
diff --git a/bin/hmmconvert b/bin/hmmconvert
diff --git a/bin/hmmemit b/bin/hmmemit
diff --git a/bin/hmmfetch b/bin/hmmfetch
diff --git a/bin/hmmlogo b/bin/hmmlogo
diff --git a/bin/hmmpgmd b/bin/hmmpgmd
diff --git a/bin/hmmpress b/bin/hmmpress
diff --git a/bin/hmmscan b/bin/hmmscan
diff --git a/bin/hmmsearch b/bin/hmmsearch
diff --git a/bin/hmmsim b/bin/hmmsim
diff --git a/bin/hmmstat b/bin/hmmstat
diff --git a/bin/jackhmmer b/bin/jackhmmer
diff --git a/bin/makehmmerdb b/bin/makehmmerdb
diff --git a/bin/muscle b/bin/muscle
diff --git a/bin/nhmmer b/bin/nhmmer
diff --git a/bin/nhmmscan b/bin/nhmmscan
diff --git a/bin/phmmer b/bin/phmmer