Skip to content

Instantly share code, notes, and snippets.

View rec3141's full-sized avatar

R. Eric Collins rec3141

  • University of Alaska Fairbanks
View GitHub Profile
@rec3141
rec3141 / get-kos.sh
Created April 24, 2018 01:17
extract specific HMMs from FOAM
#!/bin/bash
Search=$1
grep -i $Search FOAM-onto_rel1.tsv > $Search.ont
grep -i $Search FOAM-onto_rel1.tsv | cut -f5 > $Search.ko
wc -l $Search.ko
rm $Search.acc
grep -A1 -f $Search.ko ./../FOAM/FOAM-hmm_rel1a_switch.hmm | grep NAME | tr -s ' ' | cut -f2 -d' ' | sort -u > $Search.acc
wc -l $Search.acc
@rec3141
rec3141 / foamcount.sh
Created April 26, 2018 22:40
count hits from FOAM hmmsearch
for FILE in *.tsv; do
LI=`grep -v '^#' $FILE | wc -l`;
SEQS=`gunzip -c ./../NPRB/$F1.faa.gz | grep -c '^>'`;
BP=`gunzip -c ./../NPRB/$F1.faa.gz | grep -v '^>' | wc -b`;
echo -e "$FILE\t$LI" >> target_hits.txt;
done
rm target_hits2.txt;
for FILE in *.tsv; do
ps -axeo pid,etime,command | awk '$2~/^0[1234]-/ && $0~/COMMAND/ { print $1 }' | xargs kill