R/extrac_quality.R

Defines functions extrac_quality

extrac_quality = function(sr, index=1:100,rev=FALSE) {
    l=as(PhredQuality(quality(fq)),"IntegerList")
    hist(sapply(l,median),100)


    mat = as.matrix(l)
    boxplot(mat[,1:1000])

}
czhu/R_nanopore documentation built on Dec. 19, 2021, 7:10 p.m.