Closed AishaAJ closed 6 years ago
Renaming my chromosomes to "1, 2, 3....X" fixed this.
Now all chromosome format names are supported.
I am trying to run circleseq.py (downloaded on June 2018), but still getting the same error as above. My reference hg38 and mapped sam/bam files use chr1, chr2, ...chrX, instead of 1, 2, ...X. Should I change all the chromosome names on bam/sam files? Please let me know how the chromosome format names are supported.
Thanks for developing the nice tool, by the way.
[07/19 02:03:24PM][INFO][alignReads] Sorting for control_1383D2_CS02_Alu by name complete. [07/19 02:03:24PM][INFO][circleseq] Finished merging and aligning reads. [07/19 02:03:24PM][INFO][circleseq] Identifying off-target cleavage sites. [07/19 02:03:24PM][INFO][circleseq] Window: 3, MAPQ: 50, Gap: 3, Start 1, Mismatches 6, Search_Radius 20 Writing counts to ./Experimental_data/180717_CIRCLE-seq-analysis/identified/1383D2_CS02_Alu_count.txt Tabulate nuclease merged start positions. 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0ATabulate control merged start positions. 0.1A0.2A0.3A0.4A0.5A0.6A0.7A0.8A0.9A1.0A1.1A1.2A1.3A1.4A1.5A1.6A1.7A1.8A1.9A2.0A2.1A2.2A2.3A2.4A2.5A2.6A2.7A2.8A2.9A3.0A3.1A3.2A3.3A3.4A3.5A3.663.76[]A[] [07/19 02:05:38PM][ERROR][circleseq] Error identifying off-target cleavage site. [07/19 02:05:38PM][ERROR][circleseq] Traceback (most recent call last): File "./circleseq/circleseq/circleseq.py", line 168, in findCleavageSites identified_sites_file, self.all_chromosomes, merged=self.merged_analysis) File "./circleseq/circleseq/findCleavageSites.py", line 611, in compare ecdf_pos = ECDF(bg_position) File "/anaconda2/lib/python2.7/site-packages/statsmodels/distributions/empirical_distribution.py", line 140, in init y = np.linspace(1./nobs,1,nobs) ZeroDivisionError: float division by zero
I'm getting this error.