question about annotObj in DESeqDataSetFromSlidingWindows

[LoneKnightz]

Hi,

I followed the tutorial from Htseq-clip and generated everything.

But the sliding windowed annotation file that Htseq-clip created is BED6 format.

It cannot be used for DESeqDataSetFromSlidingWindows as annotObj.

May I ask do you have any idea that can convert annotation from Htseq-clip to DESeqDataSetFromSlidingWindows recoginized annotObj ?

[Distue]

Hi, Please use the htseq-clip mapToId function to convert the bed style annotation to a format which can be parsed by DEWSeq. The flow diagram at https://htseq-clip.readthedocs.io/en/latest/ might help understanding the workflow. The Bed6 format is primarily useful for visual inspection and counting, this additional step speeds up reading in large annotation files in R. We will make sure to add an error and explanation to DEWSeq. Thank you for bringing this up. If you need anything else. Please feel free to contact us. Best, Tom LoneKnightz ***@***.***> schrieb am Di., 29. Nov. 2022, 01:39:

…

Hi, I followed the tutorial from Htseq-clip and generated everything. But the sliding windowed annotation file that Htseq-clip created is BED6 format. It cannot be used for DESeqDataSetFromSlidingWindows as annotObj. May I ask do you have any idea that can convert annotation from Htseq-clip to DESeqDataSetFromSlidingWindows recoginized annotObj ? — Reply to this email directly, view it on GitHub <#6>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AAG2FZF4KBIRBNJAV6AOGU3WKVGFLANCNFSM6AAAAAASN5A73M> . You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>

[LoneKnightz]

Thanks for reply. It really helped and my question got solved.

I followed the old flow diagram and I used wrong output file from Htseq-clip.