【新手】怎么用perl写一个脚本比较两个基因文件里相同的语段再输出到新的文件里

酸奶盒

回复 40# moperyblue

awk -F'[\t]' 'NR==FNR{a[$2$3$1]=1;next}a[$2$3$4]'  GaA_GbA.txt GbA_GaA.txt > C
你给的这个表达式我可以用的 谢谢啦
我刚才说的是另外一个回答的写的那个脚本。。
不过awk看起来好难= =不知道我要学多久

。。还有四个需要匹配。。厚颜拜托OwO

moperyblue

回复 41# 酸奶盒

1. ==> GbA_GbD <==
   
2. gnl|AD2_NBI.cds|GbC30369598.1.0 AD2_NBI.cds:Gbscaffold19070.25.0        99.663  297     1       0       1       297     76      372     2.61e-154       544
   
3. 
   
4. ==> GbD_GbA <==
   
5. gnl|AD2_NBI.cds|GbC30374398.1.0 AD2_NBI.cds:Gbscaffold3437.17.0 99.822  562     1       0       1       562     1       562     0.0     1033

复制代码

这两个文件以tab分隔 是不是1和2列相同就算相同?

酸奶盒

回复 42# moperyblue

嗯，是的这两个文件是同一个物种 不同染色体组之间的相互比较

所有的相互比较的文件都是只要第一列与第二列相同就需要输出出来

moperyblue

回复 43# 酸奶盒

1. #GbD_GbA
   
2. awk 'NR==FNR{a[$1$2]=1;next}a[$1$2]' GbA_GbD GbD_GbA > GbD_GbA_C
   
3. 
   
4. #GhA_GaA
   
5. awk 'NR==FNR{a[$1$2]=1;next}a[$2$1]'  GaA_GhA GhA_GaA > GhA_GaA_C
   
6. 
   
7. #GhD_GhA
   
8. awk 'NR==FNR{a[$1$2]=1;next}a[$2$1]' GhA_GhD GhD_GhA > GhD_GhA_C
   
9. 
   
10. #GhD_GrD
    
11. awk 'NR==FNR{a[$1$2]=1;next}a[$2$1]' GhD_GrD GrD_GhD > GhD_GrD_C

复制代码

酸奶盒

回复 44# moperyblue

第一个没有输出出来orz刚才我话好像有点误解
我的意思是这两个文件也是需要GbA_GbD第一列等于GbD_GbA第二列同时GbA_GbD的第二列等于GbD_GbA的第一列才输出

剩下的都可以有

moperyblue

回复 45# 酸奶盒


从你提供的数据 文件GbA_GbD 与 文件GbD_GbA 它们第一和第二列的数据格式是一样的 你再仔细看一下 是不是

酸奶盒

回复 46# moperyblue

抱歉。。。QAQ是我没有仔细看。。。。啊啊啊orz

酸奶盒

回复 46# moperyblue
它们的数据格式是一样的。。但是代表的基因是不一样的，是对的、
因为他们是从同一个数据库里分开的，所以格式是一样的。。。

moperyblue

回复 48# 酸奶盒


按这个格式 GbA_GbD、GbD_GbA这两个文件没有找到“交集”的数据  正常吗?

酸奶盒

回复 49# moperyblue

抱歉，你的意思是说这两个文件是数据是一样的？O O