python3で指定行と特定行数先の抽出

前提・実現したいこと

ゲノム解析で大きなblastデータ(.txt)を扱っています。
例えば以下のようなデータがあります。

Query= utg000001

Length=2047151
Score E
Sequences producing significant alignments: (Bits) Value

emb|LK935835.1| Babesia divergens genome assembly 454hybrid_PBje... 3546 0.0
gb|GU447229.1| Babesia gibsoni strain TWN2 thrombospondin-relate... 3384 0.0
gb|GU447230.1| Babesia gibsoni strain TWN3 thrombospondin-relate... 3373 0.0
gb|GU447231.1| Babesia gibsoni strain TWN4 thrombospondin-relate... 3277 0.0
・
・
・

ここで"Query="を含む行と”emb|LK935835.1| Babesia divergens genome assembly 454hybrid_PBje..."の2行のみを抽出したいです。

"Query="を含む行を1行目として空行を含めて1,7行目のみを抽出すにはどうしたらいいでしょうか。

発生している問題・エラーメッセージ

エラーメッセージ

該当のソースコード

python3
for line in open("blastn.txt"):
if "Query=" in line:
print(line[:-1])

試したこと

上のコードでQuery=を含む行の抽出はできるのですが、そこから7行目の指定の仕方がわかりません。
よろしくお願いします。

補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

行動規範の内容に同意します

回答1件

ベストアンサー

直接的な方法ならば、例えばこんなふうに。

Python
1with open(r'blastn.txt') as fin:
2    for line in fin:
3        if 'Query=' in line:   # line.startswith('Query=') の方が良い
4            print(line[:-1])
5            break
6            
7    for _ in range(5):
8        next(fin)
9        
10    print(next(fin)[:-1])

ちょっと捻っても書けます。参考までに。

Python
1import collections
2
3queue = collections.deque(maxlen=7)
4for _ in range(queue.maxlen):
5    queue.append('')
6    
7with open(r'blastn.txt') as fin:
8    for line in map(str.rstrip, fin):
9        queue.append(line)
10       
11        if queue[0].startswith('Query='):
12            break
13    else:
14        raise ValueError
15        
16q, *_, emb = queue
17print(q)
18print(emb)

コメントを受けて

もし"Query="が複数ある場合でそのすべてを同じように抽出するときはどう書けばいいでしょうか。

Python
1with open(r'blastn.txt') as fin:
2    while True:
3        for line in fin:
4            if 'Query=' in line:   # line.startswith('Query=') の方が良い
5                print(line[:-1])
6                break
7        else:
8            break
9
10        for _ in range(5):
11            next(fin)
12
13        print(next(fin)[:-1])