FeatureAnnotationSetDev/illuminaHumanWG-6_v3.R

# Create FAS called "HumanWG-6_v3" with note "made from Illumina manifest 2018".

manifest <- data.table::fread("FeatureAnnotationSetDev/HumanWG-6_V3_0_R3_11282955_A.txt",
                              skip = 8, fill = T)
manifest <- manifest[ !is.na(manifest$Entrez_Gene_ID), ]
manifest <- manifest[ , colnames(manifest) %in% c("Symbol", "Probe_Id"), with=F]
colnames(manifest) <- c("Gene_Symbol", "Probe_Id")
write.table(manifest,
            file = "FeatureAnnotationSetDev/IlluminaHumanWG-6_v3.tsv",
            quote = FALSE,
            sep = "\t",
            row.names = FALSE)
ehfhcrc/UpdateAnno documentation built on May 16, 2019, 12:16 a.m.