@Istvan Albert I'm new to Biopython and trying to figure it out, I was able to print out what you suggested above. I have three gb records, genome.gb, bacteria1.gb, bacteria2.gb that have several entries inside them and I want to sort them just the records themselves by type?
Split multple gbk records by molecule type
I have multiple genbank gz files that I'm trying to sort by the molecule_type in biopython. For it to grab the protein gz files and the nucleotides in a separate entry.
• 1,479 views
•
link
updated
by
Istvan Albert
•
written
by
AveryB •
•
written
by
AveryB •
1 answer
Here is a recent post on how to extract the qualifiers of the source feature.
As you can see there the mol_type is a key in the same ordered dictionary
How can I print and write the strain /isolate/voucher number of a SeqRecord objec in biopython?
to open a gz file, use gzopen.open
• 0 views
•
link
written
by
Istvan Albert
Make a list that contains all records, then use the sort function with a key as explained here:
https://docs.python.org/3/howto/sorting.html
then you'll have a sorted list
• 0 views
•
link
written
by
Istvan Albert
Log in to answer this question.
More posts like this
-
BioPython to check multiple docking results
written by Jafs •
-
Looking for matching names in fasta file and genbank records
written by AveryB •
-
Different samples order in .ped and phynotype data files
written by Denis
-
Pulling out gene names from a FASTA file with a list of IDs
written by jvire1 •
-
Extracting a Genomic Region From Nucleotide Accession
written by 13dsc •
-
How to write FULL genbank files with Biopython SeqIO module?
written by russianconcussion •
-
Nucleotide CDS from RefSeq?
written by Prohan
-
How To Differentiate Files With One Record From Files With Multiple Records?
written by Cako Fuentes •
-
Genbank Multiple Cds For One Gene
written by matija.sosic •
-
How To Separate Protein Accessions From Nucleotide Accessions In The Livelists
written by rsingh2083 •
