This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to identify gaps in a genome?

Hi, I assembled the genome using PacBio HiFi reads and Hi-C reads. I want to know how can I identify the gaps in genome. Is there any good tool for this? I shall be grateful to you.

awk seqkit assembly fasta genome
• 1,657 views • link
updated by colindaven colindaven • written by rj.rezwan • rj.rezwan

Find related genomes (preferably finished) from NCBI and align your own to get some idea of coverage. Run BUSCO on your genome to see if you have captured most of the gene content expected to be present.

• 0 views • link
written by GenoMax GenoMax

GenoMax I have higher BUSCO than the previously published genome but I am looking to find the "NNNNNN" regions in my genome that exhibit the gaps.

• 0 views • link
written by rj.rezwan • rj.rezwan

1 answer

You can use seqtk to find gaps in a fasta

https://github.com/lh3/seqtk

seqtk gap
Usage: seqtk gap [-l 50] <in.fa>
• 0 views • link
written by colindaven colindaven

Log in to answer this question.

More posts like this