import pandas as pd rpk = pd.read_table(snakemake.input.rpk_file, index_col="gene") tpm = 1000000 * rpk / rpk.sum() tpm.to_csv(snakemake.output.tpm_file, sep="\t")