Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • SnpEff matching policy

    Hi there,

    I'm digging a little bit into SnpEff and everything works fine, but I have a doubt for which I didn't find an answer in any of the SnpEff papers/documentation.
    I have a variant (rs4372192 ) in the intronic region of one gene (SAMD11) and it returns the annotations for all Ensembl non protein-coding transcripts in an interval of +/-5KB from this variant. Everything fine up to here.
    But, not all the protein-coding transcripts for the gene (SAMD11) are matched, only one. Why is it not returning all protein-coding transcripts? Anybody knows if this is configurable? But, in any case is there any reason I should know for this working.

    SnpEff output:
    Code:
    EFF=DOWNSTREAM(MODIFIER||||NOC2L|processed_transcript|CODING|ENST00000327044|),DOWNSTREAM(MODIFIER||||NOC2L|processed_transcript|CODING|ENST00000477976|),DOWNSTREAM(MODIFIER||||NOC2L|processed_transcript|CODING|ENST00000483767|),DOWNSTREAM(MODIFIER||||NOC2L|processed_transcript|CODING|ENST00000496938|),DOWNSTREAM(MODIFIER||||SAMD11|protein_coding|CODING|ENST00000420190|),INTRON(MODIFIER||||SAMD11|protein_coding|CODING|ENST00000342066|),TRANSCRIPT(MODIFIER||||SAMD11|protein_coding|CODING|ENST00000474461|),TRANSCRIPT(MODIFIER||||SAMD11|protein_coding|CODING|ENST00000478729|),UPSTREAM(MODIFIER||||SAMD11|protein_coding|CODING|ENST00000464948|),UPSTREAM(MODIFIER||||SAMD11|protein_coding|CODING|ENST00000466827|)
    Ensembl visualization:
    http://www.ensembl.org/Homo_sapiens/...01cce786943248


    Cheers,

    Pablo.

  • #2
    I will answer myself
    A clue for this might be that the non-matched transcripts are marked as novel in Ensemble, instead of known. So, there are two alternatives: first, that the database SnpEff is using does not include these transcripts; or second that they are filtered out as they are less confident.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Exploring the Dynamics of the Tumor Microenvironment
      by seqadmin




      The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
      07-08-2024, 03:19 PM
    • seqadmin
      Exploring Human Diversity Through Large-Scale Omics
      by seqadmin


      In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
      06-25-2024, 06:43 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 07-10-2024, 07:30 AM
    0 responses
    26 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-03-2024, 09:45 AM
    0 responses
    201 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-03-2024, 08:54 AM
    0 responses
    212 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-02-2024, 03:00 PM
    0 responses
    193 views
    0 likes
    Last Post seqadmin  
    Working...
    X