diff --git a/README.md b/README.md
index 6792ee4..341a1ca 100644
--- a/README.md
+++ b/README.md
@@ -423,3 +423,26 @@ If the bam file has an MI tag, then the reads will be treated as linked reads.
The plots will be similar to short read plots, but all alignments with the same MI is plotted at the same height according to alignment with the largest gap in the group. A green line connects all alignments in a group.
+
+#### Interchromosomal gene fusion
+
+Interchromosomal variants require two location flag sets (-c, -s, and -e). For example, the BCR-ABL1 fusion from RNA-seq of the K562 chronic myeloid leukemia cell line is shown.
+
+```
+samplot plot \
+ -b test/data/k562.rna.starfusion.sort.region.bam \
+ -n K562.BCR--ABL1 \
+ -c chr22 \
+ -s 23632600 \
+ -e 133729451 \
+ -c chr9 \
+ -s 23632600 \
+ -e 133729451 \
+ -t BND \
+ -A test/data/bcr.abl1.bed.gz \
+ -o k562.bcr.abl1.png
+```
+
+
+
+
diff --git a/doc/imgs/k562.bcr.abl1.png b/doc/imgs/k562.bcr.abl1.png
new file mode 100644
index 0000000..48fdd54
Binary files /dev/null and b/doc/imgs/k562.bcr.abl1.png differ
diff --git a/test/data/bcr.abl1.bed.gz b/test/data/bcr.abl1.bed.gz
new file mode 100644
index 0000000..c09020c
Binary files /dev/null and b/test/data/bcr.abl1.bed.gz differ
diff --git a/test/data/bcr.abl1.bed.gz.tbi b/test/data/bcr.abl1.bed.gz.tbi
new file mode 100644
index 0000000..e2c2606
Binary files /dev/null and b/test/data/bcr.abl1.bed.gz.tbi differ
diff --git a/test/data/k562.rna.starfusion.sort.region.bam b/test/data/k562.rna.starfusion.sort.region.bam
new file mode 100644
index 0000000..d6f72bd
Binary files /dev/null and b/test/data/k562.rna.starfusion.sort.region.bam differ
diff --git a/test/data/k562.rna.starfusion.sort.region.bam.bai b/test/data/k562.rna.starfusion.sort.region.bam.bai
new file mode 100644
index 0000000..6d8f610
Binary files /dev/null and b/test/data/k562.rna.starfusion.sort.region.bam.bai differ