How to know the positions of contigs in a reference genome
2
1
Entering edit mode
3.4 years ago
Ashi ▴ 20

Hi I am quite new to Bioinformatics and this site has always helped me. I have generated contigs (FASTA file) of E coli from short read assemblers like PaKman, Minia. I want to know the positions of these contigs in the reference genome now. Is there any tool to do that? I have seen the earlier posts here saying that blast can give the positions but I could not get the command line to do so. I know quast gives you the position of the misaligned contig positions but I want to know the positions of the all the contigs that I have got from the short read assembler.

Thanks for all the help.

mapping positions alignment reference genome contigs • 906 views
ADD COMMENT
2
Entering edit mode
3.4 years ago
GenoMax 147k

You can simply align the contigs to E. coli genome. minimap2 or lastz would be good options.

ADD COMMENT
0
Entering edit mode
3.4 years ago
Jordi ▴ 60

Another option that works well for aligning FASTA files is MAFFT (https://mafft.cbrc.jp/alignment/software/).

ADD COMMENT

Login before adding your answer.

Traffic: 2732 users visited in the last hour
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