if(require(pepDat)){
data(pep_hxb2)
hxb2_seq <- metadata(pep_hxb2)$sequence
st<-ProteinSequenceTrack(sequence=hxb2_seq, name="env")
# Plotting amino acids
plotTracks(st, to = 20)
# When the range becomes wider, only coloured squares are displayed
plotTracks(st, to = 100)
# When overplotting, a single line will mark the ProteinSequenceTrack
plotTracks(st)
}
Run the code above in your browser using DataLab