DNAsty scheme

1 Selection of invariant genes / expression interspersed in lowDim space (scanpy)

2 BEDfile from selected genes

3 Remove SNP from dbSNP above f allelic frequency (pybedtools)

4 split gene into point positions

5 feed the lociList to SCanSNP main walker

selection biallelic only?? :S

6 ADATA layer addition (with locus-level information)

adata.layer["Alt"]

adata.layer["Ref"]