Stack Exchange Network

Stack Exchange network consists of 174 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [python]

python is a programming language, widely used in bioinformatics

2
votes
1answer
53 views

Drug with as many targets as possible through a system biology graph?

I have 4 pandas data frames containing the following data: affinity_score1: affinity between each couple (Drug-Target) measured using tool1 affinity_score2: affinity between each couple (Drug-Target)...
2
votes
0answers
30 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)" ...
1
vote
1answer
85 views

python does not quit when the input file is too big

I have a script that extracts the set of unique kmer pairs from a kmer dump file (& then using them downstream to guess ploidy). The script seems to work correctly (it's reporting a meaningful ...
1
vote
0answers
47 views

How to plot a 2D RNA structure in python using dot-bracket input?

I need some help with forgi library for visualisation of RNA secondary structure. I have an RNA sequence and dot-bracket notation of the secondary structure and I would like to plot 2D graph like ...
5
votes
1answer
93 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: <...
2
votes
1answer
70 views

Snakemake: Cannot find first rule?

Back again with another snakemake query. This time I decided to port my read cleaning and alignment pipeline to snakemake. Repeating the steps from previous question and trying not to make any typo ...
2
votes
1answer
15 views

How to remove frame on RDKit figures?

RDKit can plot molecules, thought the structures are surrounded by an ugly frame with ticks. How can I turn that off and plot the molecules without that frame? Note, the ...
1
vote
1answer
11 views

Rdkit.Chem.Draw.savefig() creates empty figures

I can draw a molecule in a jupyter notebook with RDKit. Though, when I use the savefig() method only a white canvas is saved. No molecule is saved. ...
1
vote
1answer
53 views

How to separate a microbial strain of MTB resistant and susceptible drugwise?

Im trying to classify my list into three type: XDR - A strain must be resistant to: ...
0
votes
1answer
90 views

How to get PSSM of protein by protein sequence with python? [closed]

I didn't find way that convert protein sequences in FASTA to PSSM by python, and how to write program in python to find PSSM for our data set ? Maybe better to question how to find PSSM from NCBI ...
4
votes
2answers
50 views

Group genes by functional categories suming expression values

Using the count of rpkm values from genes in a metagenome sample, I want to group these genes into established categories (for example KEGG or COG). For each sample, my goal is to determine which ...
0
votes
0answers
26 views

Python module bamtags

Does anyone know this python module, bamtags? I'm using a pipeline (Dovetail Genomics, HiRise) that imports it, but it seems it doesn't exist. It's impossible to install it and there is no ...
3
votes
1answer
48 views

RSVSim insertions from chr1 to chr2

edit explaining python tag I would still rather have a solution based on the RSVSim package in R, but while waiting for someone who might have an answer I wanted to look at other solutions as well. ...
2
votes
0answers
101 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
47 views

Store contigs and genome substrings based on mismatches

I have a string, a list of substrings, and a list of positions where substrings are found in string. ...
4
votes
0answers
63 views

How to plot character state changes from a presence-absence matrix on to a phylogeny

I am trying to assign character state changes from a presence-absence matrix to a phylogeny. I have tried assigning each character to its leaf node, and then if the leaf node's sister has the same ...
1
vote
1answer
43 views

AttributeError: module 'scater' has no attribute 'normalize'

I am using rpy2 to run normalization using scran package. After calling computeSumFactors I ...
4
votes
1answer
57 views

python package for NNI neighbors

I am working on protein sequence data files for reconstructing phylogenetic tree and I need to generate all NNI-neighbours of the tree (two trees are NNI-neighbours if one can be transformed into ...
3
votes
1answer
115 views

How to iterate protein sequences using amino acids?

I was working on certain program using python and I have been using a protein aligned sequence file in two formats, phylip (.phy) and clustal (.aln). Example clustal file: ...
2
votes
1answer
68 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. ...
3
votes
1answer
136 views

How to automate NCBI genome download

I need to download all the completely assembled cyanobacterial genome's GenBank file(.gbff) from NCBI(RefSeq or INSDC ftp data). For this I think, the steps are: Need to find the completely ...
0
votes
1answer
59 views

Check if a string is ensembl gene id or not in python

I have many scRNAseq datasets saved as hdf5 files, and when I open them keys sometimes are confused, so I need a way of ...
-1
votes
1answer
42 views

Convert TPM-normalized matrices back to UMI in python

I want to process a plenty of scRNAseq datasets in python, and I want to run TMM ...
5
votes
2answers
78 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: ...
4
votes
1answer
104 views

Using pysam with cython: htslib/kstring.h not found

I'm trying to learn to use cython to speed up some code based on pysam. My issue is not strictly speaking about bioinformatics, but rather about building tools using a bioinformatics library. I hope ...
-2
votes
1answer
82 views

Python API for working with ENSEMBL genomes

Dear bioinformaticians, I could find two python APIs to mine ENSEMBL genome databases. There is a pyensembl which is more popular and then there is ...
1
vote
1answer
40 views

How to efficiently extract coordinates from a `pybedtools.Interval`?

I'm using pybedtools to process data, and at some point, I loop over a BedTool object and use unpacking to extract coordinates from the ...
0
votes
1answer
38 views

How to calculate Moreau-Broto autocorrelation, Moran autocorrelation, Geary autocorrelation from protein (amino acid) sequences?

I want to code Moreau-Broto autocorrelation, Moran autocorrelation, Geary autocorrelation in Python 3.6 from the scratch for amino acid sequences without using any package. Unfortunately, I have not ...
2
votes
1answer
75 views

Plotting Ramachandran Plot from more than one PDB file

I want to plot a Ramachandran Plot for around 5000 PDB files in python that is saved in a folder in my system. I am searching for a way to make it without accessing each file by hard-coding each file ...
2
votes
1answer
124 views

pysam pileup: what reads appear in the pileup?

Cross-posted on Biostars (with no answers currently). I hope that's OK. I thought that iff a read matches the reference at a position, it would appear in the pileup column for that position. But ...
2
votes
1answer
65 views

Download proteomes from NCBI based only on binomial names

I'm trying to download around 2000 proteomes from NCBI, but all I have is the genus and species e.g: lophodermium_seditiosum trichoderma_harzianum Does anyone know ...
0
votes
2answers
136 views

Running python script from terminal: umap issue

Some time ago I asked the question: RunUMAP in Seurat not working: module 'umap' has no attribute 'UMAP' It was solved by following the answer in it, however, it is still not ...
3
votes
1answer
60 views

Issues installing Dendropy

I want to perform phylogenetic analysis using Dendropy. I have installed dendropy using command sudo pip install -U dendropy but when I import it in python3.5.1 it ...
6
votes
1answer
69 views

Generating random BED intervals given constraints

Problem is to generate a random BED interval given the following constraints: minimum start maximum end fixed length maximum ...
3
votes
1answer
91 views

How to extract the protein fasta file from a genbank file?

Let's say I have a genbank file, e.g. the one downloaded from here, that contains entries as ...
14
votes
1answer
229 views

Biopython Phylogenetic Tree replace branch tip labels by sequence logos

Having recently constructed a lot of phylogenetic trees with the module TreeConstruction from Phylo package from Biopython, I've been asked to replace the branch tip labels by the corresponding ...
4
votes
1answer
102 views

Finding proteins in a sequence

I have a file with protein sequences, all on a single line, like this: ...
5
votes
2answers
121 views

Is there a Python/R package with the ability to convert an alignment and reference into a CIGAR?

I'm writing a python function from scratch to do this, but I feel like this must exist in some standard bioinformatics library already. In principle, this is a simply regex operation which many must ...
2
votes
2answers
57 views

counting individuals with the present gene and absent gene for each population

I have a folder containing 500 files from different individuals: ...
6
votes
1answer
92 views

How to convert the given mathematical computation (on biological problem) to mathematical fomula, equation?

I have crossposted this question in maths StackExchange. The problem is dominantly mathematical (this question) but the application of the problem is mainly biological. Hoping that people in this ...
1
vote
1answer
159 views

I need to use my code with fasta format. Can anybody help?

I need some help. I have written a program that simulates protein cleavage but I need to modify it so that it produces fasta formatted output. here is all of my sequence in .fas file ...
7
votes
1answer
244 views

Convert R RNA-seq data object to a Python object

I have done some work in R and would like to try a Python tool. What is a good way to import the data (and its annotations etc) as a Python object? I am particularly interested in converting a ...
3
votes
2answers
71 views

Pepsin digest (cleavage) does not work using RE?

Aim is to code the theoretical peptic cleavage of protein sequences in Python. The cleavage (cutting) rule for pepsin is: 1234^567. This is a position, ...
5
votes
1answer
101 views

Evaluate clusters of individuals by using their sequence data

For a dataset of several hundred individuals, I applied a hierarchical clustering to generate clusters based on a functional trait that sets them apart. My task is now to evaluate if these clusters ...
5
votes
1answer
451 views

Extracting the CIGAR string from a BAM via Python?

Is there a standard method in Python to extract a CIGAR string from the BAM? There are great libraries which parse the CIGAR, e.g. https://pypi.python.org/pypi/cigar/0.1 ...
7
votes
1answer
152 views

How to get the count of each kmer past 255 using khmer

I have a Fastq file and I want to get the exact count of each possible kmer from this file. On a previous post called How to use khmer to count k-mers? Daniel Standage proposed a custom script based ...
5
votes
3answers
187 views

Generating the reconstructed alignment from BAM

I have a (small) BAM file with CIGAR and MD fields. Question 1: What tools exists in Python and/or R to reconstruct the alignment between the reference and the read in a BAM? Given that this is a ...
2
votes
0answers
49 views

Post ipyrad filtering of SNP loci from GBS, based on SNP quality

I want to use SNPs produced by ipyrad analysis of GBS data in RAxML to examine the monophyly of a highly-variable focal species and its phylogenetic relationships with about 40 congeners. An initial ...
12
votes
7answers
3k views

What is the fastest way to get the reverse complement of a DNA sequence in python?

I am writing a python script that requires a reverse complement function to be called on DNA strings of length 1 through around length 30. Line profiling programs indicate that my functions spend a ...