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

Absolute read depth calculation #243

Open
wants to merge 7 commits into
base: main
Choose a base branch
from

Conversation

AustinHartman
Copy link

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.

@@ -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)
Copy link
Author

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.

@AustinHartman
Copy link
Author

Depth in the log will appear as
When long reads and short reads are available:

Calculating read depth (2020-07-30 20:36:55)
--------------------------------------------
Short read depth:   54.48
Long read depth:    22.88
Total read depth:   77.36

When only one type (long or short reads) are available:

Calculating read depth (2020-07-30 20:36:55)
--------------------------------------------
Short read depth:   54.48

Base automatically changed from master to main March 15, 2021 02:55
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

Successfully merging this pull request may close these issues.

1 participant