Comments (6)
Hi Maria,
Yes the -xmfa_file true
allows you to use a xmfa file as input which can contain multiple (unconcatenated) genes. R/theta is estimated for the whole set, not for each gene separately, but you could compute the number of recombination events that hit each gene from the list of all events contained in the output file with suffix importation_status.txt
Best wishes,
Xavier
from clonalframeml.
Hello,
Yes you can use ClonalFrameML to analyse alignments of core genes. You need to prepare a xmfa file where each core gene sequence is a multi-fasta block separated by lines with just the '=' sign. So if you have each of your core genes as separate fasta files, you just need to add a '=' at the end of each file and concatenate these files to produce your xmfa input.
Best wishes,
Xavier
from clonalframeml.
Thank you so much for your reply. Just to clarify, I need to create separate multi-fasta files for each of the core genes (so 4059 multi-fasta blocks) and then concatenate them together with an '=' separating them.
As it stands my MSA is in the standard format (
Sample1
......
Sample2
.....)
and does not specify where genes start and end.
from clonalframeml.
Yes, that's right. You will need to find out where genes start and end, because otherwise the software will assume that they occur next to each other which affects linkage and recombination detection.
from clonalframeml.
One can only use the concatenation of genes under the xmfa format with the option "-xmfa_file true" included in the command. Right?
I am running it now, so am I going to get an estimate of R/theta for each gene? Is it unreliable to make an estimate per gene? Would you recommend LD based approaches instead?
Thank you!!
from clonalframeml.
from clonalframeml.
Related Issues (20)
- cfml_results.R error HOT 8
- xmfa format question; error: Segmentation fault (core dumped) HOT 1
- Question: xmfa format? HOT 7
- Recombination events total HOT 1
- convert_NewickTree_to_marginal_tree() HOT 2
- ERROR: incompatible with a bifurcating unrooted tree HOT 1
- node labels HOT 2
- ClonalFrameML for analyses of gene families HOT 2
- How to solve problems about "makefile", Thanks a lot! HOT 2
- input alignment for ClonalFrameML HOT 1
- segmentation fault (core dumped) HOT 1
- how to calculate the r/m value for subclade of a treeοΌ HOT 4
- R values for calculation of R*delta HOT 2
- Recombination filtered output vs importation_status.txt HOT 1
- why the r/m I calculated is so small using a subclone of subclone of klebsiella pneumoniae HOT 1
- Is it better to order genes by the reference genome for a xmfa file ??? HOT 2
- std::bad_alloc problem HOT 1
- Are there any detailed description of the output pdf from cfml_results.R? HOT 1
- ERROR: Token: ((( Found left bracket but no right bracket HOT 1
- High recombination sequence detection HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
π Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. πππ
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google β€οΈ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from clonalframeml.