R commands

 # Loading the library
 library("DESeq2")
 # Reading and preparing data
 c<-read.table("m.txt",head=TRUE,row.names = 1)
 c<-round(c)
 colData <-read.table("colData.txt",head=TRUE,row.names=1)

 # Running DESEQ
 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")
 
 # MA plot
 plotMA(res,main="MF",ylim=c(-5,5))

 # Printing differentially expressed genes
 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

Papers

https://genomebiology.biomedcentral.com/articles/10.1186/s13059-014-0550-8 https://genomebiology.biomedcentral.com/articles/10.1186/gb-2010-11-10-r106