brwnj

"""
qvality must be on the $PATH.
run on a file with pvalues in column 4:

    $ python qvality.py --column 4 input.bed > output.bed

output file will have extra columns for PEP (posterior error prob)
and q-value.

if input.bed has a header, it is retained.