Hi biostars,

I've discovered HHBLITS. Not being used to work with alignment files, I find it intricate to work with the output file HHBLITS provides, with extension 'hhr'.

Given an hhr output, I'd just like to get the sequences of the proteins that have an E-value below a threshold. Please note that I do not need the aligned sequences, with gaps and so on.

I'm using python most of the time, but anything will do!

Thanks.