I downloaded and imported gnomad.exomes.r2.1.sites.vcf.bgz and saved it as MatrixTable:
mt = hl.read_matrix_table('/data/gnomad.exomes.mt')
Now, I’d like to filter vars with MAF < 10e-5 by non_neuro_AF_nfe and print out rsIDs and MAFs of those vars that pass that threshold.
I tried several R-minded options, but none of them really worked:
mtf = mt.filter_cols(mt.info.non_neuro_AF_nfe < 0.00001)
Can you help please?