-ignore_user_sites,how to treat these region? about clonalframeml

[xavierdidelot]

Dear YangChao,

Ignoring a region via the "-ignore_user_sites" option is not the same thing as removing the region from your input file. The difference is that if you removed the region from your input file, the two flanking regions would become attached in the input, so that linkage would considered high between them and recombination events would be likely to stretch across them. When you ignore a region using "ignore_user_sites", the linkage between the flanks is kept consistent with the length of the ignored region. So for example if the ignore region is very large, linkage will be low between the flanks and a recombination event will not be allowed to stretch across the ignored region. So the results will be more biologically meaningful. The output files are also easier to interpret, since the location of sites in the output is the same as the location of sites in the input, no matter which sites were ignored.

If you work with Mauve, I would recommend first to extract the core regions using the stripSubsetLCBs utility, and then use the resulting xmfa file as an xmfa input into ClonalFrameML, which requires to use the "-xmfa_file true" option but not the "ignore_user_sites" option.

Best wishes,
Xavier

[YangChaoCN]

Thank you so much for your help!