This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Collect ITS sequences from bacterial genomes

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
• 1,429 views • link
updated by Mensur Dlakic Mensur Dlakic • written by star715 • star715

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.

• 0 views • link
written by h.mon h.mon

0 answers

No answers yet.

Log in to answer this question.

More posts like this