shell还没入门,确切的说是还没学,用python3写了个小脚本统计23andme、gsa、wegene各染色体位点. 代码如下:
23andme:
gsa的同上,稍做一两处调整。 wegene的没有源数据,找的生信菜鸟团jimmy的shell
` cut -f 2 jimmy_wegene.txt | unique -c |grep -v “^#” ``
结果分析
从图中可以看出,毕竟23andme是最早做的,位点最多,但是,现在来讲,大约也就上千个点可以被解读,所以意义不大。但是,国内公司的还是引入了中国人的特征位点的,更具有人群优势,适合国人的基因检测,现阶段,大概也就这样了,不知何时可以每个人做全外或全基因检测。