Questions tagged [biopython]

Biopython is a set of freely available tools for biological computation written in Python.

Filter by
Sorted by
Tagged with
0
votes
0answers
16 views

Can I pull a jackhmmer output file into ClustalO for pairwise alignments?

I'm still working on figuring out how to send ClustalO multiple jobs. I have a hmmer output and was looking to do a mixture of pairwise and local alignments on a subset of these. Is it possible to ...
1
vote
1answer
65 views

Can I automate CLUSTALO and output alignment sequence identity?

I've detected homology between targets of ligands in drugbank and proteins in the proteome of a pathogen. I've parsed the output very rudimentary and calculated my query coverage. This exists in an ...
2
votes
1answer
52 views

Possibility to save output blastn table in memory using biopython

Is there a possibility, using biopython to save output table of blastn in memory, not in file on hard drive, for process it using pandas for example and then delete it from memory? I mean, is there a ...
1
vote
0answers
14 views

Indel quality from mpileup

Could someone please confirm that I understand this correctly? I have the res and qual columns from mpileup and I would like to match them to get the qual per base. It seems that the indel initiation (...
1
vote
0answers
15 views

Visualizing cell growth

I model the following events: a rod-shaped cell with center $x$ grows symmetrically until it reaches a maximum length/age, it divides into two identical cells, and the process continues. At division, ...
2
votes
1answer
21 views

Tools for modelling and visualizing growth of cells

I am looking for a program/GUI app/package that would help me do simulations for cell growth. I have a microscopic mathematical model, and the scenario is basically the following: I start with a rod-...
1
vote
1answer
21 views

Cutting the sequence into several sequences with the information of a dataframe

I have a fasta file with several 120-concatenation protein sequences. I also have a data frame with 120 names of proteins in column 1 and their length in column 2. Using this dataframe I want to ...
0
votes
1answer
24 views

how to generate structure file from blast results of an input sequence

I have tried to incorporate blast in my application using biopython.i could extract the blast results but now i want to extract the structure file corresponding matches. My code is : ...
0
votes
0answers
46 views

Parsing a GenBank file with multiple gene entries

Basically a GenBank file consists of gene entries (announced by 'gene') followed by its corresponding 'CDS' entry (only one per gene) like the two shown here below. I would like to extract part of the ...
0
votes
0answers
18 views

What values for database keyword can be used with Bio.Blast.NCBIWWW.qblast?

I want to make a Blast search with Bio.Blast.NCBIWWW.qblast for a protein sequence. from Bio.Blast.NCBIWWW import qblast qblast("blastp", "nt", sequence=seq) The ...
1
vote
2answers
53 views

How to create a .phy file for constructing a phylogenetic tree?

I'm new here. I need to construct a phylogenetic tree from 16s Sequence data (merged 151 bp) and color the branches/nodes based on my experimental conditions. From this blog, it looks like I need a ....
1
vote
2answers
40 views

No output when trying to obtain protein sequence from PDB file

I am trying to obtain the sequence from PDB file, but no output is given when I run the following code : ...
2
votes
2answers
87 views

How to translate amino acid sequences to Nucleotide sequences

I want to convert a list of fasta ( protein sequences) in a .text file into corresponding nucleotide sequences. A Google search gives me result of DNA to protein conversion but not vice versa. Also, I ...
1
vote
1answer
38 views

renaming IDs in gff3 with BCBio.GFF

I wrote a script which should changes IDs in a GFF3 file. Unfortunately, the below script has two problems. It attaches the new ID to Parent which leads that ...
2
votes
1answer
40 views

Infer the new (Ban et. al) ribosomal nomenlature (ex. uL53 ) from the ribosomal protein's sequence

I'm trying to programmatically construct a name for each protein subchain in any ribosome from Uniprot in accordance with Ban et. al's 2014 proposal (excerpt given below) using PDB's and Uniprot's ...
1
vote
1answer
33 views

Unable to read all alignments in stockholm format : UnicodeDecodeError: 'utf-8'

I have a file in Stockholm format and I am trying to read it using the Biopython package. The code is as follows : ...
0
votes
1answer
46 views

Given a .gb file and a locus - how to get relevant annotations in Python?

Given a .gb file and a specific locus in the genome - how can I retrieve the relevant annotations in Python (i.e., annotations that include that locus)? I could ...
3
votes
1answer
68 views

protein secondary structure matching/comparison scoring

I want to compare two secondary structures of aligned proteins. I do not know exactly how to do it well. Example: ...
0
votes
2answers
49 views

What kind of analysis is practically done on GSE data files?

I have a GSE data file in csv file format containing fields such as: ID, adj.P.Val, P.Value, t, B, logFC, Gene.symbol, Gene.title. In which adj.P.Val, P.Value, t, B, logFC fields being numeric. What ...
1
vote
1answer
111 views

Is there any tool for fuzzy sequence matching?

Note: this question was also asked on biostars are you aware of any tool that is able to perform error-tolerant pattern-matching search on protein FASTA files? For example, I want to know, which ...
0
votes
2answers
70 views

read and matching pattern with python

I have read the content of a text file into pandas and needed some help matching the pattern. Here is the pattern, where a can any number greater than zero and <...
2
votes
0answers
23 views

How to load output from the kegg API in biopython into a pandas dataframe?

Biopython provides a (bit unintuitive) API to access KEGG. I am trying to make use of it, but the output is quite unhandy as a string. What is the best way to parse the data into a ...
1
vote
1answer
37 views

Biopython Genbank writer not splitting long lines

I am parsing a csv file of annotated sequences and using Biopython to generate Genbank files for each. I want to add annotations of the sequence features. My output file shows features listed without ...
2
votes
1answer
42 views

Parsing specific features from Genbank by label?

I'm trying to parse a Genbank file to find a specific feature. I can pull it out if I know the feature type (e.g. repeat_region). E.g. if I'm looking for this feature: ...
1
vote
1answer
65 views

From ensembl id to protein structure

Forgive me if this question is too trivial. I have a list of ~700 Ensembl ID. I need to extract the protein structures of each Ensembl id. How can I do this? Is there any script(python, R) so that ...
1
vote
1answer
59 views

I have a list of protein targets. How to get structural data (atomic coordinate) from PDB?

I have a list of around 500 protein targets (human). I only have their protein symbols. As an example: EGFR AR FOXA1 CXCR3 ... I want to get structural data (atomic coordinate) from PDB (Protein ...
1
vote
2answers
79 views

Biopython Failure to break out of loop

After qblast, I wanted to print only details from the first 3 BLAST output. Yet, the code below printed the details of every BLAST output. Why can't I break out of the loop to achieve my desired ...
1
vote
2answers
87 views

How to efficiently get human gene names from NCBI based on a large list of SNPs

I found a good answer related to my question here: How to get a list of genes corresponding to the list of SNPs (rs ids)? But it says about small number of SNPs. I want to get gene names based on ...
3
votes
1answer
70 views

How to select only RNA with Hetero atoms from pdb file with python?

I'm trying to separate RNA from protein in a complex protein/RNA PDB file and I want all RNA info with the hetero atoms in between the bases BUT without H20 etc. In short I want RNA part of pdb file ...
0
votes
0answers
42 views

How to compare multiple genetic profiles in Python?

I have hundreds of vcf file where each vcf file contains genome profile for a tissue. A portion of the vcf file is as follows: [ I can read each vcf file into a dataframe. So it would be hundreds of ...
4
votes
1answer
159 views

HDF5 and BioSQL solutions

I'm looking at better database/storage solutions for NCBI virus data, with all attributes particularly year and country of isolation, together with structural data, possible antibody data, T-cell data ...
1
vote
1answer
191 views

What do negative dN or dS mean in cal_dn_ds in Biopython?

I'm noobie in bioinformatics. I use cal_dn_ds from Bio.codonalign.codonseq with method 'NG86'. Sometimes I get dS=-1. What does ...
3
votes
1answer
100 views

How to learn/get started bioinformatics? [closed]

I know this is a broad question and I'm sorry if this is not the right place to post (this is my first post ever), but is there some sort of roadmap for bioinformatics? I'm an undergraduate student ...
1
vote
1answer
46 views

Retrieve ID ligand from PDB file

I have thousands of PDB files and I need to extract the following information from each of them: an ID or the name of the small molecule (ligand) present in the complex. Is there a way to do that ...
1
vote
1answer
87 views

How to get Uniprot and swissprot cross references using Biopython

I am struggling to use biopython to gather cross references such as GO annotation from Uniprot/swissprot text files ("DR" in the .txt files) using biopython. From ...
1
vote
1answer
46 views

How to allow “half windows” in biopython?

I am using Biopython to examine protein sequences including the Kyte and Doolittle hydrophobicity. However, for short sequences, this returns an empty value; the value is not allowed if the window ...
2
votes
0answers
73 views

How to get bootstrap support of phylogenetic tree?

I have used multiple sequences aligned file of protein data to generate a maximum parsimony tree, then I used the "bootstrap_trees(msa, times, tree_constructor)" ...
7
votes
2answers
156 views

Truncating branch length values of Phylogenetic tree with biopython

I have been using biopython 1.72 to display my phylogenetic tree files. 1)Using the function 'Phylo.draw(pars_tree, branch_labels=lambda c: c.branch_length)' to ...
5
votes
1answer
132 views

How to find all variable-length seqs with an exact 5' and 3' match in a FASTA file

Context I am interested in finding all of the promotors specific to a particular sigma factor. I have identified the -35 and -10 sites from the literature, bold denotes -10, -35, binding sites: <...
7
votes
3answers
308 views

Importing GFF file with Biopython

Is there a way to import a GFF file for an organism with Biopython in the same way you can for a Genbank file? For example, ...
2
votes
0answers
680 views

Searching for start and stop codons for protein sequencing of contigs

I need to convert contigs into their respective protein sequences given a reference genome (i.e. I need to take a substring, whose position is already known on the string, and I need to locate the ...
2
votes
1answer
127 views

Counting the characters in a row using python

I am using python 3 and biopython 1.72. I have been using a protein sequence file in clustal format. ...
6
votes
1answer
92 views

Getting all NNI trees of a parsimony tree

I have an aligned protein sequence file which I have been using for reconstructing a parsimonious tree. I am currently using NNITreeSearcher._get_neighbors method from Biopython 1.72 but it's way to ...
5
votes
2answers
113 views

Counting letters in phylip alignment columns with Biopython

I have been using python 3.6 and biopython 1.72 to work with protein data files. I am using a protein sequence file (phylip format), for example: ...
2
votes
2answers
194 views

Get a list of PDB ids from a list of sequence protein number?

I am working on a project about phosphorylation, and I am creating a database of some proteins. I have a list of protein accession ids for which I need to get the corresponding PDB file for each one ...
4
votes
1answer
208 views

Get results of keyword search on Pfam via python script

I'm interested in all proteins that are in any way associated with Danio rerio. I decided to look them up at Pfam data base and when I just make a keyword search, I get a a nice list which looks like ...
3
votes
0answers
98 views

Glocal\semi-local\Hybride Globale-Local alignment with Python

I was looking for a simple way to do a glocal alignment. The case I have is I have a small sequence which should be find in a bigger one, thus typically a glocal alignment. Also I can not Install ...
0
votes
1answer
100 views

Find a map/correspondence between two versions of a genome

I am working with two versions of the C. elegans genome. I am finding interesting regions (specifically, tRNA genes) in version 1 and then I would like to know if version 2 also has a tRNA gene in ...
3
votes
1answer
126 views

How to assess the quality of assembled .fasta genome files?

I have assembled 3 .fasta files from contigs infastq format of 3 different Homo sapiens. I would like to see if the assembled ...
2
votes
2answers
330 views

Retrieving NCBI Taxa IDs from refseq or GenBank assembly accession

I have about 10,000 genome files all named by either refseq or genbank accession number, do you know if it's possible to convert these numbers to the corresponding NCBI taxon ID or species? for ...