Dear all I have 2 files with 3787 gene sequences of one genome and another file containing 3762 gene sequences of another bacterial genome which are highly similar so I have decided to do MEGABLAST. (these genomes having only sequencing errors). I want starting to end exact match to query sequences and parse output. please give suggestions. and what should be the significant e value in this case and how to parse exact matched sequences based on MEGABLAST parameters. Is there any way or I have to look all sequence alignment manually.
•
written
by
sharmatina189059 •
1 answer
Use -outfmt 6 'std qlen', then you will have everything you need. You are looking for lines with pident = 100 and qend = qlen. I think that should do it. e-value should be low, but doesn't really matter as you're looking for full-length matches of gene sequences and not trying to determine whether a pattern that you found is likely to be random or not.
Log in to answer this question.
-
NCBI MEGABLAST versus UCSC BLAT
written by a615ebfb •
-
Best way to BLAST orthogroups?
written by MJ Mantas •
-
One gene shows alignment with core sequences as well as with accessory gene sequences at 50% identi…
written by sharmatina189059 •
-
Finding short, highly similar matches to long sequences using BLASTN
written by CC •
-
automate command using bash or perl script?
written by sharmatina189059 •
-
Is there a tool for making a synteny blocks of a gene with its downstream and upstream genes using …
written by sharmatina189059 •
-
Pathway annotation using KEGG
written by sharmatina189059 •
-
Pathway annotation using sequence alignment
written by sharmatina189059 •
-
Query regarding BLASTKOALA (KEGG) for pathway annotation
written by sharmatina189059 •
-
Which phylogenetic method is best for single gene tree?
written by wanderingstefan •
