Announcement

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

  • brentp
    replied
    Thanks Dameon, do you know of an example of this in Plink or Beagle?

    Leave a comment:


  • Dameon
    replied
    I would use GATK available from the Broad's Genome Analysis Toolkit and indel realign, recalibrate, and call and filter genotypes following the Broad's recommended best practices. Then I would take the finished and annotated vcf file of just the SNPs (-glm of SNP in the Unified Genotyper) and use vcftools to convert the vcf file into Plink or Beagle format. Plink or Beagle would be better tools to conduct your association testing rather than try to perform a series of intersects. If you are familiar with SAS and JMPGenomics, I believe the latest version of JMPGenomics can readily import vcf format SNP files and SAS is an excellent platform to conduct any number of test you wish to perform. Good luck.

    Leave a comment:


  • best-practices for assigning read-groups for snp-caller

    Hi, We have samples from 4 individuals from 2 tissues each == 8 total:

    tissue is blood/lung
    individuals are control/disease

    so, let's (ignore that 2 samples is too few and) say
    lc1,lc2,bc1,bc2,ld1,ld2,bd1,bd2
    where l == lung, b==blood, c == control d == disease, and the last number is the individual of that disease category.

    If we want to call lung-specific, disease-related SNPs, is it better to send to the caller with ld2,ld1 vs all others by giving ld1,ld1 a read-group of disease and the rest a read-group of control?
    Or, is it better to just call snps where each of the 8 groups has it's own read-group and pull out the SNPs by intersecting post-hoc?

    I'm not asking wether it's better to send in subsets, but how to divvy up the subsets.
    thanks.

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, 02-21-2024, 08:52 AM
0 responses
30 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-20-2024, 08:57 AM
0 responses
18 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-14-2024, 09:19 AM
0 responses
52 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-12-2024, 03:37 PM
0 responses
452 views
0 likes
Last Post seqadmin  
Working...
X