cmdoret / tinycov

Command line tool to plot genomic coverage from a BAM file
Other
13 stars 5 forks source link

feature request #7

Open parlar opened 1 year ago

parlar commented 1 year ago

Thank you for creating tinycov, I really like the concept!!

I notice some cosmetic issues with the output but I'm sure you're aware.

But one perhaps useful thing: Would it be possible allow import of cov data from another source, such as d4 or bedgraph? Then something like mosdepth could be used to get the coverage data, which might be quicker than using pysam?

Good work!

parlar commented 1 year ago

just a note:

Played around with d4. it's very quick.

import sys
import fire
from pyd4 import D4File

def main(d4_file: str):
    file = D4File(d4_file)

    down_sampled_chr1 = file.resample("chr1", bin_size = 10000)
    for v in down_sampled_chr1[0]:
        print(v)

if __name__ == '__main__':
    fire.Fire(main)
cmdoret commented 1 year ago

Hi @parlar thanks for the suggestion! I was not aware of d4 but that looks useful. It would make sense to allow d4/bedgraph inputs, but the code will need a little refactor to disentangle ingestion from plotting.

I'll look into this :)