Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • aeonsim
    started a topic Using AlleleIndices with HTSJDK

    Using AlleleIndices with HTSJDK

    Hey, anyone here experienced with HTSJDK?

    I'm currently using it to read VCF files (bi-allelic microarray derived) and grab all the genotypes for all the individuals then transform them into a Byte Array, storing each genotype as the sum of it's indices (so ./. > -2, 0/0 > 0, 0/1 > 1, 1/1 > 2 etc). Using other language (cyvcf2 etc) libraries they often natively provide genotypes as indices. However, with HTSJDK when requesting genotypes from a VariantContext with getGenotypes() you get the genotypes in base form (something like [A*, C], [A*, A*], [C, C] where each allele is it's own object). Currently I'm using getAlleleIndices() from VariantContext to convert these Genotype objects back into an Array of indices which I can then sum. This is a bit slower than I'd ideally like, as such is there a better or faster way of doing this rather than using getAlleleIndices() to transform each Genotype into it's Indices?

Latest Articles

Collapse

  • seqadmin
    Multiomics Techniques Advancing Disease Research
    by seqadmin


    New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

    A major leap in the field has
    ...
    02-08-2024, 06:33 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Today, 04:11 PM
0 responses
14 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-21-2024, 08:52 AM
0 responses
34 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-20-2024, 08:57 AM
0 responses
24 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-14-2024, 09:19 AM
0 responses
56 views
0 likes
Last Post seqadmin  
Working...
X