Skip to content

Instantly share code, notes, and snippets.

@mdshw5
Created January 7, 2018 04:02
Show Gist options
  • Save mdshw5/ad425b28807dc55744434ccf34a19e3f to your computer and use it in GitHub Desktop.
Save mdshw5/ad425b28807dc55744434ccf34a19e3f to your computer and use it in GitHub Desktop.
Minimal example for mdshw5/pyfaidx #125
from Bio import bgzf
with bgzf.open('Homo_sapiens.GRCh38.dna.primary_assembly.fa.gz') as fasta:
fasta.tell()
for line in fasta:
if line.startswith('>'):
fasta.tell()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment