Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • nullalleles
    Junior Member
    • Jun 2011
    • 5

    GATK UnifiedGenotyper SNP calling with Agilent 50Mb target enrichment

    I am trying (and failing) to call SNPs using a reference file (.BED) specifically for the Agilent 50Mb Sureselect enrichment sequencing. I am using the GATK UnifiedGenotyper with the following commands. Any ideas? Thanks!

    java -jar ~/Sting/dist/GenomeAnalysisTK.jar \
    -R ~/Documents/refFASTA/SureSelect_All_Exon_50mb.bed \
    -l INFO \
    -T UnifiedGenotyper \
    -I ~/Documents/next_gen/sorted.N.dupesremoved.recali.bam \
    -o ~/Documents/SNP/N_AgilentSS_50MB_gatk.vcf \
    -stand_call_conf 50.0 \
    -stand_emit_conf 10.0 \
    -dcov 150
  • nullalleles
    Junior Member
    • Jun 2011
    • 5

    #2
    solved

    The above code should have contained the argument "-L" for the interval list (bed) file and also a "-R" for the reference fasta file!

    Comment

    Latest Articles

    Collapse

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, Today, 06:09 AM
    0 responses
    11 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-09-2026, 11:58 AM
    0 responses
    33 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-05-2026, 10:09 AM
    0 responses
    38 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-04-2026, 08:59 AM
    0 responses
    43 views
    0 reactions
    Last Post SEQadmin2  
    Working...