Labeling features in a sequence from a feature list?
Hello,
I'm just starting out with Biopython. I have a table of DNA sequence features which I want to use to annotate unlabelled genbank sequence files. Is there an easy way to do this, or should I just code up something custom for myself? I feel like someone must have tried this.
Thanks!
• 1,452 views
•
link
written
by
npatel89 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Seqrecord object to only a single line sequence (without ID)
written by voltosj •
-
ValueError: Need a Nucleotide or Protein alphabet. Trying to write a Genbank file and getting error
written by Fishy_DNA •
-
Genbank to fasta java
written by orelsmail •
-
Create a WebLogo with gapped sequences using Biopython
written by makuba •
-
Download fungal genomes
written by d.pinili •
-
Biopython drops GenBank LOCUS data?
written by jamesfowkes •
-
Querying the UCSC Genome Browser for DNA sequences (large collection of data)
written by mshumph2 •
-
How To Differentiate Files With One Record From Files With Multiple Records?
written by Cako Fuentes •
-
Concatenate Sequences With Bioperl
written by maczerek •
-
Error Translating Genbank Cds Using Biopython
written by Tim
What does your table of features look like? Sequences or co-ordinates etc? You need to show us a sample of what what your input data looks like.
Example feature list:
Feature1 ATTAACCCTCACTAAAGGGA
Feature2 TGTAAAACGACGGCCAGT
Feature3 CAGGAAACAGCTATGACCATG
Feature4 TAATACGACTCACTATAGGG
Feature5 ATTTAGGTGACACTATAG
They are just sequences and I want to find all feature coordinates that may or may not exist in an unlabelled plasmid sequence file.
The cloning software SnapGene has this feature. I can't recall if it's available in the free version, but give it a try.