INNER CODE UNIT · Python
concurrent_mode
PharmGKB/PharmCAT · preprocessor/pcat/preprocess.py:135
concurrent_mode=False, max_processes=1, verbose: int = 0) -> Path:
# shrink input VCF down to PGx allele defining regions and selected samples
# standardize chromosome names to <chr##>
pgx_region_vcf: Path = util.extract_pgx_regions(vcf_files, samples,
output_dir, output_basename,
regions_to_retain,
concurrent_mode=concurrent_mode, max_processes=max_processes,
verbose=verbose)
# normalize the input VCF
normalized_vcf = util.normalize_vcf(reference_genome, pgx_region_vcf, output_dir, output_basename, verbose=verbose)
# extract the specific PGx genetic variants in the reference PGx VCF
# this step also generates a report of missing PGx positions in the input VCF
pgx_variants_vcf: Path = util.extract_pgx_variants(pharmcat_positions_vcf, reference_genome, normalized_vcf,
output_dir, output_basename,
absent_to_ref=absent_to_ref,
unspecified_to_ref=unspecified_to_ref,