step 1.
cat wy-KO-18d-1.R1.clean.trim.fastq| awk '{if(NR%4==2) print length($1)}' | sort -n | uniq -c > wy-KO-18d-1.R1.length.txt
step 2.
利用R作图
reads<-read.csv(file="wy-het-10d-1.R1.length.txt", sep="", header=FALSE)
KOreads<-read.csv(file="wy-KO-18d-1.R1.length.txt", sep="", header=FALSE)
plot (KOreads$V2,KOreads$V1,type="l",xlab="read length",ylab="fraction",col="blue",lwd=2,main = "18d-KOvsHET reads length distribution")
lines (reads$V2,reads$V1,type="l",xlab="read length",ylab="fraction",col="red",lwd=2)
legend(40,0.5,c("KO","HET"),col=c("blue","red"),text.col = c("blue","red"))
网友评论