SigProfilerAssignment 0.2.1__tar.gz → 0.2.3__tar.gz

{sigprofilerassignment-0.2.1/SigProfilerAssignment.egg-info → sigprofilerassignment-0.2.3}/PKG-INFO

@@ -1,6 +1,6 @@
-Metadata-Version: 2.2
+Metadata-Version: 2.4
 Name: SigProfilerAssignment
-Version: 0.2.1
+Version: 0.2.3
 Summary: Mutational signatures attribution and decomposition tool
 Home-page: https://github.com/AlexandrovLab/SigProfilerAssignment.git
 Author: Raviteja Vangara
@@ -17,7 +17,7 @@ Requires-Dist: sigProfilerPlotting>=1.4.0
 Requires-Dist: reportlab>=3.5.42
 Requires-Dist: pypdf>=5.0.0
 Requires-Dist: alive_progress>=2.4.1
-Requires-Dist: PyMuPDF>=1.21.0
+Requires-Dist: pdf2image>=1.16.0
 Provides-Extra: tests
 Requires-Dist: pytest; extra == "tests"
 Dynamic: author
@@ -26,6 +26,7 @@ Dynamic: description
 Dynamic: description-content-type
 Dynamic: home-page
 Dynamic: license
+Dynamic: license-file
 Dynamic: provides-extra
 Dynamic: requires-dist
 Dynamic: requires-python
@@ -67,6 +68,12 @@ $ python
 from SigProfilerMatrixGenerator import install as genInstall
 genInstall.install('GRCh37')
 ```
+
+If you plan to use `sample_reconstruction_plots='png'` or `'both'`, the external `poppler` binary is required. You can install it using one of the following commands:
+
+- For Conda-based environments:
+  `conda install -c conda-forge poppler`
+
 ## <a name="running"></a> Running
 
 Assignment of known mutational signatures to individual samples is performed using the `cosmic_fit` function. Input samples are provided using the `samples` parameter in the form of mutation calling files (VCFs, MAFs, or simple text files), segmentation files or mutational matrices. COSMIC mutational signatures v3.4 are used as the default reference signatures, although previous COSMIC versions and custom signature databases are also supported using the `cosmic_version` and `signature_database` parameters. Results will be found in the folder specified in the `output` parameter.
@@ -99,7 +106,7 @@ Analyze.cosmic_fit(samples, output, input_type="matrix", context_type="96",
 | export_probabilities | Boolean | Defines if the probability matrix per mutational context for all samples is created. The default value is True. |
 | export_probabilities_per_mutation | Boolean | Defines if the probability matrices per mutation for all samples are created. Only available when `input_type` is "vcf". The default value is False. |
 | make_plots | Boolean | Toggle on and off for making and saving plots. The default value is True. |
-| sample_reconstruction_plots | String | Select the output format for sample reconstruction plots. Valid inputs are {'pdf', 'png', 'both', None}. The default value is None. |
+| sample_reconstruction_plots | String | Select the output format for sample reconstruction plots. Valid inputs are {'pdf', 'png', 'both', 'none'}. The default value is 'none'. If set to 'png' or 'both', the external binary `poppler` must be installed. Install via `conda install -c conda-forge poppler` or `brew install poppler` on macOS. |
 | verbose | Boolean | Prints detailed statements. The default value is False. |
 | volume | String | Path to SigProfilerAssignment volumes. Used for Docker/Singularity. Environmental variable "SIGPROFILERASSIGNMENT_VOLUME" takes precedence. Default value is None. |
 

{sigprofilerassignment-0.2.1 → sigprofilerassignment-0.2.3}/README.md

@@ -34,6 +34,12 @@ $ python
 from SigProfilerMatrixGenerator import install as genInstall
 genInstall.install('GRCh37')
 ```
+
+If you plan to use `sample_reconstruction_plots='png'` or `'both'`, the external `poppler` binary is required. You can install it using one of the following commands:
+
+- For Conda-based environments:
+  `conda install -c conda-forge poppler`
+
 ## <a name="running"></a> Running
 
 Assignment of known mutational signatures to individual samples is performed using the `cosmic_fit` function. Input samples are provided using the `samples` parameter in the form of mutation calling files (VCFs, MAFs, or simple text files), segmentation files or mutational matrices. COSMIC mutational signatures v3.4 are used as the default reference signatures, although previous COSMIC versions and custom signature databases are also supported using the `cosmic_version` and `signature_database` parameters. Results will be found in the folder specified in the `output` parameter.
@@ -66,7 +72,7 @@ Analyze.cosmic_fit(samples, output, input_type="matrix", context_type="96",
 | export_probabilities | Boolean | Defines if the probability matrix per mutational context for all samples is created. The default value is True. |
 | export_probabilities_per_mutation | Boolean | Defines if the probability matrices per mutation for all samples are created. Only available when `input_type` is "vcf". The default value is False. |
 | make_plots | Boolean | Toggle on and off for making and saving plots. The default value is True. |
-| sample_reconstruction_plots | String | Select the output format for sample reconstruction plots. Valid inputs are {'pdf', 'png', 'both', None}. The default value is None. |
+| sample_reconstruction_plots | String | Select the output format for sample reconstruction plots. Valid inputs are {'pdf', 'png', 'both', 'none'}. The default value is 'none'. If set to 'png' or 'both', the external binary `poppler` must be installed. Install via `conda install -c conda-forge poppler` or `brew install poppler` on macOS. |
 | verbose | Boolean | Prints detailed statements. The default value is False. |
 | volume | String | Path to SigProfilerAssignment volumes. Used for Docker/Singularity. Environmental variable "SIGPROFILERASSIGNMENT_VOLUME" takes precedence. Default value is None. |
 

{sigprofilerassignment-0.2.1 → sigprofilerassignment-0.2.3}/SigProfilerAssignment/controllers/cli_controller.py

@@ -135,6 +135,19 @@ def parse_arguments_common(args: List[str], description: str) -> argparse.Namesp
         default=None,
         help="User specified directory for saving/loading template files. Note: The environment variable SIGPROFILERASSIGNMENT_VOLUME takes precedence over this parameter.",
     )
+    parser.add_argument(
+        "--sample_reconstruction_plots",
+        type=str.lower,
+        choices=["none", "pdf", "both", "png"],
+        default="none",
+        help=(
+            "Output format for sample reconstruction plots. "
+            "Options: 'none' (default, disables plotting), "
+            "'pdf' (generate only PDF), "
+            "'both' (PDF + PNG), or "
+            "'png' (PNG only, PDF removed)."
+        )
+    )
 
     return parser.parse_args(args)
 
@@ -226,5 +239,5 @@ class CliController:
             context_type=parsed_args.context_type,
             export_probabilities=parsed_args.export_probabilities,
             export_probabilities_per_mutation=parsed_args.export_probabilities_per_mutation,
-            sample_reconstruction_plots=False,
+            sample_reconstruction_plots=parsed_args.sample_reconstruction_plots,
         )

{sigprofilerassignment-0.2.1 → sigprofilerassignment-0.2.3}/SigProfilerAssignment/decompose_subroutines.py

@@ -822,154 +822,234 @@ def make_final_solution(
         check_rule_negatives = []
         check_rule_penalty = 1.0
     exposureAvg = np.zeros([processAvg.shape[1], allgenomes.shape[1]])
+    
     if cosmic_sigs == True:
         denovo_exposureAvg = denovo_exposureAvg.T
-        with alive_bar(allgenomes.shape[1]) as bar:
+        # with alive_bar(allgenomes.shape[1]) as bar:
             # print("\n")
-            for r in range(allgenomes.shape[1]):
-                # print("Analyzing Sample => " , str(r+1))
-                bar()
-                if verbose == True:
-                    print(
-                        "\n\n\n\n\n                                        ################ Sample "
-                        + str(r + 1)
-                        + " #################"
-                    )
-
-                # Record information to lognote
-                lognote = open(
-                    layer_directory
-                    + "/Solution_Stats/"
-                    + solution_prefix
-                    + "_Signature_Assignment_log.txt",
-                    "a",
-                )
-                lognote.write(
-                    "\n\n\n\n\n                    ################ Sample "
+        for r in range(allgenomes.shape[1]):
+            # print("Analyzing Sample => " , str(r+1))
+            # bar()
+            if verbose == True:
+                print(
+                    "\n\n\n\n\n                                        ################ Sample "
                     + str(r + 1)
-                    + " #################\n"
+                    + " #################"
                 )
 
-                sample_exposure = np.array(denovo_exposureAvg.iloc[:, r])
-
-                init_sig_idx = np.nonzero(sample_exposure)[0]
-                init_sigs = denovo_exposureAvg.index[init_sig_idx]
+            # Record information to lognote
+            lognote = open(
+                layer_directory
+                + "/Solution_Stats/"
+                + solution_prefix
+                + "_Signature_Assignment_log.txt",
+                "a",
+            )
+            lognote.write(
+                "\n\n\n\n\n                    ################ Sample "
+                + str(r + 1)
+                + " #################\n"
+            )
 
-                init_decomposed_sigs = []
-                for de_novo_sig in init_sigs:
-                    init_decomposed_sigs = union(
-                        init_decomposed_sigs, list(attribution[de_novo_sig])
-                    )
+            sample_exposure = np.array(denovo_exposureAvg.iloc[:, r])
 
-                # print(init_decomposed_sigs)
-                init_decomposed_sigs_idx = get_indeces(allsigids, init_decomposed_sigs)
-                init_decomposed_sigs_idx.sort()
-                init_decomposed_sigs_idx = list(
-                    set().union(init_decomposed_sigs_idx, background_sigs)
-                )
-                # print(init_decomposed_sigs_idx)
+            init_sig_idx = np.nonzero(sample_exposure)[0]
+            init_sigs = denovo_exposureAvg.index[init_sig_idx]
 
-                # get the indices of the background sigs in the initial signatures
-                background_sig_idx = get_indeces(
-                    init_decomposed_sigs_idx, background_sigs
+            init_decomposed_sigs = []
+            for de_novo_sig in init_sigs:
+                init_decomposed_sigs = union(
+                    init_decomposed_sigs, list(attribution[de_novo_sig])
                 )
 
-                fit_signatures = processAvg[:, init_decomposed_sigs_idx]
-                # fit signatures
-                newExposure, newSimilarity = ss.fit_signatures(
-                    fit_signatures, allgenomes[:, r]
-                )
+            # print(init_decomposed_sigs)
+            init_decomposed_sigs_idx = get_indeces(allsigids, init_decomposed_sigs)
+            init_decomposed_sigs_idx.sort()
+            init_decomposed_sigs_idx = list(
+                set().union(init_decomposed_sigs_idx, background_sigs)
+            )
+            # print(init_decomposed_sigs_idx)
 
-                # create the exposureAvg vector
-                # print(init_decomposed_sigs_idx)
-                # print(newExposure)
-                for nonzero_idx, nozero_exp in zip(
-                    init_decomposed_sigs_idx, newExposure
-                ):
-                    exposureAvg[nonzero_idx, r] = nozero_exp
+            # get the indices of the background sigs in the initial signatures
+            background_sig_idx = get_indeces(
+                init_decomposed_sigs_idx, background_sigs
+            )
 
-                if pcawg_rule == True:
-                    maxmutation = np.sum(allgenomes[:, r])
-                    (
-                        exposureAvg[:, r],
-                        remove_distance,
-                        _,
-                    ) = ss.remove_all_single_signatures(
-                        processAvg,
-                        exposureAvg[:, r],
-                        allgenomes[:, r],
-                        metric="l2",
-                        verbose=False,
-                        cutoff=0.02,
-                    )
-                    # get the maximum value of the new Exposure
-                    maxcoef = max(list(exposureAvg[:, r]))
-                    idxmaxcoef = list(exposureAvg[:, r]).index(maxcoef)
-
-                    exposureAvg[:, r] = np.round(exposureAvg[:, r])
-
-                    # We may need to tweak the maximum value of the new exposure to keep the total number of mutation equal to the original mutations in a genome
-                    if np.sum(exposureAvg[:, r]) != maxmutation:
-                        exposureAvg[:, r][idxmaxcoef] = (
-                            round(exposureAvg[:, r][idxmaxcoef])
-                            + maxmutation
-                            - sum(exposureAvg[:, r])
-                        )
-                    # print(exposureAvg[:, r])
-                    # print("\n")
+            fit_signatures = processAvg[:, init_decomposed_sigs_idx]
+            # fit signatures
+            newExposure, newSimilarity = ss.fit_signatures(
+                fit_signatures, allgenomes[:, r]
+            )
 
-                else:
-                    if verbose == True:
-                        print(
-                            "############################# Initial Composition #################################### "
-                        )
-                        print(pd.DataFrame(exposureAvg[:, r], index=allsigids).T)
-                        print("L2%: ", newSimilarity)
+            # create the exposureAvg vector
+            # print(init_decomposed_sigs_idx)
+            # print(newExposure)
+            for nonzero_idx, nozero_exp in zip(
+                init_decomposed_sigs_idx, newExposure
+            ):
+                exposureAvg[nonzero_idx, r] = nozero_exp
+            
+            if pcawg_rule == True:
+                maxmutation = np.sum(allgenomes[:, r])
+                (
+                    exposureAvg[:, r],
+                    remove_distance,
+                    _,
+                ) = ss.remove_all_single_signatures(
+                    processAvg,
+                    exposureAvg[:, r],
+                    allgenomes[:, r],
+                    metric="l2",
+                    verbose=False,
+                    cutoff=0.02,
+                )
+                # get the maximum value of the new Exposure
+                # maxcoef = max(list(exposureAvg[:, r]))
+                # idxmaxcoef = list(exposureAvg[:, r]).index(maxcoef)
+
+                # exposureAvg[:, r] = np.round(exposureAvg[:, r])
+                exposureAvg[:, r] = ss.roundConserveSum(exposureAvg[:, r])
+
+                # We may need to tweak the maximum value of the new exposure to keep the total number of mutation equal to the original mutations in a genome
+                # if np.sum(exposureAvg[:, r]) != maxmutation:
+                #     exposureAvg[:, r][idxmaxcoef] = (
+                #         round(exposureAvg[:, r][idxmaxcoef])
+                #         + maxmutation
+                #         - sum(exposureAvg[:, r])
+                #     )
+                # print(exposureAvg[:, r])
+                # print("\n")
+                if np.sum(exposureAvg[:, r])!=round(maxmutation): raise ValueError("Mutation count not conserved")
 
-                    lognote.write(
-                        "############################# Initial Composition ####################################\n"
+            else:
+                if verbose == True:
+                    print(
+                        "############################# Initial Composition #################################### "
                     )
-                    exposures = pd.DataFrame(exposureAvg[:, r], index=allsigids).T
-                    lognote.write(
-                        "{}\n".format(
-                            exposures.iloc[:, exposures.to_numpy().nonzero()[1]]
-                        )
+                    print(pd.DataFrame(exposureAvg[:, r], index=allsigids).T)
+                    print("L2%: ", newSimilarity)
+
+                lognote.write(
+                    "############################# Initial Composition ####################################\n"
+                )
+                exposures = pd.DataFrame(exposureAvg[:, r], index=allsigids).T
+                lognote.write(
+                    "{}\n".format(
+                        exposures.iloc[:, exposures.to_numpy().nonzero()[1]]
                     )
-                    lognote.write(
-                        "L2 Error %: {}\nCosine Similarity: {}\n".format(
-                            round(newSimilarity, 2),
-                            round(
-                                cos_sim(
-                                    allgenomes[:, r],
-                                    np.dot(processAvg, exposureAvg[:, r]),
-                                ),
-                                2,
+                )
+                lognote.write(
+                    "L2 Error %: {}\nCosine Similarity: {}\n".format(
+                        round(newSimilarity, 2),
+                        round(
+                            cos_sim(
+                                allgenomes[:, r],
+                                np.dot(processAvg, exposureAvg[:, r]),
                             ),
-                        )
+                            2,
+                        ),
+                    )
+                )
+                # remove signatures
+                
+                (
+                    exposureAvg[:, r],
+                    L2dist,
+                    cosine_sim,
+                ) = ss.remove_all_single_signatures(
+                    processAvg,
+                    exposureAvg[:, r],
+                    allgenomes[:, r],
+                    metric="l2",
+                    solver="nnls",
+                    cutoff=initial_remove_penalty,
+                    background_sigs=[],
+                    verbose=False,
+                )
+                if verbose == True:
+                    print(
+                        "############################## Composition After Initial Remove ############################### "
+                    )
+                    print(pd.DataFrame(exposureAvg[:, r], index=allsigids).T)
+                    print("L2%: ", L2dist)
+                lognote.write(
+                    "############################## Composition After Initial Remove ###############################\n"
+                )
+                exposures = pd.DataFrame(exposureAvg[:, r], index=allsigids).T
+                lognote.write(
+                    "{}\n".format(
+                        exposures.iloc[:, exposures.to_numpy().nonzero()[1]]
+                    )
+                )
+                lognote.write(
+                    "L2 Error %: {}\nCosine Similarity: {}\n".format(
+                        round(L2dist, 2), round(cosine_sim, 2)
                     )
-                    # remove signatures
+                )
+                lognote.write(
+                    "\n############################## Performing Add-Remove Step ##############################\n"
+                )
+                # Close the Lognote file
+                lognote.close()
+
+                init_add_sig_idx = list(
+                    set().union(
+                        list(np.nonzero(exposureAvg[:, r])[0]), background_sigs
+                    )
+                )
+                # print(init_add_sig_idx)
+
+                # get the background_sig_idx for the add_remove function only for the decomposed solution:
+                if background_sigs != 0:  # in the decomposed solution only
+                    background_sig_idx = get_indeces(allsigids, ["SBS1", "SBS5"])
+
+                # if the there is no other signatures to be added on top the existing signatures
+                try:
                     (
+                        _,
                         exposureAvg[:, r],
                         L2dist,
-                        cosine_sim,
-                    ) = ss.remove_all_single_signatures(
+                        similarity,
+                        kldiv,
+                        correlation,
+                        cosine_similarity_with_four_signatures,
+                    ) = ss.add_remove_signatures(
                         processAvg,
-                        exposureAvg[:, r],
                         allgenomes[:, r],
                         metric="l2",
                         solver="nnls",
-                        cutoff=initial_remove_penalty,
-                        background_sigs=[],
+                        background_sigs=init_add_sig_idx,
+                        permanent_sigs=background_sig_idx,
+                        candidate_sigs="all",
+                        allsigids=allsigids,
+                        add_penalty=add_penalty,
+                        remove_penalty=remove_penalty,
+                        check_rule_negatives=check_rule_negatives,
+                        checkrule_penalty=check_rule_penalty,
+                        connected_sigs=connected_sigs,
+                        directory=layer_directory
+                        + "/Solution_Stats/"
+                        + solution_prefix
+                        + "_Signature_Assignment_log.txt",
                         verbose=False,
                     )
+
                     if verbose == True:
                         print(
-                            "############################## Composition After Initial Remove ############################### "
+                            "####################################### Composition After Add-Remove #######################################\n"
                         )
-                        print(pd.DataFrame(exposureAvg[:, r], index=allsigids).T)
+                        print(exposureAvg[:, r])
                         print("L2%: ", L2dist)
+                    # Recond the information in the log file
+                    lognote = open(
+                        layer_directory
+                        + "/Solution_Stats/"
+                        + solution_prefix
+                        + "_Signature_Assignment_log.txt",
+                        "a",
+                    )
                     lognote.write(
-                        "############################## Composition After Initial Remove ###############################\n"
+                        "####################################### Composition After Add-Remove #######################################\n"
                     )
                     exposures = pd.DataFrame(exposureAvg[:, r], index=allsigids).T
                     lognote.write(
@@ -979,88 +1059,12 @@ def make_final_solution(
                     )
                     lognote.write(
                         "L2 Error %: {}\nCosine Similarity: {}\n".format(
-                            round(L2dist, 2), round(cosine_sim, 2)
+                            round(L2dist, 2), round(similarity, 2)
                         )
                     )
-                    lognote.write(
-                        "\n############################## Performing Add-Remove Step ##############################\n"
-                    )
-                    # Close the Lognote file
                     lognote.close()
-
-                    init_add_sig_idx = list(
-                        set().union(
-                            list(np.nonzero(exposureAvg[:, r])[0]), background_sigs
-                        )
-                    )
-                    # print(init_add_sig_idx)
-
-                    # get the background_sig_idx for the add_remove function only for the decomposed solution:
-                    if background_sigs != 0:  # in the decomposed solution only
-                        background_sig_idx = get_indeces(allsigids, ["SBS1", "SBS5"])
-
-                    # if the there is no other signatures to be added on top the existing signatures
-                    try:
-                        (
-                            _,
-                            exposureAvg[:, r],
-                            L2dist,
-                            similarity,
-                            kldiv,
-                            correlation,
-                            cosine_similarity_with_four_signatures,
-                        ) = ss.add_remove_signatures(
-                            processAvg,
-                            allgenomes[:, r],
-                            metric="l2",
-                            solver="nnls",
-                            background_sigs=init_add_sig_idx,
-                            permanent_sigs=background_sig_idx,
-                            candidate_sigs="all",
-                            allsigids=allsigids,
-                            add_penalty=add_penalty,
-                            remove_penalty=remove_penalty,
-                            check_rule_negatives=check_rule_negatives,
-                            checkrule_penalty=check_rule_penalty,
-                            connected_sigs=connected_sigs,
-                            directory=layer_directory
-                            + "/Solution_Stats/"
-                            + solution_prefix
-                            + "_Signature_Assignment_log.txt",
-                            verbose=False,
-                        )
-
-                        if verbose == True:
-                            print(
-                                "####################################### Composition After Add-Remove #######################################\n"
-                            )
-                            print(exposureAvg[:, r])
-                            print("L2%: ", L2dist)
-                        # Recond the information in the log file
-                        lognote = open(
-                            layer_directory
-                            + "/Solution_Stats/"
-                            + solution_prefix
-                            + "_Signature_Assignment_log.txt",
-                            "a",
-                        )
-                        lognote.write(
-                            "####################################### Composition After Add-Remove #######################################\n"
-                        )
-                        exposures = pd.DataFrame(exposureAvg[:, r], index=allsigids).T
-                        lognote.write(
-                            "{}\n".format(
-                                exposures.iloc[:, exposures.to_numpy().nonzero()[1]]
-                            )
-                        )
-                        lognote.write(
-                            "L2 Error %: {}\nCosine Similarity: {}\n".format(
-                                round(L2dist, 2), round(similarity, 2)
-                            )
-                        )
-                        lognote.close()
-                    except:
-                        pass
+                except:
+                    pass
 
     else:
         # when refilt de_novo_signatures

{sigprofilerassignment-0.2.1 → sigprofilerassignment-0.2.3}/SigProfilerAssignment/decomposition.py

@@ -25,7 +25,7 @@ from sigProfilerPlotting import sigProfilerPlotting as sigPlot
 import sigProfilerPlotting
 import os, sys
 from pypdf import PdfWriter, PdfReader
-import fitz
+from pdf2image import convert_from_path
 import time
 from pathlib import Path
 
@@ -60,21 +60,32 @@ def get_storage_dir(volume=None):
 
 
 def convert_PDF_to_PNG(input_file_name, output_directory, page_names):
-    pdf_doc = fitz.open(input_file_name)
-    zoom = 3
-    magnify = fitz.Matrix(zoom, zoom)
+    """
+    Converts each page of the PDF to a PNG image with names from page_names.
+    Requires the 'pdf2image' Python package and the 'poppler' binary.
+
+    Parameters:
+    - input_file_name (str): Path to the input PDF file.
+    - output_directory (str): Directory where PNGs will be saved.
+    - page_names (List[str]): List of names (without extensions) to name PNGs.
+    """
+    if not os.path.exists(output_directory):
+        os.makedirs(output_directory)
+
+    # Convert PDF pages to PIL images
+    try:
+        images = convert_from_path(input_file_name, dpi=300)
+    except Exception as e:
+        raise RuntimeError(f"Error converting PDF to images: {e}")
 
-    if pdf_doc.page_count != len(page_names):
+    if len(images) != len(page_names):
         raise ValueError(
             "Error: The number of samples and number of plots do not match."
         )
-    if not os.path.exists(output_directory):
-        os.makedirs(output_directory)
 
-    for sample_name, page in zip(page_names, pdf_doc):
-        pix = page.get_pixmap(matrix=magnify)
-        out_file_name = os.path.join(output_directory, sample_name + ".png")
-        pix.save(out_file_name)
+    for image, sample_name in zip(images, page_names):
+        output_path = os.path.join(output_directory, sample_name + ".png")
+        image.save(output_path, "PNG")
 
 
 # Create sample reconstruction plots
@@ -263,6 +274,7 @@ def spa_analyze(
         genome_build = A string. The reference genome build. List of supported genomes: "GRCh37", "GRCh38", "mm9", "mm10" and "rn6". The default value is "GRCh37". If the selected genome is not in the supported list, the default genome will be used.
         verbose = Boolean. Prints statements. Default value is False.
         exome = Boolean. Defines if the exome renormalized signatures will be used. The default value is False.
+        sample_reconstruction_plots (str): Select output format for sample reconstruction plots. Valid options are {'pdf', 'png', 'both', 'none'}. Default is 'none'.
 
     Values:
         The files below will be generated in the output folder.
@@ -951,8 +963,10 @@ def spa_analyze(
             cosmic_sig_ref = processAvg.copy(deep=True)
             cosmic_sig_ref.reset_index(inplace=True)
         else:
+            
             try:
                 processAvg = pd.read_csv(signature_database, sep="\t", index_col=0)
+                
             except:
                 sys.exit(
                     "Something is wrong with the format of signature database, Pass a text file of signatures in the format of COSMIC sig database"
@@ -1061,7 +1075,8 @@ def spa_analyze(
     recon_output_types = ["png", "pdf", "both"]
     # Generate sample reconstruction plots
     if (
-        sample_reconstruction_plots in recon_output_types
+        isinstance(sample_reconstruction_plots, str)
+        and sample_reconstruction_plots.lower() in recon_output_types
         and mutation_type == "96"
         and signature_database is None
     ):