Skip to content

Instantly share code, notes, and snippets.

@tbrittoborges
Created January 10, 2025 14:42
Show Gist options
  • Save tbrittoborges/9033a7f80145077755f0a1f3f1896042 to your computer and use it in GitHub Desktop.
Save tbrittoborges/9033a7f80145077755f0a1f3f1896042 to your computer and use it in GitHub Desktop.
converting seqname among different database styles
if (!require("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("VariantAnnotation")
library("VariantAnnotation")
vcf <- readVcf("/home/aschuhegger/genomics/mills_and_1000G.indels.vcf.gz")
vcf@rowRanges
seqlevelsStyle(vcf@rowRanges) <- 'Ensembl'
vcf@rowRanges
@tbrittoborges
Copy link
Author

tbrittoborges commented Jan 10, 2025

R is available with module load R

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment