Use vcf_predict.py strand convention for RC averaging

README.md

@@ -40,6 +40,11 @@ model directly: they select the ensemble that did not train on your query's
 chromosome, and they add the MPRA vector context and average over both strands.
 Skipping either step returns plausible-looking but wrong numbers instead of an error.
 
+`predict` follows `vcf_predict.py` from the upstream code base, which generated the
+published predictions: the reverse strand is the reverse complement of the 200 bp
+insert placed back in the forward-orientation vector, matching the assay, rather
+than a reverse complement of the whole 600 bp construct.
+
 MPAC covers autosomes only; `from_pretrained` raises on chrX, chrY and anything else
 with no held-out fold.