no read counts per allele while using gatk ASEReadCounter
0
0
Entering edit mode
8 months ago
Sara ▴ 220

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?

GATK • 177 views
ADD COMMENT

Login before adding your answer.

Traffic: 782 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6