Align EffectAllele to always be the reference genome allele
align_to_ref.Rdalign_to_ref() will:
Filter any variants where EffectAllele or OtherAllele is not the reference genome allele
Flip EffectAllele such that it is always the reference genome allele
If OtherAllele is the reference genome allele, direction of B,Z and EAF are flipped.
Usage
align_to_ref(dset, ref = c("REF_38", "REF_37"))Arguments
- dset
object created by
arrow::open_dataset()ordplyr::tibble()- ref
Reference genome allele to align to. Varies in ~0.5% of locations
Value
a dplyr::tibble() or arrow query depending on whether dset is a tibble or arrow dataset