I'm mapping 36bp single ended reads to C.elegans cDNA reference by bwa, trying to make variant calling after getting mpileup files. The problem is the...