Tweet
Hello, is there an easy way to retrieve the sequences from the hits in the blastn output? I just need the sequences without the gap insertions "-". Thank you.
-
-
Look at the "blastdbcmd" program examples on this page. http://www.ncbi.nlm.nih.gov/books/NB...BLAST_search_r .. if you are looking to retrieve the subject sequences.Last edited by GenoMax; 06-19-2013, 08:06 AM. -
Thank you GenoMax. I tried using blastdbcmd but for some reason it doesn't seem to be working as described on that page.
I have downloaded the Gallus 4 genome from NCBI and ran makeblastdb using blast+ 2.2.25. When I run:
blastdbcmd -db gga_ref_Gallus_gallus-4.0_all.db -entry all | grep ">" | head
I get:
>gnl|BL_ORD_ID|0 gi|358485511|ref|NC_006088.3| Gallus gallus breed Red Jungle fowl, inbred line UCD001 chromosome 1, Gallus_gallus-4.0, whole genome shotgun sequence
>gnl|BL_ORD_ID|1 gi|358485510|ref|NC_006089.3| Gallus gallus breed Red Jungle fowl, inbred line UCD001 chromosome 2, Gallus_gallus-4.0, whole genome shotgun sequence
>gnl|BL_ORD_ID|2 gi|358485509|ref|NC_006090.3| Gallus gallus breed Red Jungle fowl, inbred line UCD001 chromosome 3, Gallus_gallus-4.0, whole genome shotgun sequence
>gnl|BL_ORD_ID|3 gi|358485508|ref|NC_006091.3| Gallus gallus breed Red Jungle fowl, inbred line UCD001 chromosome 4, Gallus_gallus-4.0, whole genome shotgun sequence
>gnl|BL_ORD_ID|4 gi|358485507|ref|NC_006092.3| Gallus gallus breed Red Jungle fowl, inbred line UCD001 chromosome 5, Gallus_gallus-4.0, whole genome shotgun sequence
>gnl|BL_ORD_ID|5 gi|358485506|ref|NC_006093.3| Gallus gallus breed Red Jungle fowl, inbred line UCD001 chromosome 6, Gallus_gallus-4.0, whole genome shotgun sequence
>gnl|BL_ORD_ID|6 gi|358485505|ref|NC_006094.3| Gallus gallus breed Red Jungle fowl, inbred line UCD001 chromosome 7, Gallus_gallus-4.0, whole genome shotgun sequence
>gnl|BL_ORD_ID|7 gi|358485504|ref|NC_006095.3| Gallus gallus breed Red Jungle fowl, inbred line UCD001 chromosome 8, Gallus_gallus-4.0, whole genome shotgun sequence
>gnl|BL_ORD_ID|8 gi|358485503|ref|NC_006096.3| Gallus gallus breed Red Jungle fowl, inbred line UCD001 chromosome 9, Gallus_gallus-4.0, whole genome shotgun sequence
>gnl|BL_ORD_ID|9 gi|358485502|ref|NC_006097.3| Gallus gallus breed Red Jungle fowl, inbred line UCD001 chromosome 10, Gallus_gallus-4.0, whole genome shotgun sequence
But when I try running:
blastdbcmd -db gga_ref_Gallus_gallus-4.0_all.db -entry 358485511
I get:
Error: 358485511: OID not found
BLAST query/options error: Entry not found in BLAST database
The entry seems to be there so why can't blastdbcmd recognize it when I try to query it?Comment
-
Did you run makeblastdb with '-parse_deflines' ?savetherhino.orgComment
-
Ok I just ran makeblastdb with the -parse_seqids option and it works now. What is the difference between parse_seqids and parse_deflines?Comment
-
Sorry, -parse_seqids is the flag I was thinking..Originally posted by rkizen View Postsavetherhino.orgComment
-
See http://blastedbio.blogspot.jp/2012/1...cbi-blast.html for my thoughts on the difference between the NCBI style IDs and free-format user provided identifiers (and the problems using them with blastdbcmd).Comment
-
Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.
The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
...Yesterday, 10:05 AM -
With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.
Introduction
Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...05-22-2026, 06:42 AM -
Studying ecosystems means dealing with complex, multi-species communities that are hard to observe at scale. This complexity, however, hides many important questions to be answered, from how biogeochemical cycles work and how climate change can affect species distribution to how conservation strategies can work best.
Genomics, particularly since the expansion of NGS, has transformed ecosystem ecology. By sequencing environmental DNA, we can now assess biodiversity without direct...05-06-2026, 09:04 AM
Comment