In 2D you can just use pylab's triplot
from dolfin import *
from mshr import *
from pylab import triplot, show, hold, plot
c1 = Circle(Point(0.,0.),0.5)
r1 = Rectangle(Point(0.,0.),Point(1.,1.))
mesh = generate_mesh(r1-c1,15)
triplot(mesh.coordinates()[:,0],mesh.coordinates()[:,1],\
triangles=mesh.cells(),color='k')
hold('on'); plot([0.2,0.8],[0.9,0.2],'bo-',\
linewidth=2,markersize=16) ;hold('off')
show()