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

Add a QC command that calculates basic statistics on a ft bam #59

Merged
merged 30 commits into from
Aug 22, 2024
Merged

Conversation

mrvollger
Copy link
Member

@mrvollger mrvollger commented Aug 20, 2024

The idea is to create text histograms that can be used as inputs to plotting programs.

The current output idea is a text file that looks like:

statistic       value   count
fiber_lengths   1000     58
fiber_lengths   1001     100
...
read_quality   0.98   1001
...
m6a_ratio    0.06    1002
...

I'm adding a new stat that is the number of m6as in MSPs of specific sizes, whether that MSP is a FIRE element, and the number of times that was seen.

statistic       value   count
...
m6a_per_msp_size        1,1,false       4451629
m6a_per_msp_size        2,2,false       123829
m6a_per_msp_size        2,3,false       347430
m6a_per_msp_size        2,4,false       306585
m6a_per_msp_size        2,5,false       189577
m6a_per_msp_size        2,6,false       182797
m6a_per_msp_size        2,7,false       169225
...

This output is added when the option --m6a-per-msp is used.

@mrvollger mrvollger merged commit 50f1012 into main Aug 22, 2024
9 checks passed
@mrvollger mrvollger deleted the qc branch August 22, 2024 14:42
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