gwas2vcf icon indicating copy to clipboard operation
gwas2vcf copied to clipboard

check for the chrom names in the input vs dbSNP & reference fasta

Open darked89 opened this issue 4 years ago • 1 comments

Hello,

Would it be possible to add an extra check for chr name in the gwas/TSV input not present in the genomic fasta and/or dbSNP file? I am aware that this is handled by the alias file/switch, so it may be redundant. Nevertheless I was quite puzzled by:

 Variants discarded during harmonisation: 537691

Getting a line stating: "WARNING: chromosome(s) 23 and barfoo not present in Fasta: path_to_fasta_file" will make things easier to debug.

Hope it helps,

Darek Kedra

darked89 avatar Sep 07 '21 12:09 darked89

great suggestion, thanks @darked89. I will add soon

mcgml avatar Sep 13 '21 16:09 mcgml