R commands

 library("DESeq2")
 c<-read.table("m.txt",head=TRUE,row.names = 1)
 colData <-read.table("colData.txt",head=TRUE,row.names=1)
 c<-round(c)
 dds<-DESeqDataSetFromMatrix(countData=c,colData=colData,design= ~ condition)
 dds <-DESeq(dds)
 res  <- results(dds)
 resOrdered <-res[order(res$pvalue),]
 write.csv(as.data.frame(resOrdered),file="set0.csv")
 
 
 plotMA(res,main="MF",ylim=c(-5,5))
 head(resOrdered)
 
 resSig <- subset(resOrdered,padj<0.1)

References

Manuals

https://www.bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.pdf https://www.bioconductor.org/packages/devel/bioc/vignettes/DESeq/inst/doc/DESeq.pdf

Paper

https://genomebiology.biomedcentral.com/articles/10.1186/s13059-014-0550-8