Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

In readthedocs, there are some false command, and what is the forward and reverse barcode? #15

Open
khsjh opened this issue Jul 24, 2023 · 1 comment

Comments

@khsjh
Copy link

khsjh commented Jul 24, 2023

While running scHiCExplorer with tutorial, I have some questions about scHiCExplorer tools such as Demultiplex and BuildMatrix.
(The tutorial means Analysis of single-cell Hi-C data index)

At first, in case of hicBuildMatrix, three arguments should be used, restrictionSequence, restrictionCutFile and danglingSequence

But in the tutorial, there are no required options in example command.
So I tried that command including these options and it work well.
(When I used this command without three required options, some error message were occurred)
Is that right the example command is false?
also some command is skipped the space between option and argument, such as -n1 -P1-I ...
Creation of Hi-C interaction matrices index is the example.

And what is the forward barcode and reverse barcode?
I think the one is barcode sequence and the other is the i5 sequence.
Is that right?

@joachimwolff
Copy link
Owner

Hi,

Even though it has been a while since you opened the issue, I would still like to reply. It may help.

Concerning hicBuildMatrix: I think the issue arises because we made the input on HiCExplorer more stringent, but did not update this in scHiCExplorer documentation/tutorials. So yes, you have to give restrictionSequence, restrictionCutFIle and the danglingSequence to build the matrix.

Forward and reverse barcode:
Hi-C data is sequenced paired-end, but the forward and reverse strand are mapped individually as single-end. In the data used in the tutorial, the original authors provide the data in the fastq files as follows: 1. line is the forward read, second line the barcode for the forward read, third line the bar code for the reverse read and fourth line the reverse read.

That we have forward and reverse barcodes is provided in the README of the Nagano et al 2017 study: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE94489 / https://www.ncbi.nlm.nih.gov/geo/download/?acc=GSE94489&format=file&file=GSE94489%5FREADME%2Etxt

I hope this clarifies your questions.

Thanks and best,

Joachim

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants