怎么对多个文件进行匹配

yestreenstars

回复 10# blackold
嗯，我刚回复完也看到了~

   

blackold

@yestreenstars 谢谢提醒。

huang6894

@blackold@yestreenstars谢谢两位大师

huang6894

回复 8# yestreenstars


    辰哥，麻烦帮忙看一下，我哪里错了~

1. => cat XAEDAAAPEI-56_2.fq
   
2. FCC36J1ACXX:1:1101:1346:13037#GATTAGTG/2        90      79      89      iPE-3+  34      0       4       11      2
   
3. FCC36J1ACXX:1:1101:1346:13020#GATTAGTG/2       90      79      89      iPE-3+  34      0       4      11      2
   
4. 
   
5. => cat 56_2.fq
   
6. @FCC36J1ACXX:1:1101:1346:13037#GATTAGTG/2
   
7. TATAAAATCTGTCTCTACTCGGTTTTAAATTTATTTATTTTTTATAAATAATAGTTGTACATATTTGGGGCTGTCTCTTATACACATCT
   
8. +
   
9. bbbeeeeegggggiiiiiihifghiiiiiiihiiifhiiiiifghhfe____ecghgfhigggggggeecccccccccccccdccccccd
   
10. @FCC36J1ACXX:1:1101:13880:10344#GATTAGTG/2
    
11. GTGCNACATCTTAACCTCAAGCACAGGAAGCAGTATAAAATCTGTCTCTTATACACATCTCTGAGCGGGCTGGCAAGGCAGACCGGGAGG
    
12. +
    
13. Z^_cBQ\`bcccedd[efg`edfgd_defhchhbbbffghbf[cO^^_e`gghhaace_fgdcfcb`gdcb``ab_^aaa^aXQXRTXX[
    
14. @FCC36J1ACXX:1:1101:1250:2201#GGAGGAAG/2
    
15. GGGGACGAGACCTCAACTCGGTTTTAAATTTATTTATTTTTTATAAATAATAGTTGTACATATTTGGGGCTGTCTCTTATACACATCTCT
    
16. +
    
17. bbbeeeeegggggiiiiiihifghiiiiiiihiiifhiiiiifghhfe____ecghgfhigggggggeecccccccccccccdccccccd
    
18. 
    
19. => cat do.txt
    
20. GCNACATCTT
    
21. ATAAAATCTGTGTGC
    
22. TATAGACGAGACCTCA
    
23. GGGGACGATATA
    
24. GGTATAGG
    
25. 
    
26. => awk 'ARGIND==1{a["@"$1]=$7" "$8}ARGIND==2{if(/@/){t=$0;s[t]=$0;if(a[$1])k=1;next}if(k){split(a[t],b);m[t]=substr($0,b[1]+1,b[2]-b[1]+1);k=0}s[t]=s[t]"\n"$0}ARGIND==3{for(i in m)if($0~"^"m[i]||$0~m[i]"$"){print s[i]}else{print "erro"}}' XAEDAAAPEI-56_2.fq 56_2.fq do.txt
    
27. erro
    
28. erro
    
29. erro
    
30. erro
    
31. erro
    
32. 
    

复制代码

按道理不是应该输出：

1. 
   
2. @FCC36J1ACXX:1:1101:1346:13037#GATTAGTG/2
   
3. TATAAAATCTGTCTCTACTCGGTTTTAAATTTATTTATTTTTTATAAATAATAGTTGTACATATTTGGGGCTGTCTCTTATACACATCT
   
4. +
   
5. bbbeeeeegggggiiiiiihifghiiiiiiihiiifhiiiiifghhfe____ecghgfhigggggggeecccccccccccccdccccccd
   
6. @FCC36J1ACXX:1:1101:1346:13037#GATTAGTG/2
   
7. TATAAAATCTGTCTCTACTCGGTTTTAAATTTATTTATTTTTTATAAATAATAGTTGTACATATTTGGGGCTGTCTCTTATACACATCT
   
8. +
   
9. bbbeeeeegggggiiiiiihifghiiiiiiihiiifhiiiiifghhfe____ecghgfhigggggggeecccccccccccccdccccccd
   

复制代码

吗？

yestreenstars

回复 14# huang6894
XAEDAAAPEI-56_2.fq和56_2.fq关联的是FCC36J1ACXX:1:1101:1346:13037#GATTAGTG/2这个key，而这个key对应的$7和$8是0和4，那么对应的起始位置和终止位置应该是第1位和第5位，在56_2.fq里对应的字符串是TATAA，而这个字符串在do.txt中并没有出现过。

   

huang6894

回复 15# yestreenstars


    啊哈~又粗心了~谢谢你