File: biopython_fastq_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 (9 lines) | stat: -rw-r--r-- 175 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
 
import sys
from Bio import SeqIO
fq = SeqIO.index_db("{}.db".format(sys.argv[2]))

with open(sys.argv[1]) as fh:
	for line in fh:
		name = line.strip()
		print(fq[name].seq)