Question: Aligning Draft to Reference genomes, rearranging contigs?
0
gravatar for YOUSEUFS
2.1 years ago by
YOUSEUFS10
YOUSEUFS10 wrote:

Hello Bioinformagicians!

My end goal is to map some RNA-seq data to a genome iv'e been given but I do not know if this genome is assembled in the right order. In order to analyse this, I'm attempting to align a draft genome sequence (containing contigs and scaffolds) to a closely related reference genome using nucmer. When I plot the output.delta file using mummerplot I get this graph:

ref x-axis, draft y-axis

The genomes seem to align very well except that they are shifted slightly with the initial missing piece appearing to the right of the plot. I'm assuming that the start of the sequences of each genome is different, does this matter if I'm to compare the genome? Do they need to start at the same sequence? If so how do I go about rearranging this genome?

Many thanks!

P.S I hope the image link works!

alignment mummer assembly • 965 views
ADD COMMENTlink modified 2.1 years ago • written 2.1 years ago by YOUSEUFS10

Why would the order of your contigs matter in the context of your end goal?

ADD REPLYlink modified 2.1 years ago • written 2.1 years ago by 5heikki8.6k

I don't know to be honest, I'm just checking from more experienced people than myself to see if it would =)

ADD REPLYlink modified 2.1 years ago • written 2.1 years ago by YOUSEUFS10

Maybe if you were scaffolding your contigs, and assuming that it worked, then the number of mapped reads could increase a little bit (assuming your assembly is good otherwise). IRL, scaffolding without long mate pairs or other such data is pointless though..

ADD REPLYlink modified 2.1 years ago • written 2.1 years ago by 5heikki8.6k
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: 1722 users visited in the last hour