Running Mauve on circular genomes
Hi!
There are a batch of circular genomes i have to align with Mauve. The issue is that the genomes start at different positions.
I'm wondering, if it's an important thing i should manage before the running of Mauve software ? I.e. should i somehow change the start coordinates in all the genomes so that these turn to be identical before the analysis?
Thanks!
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Can the DESeq2 object for paired data be set with design= ~ batch + condition?
written by bioinfo •
-
Aligning sequences of different lengths with MEGAX CLUSTALW
written by itsmarina.n •
-
Installation of libraries as non-root user in server without internet connection
written by RiccPicc •
-
Whole-genome alignment of two or more bacterial genomes - find structural variants
written by Tim
-
STAR parameters in plants genomes
written by claudiorivero92 •
-
StringTie with Ensembl annotation: Fuzzy exon identification
written by jonasmst
-
NGS hybrid reads
written by Max Ivon
-
how to find SNP positions (for non-bioinformaticians)
written by CrazyB
-
Sequence Haplotypes From 1Kg Project Variants
written by Zev.Kronenberg
-
Aligning Mitochondrial Dna
written by Ehamberg •
Were you able to solve the issue?
Hi! To get a nicer image representation it's important to set up the same start positions for all aligned genomes. I did it in a batch mode in
Geneioussoftware. Unfortunately i don't know any free software for that.I am working on bacterial genome. I assembled the genome with canu and polished with pilon. The genomes start at different positions. Do you have any idea if the structural variants detected by mauve will be correct in this case if I call variants without changing the arrangement of genome?
I think you will be able to figure out, which discrepancies are caused by the genome rearrangements and which ones just because of the different start positions by manual alignment inspection. E.g. i'd expect to see the same LCB at the somewhere in the beginning of some genome sequences, in the middle and even in the end for other ones just because of the differences in the residues numbering.
Yes, I could figure it out after aligning with mauve (de novo genome assembly of bacterial genome) I am just concerned if calling structural variants without altering the genome would be fine.
I'm sorry, but i ran
Mauveonly for the whole genome alignments and comparisons, not for the variant calling. So, now i know nothing about how the algorithm actually works. You can open a separate post for that or try to get it from the software documentation.