How to map contigs to the reference genome and know the start and end position?
3
3
Entering edit mode
7.0 years ago
xupbuy ▴ 30

Hi experts,

I want to map my assembled contigs to the reference genome and want to know where every mapped contig starts and ends on the reference genome. And after that, I want to plot it to see where on the genome have more contigs so that I can know the enrichment region.

Is there a way to do this?

Thank you very much.
genome sequencing alignment • 14k views
ADD COMMENT •
0
Entering edit mode

What do you mean assembled contains? Have you made a genome assembly and do you want to map your genome assembly to a reference genome? If so, you can run blast using your genome assembly as query.

ADD REPLY •
0
Entering edit mode

Thank you very much for your reply. After I run Blast, is there a way to make a plot? like this: contig mapping to reference plot

ADD REPLY •
0
Entering edit mode

I don't know why you want to make a plot as you said, but if you want to show a pairwise alignment between contigs and reference assembly, you can have a look mauve tool (http://darlinglab.org/mauve/user-guide/viewer.html)

ADD REPLY •
0
Entering edit mode

I found this tutorial, which can help you.

https://taylorreiter.github.io/2017-08-10-Visualizing-NUCmer-Output/

ADD REPLY •
0
Entering edit mode

I have updated this tutorial, and it now lives at a different link: https://taylorreiter.github.io/2019-05-11-Visualizing-NUCmer-Output/

If you would like to see the original, the source code now lives here: https://github.com/taylorreiter/taylorreiter.github.io/blob/master/_Rmd/viz-nucmer-output.Rmd

ADD REPLY •
3
Entering edit mode
7.0 years ago

I think you could use minimap2 for this.
ADD COMMENT •
0
Entering edit mode

Thank you. minimap2 is like Bowtie, for reads mapping to reference. I want to map contigs to reference and make a plot like this:

ADD REPLY •
3
Entering edit mode

These tools map given nucleotide sequences to a reference, no matter where they come from. The result will be a SAM or BAM file, which contains the respective chromosome, the start position and the length of the aligned sequence. From this, you can filter out the necessary information towards start and end.

ADD REPLY •
0
Entering edit mode

The contigs are in FASTA format, and have no quality score to go into FASTQ.

ADD REPLY •
1
Entering edit mode

You don't need FASTQ or quality scores for alignment, FASTA would do just fine.

ADD REPLY •

Login before adding your answer.

Similar Posts
Loading Similar Posts
Traffic: 2201 users visited in the last hour
Content Search
Users
Tags
Badges
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6