You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Currently, uniwig can take an input file of either: bed, narrowPeak, bam
and create an output of: wig npy bedGraph bw (via an intermediate bedGraph conversion)
Some items to accomplish for this task:
add bed output support if given a non-bed-like file as input (i.e. bam).
Original bamSitesToWig adjusts positions based on flags. Noodles does support Flags but on first glance they are named differently and must be reviewed.
I believe PEPATAC pipeline only needs starts counts? So add ability to only output any counts the user wishes instead of always outputting all three (starts,ends, core).
Scaling is used in bamSitesToWig.py. Currently, uniwig does not scale counts. Should it?
if associated .bai does not exist, have uniwig create it (using noodles crate)
The text was updated successfully, but these errors were encountered:
donaldcampbelljr
changed the title
Add ability for rust uniwig to create wiggle files from .bam files
Add ability for rust uniwig to create output files from input .bam files
Oct 24, 2024
Some work accomplished with PR #40.
We would like to use this code as a drop in replacement for bamSitesToWig.py from PEPATAC: https://github.com/databio/pepatac/blob/master/tools/bamSitesToWig.py
bamSitesToWig.py creates three files as output:
Currently, uniwig can take an input file of either:
bed
,narrowPeak
,bam
and create an output of:
wig
npy
bedGraph
bw (via an intermediate bedGraph conversion)
Some items to accomplish for this task:
bed
output support if given a non-bed-like file as input (i.e.bam
).starts
counts? So add ability to only output any counts the user wishes instead of always outputting all three (starts,ends, core).https://github.com/databio/pepatac/blob/master/tools/cutsToWig.pl
https://github.com/databio/pepatac/blob/master/tools/smoothWig.pl
Nice to have:
.bai
does not exist, have uniwig create it (using noodles crate)The text was updated successfully, but these errors were encountered: