Collect ITS sequences from bacterial genomes
0
0
Entering edit mode
12 weeks ago
star715 ▴ 30

Hi, I want to build a phylogenetic tree with ITS sequences of assembled genomes. For this, I need to retrieve only ITS (internal transcribed region) sequences from a genome. Genome annotations provides annotations of rRNA sequences but not ITS. Is there any tool for that?

Currently, I can only think of blasting the genome against ITS sequences and recovering blast hits as ITS sequences from the genomes. Is there a more systematic or better approach or tool for doing it?

Thank you

rRNA ITS internal 16S transcribed taxonomy region • 237 views
ADD COMMENT
1
Entering edit mode

ITS is the region between adjacent 16S and 23S rRNA genes. As you have the 16S and 23S annotations, you can get the ITS intervals with some mathematics, and then use bedops / bedtools / samtools faidx to recover the fasta sequences of these regions.

ADD REPLY

Login before adding your answer.

Traffic: 1981 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