Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • accipiter
    Junior Member
    • Sep 2011
    • 4

    #1

    What is the header format for Interval List in PicardTools?

    Hi,

    I have aligned paired-end read data from from a gene resequencing project, and am attempting to run the CalculateHsMetrics tool from PicardTools suite to calculate various metrics, and have a question re interval lists.

    In order to derive % on-target, I need to define locations of baits, and target intervals. This is not a problem, but I cant figure out what the header should look like. The main reference often referred to re formatting of the interval list header is this: http://picard.sourceforge.net/javado...ervalList.html

    ... but I was looking for an actual example, and/or how to generate a header that will work with my processed bam file.

    Thanks in advance,

    Dave
  • accipiter
    Junior Member
    • Sep 2011
    • 4

    #2
    bump.... Anyone?

    Comment

    • lek2k
      Member
      • Aug 2011
      • 32

      #3
      There is an example on the GATK website
      http://www.broadinstitute.org/gsa/wi...line_arguments

      Here's a direct link to the file if you don't wanna read all the text
      ftp://ftp.broadinstitute.org/pub/gsa....interval_list

      Comment

      • LiLin
        Member
        • May 2011
        • 15

        #4
        such as:
        chr1:726340-727101
        chr1:121352000-121352614

        Comment

        • accipiter
          Junior Member
          • Sep 2011
          • 4

          #5
          Originally posted by lek2k View Post
          There is an example on the GATK website
          http://www.broadinstitute.org/gsa/wi...line_arguments

          Here's a direct link to the file if you don't wanna read all the text
          ftp://ftp.broadinstitute.org/pub/gsa....interval_list
          Thanks, I'm assuming that the header is specific to my reference genome.
          Probably a basic question, but how do I go about generating it so that it matches the contig locations in my reference genome specifically?

          thanks

          Comment

          • manducasexta
            Member
            • Mar 2009
            • 12

            #6
            For the possible benefit of future searchers: I found the properly formatted header in the .dict file corresponding to the genome I am using. The .dict file was generated by GATK and can be generated in isolation as described at this link referred to above.

            Comment

            Latest Articles

            Collapse

            • SEQadmin2
              Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
              by SEQadmin2



              Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

              Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
              ...
              10-14-2026, 09:55 AM
            • SEQadmin2
              New Genomics Technologies Take Aim at Long-Standing Limits
              by SEQadmin2


              Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

              We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
              ...
              09-28-2026, 10:25 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, 10-22-2026, 11:13 AM
            0 responses
            26 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-22-2026, 10:48 AM
            0 responses
            20 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-15-2026, 10:52 AM
            0 responses
            36 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-14-2026, 05:39 AM
            0 responses
            43 views
            0 reactions
            Last Post SEQadmin2  
            Working...