nloyfer / wgbs_tools

tools for working with Bisulfite Sequencing data while preserving reads intrinsic dependencies
Other
134 stars 37 forks source link

Intersect (wgbstools view -L Atlas.U25.l4.hg19.bed Left_Ventricle_STL001.IGF2.pat.gz) #68

Closed AzlanNI closed 7 months ago

AzlanNI commented 7 months ago

Test1.zip Hello everyone!

I am currently trying out stuff and want to use my ONT Bam files as input for the uxm deconvolution algorithm. I was trying to generate the sliced pat files which only have entries which are present in the Atlas.U25.l4.hg19.tsv. According to the Tutorial i used

wgbstools view -L Atlas.U25.l4.hg19.tsv Left_Ventricle_STL001.IGF2.pat.gz

I got

So i tried to adjust the tsv file to have the 5 columns 'chr' 'start' 'end' 'startCpG' and 'endCpG' in the bedfile only and ensured that the file is tab-delimited.

wgbstools view -L Atlas.U25.l4.hg19.bed Left_Ventricle_STL001.IGF2.pat.gz So i still got the same Error. And now i was trying too look for the basecode to understand why this Error is occurring but i did not find the Error Message in the view.py or cview.cpp files. Maybe someone could help me out.

The Left_Ventricle file i got from wgbs_tools/tutorials and the Atlas i got from /UXM_deconv/supplemental. I was just trying how how to filter the pat file to only have the entries which are included in the Atlas for my later usage. In the Zip file i just uploaded my files if anyone wants to have a look at them.

I removed the headers from the .bed file and it created a Intersected document which is empty. Is it cause there are no overlaps between these both files i am using ?

I am thankful for any advice! kind regards, Azlan

AzlanNI commented 7 months ago

can be closed!