發(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í)際操作過程中...
-
重測(cè)序-篩去部分樣品并統(tǒng)一編號(hào)
篩取vcf中的樣品 samplelist_vcf.txt中按列排樣本名即可 重命名，再轉(zhuǎn)為bed文件 rename_snp.txt是一個(gè)兩列文件...
-
一個(gè)練習(xí)題
15-1 定義一個(gè)新的名為Circle的類表示圓形啃炸，它的屬性有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文件中提取堿基過濾信息（python實(shí)現(xiàn)）
先準(zhǔn)備文件路徑 有些代碼是可有可無的肿嘲，忘了刪了
-
基于barcode拆分二代序列
為了節(jié)約成本，研究者們通常會(huì)把多個(gè)樣品混在一個(gè)文庫筑公，并給不同樣品加上一段 Barcode 序列雳窟。在后續(xù)的生物信息分析中，根據(jù) Barcode 序...
-
二代群體遺傳與重測(cè)序（下）五.種群歷史與基因交流 書接上回http://www.reibang.com/p/0147afc3334c[https://www.jians...