BLAST tabular output

[emc2rae]

Hi! I'm having trouble getting the output I need. The "sseqid" option seems to behave differently depending on what other options are included. Is there any way to get both the ID number AND the accession in the same output (in any format)?

--outfmt 6 sseqid evalue bitscore
1334999579 0.0 1567

--outfmt 6 sseqid stitle evalue bitscore
WP_103141275.1 WP_103141275.1 DNA gyrase subunit A [Bacillus cereus] 0.0 1567

[bbuchfink]

Hi, can you tell me what diamond version and what database format you are using?

[emc2rae]

Hi! I'm using diamond v2.0.13.151 and version 4 of the Blast database

…

On Dec 20, 2021, at 11:07 AM, Benjamin Buchfink ***@***.***> wrote: Hi, can you tell me what diamond version and what database format you are using? — Reply to this email directly, view it on GitHub <#539 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AUAUVFBEBLTRGN6ZMNKFE7LUR5IF7ANCNFSM5KJYHUIA>. Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>. You are receiving this because you authored the thread.

[bbuchfink]

Are you using a self-made blast database? If so, you need to use -parse_seqids with makeblastdb to get the accessions. Getting both the accession and the number is not supported at the moment.