...
Code Block | ||
---|---|---|
| ||
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") |
...