Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • how to validate SNPs and Indels after assembly?

    Hi all,
    We have observed some SNP sites and Indels by aligning our contigs to the reference genome. We want to distinguish real variants from those due to sequencing errors.

    I have found several tools for variants detection but all of them work with the raw reads. I think a lot of information would be lost if assembled contigs used for this task... However, we performed a de novo assembly of the 454 reads due to the design of this specific project. It would be great if we could estimate the "confidence" of the variants we found after assembly. Is there any variant validation method working with contigs?

    Or can I verify this by an independent analysis performed on the raw reads? I am afraid the variants identified before and after assembly would be somewhat inconsistent...

    And, is there significant improvement if the variants are identified by using raw reads, rather than assembled contigs?

    Thanks in advance.
    Last edited by sulicon; 02-23-2011, 11:37 AM.

  • #2
    Hi sulicon,

    Are you mapping contigs to the reference seq using gsMapper?

    if your contigs are larger than 2000bp gsMapper will not consider them for mapping to the ref seq, as i got that error during aligning contigs with gsMapper, so it is good to map reads to ref seq.

    Also the gsMapper outputs HCDiff.txt containing high confidence SNP sites and INDELS.
    Last edited by ketan_bnf; 02-22-2011, 09:00 PM.

    Comment


    • #3
      Originally posted by ketan_bnf View Post
      Hi sulicon,

      Are you mapping contigs to the reference seq using gsMapper?

      if your contigs are larger than 2000bp gsMapper will not consider them for mapping to the ref seq, as i got that error during aligning contigs with gsMapper, so it is good to map reads to ref seq.

      Also the gsMapper outputs HCDiff.txt containing high confidence SNP sites and INDELS.
      No. I used gsAssembler for the de novo assembly. Then I BLATed the contigs against reference genome to get the structure of the genes. I observed some SNPs/Indels in the alignments, but had no confidence about the results...

      Comment


      • #4
        If you want to find SNPs, you should map reads to ref seq using gsMapper or map contigs to ref seq using BWA http://bio-bwa.sourceforge.net/, get output in sam, extarct SNPs using SAMTools, Magicviewer.

        You can also further annote that SNPs using variant effect predictor on EnsEMBL.

        Comment


        • #5
          Originally posted by sulicon View Post
          No. I used gsAssembler for the de novo assembly. Then I BLATed the contigs against reference genome to get the structure of the genes. I observed some SNPs/Indels in the alignments, but had no confidence about the results...
          Why don't you just map your reads with gsMapper on your refseq? As you already have seen some SNP from BLATing your contigs, you know where to look for them ...

          Sven

          Comment


          • #6
            If you can generate a SAM/BAM file from your alignment, then the various SNP callers which work on that format should allow you to estimate confidence in the calls.

            However, they will be relying on the quality scores generated by the base caller. I've recently run into a situation on another platform (SOLiD) in another setting (RNA-Seq) in which systematic errors were reinforced, and so some of my very confident calls from the SNP caller were bogus. In the end, nothing beats verifying at least a sample of your variant calls experimentally -- which is how I discovered the trouble in my data.

            Comment


            • #7
              Thanks Ketan and Sven. I have already assembled the contigs by newbler and performed a lot of subsequent analysis. It's better if I needn't to assembled the reads again... Maybe I have to map the reads to reference seq, just for the purpose of SNP detection.

              Comment


              • #8
                @krobison

                Thanks. Could I generate SAM/BAM files from BLAT alignment between the contigs and human genome?

                I think the there would be some information lost if I worked on the contigs, instead of reads. However, compared with aligning raw reads, I guess the de novo assembler has already considered the alignment between reads, and I could provide a quality file for the contigs. But I don't whether the SNP callers could realize sequencing error rate would rise in homopolymer regions. I will have a try.

                We would perform some experiments for variants validation if I could find interesting candidates.

                Comment


                • #9
                  I would definitely just use gsMapper to align to a reference genome. There are some nice output files with confidence, I think they are called HCDiff.txt or similar.
                  We have been doing this and the validation quite a lot of late and the 454 data is very nice for SNP calling, even at low coverages which has really surprised us after fun with Illumina-predicted SNPs at low coverage.

                  Comment


                  • #10
                    Just a word about the Ensembl variant effect predictor- if you have chromosome and base pair positions (on the reference assembly), you can enter in any alleles found at that position as your input. The output will let you know any dbSNP IDs that map to the same position. In this way, you can see if there is a known dbSNP ID for the allele/alternate nucleotide you have found. More is here:



                    It also accepts VCF format:

                    1000genomes.org is your first and best source for all of the information you’re looking for. From general topics to more of what you would expect to find here, 1000genomes.org has it all. We hope you find what you are searching for!


                    The tool itself is here- both a script, and a web interface:



                    Hope that helps.

                    Comment

                    Latest Articles

                    Collapse

                    • seqadmin
                      Best Practices for Single-Cell Sequencing Analysis
                      by seqadmin



                      While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
                      06-06-2024, 07:15 AM
                    • seqadmin
                      Latest Developments in Precision Medicine
                      by seqadmin



                      Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

                      Somatic Genomics
                      “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
                      05-24-2024, 01:16 PM

                    ad_right_rmr

                    Collapse

                    News

                    Collapse

                    Topics Statistics Last Post
                    Started by seqadmin, 06-17-2024, 06:54 AM
                    0 responses
                    10 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 06-14-2024, 07:24 AM
                    0 responses
                    21 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 06-13-2024, 08:58 AM
                    0 responses
                    17 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 06-12-2024, 02:20 PM
                    0 responses
                    20 views
                    0 likes
                    Last Post seqadmin  
                    Working...
                    X