voutcn / megahit

Ultra-fast and memory-efficient (meta-)genome assembler
http://www.ncbi.nlm.nih.gov/pubmed/25609793
GNU General Public License v3.0
590 stars 135 forks source link

Megahit for meta-transcriptome assembly #269

Open ARW-UBT opened 4 years ago

ARW-UBT commented 4 years ago

Hello,

I found a few contradicting posts in bioinformatics forums whether Megahit could be used or must not be used for metatranscriptome assemblies from microbiome meta-RNA-seq data. I have it a try but realized that the contigs, that several extrem long contigs were produced (e.g. over 100k bases).

I think, there are no transcripts of this lenth in bacteria. So, what is you opinion on using Megahit for metatranscriptome analyses? Do you have special options that should be set in order to tell Megahit not to make assemblies as lon a possible?

WJT0925 commented 2 years ago

I wonder if the eukaryotic metatranscriptome will not be affected by the overlapping peak, because the transcriptome lacks intron and intergenic region sequence, and the length of each gene normally assembled according to the overlapping peak is cut off, right?

brauliovill commented 1 year ago

Did you ever figure out why your contigs were so long (>100k)? I'm getting the same issue

ARW-UBT commented 1 year ago

No, I did not figure out the reason for this, and I have not further inspected the content of the long assemblies. This would probably the next step to see if these long contigs are of 'biological' origin, or whether these are assembly artifacts. Insufficient removal of library adapters could cause such troubles...