Maize Mo17 SNPs

Abstract

Genome resequencing of Mo17 was done as part of the bioMAP project (REF). 477 million 100bp paired-end reads were generated for Mo17 giving an average of 95x coverage. Reads were first trimmed by Trimmomatic (Bolger et al. 2014) and mapped to the maize B73 genome AGPv4 (Jiao et al. 2017) using BWA-MEM (Li and Durbin 2010). PCR duplicates were marked and removed using GATK (McKenna et al. 2010). Variants were called by GATK haplotypecaller and filtered using different filters for SNPs: (QD > 2, FS < 60, MQ > 40, MQRankSum > -12.5, ReadPosRankSum > -8, SOR < 4) and for InDels: (QD > 2, FS < 200, ReadPosRankSum > -20, SOR < 10). Moreover, variants located in regions with unusually high coverage (DP > Mean + 2*SD) and heterozygous calls (GT == '0/1') were also removed. The final variant file containing 8.04 million variants with 164 thousand CDS variants was deposited in DRUM.