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

What happens if you use a smaller Seed? #64

Open
newpest opened this issue Aug 15, 2024 · 3 comments
Open

What happens if you use a smaller Seed? #64

newpest opened this issue Aug 15, 2024 · 3 comments

Comments

@newpest
Copy link

newpest commented Aug 15, 2024

If I manually split a 300kb seed into two 100kb seeds, what is the difference in the AA running results?

@jluebeck
Copy link
Member

Answer really depends on what the seed represents. In many cases, those two seeds will become joined again into one AA amplicon because they share continuously high CN and probably SVs that link them. Feel free to try it out and see for yourself. :)

If you're splitting seeds because of long runtimes, you can also consider trying to set --AA_insert_sdevs 9.

@newpest
Copy link
Author

newpest commented Aug 15, 2024

The seeds generated by CNV detection by Cnvkit seem to have higher resolution than the seeds obtained by readDepth
Seeds by ReadDepth:

8       130540901       131268201       5.796678958040257
8       132769401       133678001       7.199149293492551

Seeds by cnvkit:

8       130336500       133681505       6.785932467648188

Amplicon detected by ReadDepth Seeds has a long low-copy fragment compared to Cnvkit, which ultimately affects the classification results of AC.

@jluebeck
Copy link
Member

Interesting, thanks - if this is a public sample which you can point me to I'm happy to dig in further.

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

2 participants