search

jbloomlab / dms_tools2

software for the analysis and visualization of deep mutational scanning data
GNU General Public License v3.0
29 stars 20 forks source link

Problem with error conrtols #60

Closed gsifeng closed 2 years ago

gsifeng commented 2 years ago

Hi there,

This program is super useful and amazing! Just got one quick question. When I included --err using sequencing data of a wt gene in the diffsel program, it always terminates and shows "err counts of 0 for wildtype". Diffsel ran fine when I took the --err option out. What might have caused this issue and how could I solve it?

Thank you very much for your help!

jbloom commented 2 years ago

I think this means there is some problem with your wildtype sequencing error control? Is there a site in the wildtype error control where you don't have any sequencing counts even for the wildtype identity? That would cause this error.

In any case, you are probably fine running without the --err option?

gsifeng commented 2 years ago

Got it. Thank you!