This is to limit clutter from old analyses.

Duplicates were due to `protein_coding_*` biotypes. Those were removed
in join_all_counts.

Also introduced minor modification recently made for other pipelines.

Currently only the folds computed by DESeq2 are available for RNA-seq.

Made paired scatterplot plotting function part of libhts as
`plot_paired_scatters`. Modifications in small RNA-seq pipeline not
tested.

The genome is now configurable. This lead to large modifications in the
config file.

Output organization simplified.

Factorized linking to raw data.

It turns out that using `include:` can fail when there are symlinks in
the path to the snakefile.

This makes the DE results easier to compare with those from small
RNA-seq.

For RNA-seq and GRO-seq, log and data directories are now subdirectories
of output directories. local_annot_dir also moves as subdirectory of
GRO-seq output.

This seems to happen on data matrices with only 2 rows left after
filtering.

The goal is to be able to factorize the rule in a wrapper.

Started to apply to RNA-seq and GRO-seq technique used in small RNA-seq.

I should check to be sure that this is the correct thing to do. The
other possibility would be to sum using agg:
https://stackoverflow.com/a/35403903/1878788

lfcMLE may bring problems with low-expressed genes.

Also:
- Boxplots with labels turned 90 degrees
- do_deseq2 returns size factors

It currently fails when libraries have too many rRNA reads:
https://github.com/arq5x/bedtools2/issues/565

Also added lfc distribution and improved MA plot.