Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Breakdancer output file question

    Hi all,

    I am using breakdancer to identify SVs on some whole genome samples. This is the first time that I am using the program and am a little confused on the output. The documentation hasn't been that helpful to me either.

    Here is a line of the output:

    22 28948681 11570+11437- 22 32903818 150360+152913- ITX 2117 5 2 WholeGenome_MatePair|2 0.00 BreakDancerMax-1.0r112 |o22

    What I am confused about is columns 3 and 6 (bolded). I can't quite figure out what it is actually telling me.


    Thanks for any clarification

  • #2
    Originally posted by lre1234 View Post
    Hi all,

    I am using breakdancer to identify SVs on some whole genome samples. This is the first time that I am using the program and am a little confused on the output. The documentation hasn't been that helpful to me either.

    Here is a line of the output:

    22 28948681 11570+11437- 22 32903818 150360+152913- ITX 2117 5 2 WholeGenome_MatePair|2 0.00 BreakDancerMax-1.0r112 |o22

    What I am confused about is columns 3 and 6 (bolded). I can't quite figure out what it is actually telling me.


    Thanks for any clarification
    I think its the number of reads supporting the SV on each strand.

    SV detection from paired end reads mapping. Contribute to kenchen/breakdancer development by creating an account on GitHub.

    Comment


    • #3
      Vivek is right. For each coordinate that BreakDancer reports in that output line, the strand of every abnormal read at that location is reported. Thus, in your example, the region around the first coordiante had 11570 abnormally mapped reads on the forward strand and 11437 on the reverse. Note that this is way more than the number of reads reported as supporting the SV call. This is likely not a good call.

      Comment


      • #4
        Originally posted by vivek_ View Post
        I think its the number of reads supporting the SV on each strand.

        https://github.com/kenchen/breakdancer
        Hi,

        as I have read on the documentation you posted, the output should be constituted with 14 columns. But in my case, I have 17 !
        I think that I have :
        1. Chromosome 1
        2. Position 1
        3. Orientation 1
        4. Chromosome 2
        5. Position 2
        6. Orientation 2
        7. Type of a SV
        8. Size of a SV
        9. Confidence Score
        10. Total number of supporting read pairs
        11. Total number of supporting read pairs from each map file

        But starting from here I'm not so sure
        12. Estimated allele frequency

        And this following I don't have:
        13. Software version
        14. The run parameters

        My lines look like this:
        Code:
        #Chr1   Pos1    Orientation1    Chr2    Pos2    Orientation2    Type    Size    Score   num_Reads       num_Reads_lib   Sample_Barke_reheadered.bam     Sample_Bowman_reheadered.bam    Sample_Bulbosum_reheadered.bam  Sample_FT11_reheadered.bam      Sample_Igri_reheadered.bam      Sample_Pubiflorum_reheadered.bam
        chr1H_part1     13958   241+168-        chr1H_part1     14031   241+168-        ITX     -212    34      3       ../barleyBAMs/Sample_Barke_reheadered.bam|3     NA      NA      NA      NA      NA      NA
        chr1H_part1     22659   14+28-  chr1H_part1     22826   14+28-  ITX     -147    89      6       ../barleyBAMs/Sample_Bulbosum_reheadered.bam|6  NA      NA      NA      NA      NA      NA
        chr1H_part1     51554   3+2-    chr1H_part1     51654   3+2-    ITX     -157    38      2       ../barleyBAMs/Sample_Pubiflorum_reheadered.bam|2        NA      NA      NA      NA      NA      NA
        chr1H_part1     87327   42+24-  chr1H_part1     87429   42+24-  ITX     -179    99      18      ../barleyBAMs/Sample_Barke_reheadered.bam|5:../barleyBAMs/Sample_Bulbosum_reheadered.bam|12:../barleyBAMs/Sample_Pubiflorum_reheadered.bam|1    NA      NA      NA      NA      NA      NA
        chr1H_part1     89261   11+5-   chr1H_part1     89314   11+5-   ITX     -176    74      4       ../barleyBAMs/Sample_Bulbosum_reheadered.bam|4  NA      NA      NA      NA      NA      NA
        chr1H_part1     94809   25+24-  chr1H_part1     94860   25+24-  ITX     -154    99      15      ../barleyBAMs/Sample_Bulbosum_reheadered.bam|14:../barleyBAMs/Sample_Pubiflorum_reheadered.bam|1        NA      NA      NA      NA      NA      NA
        chr1H_part1     95422   4+4-    chr1H_part1     95467   4+4-    ITX     -169    55      3       ../barleyBAMs/Sample_Bulbosum_reheadered.bam|3  NA      NA      NA      NA      NA      NA
        chr1H_part1     97915   153+134-        chr1H_part1     98001   153+134-        ITX     -170    99      32      ../barleyBAMs/Sample_Barke_reheadered.bam|10:../barleyBAMs/Sample_Bulbosum_reheadered.bam|15:../barleyBAMs/Sample_Pubiflorum_reheadered.bam|7   NA      NA      NA      NA      NA      NA
        chr1H_part1     99599   14+8-   chr1H_part1     99668   14+8-   ITX     -154    41      2       ../barleyBAMs/Sample_Pubiflorum_reheadered.bam|2        NA      NA      NA      NA      NA      NA
        chr1H_part1     100374  7+3-    chr1H_part1     100449  7+3-    ITX     -164    39      3       ../barleyBAMs/Sample_Barke_reheadered.bam|1:../barleyBAMs/Sample_Bulbosum_reheadered.bam|2      NA      NA      NA      NA      NA      NA
        chr1H_part1     124801  45+81-  chr1H_part1     133975  25+2-   ITX     8473    99      17      ../barleyBAMs/Sample_Barke_reheadered.bam|11:../barleyBAMs/Sample_Igri_reheadered.bam|6 0.21    0.00    0.03    0.09    0.25    0.02
        chr1H_part1     124801  45+64-  chr1H_part1     134111  6+2-    ITX     5569    45      6       ../barleyBAMs/Sample_Barke_reheadered.bam|3:../barleyBAMs/Sample_Bulbosum_reheadered.bam|2:../barleyBAMs/Sample_Igri_reheadered.bam|1   0.21    0.00    0.03    0.10    0.25    0.02
        chr1H_part1     125588  18+2-   chr1H_part1     126418  1+15-   DEL     890     45      8       ../barleyBAMs/Sample_Barke_reheadered.bam|5:../barleyBAMs/Sample_Igri_reheadered.bam|3  0.21    0.02    0.21    0.06    0.38    0.05
        chr1H_part1     127470  8+33-   chr1H_part1     133703  9+33-   INV     5823    99      22      ../barleyBAMs/Sample_Barke_reheadered.bam|16:../barleyBAMs/Sample_Igri_reheadered.bam|6 0.16    0.00    0.01    0.06    0.18    0.00
        Do you have an idea what are the final column I have? Is it the estimated allele frequency for each of my sample?

        Thank you very much, have a nice day
        C.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Strategies for Sequencing Challenging Samples
          by seqadmin


          Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
          03-22-2024, 06:39 AM
        • seqadmin
          Techniques and Challenges in Conservation Genomics
          by seqadmin



          The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

          Avian Conservation
          Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
          03-08-2024, 10:41 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 03-27-2024, 06:37 PM
        0 responses
        13 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-27-2024, 06:07 PM
        0 responses
        11 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-22-2024, 10:03 AM
        0 responses
        53 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-21-2024, 07:32 AM
        0 responses
        69 views
        0 likes
        Last Post seqadmin  
        Working...
        X