Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Exome Sequencing

    can anyone suggest a pipeline for analysis of exome seq data?

  • #2
    Originally posted by Shruti Madhiwalla View Post
    can anyone suggest a pipeline for analysis of exome seq data?
    Are you looking for snps/mutations/indels? What kind of data do you want at the end of the analysis?

    Comment


    • #3
      I am looking for snps in the data.

      Comment


      • #4
        One possibility:
        1. align with BWA
        2. get variation with SAMTOOLS pileup


        Comment


        • #5
          For SNPs and Indels do try Novoalign as it performs quite well in terms of accuracy but is slower than BWA.

          We also have a Novoalign NGS guide doing this basic variant detection pipeline.

          Basically

          1. Align with novoalign
          2. Sort alignments
          3. Merge if you have multiple runs for the same library
          4. Remove PCR duplicates with samtools or Picard
          5. Run the samtools pileup variation caller
          6. Filter

          See the posted link for command line examples

          Comment


          • #6
            I am also currently trying to find out the best way to handle exome sequencing data (sureselect, sequenced on the solid).
            It seems to me, that most people map against the whole genome, not an exome sequence, probably to reduce false positive SNPs in the end. So my plan for a pipeline looks like this so far:

            1. Align against whole genome (with bioscope in my case)
            2. remove duplicates with picard
            3. SNP calling with bioscope, as well as samtools pileup (to compare results)
            4. filter out the SNPs in the targeted region

            So far I'm not quite sure, of the best way to filter in the last step. I'd be very grateful for some suggestions. :-)

            Comment


            • #7
              Align with bwa.

              Then i suggest you to use GATK. You can analyse multiple samples at once and it gives a robust set of calls and allows for a whole-exome specific pipeline..Cool!

              Once GATK analysis is completed, you can annotate called variations by annovar.

              That's all!

              Comment


              • #8
                We see some decent coverage in non-target regions. Has anyone looked at that? Probably its good data, when one sees more than 10x coverage of coding regions, even when not targeted by the capture kit. However, they could certainly be false!

                @jeckow, whats your experience with annovar, could you comment on its usage, time to run, efficiency, etc..
                --
                bioinfosm

                Comment


                • #9
                  regarding "out of target" sequence freebies... becareful with pseudogenes and paralogs... the capture kits will pull down things that can cross-hybridize with your baits

                  this can also mess up your intended targets...

                  anyone come up with a way to handle these in an automatic fashion? to ignore genes affected by pseudogenes / paralogs?

                  Comment


                  • #10
                    I am trying to understand the exome-capture dataset we got for human and mouse (separate projects). To begin with, I am interested in estimating:

                    1. How much of the exome is covered by at least 1 (or N) base(s). (breadth of exome coverage)
                    2. What is the depth at which each exon is covered. (Depth of exome coverage)

                    Anybody has done this kind of analysis?
                    Please suggest if there are tools that I could use for this purpose.
                    Harsha

                    Comment


                    • #11
                      This information is easy to get using samtools on a BAM file produced by any decent aligner (like bwa). The "samtools mpileup" function will report the number of reads covering each reference position (depth), and you can simply use awk and bed tools to generate your metrics.

                      I would caution you against using a minimum of 1 base coverage as your metric for coverage. For diploid sequences, one read is completely useless. For coverage, we count bases at which we can call a genotype with >99.9% confidence. This usually works out to be somewhere in the 10-20x range, depending on the bases seen.

                      Hope that helps!
                      --Nancy

                      Comment


                      • #12
                        Estimating breadth and depth of coverage

                        AFter trying a couple different approaches, I concluded that using the coverageBed script in BEDTools is the easiest way to determine the breadth and depth of coverage.

                        coverageBed -abam reads.bam b exons.bed -hist >result.txt

                        When the run completes, you need to look at the end of the result.txt file for the 2 column data to plot the histogram.

                        Comment


                        • #13
                          hey thats great!

                          From bedtools "New "per base depth feature" (-d) added to coverageBed. This reports the per base coverage (1-based) of each feature in file B based on the coverage of features found in file A. For example, this could report the per-base depth of sequencing reads (-a) across each capture target (-b). "

                          Guess will have to try it out to see what it really looks like!
                          --
                          bioinfosm

                          Comment


                          • #14
                            Hi!

                            I am also working on Exome annotation. I have 454 sequencing data, so i am planning to use MOSAIK for aligning reads to ref. chromosome and then extract data from aligned_sorted file to BAM/SAM format, going to analyze with samtools.

                            Is this the right pipeline for Exome annotation SNP finding of 454 data?
                            Does bwa supports long reads (>=200) of 454 seq data? can it handle that?

                            Thanks,

                            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, Yesterday, 06:37 PM
                            0 responses
                            11 views
                            0 likes
                            Last Post seqadmin  
                            Started by seqadmin, Yesterday, 06:07 PM
                            0 responses
                            10 views
                            0 likes
                            Last Post seqadmin  
                            Started by seqadmin, 03-22-2024, 10:03 AM
                            0 responses
                            51 views
                            0 likes
                            Last Post seqadmin  
                            Started by seqadmin, 03-21-2024, 07:32 AM
                            0 responses
                            67 views
                            0 likes
                            Last Post seqadmin  
                            Working...
                            X