msmc2 icon indicating copy to clipboard operation
msmc2 copied to clipboard

Published 20 hours ago • verified publisher icon stschiff

Haploid problem

Open cfz1998 opened this issue 2 years ago • 3 comments
trafficstars

Hi! @stschiff

Does MSMC2 work for haploid? My organisms is haploid.

Thank you for your reply! --CFZ

cfz1998 avatar Jun 29 '23 08:06 cfz1998

Yes, if it is recombining. If not, then you don't need MSMC, you can just use any phylogenetic software or build a neighbour-joining tree.

stschiff avatar Jul 04 '23 07:07 stschiff

Hi！@stschiff
Should make a “pseudo-diploid”(merge two sample genotypes) after bamCaller.py? Or just simply double each sample genotype and pick only one haplotype to msmc2? Like: msmc2 -t 40 -p 1*2+25*1+1*2 -i 100 -I 0,2,4,6,8,10,12,14 8 samples

--CFZ

cfz1998 avatar Jul 20 '23 11:07 cfz1998

And, if merged two Haploid samples to a “pseudo-diploid”. How to deal with the “mask.bed.gz” from bamCaller.py? Simply merged(Interval merging) the “mask.bed.gz” from two Haploid samples?

cfz1998 avatar Jul 20 '23 13:07 cfz1998