We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent eb6cd37 commit 2a3defdCopy full SHA for 2a3defd
capy/seq.py
@@ -1,6 +1,7 @@
1
from pyfaidx import Fasta as _Fasta, FastaNotFoundError as _FastaNotFoundError
2
import numpy as _np
3
import sys as _sys
4
+import os as _os
5
import mmap as _mmap
6
import pandas as _pd
7
@@ -61,7 +62,7 @@ def __set_reference(self, ref):
61
62
print("Could not load reference genome!")
63
64
# export public functions
-_fa = _FA();
65
+_fa = _FA(ref_fa_file = _os.environ["CAPY_REF_FA"] if "CAPY_REF_FA" in _os.environ else "/mnt/j/db/hg19/ref/hs37d5.fa");
66
genome_region = _fa._genome_region
67
chrpos2gpos = _fa._chrpos2gpos
68
gpos2chrpos = _fa._gpos2chrpos
0 commit comments