AWK script to parse info on a single sample and position from a VCF file

awk_vcf_anim_chrom_pos.sh

#!/bin/bash

# bash awk_vcf_anim_chrom_pos.sh in.vcf sample chrom pos

echo vcf: "$1"
echo sample: "$2"
echo chrom: "$3"
echo pos: "$4"

awk -v sample=$2 -v chrom="$3" -v pos=$4 \
	'/#CHROM/ {for(i=1;i<=NF;i++) if ($i~sample) f=i;next} ($1==chrom && $2==pos) {print $1 "_" $2 ": " $f " (" $9 ")"}' $1

FORMAT field could be added to the output by adding $9 to the print statement, for easier interpretation of the output

Edit: DONE