Header Leaderboard Ad

Collapse

the question from pairToBed output

Collapse

Announcement

Collapse

SEQanswers June Challenge Has Begun!

The competition has begun! We're giving away a $50 Amazon gift card to the member who answers the most questions on our site during the month. We want to encourage our community members to share their knowledge and help each other out by answering questions related to sequencing technologies, genomics, and bioinformatics. The competition is open to all members of the site, and the winner will be announced at the beginning of July. Best of luck!

For a list of the official rules, visit (https://www.seqanswers.com/forum/sit...wledge-and-win)
See more
See less
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • the question from pairToBed output

    Hi everyone,

    I was testing the pairToBed with version 2.17.0, but I was confused with the output.

    -- CASE 1 --
    Here is the input A in bedpe:
    Code:
    chrBB	100	230	chrBB	330	400	paired-read-1	60	+	-
    Here is the input B in bed6:
    Code:
    chrBB	200	250	gene1	1	+
    chrBB	300	350	gene1	1	+
    Here is the command for running pairToBed:
    Code:
    pairToBed -a A.bedpe -b B.bed6 -bedpe -type both
    Here is the outputs:
    Code:
    chrBB   100     230     chrBB   330     400     paired-read-1   60      +       -       chrBB   200     250     gene1   1       +
    chrBB   100     230     chrBB   330     400     paired-read-1   60      +       -       chrBB   300     350     gene1   1       +

    -- CASE 2 --
    Here is the input A in bedpe:
    Code:
    chrBB	100	230	chrBB	330	400	paired-read-1	60	+	-
    Here is the input B in bed6:
    Code:
    chrBB	200	350	gene1	1	+
    Here is the command for running pairToBed:
    Code:
    pairToBed -a A.bedpe -b B.bed6 -bedpe -type both
    Here is the outputs:
    Code:
    chrBB   100     230     chrBB   330     400     paired-read-1   60      +       -       chrBB   200     250     gene1   1       +
    chrBB   200     350     chrBB   330     400     paired-read-1   60      +       -       chrBB   200     350     gene1   1       +

    I would expect just one record. But both case seem to me, pairToBed would show the two records for each read in PE-read. Is that correct?

    Thanks.

Latest Articles

Collapse

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 08:56 PM
0 responses
8 views
0 likes
Last Post seqadmin  
Started by seqadmin, Yesterday, 07:33 AM
0 responses
8 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-31-2023, 07:50 AM
0 responses
4 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-26-2023, 09:22 AM
0 responses
9 views
0 likes
Last Post seqadmin  
Working...
X