search

jsh58 / Genrich

Detecting sites of genomic enrichment
MIT License
182 stars 27 forks source link

Follow Up: Genrich with reads from a tiny subset of genome #29 #37

Closed vanottee closed 4 years ago

vanottee commented 4 years ago

Dear Genrich team -

I am wondering what chromosomes I am missing to exclude the majority of the genome from my analysis? Similar to question 29, I am trying to run my analysis on only a subset of the genome (~75Kb). However, because of my bam header, I am also getting misleading data due to the size of the genome being considered during analysis. I have used both the -e and -E functions to try and eliminate every possible region of the genome (except my 75Kb) and it still shows the genome length as ~1.6x10^8 in length. I have attached a screenshot that includes the regions I am eliminating from analysis. Capture

jsh58 commented 4 years ago

At a first glance, I would guess chrm3 is the issue.

vanottee commented 4 years ago

Obvious, oops, sorry for the bother