發(fā)簡(jiǎn)信
IP屬地：北京
-
群體結(jié)構(gòu)分析的繪圖腳本，基于sns腳本 首先導(dǎo)入需要用到的包祭隔。--order和-n都是用于指定順序的秃诵，有一個(gè)參數(shù)就可以了鲁僚。其他都是必選參數(shù)蝴罪。 核心的排序函數(shù)和繪圖函數(shù) 讀入數(shù)據(jù)钠四，...
-
尋找SV.vcf文件中是否含有某個(gè)片段的插入現(xiàn)在已知這個(gè)片段在3號(hào)染色體上幔嫂，基本思路是先將vcf文件變成fa文件驴党，然后用blastn處理瘪撇。 第四列記錄的是INS的序列內(nèi)容 在實(shí)際操作過(guò)程中...
-
重測(cè)序-篩去部分樣品并統(tǒng)一編號(hào)
篩取vcf中的樣品 samplelist_vcf.txt中按列排樣本名即可 重命名，再轉(zhuǎn)為bed文件 rename_snp.txt是一個(gè)兩列文件...
-
一個(gè)練習(xí)題
15-1 定義一個(gè)新的名為Circle的類(lèi)表示圓形，它的屬性有center和radius倔既，其中center是一個(gè)Point對(duì)象恕曲，而radius是...
-
GWAS中的曼哈頓圖與QQ圖（CMplot）
輸入文件為exmma的輸出ps文件，整理為ps2文件渤涌，格式如下 調(diào)用命令為
-
基于sort和uniq的集合操作（簡(jiǎn)單明了）
sort:用于將文本文件內(nèi)容加以排序 uniq:刪除文件中重復(fù)的行码俩，得到文件中唯一的行 交集 （用sort將a.txt b.txt文件進(jìn)行排序，...
-
從質(zhì)控log文件中提取堿基過(guò)濾信息（python實(shí)現(xiàn)）
先準(zhǔn)備文件路徑 有些代碼是可有可無(wú)的歼捏，忘了刪了
-
基于barcode拆分二代序列
為了節(jié)約成本稿存，研究者們通常會(huì)把多個(gè)樣品混在一個(gè)文庫(kù)，并給不同樣品加上一段 Barcode 序列瞳秽。在后續(xù)的生物信息分析中瓣履，根據(jù) Barcode 序...
-
二代群體遺傳與重測(cè)序（下）五.種群歷史與基因交流 書(shū)接上回http://www.reibang.com/p/0147afc3334c[https://www.jians...