This is a test version of Biostars. For the public version, visit https://www.biostars.org.
SICER2 (or EPIC) genomes

Hi all,

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?

genome epic chip-seq sicer sequencing

1 answer

See https://github.com/biocore-ntnu/epic2

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.

Log in to answer this question.