- matplot(x, cbind(peA, peC, peG, peT), type = "b", lty = 1:4, pch = 0:3, main = "Position v.s. Percentage Sequence Error", xlab = "Position", ylab = "Percentage of Sequencing Error")
- legend("topleft", c("A", "C", "G", "T"), lty = 1:4, pch = 0:3)
+ matplot(x, cbind(peA, peC, peG, peT), type = "b", lty = 1:4, pch = 0:3, col = 1:4, main = "Position vs. Percentage Sequence Error", xlab = "Position", ylab = "Percentage of Sequencing Error")
+ legend("topleft", c("A", "C", "G", "T"), lty = 1:4, pch = 0:3, col = 1:4)