How to get consensus sequences from vcf file for hetero SNP?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Clustering sequences with low identity
written by agata88 88Hi all! I need to cluster sequences with 30-40% of identity threshold. Which tool would you recommend? I've tried CD-HIT, but it's not recommended for …
-
How to download FASTA protein sequences for Escherichia coli strains from Ensembl?
written by agata88 88Hi all! I would like to download FASTA protein sequences from all Escherichia coli strains. At the Ensembl Bacteria page I see that I should …
-
Qiime2 - make table.gza
written by agata88 88Hi all! I am new with Qiime2 and any kind of help would be much appreciated. I have abundance results for every taxonomy (phylum - …
-
Comparison of SNPs between groups of samples
written by agata88 88Hi all! Could you please recommend a tool that can statistically compare SNPs between two groups of vcf's? To be more clear, I have 30 …
-
How to create scaffolds from contigs?
written by agata88 88Hi all! I was using velveth and velvetg to create contigs.fa file. It includes all sequences with kmer=31, min contig length 100 and expected coverage …
-
problem in running vcf-consensus command
written by rajesh.msch 0Hi, I would like to create consensus sequence for an individual from 1000G where the sequence incorporates variants typed for this individual. Using tabix and …
-
How to delete two populations pseudo snps from a mutiple populations vcf file?
written by zhangdezhi008 2Hi all, I used three populations to call snps by SAMtools.And filtered biallelic snps by vcftools, I need to calculate between populations's fst, the original …
-
usearch version, which one use to get results faster?
written by agata88 88Hi all, I would like to parse 1.5GB database by USEARCH tool. Right now I am using 32bit version and need to divide database into …
-
Solanum tuberosum reference, where can I download it?
written by agata88 88Hi all, I would like to perform RNA-seq analysis for potato. Unfortunately I cannot download genome, trancriptome or annotation files from Plant Genome Database (ftp …
-
How to add protein ID transcript ("NP_")?
written by agata88 88Hi all, I would like to add "NP_" number to snp variants. I am using SnpEff and SnpSift. I was using command below to annotate …
How do you plan to phase your variants per allele?
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?
Per @WouterDeCoster's question, what if the first is het AT and the second is het GA - how would you want to report it?
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.
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?
I would discard those reads from analysis ...
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?
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.
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.