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

The documentation of alignmentSieve seems to be wrong. #1331

Open
ljw20180420 opened this issue Oct 11, 2024 · 0 comments
Open

The documentation of alignmentSieve seems to be wrong. #1331

ljw20180420 opened this issue Oct 11, 2024 · 0 comments

Comments

@ljw20180420
Copy link

If, instead, the protocol is strand-specific, then the first set of integers in a pair would be applied to fragments where read 1 precedes read 2, and the second set to cases where read 2 precedes read 1. In this case, the first value in each pair is applied to the end of read 1 and the second to the end of read 2. Take the following command as an example:

.. code:: bash

$ alignmentSieve -b paired_chr2L.bam \
--minFragmentLength 140 --BED \
--shift -5 3 -1 4 -o fragments.bedpe

Given that, the -5 3 set would produce the following::

     ----> read 1
                 read 2 <----

     ------------------------ fragment

-------------------------------- shifted fragment

and the -1 4 set would produce the following::

     ----> read 2
                 read 1 <----

     ------------------------ fragment

         --------------------- shifted fragment

As can be seen, such fragments are considered to be on the - strand, so negative values then shift to the left on its frame of reference (thus, to the right relative to the + strand).

The diagram for -1 4 seems to be wrong. Because -1 and 4 has opposite sign, the shift of read 1 and read 2 should also in the opposite directions (see the shift for -5 3).

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

1 participant