Entering edit mode
20 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.
why can't i use seqkit?
i got it
"seqkit locate --pattern-file seq.fasta ribogrove.fasta -o found.fasta"