Entering edit mode
10.5 years ago
alex
▴
250
Getting back an error in seqdb in pygr
no key 'chr9' in database <SequenceFileDB '/local/resources/hg18/karyotypic/hg18.fa'>
While it worked fine for chr15, what is the error mean, that it can't find the chromosome in the assembly? Thanks!
So I believe this was due to a corrupt assembly file (i.e. hg18.fa). I rebuilt it and it worked.
print all keys and see what is going on
print db.keys()