Entering edit mode
7.9 years ago
Jautis
▴
530
Hi,
I have a sample for which I am able to call very few genotype calls (virtually none) compared to other samples in the dataset. This is despite having relatively high total coverage and number of mapped reads.
All samples were digested with Msp1, but the digestion efficiency (proportion of reads starting and/or ending at a Msp1 cut site) does not differ between samples.
What are some potential causes for this and how can I examine what is happening?
Providing some additional info would be helpful. How many reads and what coverage do you mean by high? What organism are you working with? I am assuming you are doing RAD-Seq or something similar?
What software/pipeline are you using? etc