Возник возможно дурацкий вопрос к биоинформатикам.
Хочется отклассифицировать метагеномные короткие риды. Проблема в том что отклассифицировать надо не по стандартным базам а чисто по своим референсным последовательностям (среди которых есть и достаточно близкие). Кто это нормально умеет делать? Смотрим сейчас на бракен-кракен и centrifuge, а еще что разумное есть?
Выравнивание ридов на свои референсы, потом разделить .sam на части по тому, на какие референсы они (риды) попали