for i in *.mutect.vcf.gz
do
  nohup ./filter.sh $i >log.`basename $i .mutect.vcf.gz` &
done

for i in *oncefiltered.vcf.gz
do
  zcat $i |awk '/^#/ || $7~/PASS/' - >`basename $i .somatic_oncefiltered.vcf.gz`.clean.vcf
done

for i in *.clean.vcf
do
  vcf2bed < $i  >`basename $i .vcf`.bed
done

if [ -s merged.bed ];then rm merged.bed;fi
for i in *.clean.bed
do
  name=`basename $i .clean.bed`
  awk 'BEGIN{OFS="\t"}{print $1,$2-length($6)+2,$3,$6,$7,"'$name'"}' $i >>merged.bed
done
sort -V merged.bed >merged.sorted.bed

annotate_variation.pl -buildver hg19 -outfile merged --otherinfo merged.sorted.bed  /database/annotation/annovar/humandb/