Skip to content
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
30 changes: 18 additions & 12 deletions nextflow/vcf_prepper/bin/summary_stats.py
Original file line number Diff line number Diff line change
Expand Up @@ -225,18 +225,24 @@ def main(args=None):

# In case of no representative population - pick a random one
if population_name == "":
population = random.choice(population_data[species_patt])
population_name = population["name"]
freq_info_display = population["name"].replace(
"_", " "
) + population.get("version", "")

file = random.choice(population["files"])
freq_csq_fields.append(
file["short_name"]
+ "_"
+ random.choice(file["include_fields"])["fields"]["af"]
)
# If no population data for this species/assembly (i.e. it's been removed from assets/population_data.json), skip adding population fields
if species_patt not in population_data or not population_data[species_patt]:
population = None
else:
population = random.choice(population_data[species_patt])

if population:
population_name = population["name"]
freq_info_display = population["name"].replace(
"_", " "
) + population.get("version", "")

file = random.choice(population["files"])
freq_csq_fields.append(
file["short_name"]
+ "_"
+ random.choice(file["include_fields"])["fields"]["af"]
)

# add to header and write header to output vcf
if freq_info_display != "":
Expand Down
2 changes: 1 addition & 1 deletion nextflow/vcf_prepper/modules/local/process_gff.nf
Original file line number Diff line number Diff line change
Expand Up @@ -29,7 +29,7 @@ process PROCESS_GFF {
shell:
genome = meta.genome
genome_uuid = meta.genome_uuid
release_id = meta.release_id
release_id = params.release_id
version = params.version
out_dir = meta.genome_temp_dir
ini_file = params.ini_file
Expand Down
Loading