Hi,
I've been troubleshooting getting juicer to run on our LSF cluster using the example HIC003 data.
I get the following error when trying to create the inter.hic file running:
/juicer/scripts/juicebox pre -s /juicer/work/HIC003/aligned/inter.txt -g /juicer/work/HIC003/aligned/inter_hists.m
-q 1 /juicer/work/HIC003/aligned/merged_nodups.txt /juicer/work/HIC003/aligned/inter.hic /juicer/references/chrom.sizes;
Error:
Picked up _JAVA_OPTIONS: -Xmx16384m
Skipping chr1 249250621
Skipping chr2 243199373
Skipping chr3 198022430
Skipping chr4 191154276
Skipping chr5 180915260
Skipping chr6 171115067
Skipping chr7 159138663
Skipping chr8 146364022
Skipping chr9 141213431
Skipping chr10 135534747
Skipping chr11 135006516
Skipping chr12 133851895
Skipping chr13 115169878
Skipping chr14 107349540
Skipping chr15 102531392
Skipping chr16 90354753
Skipping chr17 81195210
Skipping chr18 78077248
Skipping chr19 59128983
Skipping chr20 63025520
Skipping chr21 48129895
Skipping chr22 51304566
Skipping chrX 155270560
Skipping chrY 59373566
Skipping chrM 16572
Not including fragment map
Start preprocess
Writing header
Writing body
java.lang.RuntimeException: No reads in Hi-C contact matrices. This could be because the MAPQ filter is set too high (-q) or because all reads map to the same fragment.
at juicebox.tools.utils.original.Preprocessor$MatrixZoomDataPP.mergeAndWriteBlocks(Preprocessor.java:1397)
at juicebox.tools.utils.original.Preprocessor$MatrixZoomDataPP.access$000(Preprocessor.java:1168)
at juicebox.tools.utils.original.Preprocessor.writeMatrix(Preprocessor.java:582)
at juicebox.tools.utils.original.Preprocessor.writeBody(Preprocessor.java:313)
at juicebox.tools.utils.original.Preprocessor.preprocess(Preprocessor.java:223)
at juicebox.tools.clt.old.PreProcessing.run(PreProcessing.java:98)
at juicebox.tools.HiCTools.main(HiCTools.java:77)
I noticed my merged_nodups.txt does not have the 11 column format described in the doc for juicebox pre <infile>
Example line of merged_nodups.txt:
0 chr1 10005 0 0 chr1 249240258 64394 0 100M CCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTA 0 46M1I21M32S GGTTAGGGTTAGGGTTAGGGTTAGGGTTAGGGTTAGGGTTAGGGTTGGGGTTAGGGGTAGGGGTAGGGGGTGGGGTAGGGATAGGGATAGGGGTAGGGGT M00336:181:000000000-A29H6:1:1110:11312:16525/1 M00336:181:000000000-A29H6:1:1110:11312:16525/2
Is that close to the problem?
Thanks,
Tyler