Mercurial > repos > mmseq > mmseq
view mmseq.sh @ 10:ca56c4ab0edf default tip
Incorporate samtools invocation
author | mwoodbri |
---|---|
date | Thu, 27 Oct 2011 19:01:46 -0400 |
parents | |
children |
line wrap: on
line source
#!/bin/sh PREFIX=temp FA=$(grep -P "index\t$2\t" $3/sam_fa_indices.loc | cut -f3) samtools view -bt $FA.fai -o $PREFIX.bam $1 2> /dev/null samtools sort -n $PREFIX.bam $PREFIX.namesorted 2> /dev/null bam2hits $FA $PREFIX.namesorted.bam > $PREFIX.hits 2> /dev/null mmseq $PREFIX.hits $PREFIX > /dev/null mv $PREFIX.mmseq $4 mv $PREFIX.identical.mmseq $5 mv $PREFIX.gene.mmseq $6