fedarko / strainFlye

Pipeline for analyzing (rare) mutations in metagenome-assembled genomes
BSD 3-Clause "New" or "Revised" License
8 stars 1 forks source link

Completely avoid division of values of p by 100 (at least in calling) #14

Closed fedarko closed 2 years ago

fedarko commented 2 years ago

Can keep everything in terms of integers (1 = 0.01%, 15 = 0.15%, 200 = 2%, etc.) and just adjust the p-mutation calling function to handle this as a comparison between integers.

Probably fine as is, tho.

fedarko commented 2 years ago

Done! See:

https://github.com/fedarko/strainFlye/blob/0458989ce3dd6c13c5dd99a4fde776b27d71535d/strainflye/call_utils.py#L580-L627