• 全国 [切换]
  • 二维码
    错改网

    手机WAP版

    手机也能找商机,信息同步6大终端平台!

    微信小程序

    微信公众号

    当前位置: 首页 » 行业新闻 » 热点新闻 » 正文

    OmegaPlus 检测正选择信号

    放大字体  缩小字体 发布日期:2024-12-11 04:02:29   浏览次数:31  发布人:7463****  IP:124.223.189***  评论:0
    导读

    OmegaPlus 是一个用于检测基因组数据中正选择痕迹的生物信息学工具。基于单倍型频率分布和重组事件检测选择信号,通过计算ω统计量(omega statistic)来识别选择位点。官方软件提供了几种不同的版本:OmegaPlus:顺序版本;OmegaPlus-F:细粒度并行版本,适合在在更小的任务单位上并行化;OmegaPlus-C:粗粒度并行版本,在较大的任务块上分配工作,适合任务量较大的计算

    OmegaPlus 是一个用于检测基因组数据中正选择痕迹的生物信息学工具。基于单倍型频率分布和重组事件检测选择信号,通过计算ω统计量(omega statistic)来识别选择位点。

    官方软件提供了几种不同的版本:

    OmegaPlus:顺序版本;

    OmegaPlus-F:细粒度并行版本,适合在在更小的任务单位上并行化;

    OmegaPlus-C:粗粒度并行版本,在较大的任务块上分配工作,适合任务量较大的计算场景;

    OmegaPlus-M:多粒度并行版本结合了细粒度和粗粒度的优点,通常在实际计算中性能表现最佳。

    示例:

    ./OmegaPlus-M -name TEST -input TEST.fasta -minwin 100 -maxwin 1000 -grid 10000 -threads 4

    -name: 设置输出文件的前缀,用于标记结果文件。

    -input: 指定输入文件(例如基因组序列文件)。

    -minwin 和 -maxwin: 设置扫描窗口的范围。

    -grid: 决定基因组上划分的网格数量,值越大分辨率越高,但计算量也越大。

    -threads: 指定线程数量(仅用于并行版本)

    其它可选参数:

    -name <STRING>

    指定运行名称,生成的输出文件将以该名称作为前缀。

    -input <STRING>

    输入对齐文件的名称。

    -grid <INTEGER>

    指定在对齐中计算 omega 值的网格数,数值越大分辨率越高,但计算量更大。

    -minwin <INTEGER>/-maxwin <INTEGER>

    设置用于计算 SNPs 之间连锁不平衡值的最小和最大窗口大小。

    -length <INTEGER>

    指定对齐长度,仅适用于 MS-like 和 MaCS-like 输入文件。

    -impute <N or GAP>

    启用对N或GAP符号的随机插补,将它们转换为有效字符。

    可以使用两次-impute分别处理N和GAP。

    示例:-impute N -impute GAP。

    -seed <INTEGER>

    指定随机数种子,用于插补N和GAP符号,或将 DNA 对齐转化为二进制数据。

    必须与-impute和-binary一起使用。

    -binary

    将 DNA 对齐数据转换为二进制格式。

    -threads <INTEGER>

    指定线程数,仅适用于并行版本 OmegaPlus-F、OmegaPlus-C 和 OmegaPlus-M。

    -minsnps <INTEGER>

    指定子区域内计算 omega 值的最小 SNP 数量。

    示例:-minsnps 10。

    -ld <STRING>

    指定连锁不平衡的测量类型,支持以下选项:

        RSQUARE(默认): SNPs 之间的相关系数 r²。

        D: D 统计量。

        ABSD: D 的绝对值。

        DOM: D_omega 统计量。

        ABSDOM: D_omega 的绝对值。

        ABSDOM2: D_omega 的绝对值,按频率乘积归一化。

    示例:-ld RSQUARE。

    -maf <FLOAT>

    排除次等位基因频率低于阈值的 SNP。

    -no-singletons

    从分析中排除单一出现的 SNPs。

    -b <INTEGER>

    允许每个位置在左、右窗口中 SNPs 数目差异不超过<INTEGER>,实现对最佳 omega 值的近似搜索。

    -all

    输出文件中显示更多结果信息。

    -reports

    为每个对齐生成单独的报告文件。

    -sampleList <STRING>

    用于 VCF 文件,指定需要包含的样本。

    -sampleList_out <STRING>

    为输入 VCF 文件生成样本列表。

    -mbs

    指定输入文件为 mbs 格式(将按 ms 格式处理)。

    -noSeparator

    禁止在输出中打印//分隔符,便于元处理(如在 R 中处理结果)。



    案例:




    Weng et al., 2024, Mol Biol Evol




    Taliadoros et al., 2023, PLOS Genet




    Sørensen et al., 2023, Sci

     
    (文/匿名(若涉版权问题请联系我们核实发布者) / 非法信息举报 / 删稿)
    打赏
    免责声明
    • 
    本文为昵称为 7463**** 发布的作品,本文仅代表发布者个人观点,本站未对其内容进行核实,请读者仅做参考,如若文中涉及有违公德、触犯法律的内容,一经发现,立即删除,发布者需自行承担相应责任。涉及到版权或其他问题,请及时联系我们154208694@qq.com删除,我们积极做(权利人与发布者之间的调停者)中立处理。郑重说明:不 违规举报 视为放弃权利,本站不承担任何责任!
    有个别老鼠屎以营利为目的遇到侵权情况但不联系本站或自己发布违规信息然后直接向本站索取高额赔偿等情况,本站一概以诈骗报警处理,曾经有1例诈骗分子已经绳之以法,本站本着公平公正的原则,若遇 违规举报 我们100%在3个工作日内处理!
    0相关评论
     

    (c)2008-现在 cuogai.com All Rights Reserved.