#Load the classic BLOSUM62 matrix
data(blosum62)
#Set slightly different gap open, and gap extend parameters
b62params <- new("MetricParams", smatrix=blosum62, gapOpen=-8,
gapExtension=-5)
data(TULASequences)
plot(TULASequences, params=b62params)
#Change the metric parameters to a more stringent substitution matrix
data(blosum90)
b90 <- new("MetricParams", smatrix=blosum90, gapOpen=-8,
gapExtension=-5)
plot(TULASequences, params=b90)Run the code above in your browser using DataLab