
Perform permutation based interaction analysis (used to calculate corrected p-values for MPs)
  * Perform 100 permutations
  * Test interaction with term from [./FastQTL/example/interaction.txt]
  * Random number generator is seeded with 1
  * Considering variants within 1e+06 bp of the MPs
  * Chunk processed 1 / 100

Scanning phenotype data in [./FastQTL/example/phenotypes.bed.gz]
  * 364 phenotypes

Reading phenotype data in [./FastQTL/example/phenotypes.bed.gz]
  * region = 22:17517460-17739125
  * 373 samples included
  * 5 phenotypes included

Reading genotype data in [./FastQTL/example/genotypes.vcf.gz] in VCF format
  * region = 22:16517460-18739125
  * 373 samples included
  * 8684 sites included

Reading covariates in [./FastQTL/example/covariates.txt.gz]
  * 373 samples included
  * 4 covariate(s) included

Reading interaction term in [./FastQTL/example/interaction.txt]
  * 373 samples included

Imputing missing genotypes

Imputing missing phenotypes

Initialize covariate 

Processing gene [ENSG00000237438.1]
  * Number of variants in cis = 7966
  * Beta distribution parameters = 1.0382 2287.7101
  * Progress = 20.0%

Processing gene [ENSG00000177663.8]
  * Number of variants in cis = 8165
  * Beta distribution parameters = 0.7994 2826.0721
  * Progress = 40.0%

Processing gene [ENSG00000183307.3]
  * Number of variants in cis = 8279
  * Beta distribution parameters = 0.9595 9906.5566
  * Progress = 60.0%

Processing gene [ENSG00000069998.8]
  * Number of variants in cis = 8466
  * Beta distribution parameters = 0.9606 4564.2130
  * Progress = 80.0%

Processing gene [ENSG00000093072.10]
  * Number of variants in cis = 8531
  * Beta distribution parameters = 1.0645 1772.4331
  * Progress = 100.0%

Running time: 162 seconds
