Conduct a PheWAS for a variant. The result will be p-values for the variant across all traits.

#download summary stat files for traits of interest
trait1.sumstats.gz
trait2.sumstats.gz
…
traitN.sumstats.gz

#variant of interest
chrom=11
pos=11324143

#conduct the phewas
for file in *.sumstats.gz; do 
  zcat $file | \
  awk -v OFS="\t" -v fname="$file" '{print fname, $0}'; done \
| (sed -u 1q; awk -v chrom="$chrom" -v pos="$pos" '$2 == chrom && $3 == pos') \
| sed '1 s/^\S\S*/Trait/' > phewas.out