diff --git a/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/BastyGenerateFasta.scala b/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/BastyGenerateFasta.scala index ead615fe31c0f73cf54893d42af1852364f14fe7..188e8f82492af0ed9846f39edb9dfa1824a64fe4 100644 --- a/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/BastyGenerateFasta.scala +++ b/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/BastyGenerateFasta.scala @@ -268,7 +268,7 @@ object BastyGenerateFasta extends ToolCommand { val genotype = vcfRecord.getGenotype(cmdArgs.sampleName) if (genotype == null) return fillAllele("", maxSize) - val AD = genotype.getAD + val AD = if (genotype.hasAD) genotype.getAD else Array.fill(vcfRecord.getAlleles.size())(cmdArgs.minAD) if (AD == null) return fillAllele("", maxSize) val maxADid = AD.zipWithIndex.maxBy(_._1)._2 if (AD(maxADid) < cmdArgs.minAD) return fillAllele("", maxSize)