Allele specific reads
2
0
Entering edit mode
6.5 years ago

I have SNVs from the patients. And I also have the ATAC-seq data from the same patient. Now I am looking for extracting allele specific reads using the SNV positions. I was wondering if there is any tool available for this purpose. Can we do this using samtools? Any suggestions will be highly appreciated?

SNV allele specific reads ATAC-seq • 1.7k views
ADD COMMENT
2
Entering edit mode
6.5 years ago

You could also use ASEReadCounter of GATK in addition to the tools mentioned by Devon.

ADD COMMENT
0
Entering edit mode
6.5 years ago

Any of the normal allele-specific tools (e.g., lapels and then suspenders) should work. Not that you'll need to randomly assign each heterozygous allele to one genome. The results will likely be better if you try to phase the variants first.

ADD COMMENT
0
Entering edit mode

Thanks for your suggestions. Actually I am very new to this field. What do you mean by phasing the variants. I have already masked the genome using the locations of the variants.

ADD REPLY
2
Entering edit mode

This previous question deals with phasing: What Are Phased And Unphased Genotypes?

ADD REPLY

Login before adding your answer.

Traffic: 2063 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