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

Maximum number of samples and Time to Process #2

Open
syngnathid opened this issue Sep 2, 2022 · 3 comments
Open

Maximum number of samples and Time to Process #2

syngnathid opened this issue Sep 2, 2022 · 3 comments

Comments

@syngnathid
Copy link

Hello SRY team,

Thank you for creating an useful tool. I am currently trying to use it to assemble a better genome that includes sex chromosomes. I was curious about few things that I couldn't find in the readme/documentation. Could you share some insights please?

  1. Could I use 100 samples (50 males and 50 females) when I run SRY program?
  2. What is the relationship between the number of the samples, the size/quality of the genome and time taken to identify the Y chromosome reads?

Context:
I have a genome of a non-model organism with about ~450 contigs of which 21 are chromosome level scaffolds. This genome is from a male sample and is a heterogametic system (XY). Ideally. I would like to identify and stitch the sex chromosomes if possible. The genome was assembled using Hi-Fi and Hi-C data. I have Whole Genome Resequencing data from Illumina at ~30X coverage for about 50 males and 50 females .

@caaswxb
Copy link
Owner

caaswxb commented Sep 6, 2022

Hi,
As the sample number increases, the speciality is up but slow when larger than ~5 in human (see the biorxiv paper), and the sensitivity is down yet. So for good result and shorter running time, I suggest 10 samples (10 M and 10 F). The script SRY sorts ONT or PB reads, without HiFi. You can stop the SRY when the SRY_kmer.txt is generated and then deliver it to the script SRY_contig. The HiFi sorting function will be added in SRY latter.

@Johnsonzcode
Copy link

Johnsonzcode commented Jun 4, 2024

Sorry to bother, but I ran with 8 males and 7 females with 80 cores, it's been 15 days. I want to speed up, is there any suggestions?

@caaswxb
Copy link
Owner

caaswxb commented Jul 13, 2024 via email

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