Skip to content
This repository has been archived by the owner on Sep 12, 2018. It is now read-only.

Calculate sequence read depth from GFA files if possible #118

Closed
fedarko opened this issue Dec 31, 2016 · 2 comments
Closed

Calculate sequence read depth from GFA files if possible #118

fedarko opened this issue Dec 31, 2016 · 2 comments
Assignees

Comments

@fedarko
Copy link
Owner

fedarko commented Dec 31, 2016

Figure out how to do it and do it (by just adding the depth info to the .db file).

@fedarko fedarko self-assigned this Dec 31, 2016
@fedarko
Copy link
Owner Author

fedarko commented Dec 31, 2016

Ah, so it looks like depth is calculated using the "read count" optional parameter (RC:i:[some int here]) as the integer in the RC divided by the length of the sequence. Or at least that's how Bandage does that. If a sequence doesn't have a read count parameter then its depth is assumed as 1.

Since I'm not entirely sure what the implications of this optional parameter are or of assuming a uniform read depth of 1 for nodes, I think I'm going to stick with just giving all GFA sequences a "N/A" depth in the viewer for now. I can come back to this later and reevaluate my approach, but there are more important things to focus on for now.

@fedarko
Copy link
Owner Author

fedarko commented Aug 9, 2017

This issue was moved to marbl/MetagenomeScope#23

Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Projects
None yet
Development

No branches or pull requests

1 participant