求助_根据一个文件中的一列去另一个文件中提取相应的行

hechao7021

a.txt
LOC_Os01g01160
LOC_Os01g06454
LOC_Os01g17030
LOC_Os01g53020
b.txt
5        29142204        +        CHH        CTC        LOC_Os05g50810
5        29142215        +        CHH        CAA        LOC_Os05g50810
5        29142259        -        CHH        CTT        LOC_Os01g01160
5        29145615        +        CHH        CCA        LOC_Os01g01160
5        29145616        +        CHH        CAT        LOC_Os01g01160
5        29145869        -        CHH        CTA        LOC_Os05g50810
5        29145890        -        CHH        CAT        LOC_Os01g06454
5        29145912        -        CHH        CTT        LOC_Os01g06454
5        29145913        -        CHH        CCT        LOC_Os05g50810
5        29146309        +        CHH        CTA        LOC_Os01g17030
5        29146385        +        CHH        CAA        LOC_Os01g17030
5        29146426        +        CHH        CAA        LOC_Os05g50810
6        853213        -        CHH        CAC        LOC_Os01g53020
得到：
5        29142259        -        CHH        CTT        LOC_Os01g01160
5        29145615        +        CHH        CCA        LOC_Os01g01160
5        29145616        +        CHH        CAT        LOC_Os01g01160
5        29145890        -        CHH        CAT        LOC_Os01g06454
5        29145912        -        CHH        CTT        LOC_Os01g06454
5        29146309        +        CHH        CTA        LOC_Os01g17030
5        29146385        +        CHH        CAA        LOC_Os01g17030
6        853213        -        CHH        CAC        LOC_Os01g53020
自己弄的代码不知道怎么回事，总是不行，代码是：awk 'NR==FNR{n++;a[n]=$1}NR>FNR{if ($10 in a[n]) print $0}' a.txt b.txt      这样也不行：awk 'NR==FNR{a[$1]=$1}NR>FNR{if ($10 in a) print $0}' a.txt b.txt

jzsjm1002

回复 1# hechao7021

1. [root@vh100 awk]# awk 'NR==FNR{a[$0]} NR!=FNR{if($6 in a) print $0}' a b
   
2. 5        29142259        -        CHH        CTT        LOC_Os01g01160
   
3. 5        29145615        +        CHH        CCA        LOC_Os01g01160
   
4. 5        29145616        +        CHH        CAT        LOC_Os01g01160
   
5. 5        29145890        -        CHH        CAT        LOC_Os01g06454
   
6. 5        29145912        -        CHH        CTT        LOC_Os01g06454
   
7. 5        29146309        +        CHH        CTA        LOC_Os01g17030
   
8. 5        29146385        +        CHH        CAA        LOC_Os01g17030
   
9. 6        853213        -        CHH        CAC        LOC_Os01g53020

复制代码

Jackis

awk 'NR==FNR{file[$0]};NR>FNR{if($6 in file) print}' a.txt b.txt

hechao7021

回复 3# Jackis

非常感谢！不知道为什么，如果换成两个大文件，就运行不出结果

hechao7021

回复 5# hechao7021

我明白了，是不是第一个文件不能用.txt的后缀，我把后缀去掉就能输出结果了

caca0215

for i in `awk '{print $1}' a.txt`; do cat b.txt |grep $i ; done

REFRESHMM

你们忘记了 grep -f a.txt b.txt