bmvdgeijn / WASP

WASP: allele-specific pipeline for unbiased read mapping and molecular QTL discovery
Apache License 2.0
102 stars 51 forks source link

Could not find test snp #90

Open SaideepGona opened 4 years ago

SaideepGona commented 4 years ago

In the: extract_haplotype_read_counts.py step, I keep getting errors such as the following:

WARNING: could not find test SNP at position 1:74570

I'm having trouble tracing the issue back to the source as it is complicated to follow. Do you know the likely causes of this? I'm pretty sure I made the bed-like file to specification.

bmvdgeijn commented 4 years ago

Hi Saideep,

I believe this means there is a mismatch between the SNP and genotype files and the bed-like file. Is there a SNP on chr1:74570 in your input genotypes? Another possibility is there is an off by one and your SNP file is 0-based and bed is 1-based.

On Sat, Nov 9, 2019 at 8:15 PM Saideep Gona notifications@github.com<mailto:notifications@github.com> wrote:

In the: extract_haplotype_read_counts.py step, I keep getting errors such as the following:

WARNING: could not find test SNP at position 1:74570

I'm having trouble tracing the issue back to the source as it is complicated to follow. Do you know the likely causes of this? I'm pretty sure I made the bed-like file to specification.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://github.com/bmvdgeijn/WASP/issues/90?email_source=notifications&email_token=AA2IHAR4QPLNZQKXCN5NW5DQS5OCVA5CNFSM4JLKLDFKYY3PNVWWK3TUL52HS4DFUVEXG43VMWVGG33NNVSW45C7NFSM4HYGYTFA, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AA2IHAVJWFJ7BOFT6PCTX3DQS5OCVANCNFSM4JLKLDFA.