restoring read2sigma code after wrong merge