Мария Нечеухина
unread,Apr 7, 2025, 7:11:36 PMApr 7Sign in to reply to author
Sign in to forward
You do not have permission to delete messages in this group
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to 3D Genomics
I use juicer version 2.0 and juicer_tools.version 2.09.00. I get a translation of the resulting hexagon into a cool hexagon, but I don't find the links: INFO:hicexplorer.hicConvertFormat:Conversion using hic2cool.
Backtracking (last call):
File "/home/brmsh/.conda/envs/myenv/bin/hicConvertFormat", line 7, in <module>
main()
File "/home/brmsh/.conda/envs/myenv/lib/python3.9/site-packages/hicexplorer/hicConvertFormat.py ", line 124, in main
hic2cool_convert(matrix, outgoing name, resolution)
File "/home/brmsh/.conda/envs/myenv/lib/python3.9/site-packages/hic2cool/hic2cool_utils.py ", line 859, in hic2cool_convert
used_chrs, permissions, masteridx, genome, metadata = read_header(req)
File "/home/brmsh/.conda/envs/myenv/lib/python3.9/site-packages/hic2cool/hic2cool_utils.py ", line 95, in read_header
key = readcstr(req)
File "/home/brmsh/.conda/envs/myenv/lib/python3.9/site-packages/hic2cool/hic2cool_utils.py ", line 58, in readcstr
returns buf.decode("utf-8")
UnicodeDecodeError: The 'utf-8' codec cannot decode bytes at position 0-1: invalid continuation byte
the value of max_perc_undetected is set to 50.0 based on the configuration file.
The note is made like this: enter the value "$(/tools/bin/conda3.sh )"
to activate juicy_private
/home/brmsh/tsv/juicer/scripts/juicer.sh -g dm6 -d $dir \
-p /store/bioinf/analysis/chipseq/proc_data/dm6/chr_len_sorted \
-y "/store/bioinf/analysis/chipseq/proc_data/dm6/genome_GTATAC_TTATAA_GTATAA_TTATAC.txt" \
-z /store/bioinf/analysis/chipseq/proc_data/dm6/genome.fa \
-t 40 \
--qc
conda is deactivated
cd $dir/aligned
conda activates myenv
hic convertformat -m inter_30.hic --InputFormat hic --OutputFormat cool -o inter_30.cool --resolutions 2000
conda deactivates
conda activates new
chromosight detect - number of streams 48 -Pearson=0.3 -minimum distance=8000 -maximum distance=800000 -minimum distance =2000 -zero value=10 inter_30_1000.cool inter_30_1000
echo "Passed before 2000"
conda deactivate