fix: use correct bcftools filter column for tumor in manta tumor only #1400
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
This PR:
Correct bcftools filter in tumor / normal Manta rule to use the tumor column instead of the normal (bug).
Issue was discovered during validation that the bcftools filter in the manta_tumor_normal rule was using the normal column for filtering (index 0) instead of the tumor (index 1). This is not an issue in the manta_tumor_only as the tumor there is index 0. By switching the sample in the tumor_normal manta rule to index 1 for the filter, it should work as expected.
Fixed: for any bug fixes.
Tests to verify that issue is correct:
Run bcftools on the raw somantic manta calls on tumor-normal cases and confirm that it applies filter on tumor column in sample: K (see https://docs.google.com/document/d/12x3ozQE62bk_yWNTBW_u8vZMqhpKFRMW3yizzM_mWZs/)
command:
bcftools filter --threads 2 --exclude 'SUM(FORMAT/PR[1:1]+FORMAT/SR[1:1]) < 4' --soft-filter 'newfilter' --mode '+'
Sets filter on 5 variants: PR:SR 66,0:82,0 (normal) 31,1:38,2 (tumor)
original command sets filter on 504 variants: PR:SR 52,0:61,0 (normal) 63,5:81,2 (tumor)
Review and tests: