Virstrain build sites cutoff default value.

[MTDouglas]

Hello,

I'm just curious as to what the default value being used in Virstrain_build --sites_cutoff is when no value is specified?

Thanks

[liaoherui]

Hi, thanks for using VirStrain!

By default, VirStrain applies "greedy covering algorithm" to find the minimum sites that can distinguish different viral strains automatically (more details can be found in the manuscript). So, this parameter is not used. This parameter may be needed/added when the sequencing depth or quality is low and we want more sites be covered for strain identification.