Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Code Block
titleExample deepSNV commands
login$login1$ R
...
> regions <- data.frame(chr="gi|254160123|ref|NC_012967.1|", 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")

...