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

      • SEQadmin2
        Nine Things a Sample Prep Scientist Thinks About Before Sequencing
        by SEQadmin2


        I’m not a sequencing expert. I’m a purification scientist who uses NGS to evaluate workflows my group develops. With this perspective, we think about the sample first and the NGS workflow second. The sequencer is an exceptionally honest reporter, but it can only report on what you give it, so whether you get clean, interpretable data from an NGS workflow is largely determined before you begin.


        Here are nine questions we think about, in roughly the order they matter, before...
        06-18-2026, 07:11 AM
      • SEQadmin2
        From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data
        by SEQadmin2


        Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.


        The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
        ...
        06-02-2026, 10:05 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 06-17-2026, 06:09 AM
      0 responses
      25 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 06-09-2026, 11:58 AM
      0 responses
      42 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 06-05-2026, 10:09 AM
      0 responses
      48 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 06-04-2026, 08:59 AM
      0 responses
      49 views
      0 reactions
      Last Post SEQadmin2  
      Working...