Blacklist
Blacklist copied to clipboard
Entire chromosome outputs as high signal
Entire chr29 as High Signal Region
I have created mappability file for chr29 (smallest chromosome for quicker analysis) with kmer 100 and 150, and used final uint output file of combine_umaps
However, it spits out entire chromosome as high signal region, not sure if this is because of low number of inputs I have provided (5 inputs).
Blacklist was installed using https://anaconda.org/bioconda/encode-blacklist
$ Blacklist chr29
chr29 0 43051100 High Signal Region
Directory structure below
$ tree
.
├── input
│ ├── Ss1_final_chr29.bam
│ ├── Ss1_final_chr29.bam.bai
│ ├── Ss2_final_chr29.bam
│ ├── Ss2_final_chr29.bam.bai
│ ├── Ss3_final_chr29.bam
│ ├── Ss3_final_chr29.bam.bai
│ ├── Ss4_final_chr29.bam
│ ├── Ss4_final_chr29.bam.bai
│ ├── Ss5_final_chr29.bam
│ └── Ss5_final_chr29.bam.bai
└── mappability
└── chr29.uint8.unique
Contents of uint file # combined kmer 100 and 150
$ od -t x1 mappability/chr29.uint8.unique | head
0000000 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00
*
0000140 00 00 00 00 00 00 00 96 96 96 96 96 96 96 96 96
0000160 96 96 96 96 96 96 96 96 96 96 96 96 96 96 96 96
*
0000220 96 96 96 96 96 96 96 96 96 64 64 64 64 64 64 64
0000240 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64
*
0000740 64 64 00 00 00 00 00 00 00 00 00 00 00 00 00 00
0000760 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00
Please could you suggest why Blacklist output is the entire chromosome. Thanks.