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.
3 answers
I think you could use minimap2 for this.
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.
You don't need FASTQ or quality scores for alignment, FASTA would do just fine.
Log in to answer this question.
-
Keep the sequencing depth info when mapping contigs in Minimap2
written by pauline10albert •
-
How to assess the quality of the assembly of unmapped reads ?
written by Sony Nguyen •
-
Export sequence from specific region from bam file.
written by malkomalktho •
-
How to plot the gap distribution of contigs wrt to reference genome
written by Ashi •
-
How to evaluate and check the correctness of the assembly?
written by JXYang •
-
BWA does contig orientation matter?
written by dp •
-
How to get a consensus of the contigs obtained from an assembler like spades?
written by DanielC
-
understanding contig statistics
written by mary •
-
Looking for a genome /draft genome aligner with the these features
written by novice
-
Plasmid Assembly Use Of Cisa Contig Intergrator
written by bobo •
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.
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
Plot like this:
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)
I found this tutorial, which can help you.
https://taylorreiter.github.io/2017-08-10-Visualizing-NUCmer-Output/
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