Hello,
I am trying to reproduce the results from the code you provide in the tutorial . However, I end up with an error that I am unable to resolve, at the step of designing all possible allele-specific guides. It appears that using the exact same input files that you had used there are no guides that meet the chosen criteria (whereas the output you have linked to there are 5 possible guides). Would you be able to look into it?
Thanks,
-Reuben
% python3 ./AlleleAnalyzer_bcft_RT/scripts/gen_sgRNAs.py
wtc_phased_hg19.bcf
mfn2_hg19_annots.h5
1:12040238-12073572
hg19_pams
hg19_pams/chr1.fa
test_sgrnas
SpCas9
20
-v
[2023-01-02 01:00:51,990 root:INFO ]bcftools version 1.16 running
[2023-01-02 01:00:51,991 root:INFO ]{'--bed': False,
'--cas-list': False,
'--crispor': None,
'--help': False,
'--hom': False,
'--max_indel': '5',
'--min_score': '0',
'--ref_guides': False,
'--sim': False,
'--strict': False,
'-c': False,
'-d': False,
'-r': False,
'-v': True,
'<annots_file>': 'mfn2_hg19_annots.h5',
'': 'wtc_phased_hg19.bcf',
'<cas_types>': 'SpCas9',
'<gene_vars>': None,
'<guide_length>': '20',
'': '1:12040238-12073572',
'': 'test_sgrnas',
'<pams_dir>': 'hg19_pams',
'<ref_fasta>': 'hg19_pams/chr1.fa'}
[2023-01-02 01:00:51,991 root:INFO ]Finding allele-specific guides.
[2023-01-02 01:00:52,183 root:INFO ]There are 3 heterozygous variants in this locus in this genome.
[2023-01-02 01:00:56,691 root:INFO ]Currently evaluating SpCas9.
[2023-01-02 01:00:56,709 root:INFO ]No sgRNAs meet the criteria for this locus, exiting.
Traceback (most recent call last):
File "./AlleleAnalyzer_bcft_RT/scripts/gen_sgRNAs.py", line 1715, in
main(arguments)
File "./AlleleAnalyzer_bcft_RT/scripts/gen_sgRNAs.py", line 1678, in main
out = get_allele_spec_guides(args).query('variant_position_in_guide > -1')
AttributeError: 'NoneType' object has no attribute 'query'