You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
I've got a set of 972 Mtb isolates that I'm trying to run snp-sites -c -o on, but it fails with the error Warning: No SNPs were detected so there is nothing to output. However, it works with removing the -c flag. How can I try and figure out which isolates are causing problems?
The text was updated successfully, but these errors were encountered:
This seems to happen because the -c option removes any call across all samples where any of your samples has a non-ACTG at that position. This could be due to a true deletion at that site, an area of a particular sample that has not aligned to your reference (low depth? poor quality?), or a site with some other non-ACTG character. With a large number of samples, even a single sample missing a call at a position will result in that call being removed entirely when -c is used. With a large number of samples, you're probably just filtering out all sites because at least one sample has a non-ACTG for every site.
I've got a set of 972 Mtb isolates that I'm trying to run
snp-sites -c -o
on, but it fails with the errorWarning: No SNPs were detected so there is nothing to output
. However, it works with removing the-c
flag. How can I try and figure out which isolates are causing problems?The text was updated successfully, but these errors were encountered: