-
Notifications
You must be signed in to change notification settings - Fork 10
Open
Description
I have an error that I can't seem to track down. Running straingr 1.3.3 (bwa 0.7.17) with something like:
straingr call refs_concat.fasta 1034_Staph_Mockcomm_5_acvrk.bam --hdf5-out 1034_Staph_Mockcomm_5_acvrk.out.hdf5 --summary 1034_Staph_Mockcomm_5_acvrk.out.tsv --tracks all`
errors with:
2024-02-28 17:23:46,908 - INFO:root:Loading reference refs_concat.fasta...
2024-02-28 17:23:46,934 - INFO:root:Reference length: 5300120
2024-02-28 17:23:46,942 - INFO:root:Start analyzing aligned reads...
Traceback (most recent call last):
File "/data/conlans/mymamba/envs/strainge/bin/straingr", line 11, in <module>
sys.exit(straingr_cli())
File "/data/conlans/mymamba/envs/strainge/lib/python3.7/site-packages/strainge/cli/main.py", line 110, in __call__
self.run(args)
File "/data/conlans/mymamba/envs/strainge/lib/python3.7/site-packages/strainge/cli/registry.py", line 83, in run
rc = subcommand_func(**args_dict)
File "/data/conlans/mymamba/envs/strainge/lib/python3.7/site-packages/strainge/cli/straingr.py", line 468, in __call__
call_data = caller.process(reference, sample_bam)
File "/data/conlans/mymamba/envs/strainge/lib/python3.7/site-packages/strainge/variant_caller.py", line 992, in process
self._assess_allele(call_data, scaffold, refpos, read)
File "/data/conlans/mymamba/envs/strainge/lib/python3.7/site-packages/strainge/variant_caller.py", line 1080, in _assess_allele
call_data.bad_read(scaffold, refpos)
AttributeError: 'VariantCallData' object has no attribute 'bad_read'
I've looked through the code and can't figure out where this bad_read attr is coming from. The only thing I'm doing differently from the tutorial is pointing --hdf5-out to a new hdf5 file instead of the per sample hdf5 files generated by straingst kmerize in an earlier stage of the readthedocs walkthrough, but I don't think that's it.
Reactions are currently unavailable
Metadata
Metadata
Assignees
Labels
No labels