Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
newreference: If gene is given but not found in genbank, error out
Currently, if the --gene option is used and the gene name is not found, the script will use the entire genome which may cause some silent undesired behaviors. This commit changes that such that the script will error out if the gene is not found in the GenBank file as this indicates the gene name may be misspelled or the user may be using the wrong GenBank file. If the --gene option is not used, the script will continue to process the entire genome as expected.
- Loading branch information
Thanks for adding this, @j23414! I'd move this conditional immediately after the for loop above, just to avoid trying to extract the sequence feature when you're going to exit with an error anyway.