hifiasm icon indicating copy to clipboard operation
hifiasm copied to clipboard

Parameters for plant genomes

Open kkrizanovic opened this issue 3 years ago • 4 comments

Hello,

I know this might not be a real issue, but I didn't know where else I could ask. I was wandering if anyone has any advice in using hifiasm for assembling plant genomes.

I have about 25X coverage with PacBio HiFi reads (and about 100X Hi-C reads). The genome in question is about 4Gbp, its not precisely known, and is diploid. I'm sure what the actual plant name is.

I've run hifiasm with default parameters and got 3.7 Gbp assembly with about 3100 contigs.

Does anyone have any parameters settings that should work well with plants that I could try?

kkrizanovic avatar Sep 24 '21 13:09 kkrizanovic

Do you expect a larger assembly or some other issues?

chhylp123 avatar Sep 24 '21 13:09 chhylp123

I'm not sure about the assembly size, but I would like to get larger contigs :).

kkrizanovic avatar Sep 24 '21 14:09 kkrizanovic

Your HiFi coverage is quite low. This will significantly affect the contiguity.

lh3 avatar Sep 24 '21 14:09 lh3

I often get thousands of contigs with plants, even those with n50 in tens of mB's. And in my experience, only the minus L parameter for purging has a significant effect on the n50. The best parameter to tweak is more or longer hifi data...

On Fri, 24 Sep 2021, 16:00 Krešimir Križanović, @.***> wrote:

I'm not sure about the assembly size, but I would like to get larger contigs :).

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/chhylp123/hifiasm/issues/188#issuecomment-926649069, or unsubscribe https://github.com/notifications/unsubscribe-auth/AARZCFLCDZHGLW56KFTOQ3LUDSAAFANCNFSM5EV7FYRA . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

HenrivdGeest avatar Sep 24 '21 14:09 HenrivdGeest