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
                      Recent Advances in Sequencing Analysis Tools
                      by seqadmin


                      The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
                      05-06-2024, 07:48 AM
                    • seqadmin
                      Essential Discoveries and Tools in Epitranscriptomics
                      by seqadmin




                      The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
                      04-22-2024, 07:01 AM

                    ad_right_rmr

                    Collapse

                    News

                    Collapse

                    Topics Statistics Last Post
                    Started by seqadmin, 05-14-2024, 07:03 AM
                    0 responses
                    19 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 05-10-2024, 06:35 AM
                    0 responses
                    44 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 05-09-2024, 02:46 PM
                    0 responses
                    54 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 05-07-2024, 06:57 AM
                    0 responses
                    42 views
                    0 likes
                    Last Post seqadmin  
                    Working...
                    X