Questions tagged [phylogeny]
Phylogeny is a statistical approach for reconstructing evolutionary relationships using trees
244
questions
0
votes
1
answer
50
views
Genetic relationships between Cyanobacteria: terrestrial vs aquatic
I believe one of the site members is an expert in Cyanobacteria hence the question.
I'm looking for a comprehensive 'phylogenetic tree' of the phylum Cyanobacteria and want to identify all the ...
M__♦
- 11.2k
1
vote
1
answer
59
views
Variation in 16S rRNA between assemblers - how do I know which is real?
I have a low-diversity metagenome (~11 bins > 80% completion). Out of the bins, 3 are of interest to me. None of the lower-completion bins that can be identified are from the group of interest. So ...
- 869
1
vote
1
answer
23
views
phytools add.species.to.genus could not match your species to a genus
I have a species tree to which I am trying to add species, which includes a tip "Geissospermum laeve". When I run:
...
- 115
1
vote
1
answer
30
views
Finding linkage disequilibrium when I have a kinship matrix but not pedigree
I have a bunch of samples of a non-model organism divided in two groups: parental generation and offspring. I do not have a pedigree so I do not actually know who is a parent of who, and which are the ...
- 111
2
votes
1
answer
147
views
Correlation heatmap of RNA-seq clusters all samples together leading to very low no. of DEGs
I am writing to you to take an input or may be you can provide a different perspective. I am at wits end :(
So I have nearly 200 samples. These are separated into two groups (group I treated with ...
- 227
1
vote
1
answer
50
views
Pal2nal translation of large multi-fasta files produces a codon translated file where some sequences half length of the average
I did sequence alignment of a large peptide multi-fasta (n= 4991 sequences). The peptide alignment has sequences with the same length and pal2nal went through just fine... except some of the codon ...
- 217
2
votes
1
answer
60
views
How to determine the phylogeny of a specific set of genes
My goal is to understand the phylogeny of genes involved in bacterial restriction-modification systems, in psychrophilic bacteria.
To do so here is what I have done.
Downloaded the protein sequences ...
- 163
0
votes
1
answer
82
views
How to parallelize BEAST
I'm running BEAST v1.10.4 with BEAGLE v4.0.0. My dataset has one large partition and I have access to a machine with many CPUs so I would like to optimize performance by parallelizing the process. ...
- 111
1
vote
1
answer
98
views
How to add bootstrap values to the phylogenetic tree generated by OrthoFinder?
When we run the OrthoFinder analysis tool on a group of genomes to get the orthologues shared by them one of the output files include a folder named 'Species_Tree' that contains a text file named '...
- 53
1
vote
1
answer
47
views
Evolutionary tree from SNP data
I would like to make an evolutionary tree using the SNP data from DArTseq. Are there any software or packages available for this purpose? I came across the software SNphylo, but its home page is not ...
- 11
1
vote
1
answer
61
views
How to use Gblocks for trimming single-copy gene sequences?
I have run OrthoFinder on a set of 11 genomes and got the results in a folder. From the output folder, I saved all the single-copy gene sequences in a single FASTA file and now wish to remove all gaps ...
- 53
4
votes
1
answer
173
views
How to get a GISAID account? I registered months ago, still no reply!
Inspired by amateur variant hunters, I would like to join the Pango lineage proposal community and help contribute to variant surveillance.
However, I cannot seem to get access to GISAID, the platform ...
- 43
0
votes
1
answer
118
views
Refactoring pandas using an iterator via chunksize
This question was also asked on Stack Overflow
Bioinformatics rationale eggNOG files can be very big and sump all available RAM for regular to medium sized desktops.
I am looking for advice on using ...
M__♦
- 11.2k
2
votes
1
answer
33
views
NCBI Blast databases - cannot find some databases
I'm working with a system that does local BLAST searches using some preconfigured BLAST databases. From what I see, these were updated within the last year, but I cannot find any current databases ...
- 175
0
votes
1
answer
125
views
Blast output file only shows 500 lines -outfmt 6
I had created databases of different sets of metagenome datasets - one with 6 runs, other with 48 runs, another with 100 runs, etc using the accession list for each of these datasets and makeblastdb ...
- 53
1
vote
0
answers
39
views
How to identify all descendents of a branch on a phylogeny
I'm trying to simulate the evolution of binary traits on phylogenies in R to end up with a dataset of the traits held by each tip. Starting at the root, at every timestep, new traits are added with a ...
3
votes
1
answer
106
views
What is the limit of partition incompleteness in Bayesian MCMC?
I am seeking advice on whether or not to incorporate incomplete data into my BEAST analysis.
I have HIV data consisting of two partitions of the pol gene: one that is 1.2 kb in length and present in ...
- 435
2
votes
1
answer
41
views
Parsimony-distance neighborhoods in tree space
I have the following question about the paper A parsimony-based metric for phylogenetic trees by V. Moulton and T. Wu.
In this paper, the authors define a "parsimony-distance" between ...
- 131
3
votes
1
answer
61
views
compare fasta sequences in pairs and collect metrics
I have 96 fasta files (A1, A2, A3...) from one plasmid assembly pipeline, and I have another 96 fasta files (B1, B2, B3 ...) from another plasmid assembly pipeline.
I would like to compare pair ...
- 119
2
votes
1
answer
34
views
Using ticks on the x axis or a scale bar when visualisizing a phylogram?
When visuliazing a phylogram is it preferable to show ticks on the x axis or a scale bar?
Given an example tree with branch lengths in newick format ...
- 43
0
votes
1
answer
49
views
Calculating total number of Unique gene in every genome from ROARY pangenome output
I performed a ROARY from 100 genome. I got several file including gene_presence_absence.csv from the results. I want to know how many unique gene(unshared gene, based on blast homology, same gene (...
- 15
3
votes
1
answer
157
views
Plot phylogenetic tree from list of edges
I have a dataset that I wish to convert to tree or phylo format like in the ape package, in order to plot the phylogenetic tree. It is formatted like a list of ...
- 33
2
votes
1
answer
100
views
Math on Pandas Columns
I have a pandas dataframe that reads in a PAF file from minimap2. What I would like to do is take the first 5 columns of the data from to create a BED file.
I used this to extract the first 5 columns:
...
- 804
2
votes
1
answer
65
views
Is there a graphical/interactive 16S rRNA clustering method?
I've been doing phylogenetics with large (hundreds) 16S rRNA sequences lately.
Usually I'm focusing on one order, and using a combination of trees and sequence similarity to assess stuff like 'is this ...
- 869
3
votes
1
answer
59
views
Why is COI metabarcoding not used for prokaryotes?
I'm a beginner in bioinformatics working with environmental metabarcoding, and while reading up on different marker genes, I understand that COI is preferably used for animals thanks to its higher ...
- 33
1
vote
1
answer
74
views
Phylogeny building in R from FASTA files:
Im working on building a phylogeny from scratch with downloaded FASTA sequences from GeneBank. I think Im doing alright up until the multi sequence alignment in the ...
- 111
2
votes
1
answer
43
views
Should I compensate for two clusters of samples in microarray by hierarchical clustering and MDS?
I seem to have two groups of samples based on hierarchical clustering as well as MDS plot analysis using the sample label from the tree cut of the hierarchial clustering. I am analyzing microarray ...
2
votes
1
answer
63
views
Is it possible to do homology inference across species using different kinds of NGS data?
Background: I have a list of species that I want to put through homology inference. The goal of homology inference is to investigate the evolution of a trait on a species tree. I want to use the ...
- 217
2
votes
1
answer
47
views
What do you use the CLANS program for?
I am reading this paper called Birth and death in the Terminal Complement pathway, doi : https://doi.org/10.1101/2022.04.21.489004 where they have performed CLANS analysis.
I understand that CLANS ...
- 105
2
votes
1
answer
60
views
How can I fix "Ill-defined genetic structure" Error message in Arlequin?
I am performing an AMOVA analysis via Arlequin.
At this point, I have:
collected sequences of mtDNA from the cox1 gene,
aligned the different genes using clustalX
converted the files into .arp ...
- 127
3
votes
2
answers
100
views
How can I resolve "missing name of population" message in PGDSpider?
I want to use Arlequin for AMOVA. I have my FASTA files which I want to convert into .arp files which are recognizable by arlequin. When I try to convert it using Population Genetics Data (format) ...
- 127
1
vote
1
answer
48
views
DNA, molecular descriptors
Are there any programs (preferably command line tools) for calculating molecular descriptors of DNA? I am looking for something like Chemopy or RDKit but specifically for DNA.
Thanks in advance!
- 95
2
votes
1
answer
60
views
How to make unrooted tree for Likelihood mapping result by using IQ2-tree?
I am a biologist, and I do not fully understand the tree topology of the experimental species.
I used four-taxon set (4 sequences) to identify the Four-cluster Likelihood-Mapping by using
...
- 21
1
vote
1
answer
64
views
Maximum likelihood estimation for tree construction
Could someone direct me to a resource (textbook chapter, lecture notes, online video, etc.) that demonstrates, in a mathematically or computationally rigorous way, how to use maximum likelihood ...
1
vote
1
answer
52
views
How to identify mutations in a viral genome
I have a fasta file with multiple sequences comprising 38 sequences. The length of the sequences are around 11000 bp.
How can i get changes in the genomes based in a reference genome? (aa subtitutions ...
2
votes
1
answer
182
views
ETE4 manual or improved code description
ETE3 * (Environment for Tree Exploration) is a major package and in particular Python package to navigate phylogenetic trees. It is the premier coding library/package for tree manipulation.
ETE4.0.0 ...
M__♦
- 11.2k
2
votes
1
answer
53
views
Converting Pairwise single- linkage clustering distance data to "newick" format
I have pairwise distance data for single linkage cluster and I would like to convert it to newick format.
I was unable to find a conversion algorithm.
I require newick format for visualization and ...
- 118
2
votes
2
answers
41
views
Is there any book or reference that explains evolutionary dynamics and fixation probability?
I am a master's student in computational and mathematical biology, with a bachelor's degree in pure mathematics.
I am really struggling with my computational biology course, we are supposed to study ...
3
votes
1
answer
149
views
Compare and Reorganize Fasta Headers Python
I want to compare the headers from the fasta file to file1, and if there's a match, reorganize the header and put the match first. If there's no match between fasta file and file1, look at file2 and ...
- 33
1
vote
1
answer
53
views
Using information criteria weights to create consensus tree
While I am familiar with calculating information criteria for various ml models of a multiple sequence alignment, I am not aware of how to use the derived weights to create a tree that combines ...
- 111
1
vote
4
answers
444
views
Removing particular line along with Fasta header in fasta file
I have a fasta file of 3'UTR sequences of the genes. Some of the gene do not have a 3'UTR sequence and reported with sequence unavailable but there is FASTA header.
...
- 19
1
vote
2
answers
106
views
How can I align two novel sequences to a reference genome and build a phylogenetic tree?
I have two new amino acid gene sequences that I want to align to a reference genome. Is it a good idea to have a combined phylogenetic tree? what will this type phylogenetic analysis called?
- 11
0
votes
1
answer
77
views
Single or partitioned model for supermatrix tree inference?
I am dealing with several genes in my dataset. For each of the genes, I have built gene tree, with their best estimated model. I am intended to look into the effect of the tree from concatenated ...
- 11
1
vote
1
answer
94
views
Fixing FASTA file for Local BLAST Database
I recently prefetched 157 SRA files from an NCBI BioProject using the SRAtoolkit. I then used the toolkit to download those files in FASTA format.
Each individual FASTA file looks something like this:
...
- 804
1
vote
0
answers
65
views
difficult in runing the R code for ggtree
I tried to run the below code but ended having issue. I have checked thoroughly and I didn't find any mismatch in fasta sequence
...
- 31
2
votes
1
answer
134
views
drawing multiple sequence alignment as tree in R
I am very new here. Please help me with my codes as am trying to draw multiple sequence alignment as tree suing ggtree package in R and i have following error.
...
- 31
3
votes
2
answers
131
views
How can I find the maximum percent identity between two sets of (unaligned) sequences?
As the title states, I've got two sets of unaligned amino acid sequences (~25k sequences in one, ~3k in the other). I want to find the minimum distance between a sequence in the first set and a ...
- 163
1
vote
2
answers
55
views
Improving list speed: glycosylation example
I discussed a question with @gaspanic Python/Biopython - Replace amino acid residue on MSA with "z" from a list of unaligned positions
. The issue emerged was speeding up lists in Python. ...
M__♦
- 11.2k
1
vote
1
answer
232
views
reverse translation from amino acid string to DNA strings
what is the opposite of .translate() function calls ?
I mean let's say I am given an amino acid string CYCLIC, how do I obtain all the possible combinations of DNA ...
- 141
0
votes
2
answers
781
views
RNA strings and amino acid string
Which of the following RNA strings could translate into the amino acid string PRTEIN? (Select all that apply.)
I am not sure if the following mass of amino acid is relevant to the question.
...
- 141