(517) 353 0368

  • Github
  • Twitter
  • Linkedin
  • Research Gate

mansfeld@msu.edu

  • Home

  • My Research

  • Publications

  • Software

  • Blog

  • More

    Archive

    March 2017

    February 2017

    December 2016

    November 2016

    October 2016

    May 2016

    Please reload

    Tags

    AddOrReplaceReadGroups

    Aspera

    BAM

    BQSR

    BWA

    Bootstrapping

    ENA

    FASTX

    FastAlternativeReferenceMaker

    FastQC

    GATK

    GenotypeGVCF

    HaplotypeCaller

    INDELs

    Joint Genotyping

    Nipponbare

    Picard

    PicardTools

    QTLseq

    Read Groups

    SNPs

    SRA

    SelectVariants

    Sequencing

    Trimmomatic

    VariantFiltration

    automation

    bioinformatics

    bwa

    dbSNP

    download from SRA

    fasta

    fastq

    ftp

    g.vcf

    hpcc

    parallel

    qsub

    reproducing

    rice

    samTools

    samtools

    unix

    Please reload

    Reproducing Yang et al - HaplotypeCaller

    March 8, 2017

    The re-calibration (BQSR) step didn't finish. I will use the on-the-fly capacity in HaplotypeCaller (HC) to apply the BQSR. 

     

    I am also trying to use the g.vcf pipeline for the first time. It is not really necessary but I want to see what the output is and wanted to give it a shot. I specifically want to see if it gives me allele depth in positions that are snps in one sample and not the other. This is important for the downstream analysis.

     

    The difference between a vcf and g.vcf?

     

    NOTES: "When running in `-ERC GVCF` or `-ERC BP_RESOLUTION` modes, the emitting and calling confidence thresholds are automatically set to 0. This cannot be overridden by the command line. The thresholds can be set manually to the desired levels in the next step of the workflow (GenotypeGVCFs)"

     

     Update:

    I knew I was missing something when I couldn't find how to parallelize GATK. Following the instructions and suggestions here, I added the -nct 8 option to the PrintReads function it greatly increased the speed. I then submitted a HC function (on the new bam) with -nct 4 and it is running fast too.

     

     

    Please reload

    Recent Posts

    Reproducing Yang et al. - Chr1 Joint Genotyping

    March 9, 2017

    Reproducing Yang et al - HaplotypeCaller

    March 8, 2017

    Reproducing Yang et al - BQSR

    March 6, 2017

    Please reload

    mansfeld@msu.edu

    Plant Soil Science Bldg 1066 Bogue St., East Lansing Michigan

    • Github
    • twitter
    • linkedin
    • generic-social-link

    ©2016 by Ben Mansfeld. Proudly created with Wix.com