(517) 353 0368

  • Github
  • Twitter
  • Linkedin
  • Research Gate

mansfeld@msu.edu

  • Home

  • My Research

  • Publications

  • Software

  • Blog

  • More

    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

    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

    Recent Posts

    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...

    Read More

    Reproducing Yang et al - BQSR

    March 6, 2017

    So since this paper is doing QTLseq in rice I decided to take advantage of the fact that rice has a dbSNP database.

    I want to see if I can run BQSR with a merged rice dbSNP file I received from a colleague.

    So first thing I found out was that the contigs in the two files...

    Read More

    Back for more BQSR

    November 19, 2016

    I was happy with the BQSR bootstrapping results. 

    It seemed that after two rounds of bootstrapping we reached convergence, i.e. the correct quality scores were not changed much after the second round.

    Variant quality score distributions:

    However, bootstrapping BQSR is mea...

    Read More

    An alternate reference - making the fasta

    October 27, 2016

    Ok so once the bases have been calibrated using a bootstrapped BQSR, variants have been called on calibrated data, snps and indels filtered and combined..... we can now make the Gy_genome.fasta.

    Using FastAlternativeReferenceMaker (still part of GATK).

    Edit 10/31/16

    Make...

    Read More

    An alternate reference - recalibrating bases

    October 27, 2016

    From GATK tutorial:

    "Variant calling algorithms rely heavily on the quality scores assigned to the individual base calls in each sequence read. These scores are per-base estimates of error emitted by the sequencing machines. Unfortunately the scores produced by the mach...

    Read More
    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