bcftools
bcftools copied to clipboard
How to understand merge.txt in RoH
I cannot find a readme file for RoH output and recently got puzzled with results in merged.txt. I always thought "0" means non-RoH for an individual and "1" means RoH for an individual. However, recently I ran RoH with multiple samples and found an example of a ~900,000 bp region reported as RoH in homo_range.txt but "0" in merged.txt in that region.
It is possible you encountered a bug. Can you please share the output and the *.txt.gz
files to reproduce the problem? If it is sensitive data, can you send me an email offline? The address is on my profile page. Thank you.