-
Notifications
You must be signed in to change notification settings - Fork 135
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
Absolute read depth calculation #243
base: main
Are you sure you want to change the base?
Conversation
@@ -222,6 +223,9 @@ def main(): | |||
graph.save_to_gfa(final_assembly_gfa) | |||
graph.save_to_fasta(final_assembly_fasta, min_length=args.min_fasta_length) | |||
|
|||
# Align reads to assembly.fasta and report depth | |||
get_read_depth(final_assembly_fasta, args) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Due to the reliance on assembly.fasta, read depth is calculated at the very end of the pipeline. Thus, the log header "Calculating read depth" and the depth information comes after the "Assembly complete" header. Maybe this behavior should be modified.
Depth in the log will appear as
When only one type (long or short reads) are available:
|
This doesn't take all the long to compute, and can be useful information to know. Also, databases like GenBank require read depth/coverage when assemblies are submitted.