get_merged_fasta: merge footprints and return corresponding sequence-

In jiang-junyao/IReNA: IReNA

merge footprints and return corresponding sequence-

Description

This function first merges footprints whose distances is less than 4, and then obtains the sequence of each footprint based on the reference genome, and output the sequence in FASTA format

Usage

get_merged_fasta(footprints, fastadir, distance = 4)

Arguments

footprints	footprints file, generated by dnase2tf or Hint or other related software, name of first column should be 'chr', name of second column should be 'start' and name of thrid column should be 'end'.
fastadir	character, indiating the path of reference genome
distance	numeric, indicating the cutoff for merging footprints

Value

return sequence of footprints in fasta format

Examples

fastadir='Genome/GRCm38Chr.fasta'
#merged_fasta <- get_merged_fasta(fdr005,fastadir)

jiang-junyao/IReNA documentation built on Nov. 15, 2024, 1:51 a.m.