diff --git a/app.py b/app.py
index 3942ee5..8695676 100644
--- a/app.py
+++ b/app.py
@@ -6,8 +6,7 @@
 import pandas as pd
 import streamlit as st
 from streamlit_extras.stylable_container import stylable_container
-## package
-from bridge_rna_designer.run import design_bridge_rna
+from bridgernadesigner.run import design_bridge_rna
 
 # Functions
 def get_image_as_base64(path):
diff --git a/bridge_rna_designer/__init__.py b/bridge_rna_designer/__init__.py
deleted file mode 100644
index e69de29..0000000
diff --git a/bridge_rna_designer/classes.py b/bridge_rna_designer/classes.py
deleted file mode 100644
index cef6dba..0000000
--- a/bridge_rna_designer/classes.py
+++ /dev/null
@@ -1,190 +0,0 @@
-from Bio.Seq import reverse_complement
-from bridge_rna_designer import errors
-import re
-
-class WTBridgeRNA177nt:
-    """
-    A class to represent a 177 nt bridge RNA.
-    """
-
-    # Template sequence
-    TEMPLATE = "AGTGCAGAGAAAATCGGCCAGTTTTCTCTGCCTGCAGTCCGCATGCCGTNNNNNNNNNTGGGTTCTAACCTGTNNNNNNNNNTTATGCAGCGGACTGCCTTTCTCCCAAAGTGATAAACCGGNNNNNNNNATGGACCGGTTTTCCCGGTAATCCGTNNTTNNNNNNNTGGTTTCACT"
-    DOT_STRUC_P1 = "...(((((((((((......))))))))))).((((((((((.(((............(((((....)))))..............))))).)))))))).........((((....(((.............((((((.....)))))...)...............)))..))))"
-    DOT_STRUC_P2 = "((.(((((((((((......)))))))))))))(((((((((.(((.............<(((.<>.)))>...............))))))))))))...........((((...<(((..........<.(((((((.....)))))...))....>.........)))>.))))"
-    GUIDES = ".................................................LLLLLLLCC...............RRRRRCCHH........................................lllllllc..........................rr..rrrcchh.........."
-
-    # Sequence regions
-    LTG_REGION = (49, 58)
-    RTG_REGION = (73, 80)
-    LDG_REGION = (122, 130)
-    RDG1_REGION = (160, 165)
-    RDG2_REGION = (156, 158)
-    TBL_HSG_REGION = (80, 82)
-    DBL_HSG_REGION = (165, 167)
-
-    def __init__(self):
-        self.bridge_sequence = str(self.TEMPLATE)
-        self.target = None
-        self.donor = None
-
-    def update_target(self, target: str):
-        """
-        Update the target sequence.
-        Args:
-            target: The target sequence (14 bp).
-        """
-        # Store target
-        self.target = target
-
-        # Check core of target
-        core = target[7:9]
-        if core not in ["CT", "GT"]:
-            errors.NotCTorGTCoreWarning()
-        if core not in ["CT", "GT", "AT", "TT"]:
-            errors.NotNTCoreWarning()
-
-        # Update bridge sequence
-        bridgeseq = list(self.bridge_sequence)
-        bridgeseq[self.LTG_REGION[0]:self.LTG_REGION[1]] = list(target[:9])
-        bridgeseq[self.RTG_REGION[0]:self.RTG_REGION[1]] = list(reverse_complement(target[7:]))
-        self.bridge_sequence = "".join(bridgeseq)
-
-    def update_donor(self, donor: str):
-        """
-        Update the donor sequence.
-        Args:
-            donor: The donor sequence (14 bp).
-        """
-        # Store donor
-        self.donor = donor
-
-        # Check input
-        if self.has_donor_p7C():
-            errors.DonorP7CWarning()
-
-        # Update bridge sequence
-        bridgeseq = list(self.bridge_sequence)
-        bridgeseq[self.LDG_REGION[0]:self.LDG_REGION[1]] = list(donor[:8])
-        bridgeseq[self.RDG1_REGION[0]:self.RDG1_REGION[1]] = list(reverse_complement(donor[7:-2]))
-        bridgeseq[self.RDG2_REGION[0]:self.RDG2_REGION[1]] = list(reverse_complement(donor[-2:]))
-        self.bridge_sequence = "".join(bridgeseq)
-
-    def update_hsg(self):
-        """
-        Update the HSG region of the bridge RNA.
-        """
-        # Update bridge sequence
-        bridgeseq = list(self.bridge_sequence)
-        target_p6p7, donor_p6p7 = self.get_p6p7()
-
-        # Ff P6-P7 do not match, reverse complement the donor P6-P7 and target P6-P7
-        if target_p6p7 != donor_p6p7:
-            bridgeseq[self.TBL_HSG_REGION[0]:self.TBL_HSG_REGION[1]] = list(reverse_complement(donor_p6p7))
-            bridgeseq[self.DBL_HSG_REGION[0]:self.DBL_HSG_REGION[1]] = list(reverse_complement(target_p6p7))
-        else:
-            errors.MatchingP6P7Warning()
-            bridgeseq[self.TBL_HSG_REGION[0]:self.TBL_HSG_REGION[1]] = list(reverse_complement(donor_p6p7))
-            bridgeseq[self.DBL_HSG_REGION[0]:self.DBL_HSG_REGION[1]] = list(target_p6p7)
-
-        self.bridge_sequence = "".join(bridgeseq)
-
-    @staticmethod
-    def check_target_length(target):
-        if len(target) != 14:
-            raise errors.TargetLengthError()
-
-    @staticmethod
-    def check_donor_length(donor):
-        if len(donor) != 14:
-            raise errors.DonorLengthError()
-
-    @staticmethod
-    def check_core_match(target, donor):
-        if target[7:9] != donor[7:9]:
-            raise errors.CoreMismatchError()
-
-    @staticmethod
-    def check_donor_is_dna(donor):
-        if re.match("^[ATCG]*$", donor) is None:
-            raise errors.DonorNotDNAError()
-
-    # @staticmethod
-    def check_target_is_dna(target):
-        if re.match("^[ATCG]*$", target) is None:
-            raise errors.TargetNotDNAError()
-
-    def get_p6p7(self):
-        return self.target[5:7], self.donor[5:7]
-
-    def has_donor_p7C(self):
-        if self.donor[6] == "C":
-            return True
-        else:
-            return False
-
-    def has_matching_p6p7(self):
-        if self.target[5:7] == self.donor[5:7]:
-            return True
-        else:
-            return False
-
-    def format_fasta(self, line_wrap: int=80) -> str:
-        """
-        Format the bridge RNA as a FASTA sequence.
-        Args:
-            line_wrap: The number of characters to wrap the sequence at. If set to 0, the sequence will not be wrapped.
-        Returns:
-            str: The formatted FASTA sequence.
-        """
-        out = ">BridgeRNA_tgt_{}_dnr_{}\n".format(self.target, self.donor)
-        if line_wrap > 1:
-            for i in range(0, len(self.bridge_sequence), line_wrap):
-                out += self.bridge_sequence[i:i+line_wrap] + "\n"
-        else:
-            out += self.bridge_sequence
-        return out
-
-    def format_stockholm(self, whitespaces: int=5) -> str:
-        """
-        Format the bridge RNA as a Stockholm sequence.
-        Args:
-            whitespaces: The number of whitespaces to use for leader columns.
-        Returns:
-            str: The formatted Stockholm sequence.
-        """
-        seqname = "BridgeRNA_tgt_{}_dnr_{}".format(self.target, self.donor)
-        leader_cols = len(seqname)+whitespaces
-        out = "# STOCKHOLM 1.0\n"
-        out += seqname+" "*(leader_cols-len(seqname))+self.bridge_sequence+"\n"
-        template_feat = "#=GC bRNA_template"
-        out += template_feat + " "*((leader_cols-len(template_feat))) + self.TEMPLATE + "\n"
-        guide_feat = "#=GC guides"
-        out += guide_feat + " "*((leader_cols-len(guide_feat))) + self.GUIDES + "\n"
-        ss_feat = "#=GC SS"
-        out += ss_feat + " " * ((leader_cols - len(ss_feat))) + self.DOT_STRUC_P2 + "\n"
-
-        donor_p7 = self.donor[6]
-        target_p6p7 = self.target[5:7]
-        donor_p6p7 = self.donor[5:7]
-        core = self.target[7:9]
-
-
-        warning_feat = "#=GF WARNING"
-        if donor_p7 == 'C':
-            out += warning_feat + " Donor P7 is C, this was found to be very inefficient in a screen.\n"
-        if target_p6p7 == donor_p6p7:
-            out += warning_feat + " Target P6-P7 and Donor P6-P7 match, efficiency is unclear.\n"
-        if core != 'CT' and core != 'GT':
-            out += warning_feat + " Core is not CT or GT, efficiency is unclear.\n"
-        if core != 'CT' and core != 'GT' and core != 'AT' and core != 'TT':
-            out += warning_feat + " Core does not follow the expected NT format, likely inefficient.\n"
-        out += "//"
-        return out
-
-    # Class custom print function 
-    def __str__(self):
-        out = "<WTBridgeRNA177nt>\n"
-        out += '- Programmed target: {}\n'.format(self.target)
-        out += '- Programmed donor: {}\n'.format(self.donor)
-        out += '- Bridge RNA sequence: {}'.format(self.bridge_sequence)
-        return out
\ No newline at end of file
diff --git a/bridge_rna_designer/errors.py b/bridge_rna_designer/errors.py
deleted file mode 100644
index 3373e7f..0000000
--- a/bridge_rna_designer/errors.py
+++ /dev/null
@@ -1,38 +0,0 @@
-import warnings
-
-class TargetNotDNAError(Exception):
-    def __init__(self, message="The target sequence must be DNA."):
-        self.message = message
-        super().__init__(self.message)
-
-class DonorNotDNAError(Exception):
-    def __init__(self, message="The donor sequence must be DNA."):
-        self.message = message
-        super().__init__(self.message)
-
-class TargetLengthError(Exception):
-    def __init__(self, message="The target sequence must be exactly 14 nucleotides long."):
-        self.message = message
-        super().__init__(self.message)
-
-class DonorLengthError(Exception):
-    def __init__(self, message="The donor sequence must be exactly 14 nucleotides long."):
-        self.message = message
-        super().__init__(self.message)
-
-class CoreMismatchError(Exception):
-    def __init__(self, message="That target and donor cores must be the same."):
-        self.message = message
-        super().__init__(self.message)
-
-def DonorP7CWarning():
-    warnings.warn("Donor P7 is C, this was found to be very inefficient in a screen.")
-
-def MatchingP6P7Warning():
-    warnings.warn("Target P6-P7 and Donor P6-P7 match, efficiency is unclear.")
-
-def NotCTorGTCoreWarning():
-    warnings.warn("Core is not CT or GT, efficiency is unclear.")
-
-def NotNTCoreWarning():
-    warnings.warn("Core does not follow the expected NT format, likely inefficient.")
\ No newline at end of file
diff --git a/bridge_rna_designer/main.py b/bridge_rna_designer/main.py
deleted file mode 100644
index f267441..0000000
--- a/bridge_rna_designer/main.py
+++ /dev/null
@@ -1,65 +0,0 @@
-import click
-import sys
-from bridge_rna_designer import errors, run
-from bridge_rna_designer.classes import WTBridgeRNA177nt
-
-@click.group(invoke_without_command=True)
-@click.option('--target', '-t', required=True)
-@click.option('--donor', '-d', required=True)
-@click.option('--output-format', '-of', default='stockholm', type=click.Choice(['stockholm', 'fasta']))
-@click.pass_context
-def cli(ctx, target, donor, output_format):
-    if ctx.invoked_subcommand is None:
-        ctx.invoke(default_command, target=target, donor=donor, output_format=output_format)
-
-@cli.command()
-@click.pass_context
-def default_command(ctx, target, donor, output_format):
-
-    target = target.upper()
-    donor = donor.upper()
-
-    try:
-        WTBridgeRNA177nt.check_target_length(target)
-    except errors.TargetLengthError as e:
-        print("ERROR:", e)
-        print("Exiting...")
-        sys.exit()
-
-    try:
-        WTBridgeRNA177nt.check_donor_length(donor)
-    except errors.DonorLengthError as e:
-        print("ERROR:", e)
-        print("Exiting...")
-        sys.exit()
-
-    try:
-        WTBridgeRNA177nt.check_core_match(target, donor)
-    except errors.CoreMismatchError as e:
-        print("ERROR:", e)
-        print("Exiting...")
-        sys.exit()
-
-    try:
-        WTBridgeRNA177nt.check_target_is_dna(target)
-    except errors.TargetNotDNAError as e:
-        print("ERROR:", e)
-        print("Exiting...")
-        sys.exit()
-
-    try:
-        WTBridgeRNA177nt.check_donor_is_dna(donor)
-    except errors.DonorNotDNAError as e:
-        print("ERROR:", e)
-        print("Exiting...")
-        sys.exit()
-
-    brna = run.design_bridge_rna(target, donor)
-    if output_format == "fasta":
-        print(brna.format_fasta())
-    elif output_format == "stockholm":
-        print(brna.format_stockholm())
-
-
-if __name__ == '__main__':
-    cli()
diff --git a/bridge_rna_designer/run.py b/bridge_rna_designer/run.py
deleted file mode 100644
index 5550ac5..0000000
--- a/bridge_rna_designer/run.py
+++ /dev/null
@@ -1,31 +0,0 @@
-from bridge_rna_designer.classes import WTBridgeRNA177nt
-from bridge_rna_designer import errors
-
-def design_bridge_rna(target: str, donor: str) -> WTBridgeRNA177nt:
-    """
-    Design a bridge RNA from a target and donor sequence.
-    Args:
-        target: The target sequence (14 bp).
-        donor: The donor sequence (14 bp).
-    Returns:
-        WTBridgeRNA177nt: A bridge RNA object.
-    """
-    # Convert input sequences to uppercase
-    target = target.upper()
-    donor = donor.upper()
-
-    # Check input
-    WTBridgeRNA177nt.check_target_length(target)
-    WTBridgeRNA177nt.check_donor_length(donor)
-    WTBridgeRNA177nt.check_core_match(target, donor)
-    WTBridgeRNA177nt.check_target_is_dna(target)
-    WTBridgeRNA177nt.check_donor_is_dna(donor)
-
-    # Design bridge RNA
-    brna = WTBridgeRNA177nt()
-    brna.update_target(target)
-    brna.update_donor(donor)
-    brna.update_hsg()
-
-    return brna
-
diff --git a/requirements.txt b/requirements.txt
index c86fab4..15818a7 100644
--- a/requirements.txt
+++ b/requirements.txt
@@ -4,4 +4,5 @@ python-dotenv==1.0.1
 streamlit==1.32.2
 tabulate==0.9.0
 streamlit-extras==0.4.0
+bridgernadesigner==0.0.1
 # pyperclip==1.8.2
\ No newline at end of file