I’ve been using MEGAN with LAST according to your MEGAN-LR paper. I’ve replaced NR database with another set of sequences (a subset of NT database). I appreciate the help that I’ve already gotten on this forum.
In this protocol, reads are mapped to a reference database using LAST. MAF file produced by LAST is quite large - hundreds of GB in my case. I’ve tried piping LAST output to maf2daa tool. However this also uses a lot of temporary disk space (if somewhat less then MAF file itself). For a last run it took almost 1TB, when my server ran out of disk :). I’ve been running LAST on 16 threads, not sure f that matters.
Is there a way to make the whole process use less disk?
University of Zagreb,
Faculty of Electrical Engineering and Computing