diff --git a/gtk2_ardour/midi_region_view.cc b/gtk2_ardour/midi_region_view.cc
index d76c173af1..40f267717b 100644
--- a/gtk2_ardour/midi_region_view.cc
+++ b/gtk2_ardour/midi_region_view.cc
@@ -2368,17 +2368,12 @@ MidiRegionView::cut_copy_clear (Editing::CutCopyOp op)
 MidiCutBuffer*
 MidiRegionView::selection_as_cut_buffer () const
 {
-	NoteList notes;
+	Notes notes;
 
 	for (Selection::iterator i = _selection.begin(); i != _selection.end(); ++i) {
-		notes.push_back (boost::shared_ptr<NoteType> (new NoteType (*((*i)->note().get()))));
+		notes.insert (boost::shared_ptr<NoteType> (new NoteType (*((*i)->note().get()))));
 	}
 
-	/* sort them into time order */
-
-	Evoral::Sequence<Evoral::MusicalTime>::LaterNoteComparator cmp;
-	sort (notes.begin(), notes.end(),  cmp);
-
 	MidiCutBuffer* cb = new MidiCutBuffer (trackview.session());
 	cb->set (notes);
 
@@ -2399,16 +2394,16 @@ MidiRegionView::paste (nframes64_t pos, float times, const MidiCutBuffer& mcb)
 	Evoral::MusicalTime duration;
 	Evoral::MusicalTime end_point;
 
-	duration = mcb.notes().back()->end_time() - mcb.notes().front()->time();
+	duration = (*mcb.notes().rbegin())->end_time() - (*mcb.notes().begin())->time();
 	paste_pos_beats = frames_to_beats (pos - _region->position());
-	beat_delta = mcb.notes().front()->time() - paste_pos_beats;
+	beat_delta = (*mcb.notes().begin())->time() - paste_pos_beats;
 	paste_pos_beats = 0;
 
 	_selection.clear ();
 
 	for (int n = 0; n < (int) times; ++n) {
 
-		for (NoteList::const_iterator i = mcb.notes().begin(); i != mcb.notes().end(); ++i) {
+		for (Notes::const_iterator i = mcb.notes().begin(); i != mcb.notes().end(); ++i) {
 
 			boost::shared_ptr<NoteType> copied_note (new NoteType (*((*i).get())));
 			copied_note->set_time (paste_pos_beats + copied_note->time() - beat_delta);
@@ -2516,13 +2511,13 @@ MidiRegionView::goto_previous_note ()
 }
 
 void
-MidiRegionView::selection_as_notelist (NoteList& selected)
+MidiRegionView::selection_as_notelist (Notes& selected)
 {
 	time_sort_events ();
 
 	for (Events::iterator i = _events.begin(); i != _events.end(); ++i) {
 		if ((*i)->selected()) {
-			selected.push_back ((*i)->note());
+			selected.insert ((*i)->note());
 		}
 	}
 }
diff --git a/gtk2_ardour/midi_region_view.h b/gtk2_ardour/midi_region_view.h
index 36b7c99f54..22c2b845fd 100644
--- a/gtk2_ardour/midi_region_view.h
+++ b/gtk2_ardour/midi_region_view.h
@@ -65,7 +65,7 @@ class MidiRegionView : public RegionView
 {
   public:
 	typedef Evoral::Note<Evoral::MusicalTime> NoteType;
-	typedef Evoral::Sequence<Evoral::MusicalTime>::Notes NoteList;
+	typedef Evoral::Sequence<Evoral::MusicalTime>::Notes Notes;
 
 	MidiRegionView (ArdourCanvas::Group *,
 	                RouteTimeAxisView&,
@@ -285,7 +285,7 @@ class MidiRegionView : public RegionView
 
 	void show_list_editor ();
 
-	void selection_as_notelist (NoteList& selected);
+	void selection_as_notelist (Notes& selected);
 
   protected:
 	/** Allows derived types to specify their visibility requirements
diff --git a/gtk2_ardour/midi_streamview.cc b/gtk2_ardour/midi_streamview.cc
index 01297abf89..cb5537d601 100644
--- a/gtk2_ardour/midi_streamview.cc
+++ b/gtk2_ardour/midi_streamview.cc
@@ -574,11 +574,10 @@ MidiStreamView::update_rec_regions (boost::shared_ptr<MidiModel> data, nframes_t
 
 						/* draw events */
 						MidiRegionView* mrv = (MidiRegionView*)iter->second;
-						for (size_t i=0; i < data->n_notes(); ++i) {
+						for (MidiModel::Notes::const_iterator i = data->notes().begin();
+								i != data->notes().end(); ++i) {
 
-							// FIXME: slooooooooow!
-
-							const boost::shared_ptr<MidiRegionView::NoteType> note = data->note_at(i);
+							const boost::shared_ptr<MidiRegionView::NoteType>& note = *i;
 
 							if (note->length() > 0 && note->end_time() + region->position() > start)
 								mrv->resolve_note(note->note(), note->end_time());
@@ -598,7 +597,6 @@ MidiStreamView::update_rec_regions (boost::shared_ptr<MidiModel> data, nframes_t
 							}
 
 							mrv->add_note (note, !update_range);
-
 						}
 
 						mrv->extend_active_notes();
diff --git a/libs/evoral/evoral/Sequence.hpp b/libs/evoral/evoral/Sequence.hpp
index b4e11ba808..22d81b6948 100644
--- a/libs/evoral/evoral/Sequence.hpp
+++ b/libs/evoral/evoral/Sequence.hpp
@@ -61,7 +61,7 @@ public:
 template<typename Time>
 class Sequence : virtual public ControlSet {
 public:
-	Sequence(const TypeMap& type_map, size_t size=0);
+	Sequence(const TypeMap& type_map);
 
 	void write_lock();
 	void write_unlock();
@@ -78,12 +78,8 @@ public:
 	bool writing() const { return _writing; }
 	void end_write(bool delete_stuck=false);
 
-	/** Resizes vector if necessary (NOT realtime safe) */
 	void append(const Event<Time>& ev);
 
-	inline const boost::shared_ptr< const Note<Time> > note_at(size_t i) const { return _notes[i]; }
-	inline const boost::shared_ptr< Note<Time> >       note_at(size_t i)       { return _notes[i]; }
-
 	inline size_t n_notes() const { return _notes.size(); }
 	inline bool   empty()   const { return _notes.size() == 0 && ControlSet::controls_empty(); }
 
@@ -92,6 +88,20 @@ public:
 		return a->time() < b->time();
 	}
 
+	struct NoteNumberComparator {
+		inline bool operator()(const boost::shared_ptr< const Note<Time> > a,
+		                       const boost::shared_ptr< const Note<Time> > b) const {
+			return a->note() < b->note();
+		}
+	};
+
+	struct EarlierNoteComparator {
+		inline bool operator()(const boost::shared_ptr< const Note<Time> > a,
+		                       const boost::shared_ptr< const Note<Time> > b) const {
+			return a->time() < b->time();
+		}
+	};
+
 	struct LaterNoteComparator {
 		typedef const Note<Time>* value_type;
 		inline bool operator()(const boost::shared_ptr< const Note<Time> > a,
@@ -108,7 +118,7 @@ public:
 		}
 	};
 
-	typedef std::vector< boost::shared_ptr< Note<Time> > > Notes;
+	typedef std::set<boost::shared_ptr< Note<Time> >, EarlierNoteComparator> Notes;
 	inline       Notes& notes()       { return _notes; }
 	inline const Notes& notes() const { return _notes; }
 
@@ -199,7 +209,7 @@ private:
 	Notes   _notes;
 	SysExes _sysexes;
 
-	typedef std::vector<size_t> WriteNotes;
+	typedef std::set<boost::shared_ptr< Note<Time> >, NoteNumberComparator> WriteNotes;
 	WriteNotes _write_notes[16];
 	bool       _writing;
 
diff --git a/libs/evoral/src/Sequence.cpp b/libs/evoral/src/Sequence.cpp
index 124a0acd4c..81d515afbf 100644
--- a/libs/evoral/src/Sequence.cpp
+++ b/libs/evoral/src/Sequence.cpp
@@ -427,10 +427,9 @@ Sequence<Time>::const_iterator::operator=(const const_iterator& other)
 // Sequence
 
 template<typename Time>
-Sequence<Time>::Sequence(const TypeMap& type_map, size_t size)
+Sequence<Time>::Sequence(const TypeMap& type_map)
 	: _edited(false)
 	, _type_map(type_map)
-	, _notes(size)
 	, _writing(false)
 	, _end_iter(*this, DBL_MAX)
 	, _percussive(false)
@@ -569,14 +568,13 @@ Sequence<Time>::end_write(bool delete_stuck)
 
 	if (!_percussive && delete_stuck) {
 		for (typename Notes::iterator n = _notes.begin(); n != _notes.end() ;) {
+			typename Notes::iterator next = n;
+			++next;
 			if ((*n)->length() == 0) {
 				cerr << "WARNING: Stuck note lost: " << (*n)->note() << endl;
-				n = _notes.erase(n);
-				// we have to break here because erase invalidates the iterator
-				break;
-			} else {
-				++n;
+				_notes.erase(n);
 			}
+			n = next;
 		}
 	}
 
@@ -611,7 +609,7 @@ Sequence<Time>::append(const Event<Time>& event)
 
 	const MIDIEvent<Time>& ev = (const MIDIEvent<Time>&)event;
 
-	assert(_notes.empty() || ev.time() >= _notes.back()->time());
+	assert(_notes.empty() || ev.time() >= (*_notes.rbegin())->time());
 	assert(_writing);
 
 	if (!midi_event_is_valid(ev.buffer(), ev.size())) {
@@ -678,14 +676,14 @@ Sequence<Time>::append_note_on_unlocked(uint8_t chan, Time time, uint8_t note_nu
 		_highest_note = note_num;
 
 	boost::shared_ptr< Note<Time> > new_note(new Note<Time>(chan, time, 0, note_num, velocity));
-	_notes.push_back(new_note);
+	_notes.insert(new_note);
 	if (!_percussive) {
 		DUMP(format("Sustained: Appending active note on %1% channel %2%\n")
 				% (unsigned)(uint8_t)note_num % chan);
-		_write_notes[chan].push_back(_notes.size() - 1);
+		_write_notes[chan].insert(new_note);
 	} else {
-	 	DUMP("Percussive: NOT appending active note on\n");
-	 }
+		DUMP("Percussive: NOT appending active note on\n");
+	}
 }
 
 template<typename Time>
@@ -704,19 +702,15 @@ Sequence<Time>::append_note_off_unlocked(uint8_t chan, Time time, uint8_t note_n
 		return;
 	}
 
-	/* FIXME: make _write_notes fixed size (127 noted) for speed */
-
-	/* FIXME: note off velocity for that one guy out there who actually has
-	 * keys that send it */
+	// TODO: support note off velocity
 
 	bool resolved = false;
-
-	for (WriteNotes::iterator n = _write_notes[chan].begin(); n
-			!= _write_notes[chan].end(); ++n) {
-		Note<Time>& note = *_notes[*n].get();
-		if (note.note() == note_num) {
-			assert(time >= note.time());
-			note.set_length(time - note.time());
+	for (typename WriteNotes::iterator n = _write_notes[chan].begin();
+			n != _write_notes[chan].end(); ++n) {
+		boost::shared_ptr< Note<Time> > note = *n;
+		if (note->note() == note_num) {
+			assert(time >= note->time());
+			note->set_length(time - note->time());
 			_write_notes[chan].erase(n);
 			DUMP(format("resolved note, length: %1%\n") % note.length());
 			resolved = true;
diff --git a/libs/evoral/test/SMFTest.hpp b/libs/evoral/test/SMFTest.hpp
index cc5feb684c..77ccb54572 100644
--- a/libs/evoral/test/SMFTest.hpp
+++ b/libs/evoral/test/SMFTest.hpp
@@ -60,7 +60,7 @@ class SMFTest : public CppUnit::TestFixture
         void setUp() { 
            type_map = new DummyTypeMap();
            assert(type_map);
-           seq = new MySequence<Time>(*type_map, 0);
+           seq = new MySequence<Time>(*type_map);
            assert(seq);
         }
         
diff --git a/libs/evoral/test/SequenceTest.hpp b/libs/evoral/test/SequenceTest.hpp
index 9c6f084643..e5e5015bfc 100644
--- a/libs/evoral/test/SequenceTest.hpp
+++ b/libs/evoral/test/SequenceTest.hpp
@@ -54,7 +54,7 @@ public:
 template<typename Time>
 class MySequence : public Sequence<Time> {
 public:
-	MySequence(DummyTypeMap&map, int size) : Sequence<Time>(map, size) {}
+	MySequence(DummyTypeMap&map) : Sequence<Time>(map) {}
 	
 	boost::shared_ptr<Control> control_factory(const Parameter& param) {
 		
@@ -107,7 +107,7 @@ class SequenceTest : public CppUnit::TestFixture
         void setUp (void) { 
         	type_map = new DummyTypeMap();
            	assert(type_map);
-           	seq = new MySequence<Time>(*type_map, 0);
+           	seq = new MySequence<Time>(*type_map);
            	assert(seq);
            	
            	for(int i = 0; i < 12; i++) {