# Finding specific sequences in the database

**URL:** <https://forum.biobakery.org/t/finding-specific-sequences-in-the-database/4405>\
**Category:** MetaPhlAn\
**Created:** [November 18, 2022, 7:54am UTC](https://forum.biobakery.org/t/finding-specific-sequences-in-the-database/4405 "2022-11-18T07:54:01Z")\
**Posts on this page:** 3\
**Page:** 1

<div class="post-metadata">

**Author:** ![delahondes](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/delahondes/32/1774_2.png) [@delahondes](https://forum.biobakery.org/u/delahondes)\
**Post date:** [November 18, 2022, 7:54am UTC](https://forum.biobakery.org/t/finding-specific-sequences-in-the-database/4405/1 "2022-11-18T07:54:01Z")

</div>

Dear bioBakery team, thanks for all your excellent work. We are trying to compare different classification software with metaphlan4 (which is really excellent) and in order to align the different answers of the different tools we need to assess how some specific bacterial species report in other tools and for this we need to access the reference sequence. For a lot of them, we find them in the database ([http://cmprod1.cibio.unitn.it/biobakery4/metaphlan\_databases/mpa\_vJan21\_CHOCOPhlAnSGB\_202103.tar](http://cmprod1.cibio.unitn.it/biobakery4/metaphlan_databases/mpa_vJan21_CHOCOPhlAnSGB_202103.tar)) but there are a few of them that we cannot find.

For instance  
k\_\_Bacteria|p\_\_Firmicutes|c\_\_Clostridia|o\_\_Clostridiales|f\_\_Lachnospiraceae|g\_\_Lachnospiraceae\_unclassified|s\_\_Eubacterium\_rectale|t\_\_SGB4933\_group

We cannot find either this exact string in the database, or SGB4933 or s\_\_Eubacterium\_rectale. However this string was given by metaphlan4 so there is something we do not understand, and very afraid this is a mistake from us - sorry about that.

Raynald

---

<div class="post-metadata">

**Author:** ![aitor.blancomiguez](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/aitor.blancomiguez/32/86_2.png) [@aitor.blancomiguez](https://forum.biobakery.org/u/aitor.blancomiguez)\
**Post date:** [November 21, 2022, 1:08pm UTC](https://forum.biobakery.org/t/finding-specific-sequences-in-the-database/4405/2 "2022-11-21T13:08:02Z")

</div>

Hi @delahondes  
The identifiers of the marker genes might not correspond to the SGB reported in the profiles for the SGB\_groups (they can be retrieved from one of the SGBs in the group that is not the group representative). However, you can find here ([http://cmprod1.cibio.unitn.it/biobakery4/metaphlan\_databases/mpa\_vJan21\_CHOCOPhlAnSGB\_202103\_marker\_info.txt.bz2](http://cmprod1.cibio.unitn.it/biobakery4/metaphlan_databases/mpa_vJan21_CHOCOPhlAnSGB_202103_marker_info.txt.bz2)) the association of each marker of the database to its correspondent SGB

---

<div class="post-metadata">

**Author:** ![delahondes](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.biobakery.org/delahondes/32/1774_2.png) [@delahondes](https://forum.biobakery.org/u/delahondes)\
**Post date:** [March 3, 2023, 4:40pm UTC](https://forum.biobakery.org/t/finding-specific-sequences-in-the-database/4405/4 "2023-03-03T16:40:28Z")

</div>

Hi @aitor.blancomiguez, I had not taken the time to thank you, this was really helpful.