Question: How to prep an input (nexus) file for BEAST/BEAUTi?
gravatar for hallishtw
23 months ago by
hallishtw10 wrote:

I am working with bacterial whole genome sequences. I understand I have to do a multiple alignment (for which I used parsnp on my de novo assembled genomes) and remove recombination first. I detected the recombination sites with Gubbins but did not get an output of multiple alignment where the recombination sites are removed. So now I am unsure as to how to proceed. I assume I have to somehow remove this recombination sites and convert the alignment to a nexus file for BEAUTi. Alternative ways to proceed are also appreciated. Thanks in advance!

alignment genome • 1.5k views
ADD COMMENTlink modified 22 months ago by Biostar ♦♦ 20 • written 23 months ago by hallishtw10

I am not sure if you have to remove the recombination. Isnt the recombination information useful to help you identify evolution and understand the phylodynamics?

I have created a Multiple Sequence Alignment file and then done the analysis in the past.

ADD REPLYlink written 23 months ago by sridhar56100

Hi thanks for your comment.That was also my first assumption as well. But after looking into literature, it seems like recombination is "noise" in the process. Even if you skipped the recombination filtering process, how would you convert the .xmfa (parsnp output) or converted .maf file to nexus for BEAUTi? Multiple sequence alignment of sequenecs of different length

ADD REPLYlink modified 23 months ago • written 23 months ago by hallishtw10
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 757 users visited in the last hour