Question: Mapping contigs to a reference genome
1
gravatar for SpamChop
12 months ago by
SpamChop10
Edinburgh, UK
SpamChop10 wrote:

Hi!

I have ~163000 contigs (a single FASTA file) and I want to map them to a refernce genome (also a FASTA file). Is there any way to do it? I have tried bowtie2, but could not work out its nuances.

Any help appreciated!

Thanks

alignment genome • 602 views
ADD COMMENTlink modified 11 months ago by evelyn90 • written 12 months ago by SpamChop10
2

Is this a related reference (i.e. you expect high homology)? If so you could use blat or even blast+. If these are very large contigs a program like LASTZ would also be valuable.

ADD REPLYlink written 12 months ago by genomax83k
1

Try minimap2 with preset -x:

asm5/asm10/asm20: asm-to-ref mapping, for ~0.1/1/5% sequence divergence

ADD REPLYlink modified 12 months ago • written 12 months ago by SMK1.9k
2
gravatar for h.mon
12 months ago by
h.mon29k
Brazil
h.mon29k wrote:

Bowtie2 is a short read mapper and, although it can be used to map long sequences, it probably won't be good at it, specially if the query and reference are even moderately divergent.

In addition to the suggestions by SMK and genomax , there is also LAST.

Another option to consider is QUAST, an assembly evaluation tool. Quast uses minimap2 to align one (or more) query genomes to a reference genome, and in addition to the alignment, it will provide a number of metrics comparing the query to the reference.

ADD COMMENTlink written 12 months ago by h.mon29k
1
gravatar for Corentin
12 months ago by
Corentin450
Corentin450 wrote:

Hi,

You can use D-genies, it can map two fasta to produce a dotplot (http://dgenies.toulouse.inra.fr/)

Alternatives are:

If you are more interested in ordering your contigs, then you could try "show-tiling" from Mummer

ADD COMMENTlink written 12 months ago by Corentin450
0
gravatar for evelyn
11 months ago by
evelyn90
evelyn90 wrote:

You can try nucmer from mummer:

nucmer -p output_prefix ref.fa example-contigs.fa
ADD COMMENTlink written 11 months ago by evelyn90
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: 998 users visited in the last hour