Question: Converting Mugsy Output (Maf) To Other Formats, Problems With +/- Strand And Overlapping Blocks
0
gravatar for tmortimer
6.2 years ago by
tmortimer0
Madison, WI
tmortimer0 wrote:

I am aligning bacterial draft genomes to a reference genomes using mugsy v1.2.3, and I'm using the branch of Biopython that supports MAF files to parse the output and convert to other alignment formats (fasta, nexus).

I've run into a couple problems: 1) In some of the blocks the alignment is to the negative strand of the reference sequence. I think I've solved this problem by identifying those blocks and using Biopython's reverse complement feature to change all the sequence records in that block. 2) It appears that some of the blocks overlap each other in reference to the positions in the reference sequence. This becomes a problem when I try to create a MAF index to help convert to other formats.

Is there any way to specify a particular sequence as reference with mugsy? Anyone have any tools they would recommend to deal with the overlapping blocks?

Thanks!

maf alignment • 3.1k views
ADD COMMENTlink written 6.2 years ago by tmortimer0
Please log in to add an answer.

Help
Access

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