一个超简单的admixture结果排序脚本

发表于 分类于
用的是（本人学艺不精但觉得很能装逼又实用的）python，真的简单到令人发指（？）

刚开始用admixture做种群结构分析的时候还没开始学perl和python（虽然现在也没学会多少），每次做出结果都只会按vcf默认的顺序排列作图，做完了再用ai慢慢挪动条形图。学了编程以后发现竟然就这就能搞定了：

1
2
3
4
5
6
7
8
9
10
#!/usr/bin/python3
# -*- coding: UTF-8 -*-


q = open("output.Q","r")

line = q.readlines()
print(line[65],line[66],……,line[64]) #这里用的其实还是笨办法，手动把样品顺序全部罗列一遍。理论上也可以用字符串自动比对来实现

q.close()