tissue_of_interest = "Liver"
library(here)
source("/restricted/projectnb/waxmanlab/kkarri/scRNAseq_data_integration/boilerplate.R")
#tiss = load_tissue_droplet(tissue_of_interest)
#library(scater)
library(dplyr)
library(Seurat)
library(cowplot)
#library(MAST)
require(stringr)
require(reshape2)
require(ggplot2)
require(MASS)
library(tools)
require(data.table)
library(ggfortify)
library(tidyverse)
require(dplyr)
library(miscTools)
library(caret)
library(Rtsne)
library(ggrepel)
library(reticulate)
library(hdf5r)

G172h5 <- Read10X_h5("/net/waxman-server/mnt/data/waxmanlabvm_home/kkarri/G172/10X_TCPO_premRNA_Gene/outs/molecule_info.h5",use.names = TRUE, unique.features = TRUE)
G171h5 <- Read10X_h5("/net/waxman-server/mnt/data/waxmanlabvm_home/kkarri/G171/10X_premrna_G171B_mm10_ensembl_onlyPCG_transcript_UseThis/outs/molecule_info.h5",use.names = TRUE, unique.features = TRUE)
write.table(as.matrix(G172h5), "/net/waxman-server/mnt/data/waxmanlabvm_home/kkarri/G172/analysis/demuxEM_HTodemux_analysis/G172h5", sep="\t")
write.table(as.matrix(G171h5), "/net/waxman-server/mnt/data/waxmanlabvm_home/kkarri/G172/analysis/demuxEM_HTodemux_analysis/G171h5", sep="\t")