Skip to content

Instantly share code, notes, and snippets.

@genomewalker

genomewalker/update_unk.sh

Last active October 31, 2017 16:44
Show Gist options
  • Save genomewalker/e56e167e896d01c9d6ed57968d8412f5 to your computer and use it in GitHub Desktop.
Save genomewalker/e56e167e896d01c9d6ed57968d8412f5 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
update_unk.sh
mmseqs concatdbs /PROCESSING/OSD/UNKNOWNS/MMSEQ_MPI/unkdb_update_hmp/unkprot_db \
/PROCESSING/OSD/UNKNOWNS/MMSEQ_MPI/unkdb_update_hmp/HMPSOAP_db \
/PROCESSING/OSD/UNKNOWNS/MMSEQ_MPI/unkdb_update_hmp/unkprot_db.withNewSequences
mmseqs concatdbs /PROCESSING/OSD/UNKNOWNS/MMSEQ_MPI/unkdb_update_hmp/unkprot_db_h \
/PROCESSING/OSD/UNKNOWNS/MMSEQ_MPI/unkdb_update_hmp/HMPSOAP_db_h \
/PROCESSING/OSD/UNKNOWNS/MMSEQ_MPI/unkdb_update_hmp/unkprot_db.withNewSequences_h
mawk '{print NR"\t"$1}' unkprot_db.withNewSequences_h | sed 's/\x0//g' > unkprot_db.withNewSequences.lookup
mmseqs clusterupdate /PROCESSING/OSD/UNKNOWNS/MMSEQ_MPI/unkdb_update_hmp/unkprot_db \
/PROCESSING/OSD/UNKNOWNS/MMSEQ_MPI/unkdb_update_hmp/unkprot_db.withNewSequences \
/PROCESSING/OSD/UNKNOWNS/MMSEQ_MPI/unkdb_update_hmp/unkclust30_2016_11 \
/PROCESSING/OSD/UNKNOWNS/MMSEQ_MPI/unkdb_update_hmp/newMappedSequenceDB \
/PROCESSING/OSD/UNKNOWNS/MMSEQ_MPI/unkdb_update_hmp/test_default \
/PROCESSING/OSD/UNKNOWNS/MMSEQ_MPI/unkdb_update_hmp/tmp --threads 128
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment