--similarity cutoff default about pstrain

[wshuai294]

Hi Shriram,

Thank you for your careful observation. The default cutoff I set is 0.80 actually. I will revise it.

In practice, please use the 'python3 merge.py ' to cluster the strains with different cutoff values like 0.85, 0.90, 0.95, 0.97, and so on. The script will save the results for different cutoff values. You can choose a meaningful value in your analysis.

In the script, I set --qual to filter SNP with SNP quality. Also, I set --snp_ratio to filter SNP according to the comparison of locus' depth and species' mean depth. You can use these two parameters to filter the SNP. So, I think it is not very necessary to use DP and MAPQ to filter SNP.

I think you can use 'samtools tview' or other tools to look at the reads mapping at the SNP locus. If the SNP is not reliable, you can filter SNP with stricter criteria. It is useful to test with different criteria.

Best,
Shuai