已有“OrthoFinder/Results_Jul11/Orthogroups ”
ll Orthogroups
-rw-r--r-- 1 552K Jul 11 01:15 Orthogroups.GeneCount.tsv
-rw-r--r-- 1 36K Jul 11 01:15 Orthogroups_SingleCopyOrthologues.txt
-rw-r--r-- 1 8.1M Jul 11 01:15 Orthogroups.tsv
-rw-r--r-- 1 9.3M Jul 11 01:15 Orthogroups.txt
-rw-r--r-- 1 1.6M Jul 11 01:15 Orthogroups_UnassignedGenes.tsv
#随后整理
#https://www.biostars.org/p/323537/
grep -Fwf SingleCopyOrthogroups.txt Orthogroups.txt > SingleCopyOrthologsGeneIDS.txt
#https://www.biostars.org/p/274269/
再进行提取ID
i=$(cat Orthogroups_SingleCopyOrthologues.txt | xargs )
for j in $i; do grep "$j" SingleCopyOrthologsGeneIDS.txt | awk '{for (i=2; i<=NF; i++) print $i}'> $j.txt; done;
网友评论