Hi,
I need to extract all 16S sequences from a metagenome illumina run: 12 files of ~30M reads each (after qc). I have the java based rRNASelector, which uses hmmsarch and a custom database, it seems to be working but may take several hours per file. Does anyone know if there is a command line alternative so I can script this (rRNAselector is a GUI only).
Thanks,
S.
I need to extract all 16S sequences from a metagenome illumina run: 12 files of ~30M reads each (after qc). I have the java based rRNASelector, which uses hmmsarch and a custom database, it seems to be working but may take several hours per file. Does anyone know if there is a command line alternative so I can script this (rRNAselector is a GUI only).
Thanks,
S.
Comment