This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How can I get the number of reads per variant in hisat2?

Hello,

I would like to see how many reads I have per variant. Is there any way to get this information using hisat2 or gatk or any other tool?

alignemnt variant gatk covid hisat2
• 1,237 views • link
updated by Arup Ghosh Arup Ghosh • written by Nemo • Nemo

It's called variant calling. The VCF contains that information. Given your post history it seems you are essentially going blind in your project. I suggest to search for some guided tutorials to get started. The Biostars Handbook would be one option, but it's not entirely free though.

• 0 views • link
written by ATpoint ATpoint

Thanks @ATpoint I have already the variants (vcf format) using gatk haplotypecaller. But I do not see such information like reads per variant in the resulted output file (I see depth of reads and so on but not number of reads per variant). I also looked into biostar handbook the Variant Calling section, and I couldnt find any information related to this.

• 0 views • link
written by Nemo • Nemo

0 answers

No answers yet.

Log in to answer this question.

More posts like this