This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Circular scaffolds in a genome assembly

Hi!

When I assemble PacBio HiFi reads of a Eukaryotic linear genome through HiFiasm, I get a couple of circular scaffolds. The linear scaffolds are appended with an l in the fasta header, and the circular scaffolds are appended with a c in the fasta header.

On what basis is l and c assigned to these scaffolds?

Thanks!

assembly pacbio genome hifi hifiasm
• 1,622 views • link
updated by gconcepcion gconcepcion • written by Shriram • Shriram

2 answers

Quite possible that you have assembled your organelle genome as well. You can verify that by doing a homologous search. When the assembler assembles sequence and finds support of reads to be able to circularise by overlapping the ends it adds _c in the fasta header. Not always the circular sequences are the organelles, the repeat region of your genome also comes as circular sequences because of possibility of overlaps end to end because of repeat nature.

• 0 views • link
written by bagdevi.mishra bagdevi.mishra

To elaborate slightly on bagdevi.mishra's reply which I agree with, try taking one or more of the circular sequences and BLAST'ing them (or part of them if they are too large) to see what sort of hits come up in NCBI's database: https://blast.ncbi.nlm.nih.gov/Blast.cgi

• 0 views • link
written by gconcepcion gconcepcion

Log in to answer this question.

More posts like this