Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

how to use custom genome i.e. E. Coli in epic2 ? #59

Open
Ankana001 opened this issue May 10, 2022 · 1 comment
Open

how to use custom genome i.e. E. Coli in epic2 ? #59

Ankana001 opened this issue May 10, 2022 · 1 comment

Comments

@Ankana001
Copy link

Hi, I am using epic2 for my E.coli ChIP-seq data, But it runs with hg19 genome as default, so I used custom genome fasta file which is in the same folder as the data. I am using this following command:
epic2 -t SO_4933_A_CHR1_R1_sorted.bam -c SO_4933_A_INR1_R1_sorted.bam --genome E.coli_BW25113.fasta --false-discovery-rate-cutoff 0.001 -a -o SO_4933_A_enriched_regions.csv

But it's showing me error Chromosomes removed from analysis: NZ_CP064677.1 (and 0 more.)

Traceback (most recent call last): File "/home/ibab/.local/bin/epic2", line 257, in _main(args)
File "/home/ibab/.local/lib/python3.8/site-packages/epic2/main.py", line 35, in _main
effective_genome_length, chromsizes = egl_and_chromsizes(args)
File "epic2/src/genome_info.pyx", line 325, in epic2.src.genome_info.egl_and_chromsizes
File "epic2/src/genome_info.pyx", line 267, in epic2.src.genome_info.get_effective_genome_fraction
KeyError: 'e.coli_bw25113_75.txt'

how to run this analysis with epic2 for E.coli genome?

@endrebak
Copy link
Member

endrebak commented Jun 9, 2022

You need to use the chromsizes argument where the argument is a file that looks like this:
https://github.com/biocore-ntnu/epic2/blob/master/epic2/chromsizes/balAcu1.chromsizes

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants