codes by  Ben Langmead