QTLtools-extract(1)

QTLtools-extract1

Data extraction mode.

Usage

QTLtools extract [--vcf --bed --cov] relevant_file --out output_prefix [OPTIONS]

Overview

Extract all the data from the provided files into one flat file. This will also impute all the missing data.

Options

--vcf [in.vcf|in.bcf|in.vcf.gz]

Genotypes in VCF/BCF format. At least one of --vcf --bed --cov is REQUIRED.

--bed quantifications.bed.gz

Molecular phenotype quantifications in BED format. At least one of --vcf --bed --cov is REQUIRED.

--cov covariates.txt

Covariates data. At least one of --vcf --bed --cov is REQUIRED.

--out filename_prefix

Output file name prefix. REQUIRED.

--region chr:start-end

Genomic region to be processed. E.g. chr4:12334456-16334456, or chr5.

Output files

.content.txt.gz

Flat file with samples in rows and all the extracted variables in columns.

.header.txt

Genomic locations of the extracted variables.

Examples

o

Extract genotypes, gene quantifications, and covariates from a certain region:

QTLtools extract --vcf genotypes.chr22.vcf.gz --bed genes.50percent.chr22.bed.gz --cov genes.covariates.pc50.txt.gz --region chr22:17000000-18000000 --out extract_output

Citation

Delaneau, O., Ongen, H., Brown, A. et al. A complete tool set for molecular QTL discovery and analysis. Nat Commun 8, 15452 (2017). <https://doi.org/10.1038/ncomms15452>

Releases

Unchanged in 3 releases since Debian 11 (Bullseye).

  1. 11–13

Bugs

o

Versions up to and including 1.2, suffer from a bug in reading missing genotypes in VCF/BCF files. This bug affects variants with a DS field in their genotype's FORMAT and have a missing genotype (DS field is .) in one of the samples, in which case genotypes for all the samples are set to missing, effectively removing this variant from the analyses.

Please submit bugs to <https://github.com/qtltools/qtltools>

Authors

Olivier Delaneau ([email protected]), Halit Ongen ([email protected])

From qtltools 1.3.1+dfsg-4+b3 in Debian 13 (Trixie). Parsed by carta from the page's man source.