Haplotype assembly using Riemannian trust-region method