reformat print line to get rid off unwanted spaces #606
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.
Hi, Thanks for the software !
I've noticed that the extract_snps script introduces spaces directly after the geno[1] field, at the position where the print line is split up over multiple lines. This causes issues in the assignment later on as no reads can be assigned to the alternative genotype. The pipeline fails later in the matrix generation.
Reformatting this line gets rid of the spaces and fixes the issues downstream.
I've tested with mgp_REL2021_snps.vcf.gz, corresponding to GRCm39.
Related to #572.
Kind regards,
WardDeb