This is probably a very basic question, but I am trying to use SICER (or EPIC) to analyze my ChIP-seq data. The samples are Arabidopsis, but TAIR10 is not an option for the genome. Does anyone know how I can use a genome not listed?
You either specify a hardcoded genome OR use --chromsizes & --effective-genome-fraction to indicate chromosome length and the effective (=mappable / unqiue) fraction.
https://deeptools.readthedocs.io/en/develop/content/feature/effectiveGenomeSize.html says TAIR10 has an effective size of 119481543 so I'd guess that you could use simply --chromsizes 119481543 --effective-genome-fraction 1 as a simple workaround.
--chromsizes 119481543 --effective-genome-fraction 1
Login before adding your answer.
Use of this site constitutes acceptance of our User Agreement and Privacy