miRNA_isoform <- read.table(miRNA_isoform_file, header = TRUE, stringsAsFactors = FALSE)
miRNA_isoform$unique_id <- paste(miRNA_isoform$miRNA_ID, miRNA_isoform$miRNA_region, sep = "_")

miRNA_isoform_unique_id <- unique(miRNA_isoform$unique_id)
miRNA_isoform_sample_id <- unique(miRNA_isoform$SampleId)


miRNA_isoform_reads <- matrix(data = 0, ncol = length(miRNA_isoform_sample_id), nrow = length(miRNA_isoform_unique_id))
colnames(miRNA_isoform_reads) <- miRNA_isoform_sample_id
rownames(miRNA_isoform_reads) <- miRNA_isoform_unique_id

for (i in rownames(miRNA_isoform_reads)){
  for (j in colnames(miRNA_isoform_reads)){
    miRNA_isoform_reads[i,j] <- sum(miRNA_isoform$read_count[miRNA_isoform$unique_id == i & miRNA_isoform$SampleId == j], na.rm = TRUE)
  }
}