void check_leaf_sequences(const alignment& A,int n_leaves) { vector<sequence> sequences = A.convert_to_sequences(); const alphabet& a = A.get_alphabet(); for(int i=0;i<n_leaves;i++) { sequences[i].strip_gaps(); if (not (a(sequences[i]) == a(A.seq(i)))) { cerr<<"leaf sequence "<<i<<" corrupted!\n"; cerr<<"orig: "<<A.seq(i)<<endl; cerr<<"new : "<<sequences[i]<<endl; std::abort(); } } }