Hi Peter I changed my sequence parser script to use the SeqIO module and tried your suggestion again but this time looking like coreend = corecur_seq[0]._end.position instead of corestart = corecur_seq[0].position and it works, many thanks for the suggestion Regards Liam