...
Code Block | ||
---|---|---|
| ||
login1$ R
...
> library("deepSNV")
> regions <- data.frame(chr="NC_012967", start = 1, stop=100000)
> result = deepSNV(test = "SRR032376.sorted.bam", control = "SRR030252.sorted.bam", regions=regions)
> sig_result <- summary(result, sig.level=0.05, adjust.method="BH")
> pdf("SRR032374.pdf")
> plot(mix)
> dev.off()
> write.csv(sig_result, "SRR032374")
|
...