Skip to content

Commit

Permalink
cherrypick: manta -> sv (#173)
Browse files Browse the repository at this point in the history
  • Loading branch information
pdiakumis committed Oct 3, 2024
1 parent 5286e06 commit 680a7c0
Showing 1 changed file with 6 additions and 6 deletions.
12 changes: 6 additions & 6 deletions inst/rmd/rnasum.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -1418,18 +1418,18 @@ fusions_ann <- RNAsum::fusions_annot(
)
```

```{r fusions_and_manta_data_prep, comment = NA, message=FALSE, warning=FALSE, eval = runSVsChunk}
manta_sv <- sample_data.list[["wgs"]][["manta_tsv"]] |>
RNAsum::manta_process()
tmp_sv_genes <- manta_sv$genes_all$Genes
```{r fusions_and_sv_data_prep, comment = NA, message=FALSE, warning=FALSE, eval = runSVsChunk}
sv <- sample_data.list[["wgs"]][["sv_tsv"]] |>
RNAsum::sv_process()
tmp_sv_genes <- sv$genes_all$Genes
tmp_expressed_sv_genes <- tmp_sv_genes[tmp_sv_genes %in% ref_dataset.list[[dataset]][["gene_annot_all"]]$SYMBOL]
if (length(tmp_expressed_sv_genes) == 0) {
# from now on don't evaluate SV chunks
runSVsChunk <- FALSE
}
rm(tmp_sv_genes, tmp_expressed_sv_genes)
##### Compare fusion genes called by MANTA
##### First limit MANTA output to fusions only
##### Compare fusion genes called by SV caller
##### First limit SV output to fusions only
if (runFusionChunk) {
sv_fusions <- sv[["melted_variants"]] |>
dplyr::filter(.data$is_fusion) |>
Expand Down

0 comments on commit 680a7c0

Please sign in to comment.