OmegaPlus 檢測正選擇信號

OmegaPlus?是一個用于檢測基因組數(shù)據(jù)中正選擇痕跡的生物信息學(xué)工具萧求∪糇ǎ基于單倍型頻率分布和重組事件檢測選擇信號,通過計算ω統(tǒng)計量(omega statistic)來識別選擇位點(diǎn)丙挽。

官方軟件提供了幾種不同的版本:

OmegaPlus:順序版本毕谴;

OmegaPlus-F:細(xì)粒度并行版本,適合在在更小的任務(wù)單位上并行化庭惜;

OmegaPlus-C:粗粒度并行版本,在較大的任務(wù)塊上分配工作穗酥,適合任務(wù)量較大的計算場景护赊;

OmegaPlus-M:多粒度并行版本結(jié)合了細(xì)粒度和粗粒度的優(yōu)點(diǎn),通常在實(shí)際計算中性能表現(xiàn)最佳砾跃。

示例:

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

-name: 設(shè)置輸出文件的前綴骏啰,用于標(biāo)記結(jié)果文件。

-input: 指定輸入文件(例如基因組序列文件)抽高。

-minwin 和 -maxwin: 設(shè)置掃描窗口的范圍判耕。

-grid: 決定基因組上劃分的網(wǎng)格數(shù)量,值越大分辨率越高翘骂,但計算量也越大壁熄。

-threads: 指定線程數(shù)量(僅用于并行版本)

其它可選參數(shù):

-name <STRING>

指定運(yùn)行名稱帚豪,生成的輸出文件將以該名稱作為前綴。

-input <STRING>

輸入對齊文件的名稱草丧。

-grid <INTEGER>

指定在對齊中計算 omega 值的網(wǎng)格數(shù)狸臣,數(shù)值越大分辨率越高,但計算量更大昌执。

-minwin <INTEGER>/-maxwin <INTEGER>

設(shè)置用于計算 SNPs 之間連鎖不平衡值的最小和最大窗口大小烛亦。

-length <INTEGER>

指定對齊長度,僅適用于 MS-like 和 MaCS-like 輸入文件懂拾。

-impute <N or GAP>

啟用對N或GAP符號的隨機(jī)插補(bǔ)煤禽,將它們轉(zhuǎn)換為有效字符。

可以使用兩次-impute分別處理N和GAP岖赋。

示例:-impute N -impute GAP檬果。

-seed <INTEGER>

指定隨機(jī)數(shù)種子,用于插補(bǔ)N和GAP符號贾节,或?qū)?DNA 對齊轉(zhuǎn)化為二進(jìn)制數(shù)據(jù)汁汗。

必須與-impute和-binary一起使用。

-binary

將 DNA 對齊數(shù)據(jù)轉(zhuǎn)換為二進(jìn)制格式栗涂。

-threads <INTEGER>

指定線程數(shù)知牌,僅適用于并行版本 OmegaPlus-F、OmegaPlus-C 和 OmegaPlus-M斤程。

-minsnps <INTEGER>

指定子區(qū)域內(nèi)計算 omega 值的最小 SNP 數(shù)量角寸。

示例:-minsnps 10。

-ld <STRING>

指定連鎖不平衡的測量類型忿墅,支持以下選項:

????RSQUARE(默認(rèn)): SNPs 之間的相關(guān)系數(shù) r2扁藕。

????D: D 統(tǒng)計量。

????ABSD: D 的絕對值疚脐。

????DOM: D_omega 統(tǒng)計量亿柑。

????ABSDOM: D_omega 的絕對值。

????ABSDOM2: D_omega 的絕對值棍弄,按頻率乘積歸一化望薄。

示例:-ld RSQUARE。

-maf <FLOAT>

排除次等位基因頻率低于閾值的 SNP呼畸。

-no-singletons

從分析中排除單一出現(xiàn)的 SNPs痕支。

-b <INTEGER>

允許每個位置在左、右窗口中 SNPs 數(shù)目差異不超過<INTEGER>蛮原,實(shí)現(xiàn)對最佳 omega 值的近似搜索卧须。

-all

輸出文件中顯示更多結(jié)果信息。

-reports

為每個對齊生成單獨(dú)的報告文件。

-sampleList <STRING>

用于 VCF 文件花嘶,指定需要包含的樣本笋籽。

-sampleList_out <STRING>

為輸入 VCF 文件生成樣本列表。

-mbs

指定輸入文件為 mbs 格式(將按 ms 格式處理)察绷。

-noSeparator

禁止在輸出中打印//分隔符干签,便于元處理(如在 R 中處理結(jié)果)。



案例:

Weng et al., 2024, Mol Biol Evol
Taliadoros et al., 2023, PLOS Genet
S?rensen et al., 2023, Sci
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末拆撼,一起剝皮案震驚了整個濱河市容劳,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌闸度,老刑警劉巖竭贩,帶你破解...
    沈念sama閱讀 219,539評論 6 508
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異莺禁,居然都是意外死亡留量,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,594評論 3 396
  • 文/潘曉璐 我一進(jìn)店門哟冬,熙熙樓的掌柜王于貴愁眉苦臉地迎上來楼熄,“玉大人,你說我怎么就攤上這事浩峡】善瘢” “怎么了?”我有些...
    開封第一講書人閱讀 165,871評論 0 356
  • 文/不壞的土叔 我叫張陵翰灾,是天一觀的道長缕粹。 經(jīng)常有香客問我,道長纸淮,這世上最難降的妖魔是什么平斩? 我笑而不...
    開封第一講書人閱讀 58,963評論 1 295
  • 正文 為了忘掉前任,我火速辦了婚禮咽块,結(jié)果婚禮上绘面,老公的妹妹穿的比我還像新娘。我一直安慰自己侈沪,他們只是感情好飒货,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,984評論 6 393
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著峭竣,像睡著了一般。 火紅的嫁衣襯著肌膚如雪晃虫。 梳的紋絲不亂的頭發(fā)上皆撩,一...
    開封第一講書人閱讀 51,763評論 1 307
  • 那天,我揣著相機(jī)與錄音,去河邊找鬼扛吞。 笑死呻惕,一個胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的滥比。 我是一名探鬼主播亚脆,決...
    沈念sama閱讀 40,468評論 3 420
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼盲泛!你這毒婦竟也來了濒持?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 39,357評論 0 276
  • 序言:老撾萬榮一對情侶失蹤寺滚,失蹤者是張志新(化名)和其女友劉穎柑营,沒想到半個月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體村视,經(jīng)...
    沈念sama閱讀 45,850評論 1 317
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡官套,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,002評論 3 338
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了蚁孔。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片奶赔。...
    茶點(diǎn)故事閱讀 40,144評論 1 351
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖杠氢,靈堂內(nèi)的尸體忽然破棺而出站刑,到底是詐尸還是另有隱情,我是刑警寧澤修然,帶...
    沈念sama閱讀 35,823評論 5 346
  • 正文 年R本政府宣布笛钝,位于F島的核電站,受9級特大地震影響愕宋,放射性物質(zhì)發(fā)生泄漏玻靡。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,483評論 3 331
  • 文/蒙蒙 一中贝、第九天 我趴在偏房一處隱蔽的房頂上張望囤捻。 院中可真熱鬧,春花似錦邻寿、人聲如沸蝎土。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,026評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽誊涯。三九已至,卻和暖如春蒜撮,著一層夾襖步出監(jiān)牢的瞬間暴构,已是汗流浹背跪呈。 一陣腳步聲響...
    開封第一講書人閱讀 33,150評論 1 272
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留取逾,地道東北人耗绿。 一個月前我還...
    沈念sama閱讀 48,415評論 3 373
  • 正文 我出身青樓,卻偏偏與公主長得像砾隅,于是被迫代替她去往敵國和親误阻。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,092評論 2 355

推薦閱讀更多精彩內(nèi)容