This is a test version of Biostars. For the public version, visit https://www.biostars.org.
genbank output incomplete with long feature name

I want to insert a new feature to an existing genbank record. E.g.

new_feature = SeqFeature(FeatureLocation(1, 100), type="integration_excision_site", strand=1)
gb_record.features.append(new_feature)
SeqIO.write(gb_record, gb_file, "genbank")

In the output gb file, this new feature is shown as

integration_excicomplement(1..100)

The feature name "integration_excision_site" is only half shown. I would like to ask how to handle feature with long name without renaming the feature. Thanks!

biopython seqfeature genbank seqio

Been a while since I've played with this bit of biopython, but I think you will have to add the description as a different feature field such as a product. Take a look at genbank flat files and the module object model to see which ones you can choose from.

0 answers

No answers yet.

Log in to answer this question.