Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Effective genome size T2T-CHM13 reference genome #1284

Closed
smsamin opened this issue Jan 4, 2024 · 2 comments
Closed

Effective genome size T2T-CHM13 reference genome #1284

smsamin opened this issue Jan 4, 2024 · 2 comments

Comments

@smsamin
Copy link

smsamin commented Jan 4, 2024

Hi,

I was wondering to ask if its possible to update the effective genome size table (https://deeptools.readthedocs.io/en/develop/content/feature/effectiveGenomeSize.html) for the new human reference genome T2T-CHM13 based on the read length for processed bam files.
Thanks so much in advance!

@ecroot
Copy link

ecroot commented Jan 5, 2024

I would also find this useful please - I have tried running faSize and faCount on a copy of the hs1 genome fasta file, but either due to the "complete" nature of t2t or something I'm doing wrong, I'm getting an effective genome size exactly the same as the genome size

@WardDeb
Copy link
Member

WardDeb commented Jan 10, 2024

Hi both,

The latest develop version has these tables updated (both for the non-N count as the khmer approximation), main branch will be updated when the other features are ready to implement.
Note that T2T indeed has no N's so for the first method the ESS coincides with the actual genome size.

Feel free to reopen if something is still unclear !

WardDeb

@WardDeb WardDeb closed this as completed Jan 10, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants