Detecting heterogeneity in single-cell RNA-Seq data by non-negative matrix factorization