Gene variation study through blast tab limited out put
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Generate SNP/indel annotation
written by oussriradnane 0I have a .tab file with SNPs from a fungi genome and i want to find if the SNPs are located inside genes, are there …
-
How to predict mutation without using variant callers like GATK, Mutect?
written by Mask 18I need to find out mutation in custom gene (gene of my interest). Say suppose I have 100 genes in the list I have just …
-
Need help : gene with highest number of high quality SNPs in reference genome.
written by Varshney 2Hello everyone, this is my problem: I am very new in this analysis. I have done bowtie2 for my genome sequence data and it was …
-
SNP frequencies in Fasta Alignment
written by Eva_Maria 19Hai everyone I have a set of fasta aligned file from this alignment file how i will get snps and its frequencies
-
Extraction of best alignment in blast out put
written by Eva_Maria 19Hey everyone I have about 4000 blast tab limited output, from thees output I want to extract best alignment (first line only) and want to …
-
snp and copy number of variation detection
written by Eva_Maria 19<p>In blast output file (fasta format) how to find all SNPs and copy number of variations . which tool is suits for doing this</p>
-
standalone blast tab delimited out put manipulation
written by Eva_Maria 19I generated a tab delimited output from blast. How can I fetch sequences(based on alignment) from here.
-
Gene specific region visualization in Alignment (.bam) or Variant calling file (.vcf) in IGV
written by ravi.uhdnis 22Hi, I have multiple whole genome alignment files (.bam of a human sample) and also have one variant calling file (.vcf) generated using GATK HaploTypeCaller, …
-
Bacterial RNAseq transcription start sites
written by D.a.forrest 0Hello folks I'm new to the RNAseq field and I was wondering if I could have some advice? I want to identify transcription start sites …
-
SNP to gene mapping into reference genome
written by sadsam 5<p>I have a vcf file containing a list of SNPs from a bacterial isolate which was created by comparing it with a reference strain using …
Using BLAST to call SNPs is not a great idea. BLAST is a local aligner, AFAIK not a tool designed to work with the probabilistic algorithms used to call variants.
The protocol would involve a LOT of manual work, where you'd have to look for matches covering all of the query, with as few mismatches as possible, then examine each mismatch to identify what kind of variant it is, all the while never losing sight of the bigger picture. And because this involves carefully eyeballing each sequence, tab-delimited output will not work for you.
You may want to look up possible alternatives to this protocol.