Draw regional association plots

[dariushghasemi]

We use standalone locusZoom tool to generate the regional association plots at each locus of pQTL signals.

[dariushghasemi]

The rule is now able to create locusZoom plot with a the user-defined LD pattern for the target SNP in the region.
As teh wildcard is defined on seqid, we had plenty of steps to do:

• we opened locus breaker results for each seqid,
• we took region coordinates and top SNP of the region from locus breaker results,
• meta-analysis results are based on METAL, we used tabix to extract the region of the locus from GWAS file and feed it to locusZoom standalone tool
• finally, since LZ standalone tool could not remove the temporary folders for each locus, we remove the output file path from the prefix of the oputput files.
• then, we save them in the main directory.