How do i search a FASTA database by sequence in seqkit?
1
1
Entering edit mode
12 months ago
Mahmoud Reda ▴ 10

Database: ribogrove_11.217_sequences.fasta (https://mbio.bas-net.by/cager/en/content/56-ribogrove#downloads)

Query sequence: sequence in .fasta file obtained from NCBI. (example: https://www.ncbi.nlm.nih.gov/nucleotide/NR_117620.1)

Iَ want it as: 100% identity matching but query sequence can be shorter/longer than those in the database.

seqkit rRNA genomics 16s gene • 775 views
ADD COMMENT
4
Entering edit mode
12 months ago
GenoMax 141k

You could do it using seqkit grep or locate but in this case you should use a proper search program like blat instead.

ADD COMMENT
0
Entering edit mode

why can't i use seqkit?

ADD REPLY
0
Entering edit mode

i got it

"seqkit locate --pattern-file seq.fasta ribogrove.fasta -o found.fasta"

ADD REPLY

Login before adding your answer.

Traffic: 3162 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6