Mercurial > repos > ulfschaefer > phe_samtools_mpileup

annotate phe_samtools_mpileup.sh @ 11:714a4c13d64c draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author ulfschaefer
date Mon, 21 Dec 2015 10:17:56 -0500
parents f78ed0f61f52
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
1 #!/bin/bash
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
2
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
3 echo $@
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
4
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
5 BAM=$1
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
6 shift
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
7 REF=$1
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
8 shift
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
9 OUTVCF=$1
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
10
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
11 CMD="samtools mpileup -t DP,DV,DP4,DPR,SP -Auf $REF $BAM | bcftools call -m -f GQ > $OUTVCF"
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
12
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
13 echo "CMD: "$CMD
65b5406977c9 Uploaded
ulfschaefer
parents:
diff changeset
14 eval $CMD