Skip to content

Instantly share code, notes, and snippets.

@kdauria
kdauria / read_bam.rs
Created July 22, 2023 22:10
Read BAM with noodles that has a non-compliant platform tag
use std::env;
use std::fs::File;
use std::io;
use std::io::Read;
use std::path::Path;
use std::process;
use byteorder::{LittleEndian, ReadBytesExt};
use noodles::bam;
use noodles::bgzf::{self as bgzf, VirtualPosition};
stat_smooth_func <- function(mapping = NULL, data = NULL,
geom = "smooth", position = "identity",
...,
method = "auto",
formula = y ~ x,
se = TRUE,
n = 80,
span = 0.75,
fullrange = FALSE,
level = 0.95,