Bug Description
callAltTranslation has an issue that makes it unusable in v1.4.6:
-
Tool requires at least one flag:
--selenocysteine-termination OR
--w2f-reassignment
-
Without flags: ValueError: At least one of --selenocysteine-termination and --w2f-reassignment must be given.
-
With either flag: Tool crashes after ~1h runtime:
ValueError: ''*'' is not a valid unambiguous letter for protein
at Bio/SeqUtils/__init__.py:474 in molecular_weight()
Root Cause
Both flags generate peptides containing stop codons (*), which BioPython's
molecular_weight() function cannot handle when checking --min-mw threshold.
Environment
- moPepGen: v1.4.6
- Python: 3.9
- Reference: Ensembl 114
- BioPython: (via moPepGen dependencies)
Expected Behavior
Tool should either:
- Filter out stop-codon-containing peptides before molecular weight check, OR
- Skip molecular weight validation for these specific peptide types, OR
- Make flags truly optional (don't require at least one)
Workaround
Currently none - callAltTranslation is completely unusable.
Steps to Reproduce
moPepGen callAltTranslation \
--index-dir /path/to/index \
--reference-source ENSEMBL \
--output-path output.fasta \
--selenocysteine-termination \
--cleavage-rule trypsin \
--miscleavage 2 \
--min-length 7 \
--max-length 25 \
--min-mw 500.0
Crashes after ~1 hour at peptide validation stage.
Bug Description
callAltTranslationhas an issue that makes it unusable in v1.4.6:Tool requires at least one flag:
--selenocysteine-terminationOR--w2f-reassignmentWithout flags:
ValueError: At least one of --selenocysteine-termination and --w2f-reassignment must be given.With either flag: Tool crashes after ~1h runtime:
at
Bio/SeqUtils/__init__.py:474inmolecular_weight()Root Cause
Both flags generate peptides containing stop codons (*), which BioPython's
molecular_weight()function cannot handle when checking--min-mwthreshold.Environment
Expected Behavior
Tool should either:
Workaround
Currently none -
callAltTranslationis completely unusable.Steps to Reproduce
Crashes after ~1 hour at peptide validation stage.