def printBond(bond, shift, molecule, alchemicalTransformation):
    """Generate bond line

    Parameters
    ----------
    bond : Bond Object
        Bond Object
    shift : int
        Shift produced by structural dummy atoms
    molecule : molecule object
        Molecule object
    alchemicalTransformation : bool
        True if alchemical transformation

    Returns
    -------
    bondLine : str
        Bond line data
    """


    k0 = bond.K0*836.80
    R0 = bond.R0*0.1
    ftype = 1

    line = ''

    atomAbond = molecule.atoms[bond.atomA.serialOriginal -1].serial-shift
    atomBbond = molecule.atoms[bond.atomB.serialOriginal -1].serial-shift


    if alchemicalTransformation:

        k0_B = bond.K0_B*836.80
        R0_B = bond.R0_B*0.1

        line = '%5d%5d%5d %11.4f %13.3f%12.4f%11.3f\n' % (atomAbond, atomBbond, ftype, R0, k0, R0_B, k0_B)
        # line = '%5d%5d%5d %11.4f %13.3f%12.4f%11.3f\n' % (bond.atomA.serial-shift, bond.atomB.serial-shift, ftype, R0, k0, R0_B, k0_B)

    else: line = '%5d%5d%5d %11.4f %13.3f\n' % (atomAbond, atomBbond, ftype, R0, k0)

    return line