reads length distribution

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"))