Entering edit mode
3.6 years ago
QL
•
0
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!
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.