Entering edit mode
7.1 years ago
horsedog
▴
60
Hi, I'm a newbie to blast, I have a batch of sequences with different accession number (NZ_ *)in one fasta file, they're all around 500 bp. Now I want to use them to verify if there is gene fragment that exist on the same genome has not been annotated before. Here is an example: My fasta file is like this :
>Bacteroidia_NZ_CP012938.1_1311892
ATGGAATATATATTGATATTTATTTCGGCAATCTTTGTAAACAACATCGTGTTGTCGCAGTTCCTGGGTATCTGTCCGTT
CCTGGGCGTATCCAAGAAGGTGGAAACTGCGATGGGTATGAGCGCTGCGGTTGCTTTCGTGTTGACCATCGCTACCATCG
TCACTTTCCTGATTCAGAAATTTGTACTGGATGTTTTCGGACTGGGATACTTACAGACGATTACATTTATCCTGGTGATT
GCCGGATTGGTGCAGATGGTGGAAATCATCCTGAAAAAGGTGTCTCCGGCCCTGTATCAGGCGTTGGGTGTGTTCTTGCC
>Bacteroidia_NZ_DS499672.1_245384
ATGGAATATATATTGATATTTATCTCGGCAATCTTTGTCAATAACATTGTGTTGTCGCAATTCCTGGGTATCTGTCCGTT
TCTGGGCGTATCCAAGAAAGTGGAGACAGCGTTGGGCATGTCGGCAGCAGTGGCATTCGTGCTTACTATTGCTACAATCG
TAACGTTCCTCATTCAGAAGTTTGTGCTCGATGCTTTCGGTCTGGGCTATCTTCAAACCATCACCTTTATCCTTGTGATT
GCCGCATTGGTGCAGATGGTGGAGATTATCCTGAAAAAAGTATCTCCTTCTTTATACCAGGCTTTAGGTGTGTTCTTGCC
>Bacteroidia_NZ_CP012937.1_2643295
ATGGAATATATATTGATATTTATCTCGGCAATCTTTGTAAACAACATTGTGTTGTCGCAGTTCCTGGGAATCTGTCCGTT
CCTGGGTGTTTCCAAGAAGGTGGAAACAGCGATGGGTATGAGTGCGGCAGTTGCTTTCGTGCTGACTATCGCTACCATCG
TCACATTCCTGATTCAGAAGTTTGTGTTGGACGCTTTCGGATTGGGGTACTTACAGACAATCACTTTTATTCTGGTTATT
GCCGGACTGGTGCAGATGGTAGAAATTA
Now I'd like to do blast for the whole genome with AC number NZ_CP012938.1 for the first sequence, NZ_DS499672.1 for second one, and NZ_CP012937.1 for the third one to see if there is any potential "same" gene in the same genome that has not been found before . Could anybody provide any advice?
The number after NZ_* is the starting position of CDs so just forget it. And all my sequences are bacterial sequences.