Mercurial > repos > biomonika > linkyx

view LINKYX_mpileup_wrapper.sh @ 9:695d28139f3e

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
toolshed8
author biomonika <biomonika@psu.edu>
date Tue, 09 Sep 2014 14:31:02 -0400
parents 1955f03f092e
children ea52ef8d7f3a
line wrap: on
line source

#!/bin/bash
results=$1;
bam_file=$2;
reference=$3;

sorted_bam=${bam_file}_sorted

#LINKYX_PATH set during toolshed installation

samtools mpileup -uf $reference $bam_file 2>log.txt | bcftools view -p 0.85 -cgv - 2>log.txt >$results

#samtools mpileup -uf $reference $2 2>log.txt | bcftools view -p 0.85 -cgv - 2>>log.txt| vcfutils.pl varFilter >$results 2>>log.txt