Mercurial > repos > bcclaywell > microbiome_pplacer_suite

annotate alignment-wrapper.sh @ 8:35f9b833e12b draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 2774930eebe258ecd56c8f1c5ddabf5092282ab9
author bcclaywell
date Mon, 12 Oct 2015 16:09:02 -0400
parents 221a97bf686d
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
1 #!/bin/bash
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
2
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
3 source $1
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
4
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
5 # adapted from yapp/bin/refpkg_align
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
6 ref_sto=$(taxit rp ${REFPKG} aln_sto)
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
7 profile=$(taxit rp ${REFPKG} profile)
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
8
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
9 sto=$(mktemp -u).sto
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
10
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
11 cmalign --cpu ${GALAXY_SLOTS:-4} -o "$sto" --sfile "${ALIGNED_SCORES}" --noprob --dnaout "$profile" "${INPUT_SEQS}" | grep -E '^#'
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
12
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
13 esl-alimerge --dna --outformat afa "$ref_sto" "$sto" | \
221a97bf686d Uploaded
bcclaywell
parents:
diff changeset
14 seqmagick convert --output-format fasta --dash-gap - "${ALIGNED_SEQS}"