Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • g781
    Member
    • May 2010
    • 25

    #1

    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

  • kinpict.com
    Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
    by kinpict.com
    The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
    Today, 06:37 PM
  • SEQadmin2
    16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
    by SEQadmin2





    NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
    10-30-2026, 06:46 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 11-12-2026, 07:59 AM
0 responses
27 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 11-05-2026, 09:25 AM
0 responses
35 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 11-03-2026, 11:54 AM
2 responses
56 views
0 reactions
Last Post kinpict.com  
Started by SEQadmin2, 10-28-2026, 09:57 AM
0 responses
175 views
0 reactions
Last Post SEQadmin2  
Working...