no read counts per allele while using gatk ASEReadCounter
I am using gatk ASEReadCounter to get the read counts per allele. To do so, I used the following command:
gatk ASEReadCounter -R /path_to_genome/hg38_genome/GRCh38.p13.genome.fa -I sample.sorted.bam -V sample.vcf.gz -O output.table
I am not getting any error but output.table file has only header. All the input files in the above command have index file too. I checked the vcf file and bam file and both have similar chromosome notation (with chr in the beginning of chromosome names). vcf file has only SNPs and heterozygote ones. I used GATK4.
Do you know what the problem could be?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
GATK4 ASEReadCounter Function of VCF
written by hgal4248 •
-
Removing multi-variant records from vcf file
written by Emili •
-
dbSNP with RefSeq chromosome notation
written by Tom •
-
create index for genome instead of transcriptome using Salmon
written by Sara
-
GATK issue in running Mutect2
written by jyotsana18237 •
-
HOMER Motif Analysis code not working
written by danielcgingerich •
-
Getting the db SNP VCF file for GRCh38.p13
written by nhaus •
-
GATK BaseRecalibrator Step
written by Seigfried
-
Empty output of allele specific expression analysis by GATK
written by zizigolu
-
ASEReadCounter More then one variant context at position
written by mb2subi •
