Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • binfoUser
    Member
    • Jan 2016
    • 22

    Using library htsjdk vcf files

    I want to analyse VCF files that have annotations. I used programs like VEP to annotate the file and now I want to see and manipulate these annotations that are in the file.
    I'm using a Java library htsjdk but until now I didn't find any method to manipulate one column (the annotations are added in the INFO column).

    Can someone help?
  • lindenb
    Senior Member
    • Apr 2010
    • 143

    #2
    there is currently no automatic method to parse the ouptut of VEP. you need to find the description of the INFO field for VEP from the https://samtools.github.io/htsjdk/ja...VCFHeader.html ( getInfoHeaderLine("CSQ") )

    Code:
    ##INFO=<ID=CSQ,Number=.,Type=String,Description="Consequence type as predicted by VEP. Format: Allele|Gene|Feature|Feature_type|Consequence|cDNA_position|CDS_position|Protein_position|Amino_acids|Codons|Existing_variation|DISTANCE|STRAND|SYMBOL|SYMBOL_SOURCE|HGNC_ID">
    and then, using something like "variantContext.getAttribute("CSQ")" , get the CSQ strings and split the result.

    Comment

    • binfoUser
      Member
      • Jan 2016
      • 22

      #3
      But the VCFHeader only access the file header and it gives me that line that you showed with id "CSQ". But now I want to extract the info added that is present in the data lines after the header.
      I don't understand how to use the "variantContext.getAttribute("CSQ")" that you said. Does that give me that information?

      Comment

      Latest Articles

      Collapse

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 06-09-2026, 11:58 AM
      0 responses
      15 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 06-05-2026, 10:09 AM
      0 responses
      26 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 06-04-2026, 08:59 AM
      0 responses
      37 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 06-02-2026, 12:03 PM
      0 responses
      61 views
      0 reactions
      Last Post SEQadmin2  
      Working...