gcContent {ORFik} | R Documentation |
0.5 means 50
gcContent(seqs, fa = NULL)
seqs |
a character vector of sequences, or ranges as GRangesList |
fa |
fasta index file .fai file, either path to it, or the loaded FaFile, default (NULL), only set if you give ranges as GRangesList |
a numeric vector of gc content scores
# Here we make an example from scratch seqName <- "Chromosome" ORF1 <- GRanges(seqnames = seqName, ranges = IRanges(c(1007, 1096), width = 60), strand = c("+", "+")) ORF2 <- GRanges(seqnames = seqName, ranges = IRanges(c(400, 100), width = 30), strand = c("-", "-")) ORFs <- GRangesList(tx1 = ORF1, tx2 = ORF2) # get path to FaFile for sequences faFile <- system.file("extdata/Danio_rerio_sample", "genome_dummy.fasta", package = "ORFik") gcContent(ORFs, faFile)