Added NJ negative branch calculation 

Added a function that recalculates the negative branch lengths as proposed by Kuhner & Felsenstein (1994).