#print column names
print(dimnames(annoVar)[2])

#
#DEFINE YOUR FILTERING STRATEGY HERE
#
sub = subset(annoVar,(annoVar$Qual1>=55|annoVar$Qual1==0))
sub = subset(sub, dbsnp.1000g=='')
#
# END OF FILTERING STRATEGY
#

#check which of the selected genes occur in 4 patients
t = table(sub$Sample,factor(sub$Gene))
selection = t[,which(apply(t,2,min)!=0)]

#print results
cat("This filtering leads to: ",dim(selection)[2]," genes.\n")
cat("The first 10 genes are: \n")
print(selection[,1:10])
print(dimnames(selection)[2])