2017-09-25 4 views
0
におけるFASTAファイルから最初のレコード

を削除します。これは大きなPythonスクリプトの一部です。そのseqで作業して興味深い部分を抽出したら、ファイルから削除したいと思います。最終的に、ファイルは空になりますので、フォルダから削除できます。 私はBiopythonを使用していたので、私はBiopython内でこの機能を見つけることができるかどうかを確認しようとしていましたが、fastaファイルからレコードを削除する方法についてのヒントは見つかりませんでした。 SeqIOには一種のリストがあるので、私はlist.remove(list[0])メソッドで試しましたが、うまくいきません。これについてのアイデア?私は誰かが私にスクリプトを投稿するように頼む場合、スクリプトは約200行です。おかげさまで それは、長いシーケンスも、我々はメモリにそれらのすべてを読んで回避したいのではなく、1つずつを処理する多数の配列には最適ではないのですは、私は次の形式で小さなFASTAファイルを持っているのPython

import sys 
from Bio import SeqIO 

filename = sys.argv[1] 

sequences = [] 

for not_first, record in enumerate(SeqIO.parse(filename, 'fasta')): 
    if not_first: 
     sequences.append(record) 

SeqIO.write(sequences, 'truncated-' + filename, 'fasta') 

答えて

0

は、ここで簡単な何かあなたが始めるためにです。しかし、それは始める場所です。

+0

Yepp私はこれらの用語について考えていました。私はbiopythonがこの問題を助ける何かを持っているかもしれないと思った。私はリスト(SeqIO.parse)を使用しますが、リストにあるすべての機能を持っているわけではありません。ありがとう、あなたの助けにたくさん。 – Ana

関連する問題