23andme、gsa、wegene各染色体位点统计

zd00572

Posted by zd200572 on July 20, 2017

shell还没入门,确切的说是还没学,用python3写了个小脚本统计23andme、gsa、wegene各染色体位点. 代码如下:

23andme:

gsa的同上,稍做一两处调整。 wegene的没有源数据,找的生信菜鸟团jimmy的shell

` cut -f 2 jimmy_wegene.txt | unique -c |grep -v “^#” ``

结果分析

从图中可以看出,毕竟23andme是最早做的,位点最多,但是,现在来讲,大约也就上千个点可以被解读,所以意义不大。但是,国内公司的还是引入了中国人的特征位点的,更具有人群优势,适合国人的基因检测,现阶段,大概也就这样了,不知何时可以每个人做全外或全基因检测。