File: biopython_fasta_random_access.py

package info (click to toggle)
pyfastx 2.2.0-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 1,336 kB
  • sloc: ansic: 4,820; python: 1,817; sh: 505; perl: 66; makefile: 31
file content (12 lines) | stat: -rw-r--r-- 206 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
 
import sys
from Bio import SeqIO

idfile, fafile = sys.argv[1:]

fa = SeqIO.index_db("{}.db".format(fafile))

with open(idfile) as fh:
	for line in fh:
		seqid = line.strip()
		s = fa[seqid].seq
		print(s)