How to get consensus sequences from vcf file for hetero SNP?
0
0
Entering edit mode
4.8 years ago
agata88 ▴ 820

Hi again,

How to divide vcf file into two files or call two consensus sequences from vcf file for hetero SNPs?

Best,

Agata

consensus vcf • 1.6k views
ADD COMMENT
1
Entering edit mode

How do you plan to phase your variants per allele?

ADD REPLY
0
Entering edit mode

For example I have two snps in my vcf, one is hom TT and other is het GA. I would like to get two consensus sequences where one is with TG and second TA. Is there a tool that can do something like this?

ADD REPLY
0
Entering edit mode

Per @WouterDeCoster's question, what if the first is het AT and the second is het GA - how would you want to report it?

ADD REPLY
0
Entering edit mode

hom TT, het AT and another het GA I would like to report as: TAG and TTA assuming that TAG is from one read and TTA from another.

ADD REPLY
0
Entering edit mode

So the tools should also look into the bam file to find out which reads are supporting which variant calls? That's making the story more difficult, as you can imagine. What if there is no evidence to derive phase from the bam (no reads spanning) between position 1 (GA) and position 2 (TC)? How to report that?

ADD REPLY
0
Entering edit mode

I would discard those reads from analysis ...

ADD REPLY
0
Entering edit mode

This is getting weird. What do you really want to do, I mean, what's the goal of the analysis? Why do you want to split a vcf file without biological meaning?

ADD REPLY
0
Entering edit mode

I would like to perform HLA typing ... and wanted to follow this article:

https://bmcgenomics.biomedcentral.com/articles/10.1186/1471-2164-14-355

I might be wrong, I am getting confused here. They are using two "original own perl scripts". I am trying to write that. I have amplicons for whole genes not only for HLA but also KIR. I tested few pipelines which are "ready to use" and I obtained different results :/ So, I am trying to figure out the best solution.

ADD REPLY
0
Entering edit mode

Right, so you performed long range PCR of your targets of interest, followed by NGS library prep. It would make things much more clear if you would have stated that in your original post. Have you tried asking the authors for the perl script? In my opinion it is unacceptable not to publish such an important part of their work.

If this project is something you want to continue, you might want to consider getting a MinION and sequence the longe range PCR products directly without shearing.

ADD REPLY

Login before adding your answer.

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