Skip to content

remove prefix on snarl intermediate file input and output of stoat vcf - #59

Open
Plogeur wants to merge 3 commits into
mainfrom
remove_prefix
Open

remove prefix on snarl intermediate file input and output of stoat vcf#59
Plogeur wants to merge 3 commits into
mainfrom
remove_prefix

Conversation

@Plogeur

@Plogeur Plogeur commented Apr 1, 2026

Copy link
Copy Markdown
Member

PR Log Entry

Whoever merges this PR should copy the following bullet points to the PR Log:

  • Summary of the PR
    Linked to this issue: Deal with reference path prefixes #36
    Add new option on stoat vcf that remove the snarl chromosome prefix for both case :
    • stoat vcf snarl decomposition process
    • stoat vcf genotyping process

Description

@Plogeur Plogeur linked an issue Apr 1, 2026 that may be closed by this pull request
@Plogeur
Plogeur requested a review from jmonlong April 3, 2026 15:12
Comment thread src/subcommand/vcf.cpp Outdated
// JEAN would reduce memory to write the collection while genotyping the snarls, one chr at a time, appending to the output file (or in separate chr files).
auto start_writegt_timer = std::chrono::high_resolution_clock::now();
snarl_collection.write_snarl_data_collection(*gt_writer);
snarl_collection.write_snarl_data_collection(*gt_writer, remove_prefix_str);

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This one is not needed because we should already have homogenized the path names earlier (snarl_collection.genotype_snarls_by_chr_from_vcf).

Comment thread src/subcommand/vcf.cpp
snarl_writer.reset(new StdWriter(snarls_filename));
}

// equivalent to what was done before in stoat vcf: enumerate all walks through a snarl

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

If we want to remove the prefix in the intermediate file, we should do it around here. Otherwise, another approach, would be to keep the "real" path names in those file and just remove the prefix on the fly when reading/matching the VCF records (somewhere around snarl_collection.genotype_snarls_by_chr_from_vcf).

… decomposition are done in the same time that genotyping
@jmonlong

Copy link
Copy Markdown
Member

Let me know when this is ready to review.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Deal with reference path prefixes

2 participants