Cool Bioinformatics Scripts With Python

You can use this script in two ways

• read tons of millions of P values from stdin

# python 
zcat pval.txt.gz | qqplot.py -out test -title "QQ plot on the fly"
# julia
zcat pval.txt.gz | qqplot.jl --out test --title "QQ plot on the fly"

warning : If you have 100 billion P values to process you should definitely use qqplot.jl instead of qqplot.py. The hourly processed lines of julia version is 3 billion while python is only 700 million on my server.

• use qqplot.py in your script