Questions tagged [biopython]

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

25 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
5
votes
0answers
105 views

Using dssp after chain extraction

I have a list of PDB IDs with realtive chains that are to be extracted, and the run on dssp. For the single chain extraction I tried several methods, such as: ...
5
votes
0answers
726 views

changing blast parameters using NCBIWWW module

I have found a blog post with a script that I would like to use for my current research project: link The script is incredibly fast and produces a smooth conservation plot. In the blog post, the ...
3
votes
1answer
40 views

Counting hexamers in fasta sequence and identify its structure (and interruptions)

I have a lot of fasta files, each one with thousand of reads containing the hexameric motif "CCCTCT". The hexameric motif is highly continuous in most cases but interruptions may occur. I ...
2
votes
0answers
52 views

How to extract clusters from a phylogenetic tree?

I am using Biopython Phylo and RDkit Cluster to obtain a UPGMA tree from a distance matrix of 5k x 5k entries. How can I extract clustered entry names given a number of clusters? For example, Bio....
2
votes
0answers
155 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)" ...
2
votes
1answer
406 views

Biopython - extracting query coverage from XML Blast output

I need to extract (or to calculate precisely) the query coverage from Blast output. I am using NCBIWWW.qblast to query Blast in XML format. Is it possible to know the query coverage parsing the XML ...
1
vote
0answers
15 views

Download COX1 (COI) gene via biopython using accessions for entire mitochondrial genomes

I have a list of accessions for the the entire mitochondrial genomes for big cats. I need to download the COX1 genes for each of these accessions. Here is one accession and here is a link to its COX1 ...
1
vote
0answers
27 views

AttributeError: 'FeatureDB' object has no attribute 'strand'

It appers when a gene has a negative strand then I need to reverse the sequence. Unfortunately, I got the following error: ...
1
vote
0answers
45 views

3d model of a protein via nucleotide sequence

I really appreciate @Matteo Ferla's detailed answer to my previous question. I've already tested several tools to know if a mutation can cause a modification in the function of protein or not. There ...
1
vote
0answers
20 views

Finding annotated counterpart after BLASTn with efetch (Biopython)

I am creating a pipeline for the identification of unknown transcripts. After a local BLASTn search of the transcripts, I have a large list of the respective hits with different genomes. I have the ...
1
vote
0answers
19 views

Detecting Allelic Imbalance

I am interested in detecting allelic imbalance for a credible set of SNPs. Currently, I am looking for some packages, hopefully, for Python, that can assist with this. I am aware of this R package, ...
1
vote
5answers
170 views

How to remove duplicates from a fasta file using python

I am using the following command: ...
1
vote
0answers
59 views

UTR features in gbk files

I have a gbk files for a viral genome that I have used biopython to add features of the 3' and 5' ...
1
vote
0answers
28 views

How to (efficiently) search nuccore to retrieve only annotated sequences?

I wish to efficiently search nuccore to retrieve only annotated sequences. By "annotated" I mean that the nuccore entry contains annotation information. An equivalent definition (if I understand ...
1
vote
0answers
22 views

Is there a built-in parser for the textual feature table that Entrez.efetch returns?

Entrez.efetch helpfully gives me the feature table of a nuccore entry. If I understand correctly, efetch can only return the ...
1
vote
0answers
22 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 (...
0
votes
0answers
13 views

Converting miRNA names to miRBase version IDs

I have a list of miRNAs IDs (2000-2500) that I want to find miRBase IDs for them. For example: hsa-miR-106a ---> hsa-miR-106a-5p hsa-miR-373* --> hsa-miR-373-5p hsa-miR-33 --> hsa-miR-33a-5p ...
0
votes
1answer
34 views

How can I get a list of neighboring Hydrogen atoms of an alpha-carbon?

Suppose, I want to find the neighboring Hydrogen atoms of a specific C-alpha atom. The following was my effort. However, it produces a bunch of empty lists. ...
0
votes
0answers
39 views

BioPython - Retrieve sequence records from pubmed database

I'm assigned a task to fetch sequence records from "pubmed" database. I know this database is for publications and citations, but how can I retrieve those records from it? is that even ...
0
votes
0answers
51 views

Python/Biopython - Replace amino acid residue on MSA with "z" from a list of unaligned positions

I'm trying to programmatically replace a set of amino acid residues on an MSA with a "Z" from a list of unaligned positions. Any ideas on how I could do this? Input: a list of unaligned ...
0
votes
0answers
211 views

GTF parser for Python

I found this GTF parser for Python. However, it does not work with Scallop's GTF file: ...
0
votes
0answers
26 views

Blastp MSA to the same length

To generate aligned protein sequence truncate to the same length, I use blastp in biopython for MSA. I have multiple queries and one subject sequence. However, the alignment never returns the same ...
0
votes
2answers
107 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 ...
0
votes
0answers
53 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 ...
-1
votes
1answer
90 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 ...