Questions tagged [phylogenetics]

Phylogenetics is the study of evolutionary relationships among biological entities - often species, individuals or genes.

Filter by
Sorted by
Tagged with
0
votes
1answer
51 views

How does one convert a log-likelihood substitution matrix into a probability matrix?

I am trying to calculate the substitution probabilities for amino acids (i.e. the probability of one amino acid given another). In theory, you should be able to derive them from a log-likelihood ...
0
votes
1answer
21 views

Consensus on how to select SNPs for creating a phylogeny when the full dataset is too large

I have got a quite large VCF file (about 20k samples with 160k SNPs after filtering for quality etc.) and I would like to get a phylogeny for it. However, the whole dataset's too large for my ...
4
votes
1answer
98 views

What origins for virus HIV1?

What origins for virus HIV1? This is the alignment of a small part of the POL protein of the HIV1 virus GenBank KU749412.1, with the POL protein of the Visna virus GenBank L06906.1 ...
2
votes
0answers
22 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
1answer
37 views

Plot information stored in dataframe within a tree (ggtree)

Hello to the entire Stackoverflow community! I'm writing to you because I'm currently building a phylogeny with ggtree and I have in parallel a table like this one: df : ...
1
vote
0answers
37 views

Comparison of phylogeny of SARS2 whole genome vs spike genome

I’ve built two phylogeny at minimum evolution in MEGA-X and extracted these as newick files. My next step is to compare the two phylogeny to look at if the mutations within SARS2 are more prevalent ...
1
vote
1answer
28 views

Is there a Python program that is able to show convergent and divergent evolution on a phylogenetic tree?

I was wondering if there was a Python program that detects convergent and divergent evolution on a tree. I am also curious if it is just better to look at a phylogenetic tree and conclude what is ...
1
vote
2answers
118 views

Download complete Phylogenetic Tree as JSON

I would like to integrate the phylogenetic tree into a website. Does anyone know where I can download the complete tree? Best would be a JSON, which looks something like this: ...
1
vote
1answer
63 views

How to get % similarity between strains and mutation files

I’m very new to python, and having some difficulty getting hang of some more complicated things I have multiple files which look like so: hCoV-19/Singapore/4/2020|EPI_ISL_410535|2020-02-03 hCoV-19/...
1
vote
1answer
26 views

panic: division of zero by zero or infinity by infinity error when running GODON positive selection analysis

I am using godon to find genes under positive selection with the branch-site test. However, after running for a little while I get the error "panic: division ...
1
vote
1answer
14 views

Comparing groups of 16S rRNA sequences - how to?

I have several groups of 16S rRNA sequences associated with taxonomic groups (let's call them A, B and C). The sequences within each group have a common ancestor and are on average more closely ...
1
vote
1answer
30 views

Add segregating sites to branches of a tree

I'm trying to figure out how I would plot a tree with number of segregating sites on display on the branches I'm using acctran from phangorn to plot a parsimony tree (using phydat object from fasta ...
1
vote
1answer
46 views

how to estimate phylogeny without DNA sequences?

I am new to bioinformatics, I am reading ‘Analysis of Phylogenetics Second Edition and Evolution with R’ from Emmanuel Paradis. I can create phylogeny from DNA sequences, by first calculating the ...
0
votes
1answer
27 views

Unable to use phangorn::phyDat

I am following an example from ‘Analysis of Phylogenetics Second Edition and Evolution with R’ from Emmanuel Paradis. He is doing: I am doing pretty much the same: ...
0
votes
1answer
40 views

Minimal working example for phylogenetic tree construction in R using an alignment of DNA strings and a vector of names?

I am looking at how to build a phylogenetic tree in R from aligned genetic information (dna sequences). There are multiple ...
1
vote
1answer
28 views

How to assess if this strategy is good to have some phylogenetic information?

I have several organisms and I want to calculate for each one of these their codon composition implementing a pairwise distance matrix between these codon vectors to build a phylogenetic tree. ...
2
votes
1answer
41 views

Standard for when to collapse phylogenetically uninformative nodes

I have been creating phylogenetic trees using RAxML (using Boostrap + ML) and I end up with a bestTree file, outputted in newick format. When I open up the tree in FigTree, I notice that it collapses ...
1
vote
1answer
114 views

Relaxed and sequential Phylip format conversion

I routinely use AlignIO in BioPython for manipulating alignments and in context here moving from fasta format to phylip format. Alot of phylogeny packages accept a relaxed (sequence ID) & ...
2
votes
1answer
55 views

What does instantaneous rate matrix mean?

PhyML is a tool used for maximum likelihood estimation analysis. I was running a PhyML analysis to check the rate of nucleotide change in DNA sequences and in the PhyML output, there are rate ...
0
votes
1answer
42 views

Finding out-group in unrooted tree

Im trying to find the out-group in the alignment my data files file a and file b both the case i'm getting NA when i run this ...
0
votes
1answer
60 views

Removing stop codons from alignment file for dn/ds calculation

This is one of the question which is very similar to my problem or almost same to my problem I ran with dummy data nucleotide alignment and the amino file it works perfectly fine,i can run the dn/ds ...
1
vote
1answer
27 views

Calculating Synonymous and non-synonymous

library(seqinr) library(ape) library(phangorn) alignment file data file sylvia.aln <- read.alignment("abhi_seq/all_seq_no_gal.fas", format="fasta") ...
1
vote
2answers
77 views

parsimony and maximum likelihood tree comparison in R

Maximum Likelihood analysis compares the tips or the species based on their sequence similarity while Parsimony analysis compares the characteristic features among the species. So far I did ...
4
votes
2answers
58 views

Motivation behind the neighbor-joining distance matrix recomputation

In each iteration of the neighbour-joining method of phylogenetic trees construction, after joining the nearest neighbours the (additive) distance matrix $D_{m \times m}$ is recomputed as $$Q_{ij} = (...
4
votes
4answers
257 views

Comparing phylogeny in R

So I want to compare the phylogeny created using two methods for example Maximum likelihood and maximum parsimony.Is there any way to compare the two phylogeny ? I did read about phangorn but not ...
1
vote
1answer
20 views

Phylogenetics software that can represent collapsed nodes as triangles

I am trying to represent a large (100-300+ sequences) phylogenetic tree in a way that gets my point across. I figured that reducing the number of branches (by collapsing clades which include highly ...
1
vote
1answer
57 views

phangorn fasta file read error

Im trying to read a multifasta sequence file into phangorn I get the following error ...
1
vote
1answer
36 views

After KEGG and GO analysis, how to make tables+phylogenetic trees

hope everyone is ok. I used Trinity to do a de novo transcriptome assembly, then blastp/blastx and then used Blast2GO software to do KEGG and GO analysis. So i got some txt files with header : for GO ...
0
votes
1answer
36 views

Phylogenetics analysis between diverse species

Hope everyone is well. So i was appointed with a task, and i need some guidance. Attention, the one that gave me the task may not have explained this very clearly. I will write what he exactly asked ...
1
vote
1answer
39 views

Why do we scale the rate matrix of a substitution model to make the average rate of substitution equal 1?

The matrix of transition probabilities for a substitution model Q over time t is found as follows: $e^{Qt}$. Since Q and t only show up as a product, we cannot distinguish time and rate of change. So ...
0
votes
1answer
52 views

Relationship between the instantaneous rate of change of a continuous time Markov model and the rate of a global clock model?

If continuous-time Markov models contain a parameter q that denotes the instantaneous rate of change and a global/strict clock model is the overall rate of evolution, how are these two parameters ...
3
votes
1answer
405 views

Are the conclusions in “The proximal origin of SARS-CoV-2” legit?

I don't have any background in genetics and bioinformatics, so I ask you if you think that the arguments provided in the article The proximal origin of SARS-CoV-2 by Andersen et al. are convincing. In ...
4
votes
1answer
114 views

Comparing phylogenetic models with different datasets

I'm a linguist interested in phylogenetic tree inference using language data. I'm posting here because I'm using Bayesian phylogenetic methods in my work (probably using BEAST and/or RevBayes). For ...
2
votes
1answer
40 views

How to construct the phylogenetic tree of the “High Recombination” genus streptomyces?

I need to construct the phylogenetic tree of the whole genus of streptomyces. And my goal is to find out how my target strep species distributed in the whole strep genus and obtain the distribution ...
3
votes
1answer
94 views

Coronavirus phylogeny and evolution

This new article presents the phylogeny of the coronaviruses and the placement of the new coronavirus in it. This made me wonder of what we know about the coronavirus evolution, particularly where it ...
2
votes
1answer
78 views

How many single copy gene are needed to construct the phylogenetic tree?

I ran the orthomcl program with genomes of the whole streptomyces genus(from ncbi) as input, but I only got 22 orthologous groups(streptomyces's genome size vary from 2M to 15M, which shows the great ...
1
vote
1answer
29 views

How to seperate the genus into different group based on th ANI calculation?

I've calculated the ANI of the streptomyces genus, but I'm confused about how to separate them into a different group, could anybody kindly tell me how to do that? Thanks!!!
1
vote
1answer
42 views

How to construct the phylogenetic tree of the whole genus?

I am confused about how to construct the phylogenetic tree of the streptomyces genus. The taxonomy of the streptomyces genus is mostly based on the MLST method, but now I need to use the gene data ...
-2
votes
1answer
149 views

Algorithmic investigation of viral proteins for vaccine design [closed]

I wish to investigate how a virus could be attenuated through targeted mutagenesis by a priori prediction of attenuating amino acid mutations. In particular, I am seeking bioinformatic solutions where ...
1
vote
2answers
45 views

Download sequences of isoforms

I want to collect all isoforms of all genes (Fasta/Fastq file of nucleotide and protein sequences) Arabidopsis Col-0. I am wondering if there is a straightforward way to download the file from any ...
1
vote
1answer
87 views

a package with randomization/permutation tests for genomic analyses

please would you let me know, is there a package in R/python/C that has implemented randomization/permutation tests for joint genomic analyses (an example of joint genomic analyses -- when jointly ...
0
votes
2answers
63 views

efetch nucleotide -> protein ids

I'm downloading translated viral genomes initially via Blast - which had shortcomings - and now by efetch. What I want to do is obtain a complete list of all ...
3
votes
1answer
31 views

Database with antigens and current vaccines

I am seeking a virological database comprising the following information: vaccine. antibody receptor. Preferably hosting multiple strains/samples of a given virus samples of strain. Any assistance ...
0
votes
1answer
82 views

How can i get data set of illness and not illness people?

I'm reading a virology paper where the authors made a classification for the patients as to whether they had a specific disease outcome or not against the underlying sequence ... The research focused ...
1
vote
2answers
233 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
votes
1answer
275 views

How to generate ultrametric phylogenetic trees?

I have generated rooted species tree with 1000 bootstraps by RAxML. Further, I need to generate ultrametric tree for the same. Therefore, Please suggest me how to generate ultrametric tree from the ...
2
votes
2answers
93 views

Detecting whether genome is genetically engineered?

Is there a methodology for detecting whether a genome has been tampered with by a biologist? For instance, say there has been an ebola outbreak such as in "The Andromeda Strain". Is there a way a ...
1
vote
3answers
213 views

Edit all the fasta headers using awk

kindly your comment highly appreciated have more than 4k header sequences look like: ...
1
vote
0answers
35 views

Algorithmic Information Measurement of Human Genome Databases?

What is the highest data compression ratio achieved on a population of genomes measured to include the compressing algorithm itself? The size of the compressed database approximates its algorithmic ...
28
votes
5answers
8k views

Why does the FASTA sequence for coronavirus look like DNA, not RNA?

I'm looking at a genome sequence for 2019-nCoV on NCBI. The FASTA sequence looks like this: ...