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: cnvpytor header float #1182

Merged
merged 7 commits into from
Sep 19, 2023
Merged

fix: cnvpytor header float #1182

merged 7 commits into from
Sep 19, 2023

Conversation

mathiasbio
Copy link
Collaborator

@mathiasbio mathiasbio commented Jun 21, 2023

This PR:

Attempts to fix issues:

With the solution described here:

Added:

  • sed command to cnvpytor rule to convert CN description in header from integer to float.

Tests:

WGS tumor only cases previously failed with: [E::vcf_parse_format] Invalid character '.' in 'CN' FORMAT in bcftools_quality_filter_svdb.

petmantis: previously had this error, now this error should not appear. Case has variant with float CN at:

21	16700001	CNVpytor_dup20	.	<DUP>	.	PASS	END=44200000;IMPRECISE;SVLEN=27500000;SVTYPE=DUP;pytorRD=2.2524309640529316;pytorP1=0.0;pytorP2=0.0;pytorP3=0.0;pytorP4=0.0;pytorQ0=0.0015552215961759483;pytorPN=0.0018927272727272727;pytorDG=0.0;pytorCL=rd_mean_shift	GT:CN	./1:2.25
  • CN with float value exists
  • Analysis finishes successfully

awareraptor: as described in issue (#1157) has variant with float CN at:

21      22900001        CNVpytor_dup22  .       <DUP>   .       PASS    END=44200000;IMPRECISE;SVLEN=21300000;SVTYPE=DUP;pytorRD=3.70678765025364;pytorP1=0.0;pytorP2=0.0;pytorP3=0.0;pytorP4=0.0;pytorQ0=0.0011310821079172444;pytorPN=0.002443661971830986;pytorDG=0.0;pytorCL=rd_mean_shift;set=cnvpytor;FOUNDBY=1;cnvpytor_CHROM=CNVpytor_dup22|21;cnvpytor_POS=CNVpytor_dup22|22900001;cnvpytor_QUAL=CNVpytor_dup22|.;cnvpytor_FILTERS=CNVpytor_dup22|PASS;cnvpytor_SAMPLE=CNVpytor_dup22|TUMOR|GT:./1|CN:3.71;cnvpytor_INFO=CNVpytor_dup22|END:44200000|IMPRECISE|SVLEN:21300000|SVTYPE:DUP|pytorRD:3.70678765025364|pytorP1:0.0|pytorP2:0.0|pytorP3:0.0|pytorP4:0.0|pytorQ0:0.0011310821079172444|pytorPN:0.002443661971830986|pytorDG:0.0|pytorCL:rd_mean_shift;svdb_origin=cnvpytor;SUPP_VEC=10000     GT:CN   ./1:3.71
  • CN with float value exists
  • Analysis finishes successfully

Other integrity tests:

  • SV variants uploaded to Scout are unaffected by float header change
  • SV variants can still be uploaded as usual to LoqusDB

Review and tests:

  • Tests pass
  • Code review
  • New code is executed and covered by tests, and test approve

@sonarqubecloud
Copy link

Kudos, SonarCloud Quality Gate passed!    Quality Gate passed

Bug A 0 Bugs
Vulnerability A 0 Vulnerabilities
Security Hotspot A 0 Security Hotspots
Code Smell A 2 Code Smells

No Coverage information No Coverage information
0.0% 0.0% Duplication

@codecov
Copy link

codecov bot commented Jun 21, 2023

@mathiasbio mathiasbio changed the base branch from master to develop June 21, 2023 08:06
@mathiasbio mathiasbio marked this pull request as ready for review September 12, 2023 15:30
@mathiasbio mathiasbio requested a review from a team as a code owner September 12, 2023 15:30
@mathiasbio mathiasbio linked an issue Sep 14, 2023 that may be closed by this pull request
Copy link
Contributor

@ivadym ivadym left a comment

Choose a reason for hiding this comment

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

Well done 🧁

@sonarqubecloud
Copy link

Kudos, SonarCloud Quality Gate passed!    Quality Gate passed

Bug A 0 Bugs
Vulnerability A 0 Vulnerabilities
Security Hotspot A 0 Security Hotspots
Code Smell A 0 Code Smells

No Coverage information No Coverage information
No Duplication information No Duplication information

@mathiasbio mathiasbio merged commit 9e5ad3c into develop Sep 19, 2023
@mathiasbio mathiasbio deleted the fix_cnvpytor_float branch September 19, 2023 13:16
mathiasbio added a commit that referenced this pull request Nov 4, 2023
This PR aims to change the name of the sample columns in the final VCFs to match the changes in the scout_load.yaml for sample-id, which will use LIMS-ID instead of TUMOR/NORMAL. (Clinical-Genomics/cg#2650)

Added:
- New rule to create namemap to translate sample type names in vcf header to sampleID

Changed:
- Renamed names of sample columns in final vcf (clinical.pass.vcf.gz) to sampleID

Removed:
- Removed sed command in CNVpytor rule introduced here (fix: cnvpytor header float #1182) as it has been fixed by updating CNVpytor and is no longer necessary ([ Failed Analysis ] CNVpytor float values not accepted by bcftools #1152)
@ivadym ivadym mentioned this pull request Nov 6, 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.

fix cnvpytor CN integer-float conflict
2 participants