求助：如何grep一个文件中含有关键词“OTU_”之后第八行的信息

ottchen

大家好！我刚解除Linux，花了一个下午的时间才明白Grep命令的强大，但是还有很多不懂。如下：
Query= OTU_527

Length=134


                                                                   Score     E
Sequences producing significant alignments:                       (Bits)  Value

pdb|4UQ4|1  Chain 1, Mammalian 80s Hcv-ires Initiation Complex...    248   1e-62
pdb|4W23|2  Chain 2, Structure Of The 80s Mammalian Ribosome B...    248   1e-62
ref|XM_008252668.1|  PREDICTED: Oryctolagus cuniculus serine/a...    248   1e-62
ref|NR_046233.2|  Mus musculus 45S pre-ribosomal RNA (Rn45s), ...    248   1e-62
gb|KC876030.1|  Homo sapiens clone LA13_165F6 sequence               248   1e-62
gb|KC876029.1|  Homo sapiens clone LA14_101B3 sequence               248   1e-62
gb|KC876028.1|  Homo sapiens clone LA15_25H3 sequence                248   1e-62
ref|XR_233610.1|  PREDICTED: Bos taurus uncharacterized LOC101...    248   1e-62
tpg|BK000964.3|  TPA: Mus musculus ribosomal DNA, complete rep...    248   1e-62
ref|XM_004474322.1|  PREDICTED: Dasypus novemcinctus translati...    248   1e-62
pdb|3J3D|2  Chain 2, Structure Of The H. Sapiens 40s Rrna And ...    248   1e-62
gb|JF894380.1|  Crocuta crocuta isolate CC8 18S ribosomal RNA ...    248   1e-62
gb|JQ247698.1|  Mus musculus cell-line Neuro-2a 18S ribosomal ...    248   1e-62
ref|NR_046261.1|  Sus scrofa 18S ribosomal

Query= OTU_149

Length=131


                                                                   Score     E
Sequences producing significant alignments:                       (Bits)  Value

gb|KJ763281.1|  Uncultured eukaryote clone SGUH718 18S ribosom...    215   1e-52
gb|KJ763240.1|  Uncultured eukaryote clone SGUH649 18S ribosom...    215   1e-52
gb|KJ763193.1|  Uncultured eukaryote clone SGUH584 18S ribosom...    215   1e-52
gb|KJ762994.1|  Uncultured eukaryote clone SGUH1374 18S riboso...    215   1e-52
gb|JN048122.1|  Uncultured eukaryote clone BF9_20C1E 18S ribos...    169   1e-38
emb|FR874449.1|  uncultured marine picoeukaryote partial 18S r...    169   1e-38
gb|DQ314819.1|  Uncultured marine eukaryote clone NW617.36 18S...    169   1e-38
gb|DQ314818.1|  Uncultured marine eukaryote clone NOR46.02 18S...    169   1e-38
gb|DQ314817.1|  Uncultured marine eukaryote clone NOR26.21 18S...    169   1e-38
gb|DQ369016.1|  Uncultured marine eukaryote clone UEPACAHp3 18...    163   5e-37
gb|DQ369015.1|  Uncultured marine eukaryote clone UEPACLp3 18S...    163   5e-37
gb|EU371163.1|  Uncultured marine eukaryote clone NPK2_123 18S...    150

我想要的输出结果是：
Query= OTU_527
pdb|4UQ4|1  Chain 1, Mammalian 80s Hcv-ires Initiation Complex...    248   1e-62
Query= OTU_149
gb|KJ763281.1|  Uncultured eukaryote clone SGUH718 18S ribosom...    215   1e-52

如果使用grep -A 8 OTU_ /home/*.txt 就会把中间6行的无用信息也抓出来。
所以，我的问题是，?如何grep一个文件中含有关键词“OTU_”之后指定第N行的信息。
谢谢！

jy02150795

同新手，用sed删掉中间6行？

cu_shell

grep -A 8 就是输出关键词后边的1-8行。至于指定第8行，，这个还得等大神来解决哈。。我也不会。。

qq58945591

用awk 比较方便，grep 不用管道我也不知道怎么搞了。

awk 方法
awk '/Query/{print;NR=1;next}{if(NR==9)print $0}' file
Query= OTU_527
pdb|4UQ4|1  Chain 1, Mammalian 80s Hcv-ires Initiation Complex...    248   1e-62
Query= OTU_149
gb|KJ763281.1|  Uncultured eukaryote clone SGUH718 18S ribosom...    215   1e-52

ottchen

我试了，怎么不行啊，结果是：把原文件的内容重新Print一遍，却不是你刚才显示的结果。 回复 4# qq58945591


   

ottchen

再试一次，可以了，太感谢大神您了！不过，这语句怎么看，能解释一下吗？ 回复 4# qq58945591


   

qq58945591

回复 6# ottchen


    我awk水平很水，简单点处理会，复杂的就蛋疼了。

上面命令是这样。
先找出Query 这一行，然后打印出来，再把行号设置为1，继续处理下一行，然后判断如果到了第9行的话就把第9行给打印出来。

ottchen

哈哈，很有才，简单实用。谢谢啦！

qq58945591 发表于 2014-12-14 19:44
回复 6# ottchen

ottchen

再打扰您最后一次，若我需要第八行和第十行的信息，该怎么改？谢谢！ 回复 4# qq58945591


   

lbseraph

回复 1# ottchen

如果不要Query那一行的话，加管道可以拿到第8行：
$ grep -A 8 OTU_527 <file> | tail -n 1