Recently, I’ve tried using a newer version of MEGAN6 Community Edition (MEGAN 6.7.6) to parse the BlastTAB output from Diamond’s BlastX mode, using NCBI nr from February 2017.
For one of the data sets, taxonomy algorithm malfunctioned with no error. The reported number of reads corresponded to the number of unique reads in the BlastTAB ouptput (= reads with hits), but the sum of all assigned reads and the unassigned reads was not even close to that amount. A problem must have occurred during the Naive LCA algorithm, as the “Total reads” and “Alignments” after the parsing phase drastically exceeded the “Total reads” and “Alignments” after the Naive LCA algorithm phase, as stated in the Messages box.
After a specific point, reads in the file got dropped and were not represented anywhere in the taxonomy tree. The specific point is reproducible and is included within the attached segment of the BlastX/Diamond output. The last properly parsed read is 63H9Q:00896:01124. The same file could be parsed properly with MEGAN 6.5.8. Additional details about the parameters and conditions of parsing are included below. Do you have any ideas about the problem?
Thanks for the bug report. It took me a couple of hours to figure this one out. Your file exposed a subtle bug in the code that accesses reads and alignments in an rma6 file. It threw an exception that caused the remaining reads to be ignored. I have fixed the problem, please update to 6.7.10.
D
Thank you for the swift reply!
I have just tried parsing the same file with the freshly installed 6.7.10, but the outcome is unfortunately still the same for me.
That is strange. I used the same file that I uploaded here. After processing the file on 6.7.10 and again after uninstalling and reinstalling 6.7.10, the result was still the same (Screenshots below, 6.5.8 for comparison).
For what it’s worth, we’re using a BioLinux system and MEGAN 6.7.10 was installed using the MEGAN_Community_unix_6_7_10.sh installer. It was installed in parallel in its own directory, leaving the previous versions intact.
could you please re-download the MEGAN installer. Apparently, there was a slip-up when I uploaded the latest version of MEGAN and the version number 6.7.10 was used twice (shouldn’t happen, but I was boarding a flight when I uploaded the fix to the problem that you had reported, so must have messed up). The build date on your MEGAN window title bar is March 24th, on mine it is March 25th…
D
