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.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
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
-
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
Latest Articles
Collapse
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
-
by seqadmin
Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...-
Channel: Articles
03-22-2024, 06:39 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
24 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
25 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
21 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
52 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Comment