Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • bongbimit
    Member
    • Dec 2013
    • 29

    #1

    snpeff error?

    Dear all, i am using SnpEff to annotate my VCF file but i see some error

    i use hg19 of Snpeff database ( i didn't build) to annotate my vcf file

    here is the result

    java -Xmx10g -jar snpEff.jar eff -c snpEff.config -v hg19 accepted_hits-snp-dbsnp.vcf > accepted_hits-snp-dbsnp-eff.vcf
    00:00:00.000 Reading configuration file '../../../../../tools/snpeff_snpsift/snpEff/snpEff.config'. Genome: 'hg19'
    00:00:00.302 done
    00:00:00.302 Reading database for genome version 'hg19' from file '/home/huypham/tools/snpeff_snpsift/snpEff/./data/hg19/snpEffectPredictor.bin' (this might take a while)
    00:00:07.050 done
    00:00:07.081 Building interval forest
    00:00:15.896 done.
    00:00:15.896 Genome stats :
    # Genome name : 'Homo_sapiens (USCS)'
    # Genome version : 'hg19'
    # Has protein coding info : true
    # Genes : 26346
    # Protein coding genes : 20775
    # Transcripts : 47313
    # Avg. transcripts per gene : 1.80
    # Protein coding transcripts : 38612
    # Length errors : 210 ( 0.54% )
    # STOP codons in CDS errors : 193 ( 0.50% )
    # START codon errors : 14 ( 0.04% )
    # STOP codon warnings : 0 ( 0.00% )
    # Total Errors : 239 ( 0.62% )

    # Cds : 388149
    # Exons : 459403
    # Exons with sequence : 459403
    # Exons without sequence : 0
    # Avg. exons per transcript : 9.71
    # Number of chromosomes : 93
    # Chromosomes names [sizes] :
    ...

    NEW VERSION!
    There is a new SnpEff version available:
    Version : 3.6
    Release date : 2014-04-21
    Download URL : http://sourceforge.net/projects/snpe...atest_core.zip


    how can i explain about the error which i bold?
    thank you.
  • ebioman
    Member
    • Aug 2013
    • 41

    #2
    Hello
    These stem from the database which (I guess) is in your case provided by SnpEff.
    If I recall correctly SNPeff updates from time to time their database based on ENSEMBL. Now it all depends how well the databse there is curated.
    These databases will never be flawless since many predicted protein coding genes derive still from automated annotation tools such as e.g. Augustus or SNAP.
    Since SNPeff contains as well the protein length information it will detect whether there are differences between the predicted length according to the CDS and the final protein.
    Somehow you have a few genes included which have e.g. a stop codon within the CDS or some transcripts which are incomplete (START codon errors and maybe as well Length errors). E.g. Augustus is capable to predict as well incomplete proteins (for example close to gaps or the end/start of scaffolds) and might generate a protein without a start-codon.

    Cheers

    Comment

    Latest Articles

    Collapse

    • kinpict.com
      Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
      by kinpict.com
      The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
      Today, 06:37 PM
    • SEQadmin2
      16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
      by SEQadmin2





      NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
      10-30-2026, 06:46 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 11-12-2026, 07:59 AM
    0 responses
    27 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 11-05-2026, 09:25 AM
    0 responses
    35 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 11-03-2026, 11:54 AM
    2 responses
    56 views
    0 reactions
    Last Post kinpict.com  
    Started by SEQadmin2, 10-28-2026, 09:57 AM
    0 responses
    175 views
    0 reactions
    Last Post SEQadmin2  
    Working...