Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fix Scramble #6

Merged
merged 2 commits into from
May 19, 2022
Merged

Fix Scramble #6

merged 2 commits into from
May 19, 2022

Conversation

vladsavelyev
Copy link

  1. Drop HLA contigs from VCF to avoid:
Writing to /tmp/bcftools-sort.CBfZdQ
[W::vcf_parse] Contig 'HLA-DRB1*10' is not defined in the header. (Quick workaround: index the file with tabix.)
Error encountered while parsing the input at HLA-DRB1*10:2

HLA contig names are not well standardised, and even using the same reference build, SN:HLA-DRB1*10:01:01 in the reference genome somehow turns into HLA-DRB1*10 in a VCF. I don't think anything in GATK-SV does HLA typing, so it's safe to just drop them.

  1. Add more memory to the task to avoid OOM for SCRAMble.R:
/cromwell_root/script: line 119:    19 Killed                  Rscript --vanilla $scrambleDir/cluster_analysis/bin/SCRAMble.R --out-name $clusterFile --cluster-file $clusterFile --install-dir $scrambleDir/cluster_analysis/bin --mei-refs $meiRef --ref $xDir/ref --no-vcf --eval-meis

/cromwell_root/script: line 119:    19 Killed                  Rscript --vanilla $scrambleDir/cluster_analysis/bin/SCRAMble.R --out-name $clusterFile --cluster-file $clusterFile --install-dir $scrambleDir/cluster_analysis/bin --mei-refs $meiRef --ref $xDir/ref --no-vcf --eval-meis
@vladsavelyev vladsavelyev requested a review from vivbak May 16, 2022 05:09
@vladsavelyev vladsavelyev merged commit e9ef5c2 into merge-upstream-14may May 19, 2022
MattWellie added a commit that referenced this pull request Jul 26, 2023
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.

2 participants