diff --git a/README.md b/README.md
index 430300e4154ba5b1b1e1f5eba25e8a3e8d09c402..b941e9a83abc03fefd90a315858be6f9008f66df 100644
--- a/README.md
+++ b/README.md
@@ -12,7 +12,11 @@ later be used for identifying crossovers through finding haplotype shifts (see [
 It takes the large bam file which contains all aligned DNA reads from gamete cells and 
 summarizes allele counts for the provided informative SNP markers. While counting
 the alleles, the Viterbi algorithm is executed for finding the haplotype 
-sequence for the list of SNP markers.
+sequence for the list of SNP markers. 
+
+### Publication
+
+![BioRxiv paper](https://www.biorxiv.org/content/10.1101/2022.02.10.479822v1)
 
 ![sgcocaller_fig](images/sgcocaller_fig.png)