This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to know the positions of contigs in a reference genome

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
• 1,465 views • link
updated by GenoMax GenoMax • written by Ashi • Ashi

2 answers

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

• 0 views • link
written by GenoMax GenoMax

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

• 0 views • link
written by Jordi • Jordi

Log in to answer this question.

More posts like this