I recently made a mistake in a script, and it still confuses me what exactly happened. Since part of my analysis relies on some of these details, perhaps some of you can elaborate.
I recently assembled a metagenome (metaspades) and mapped the reads back onto it (bwa). Usually, not all reads map back, because not all reads are assembled. So far so good.
However, I recently accidentally re-assembled only those unmapped reads, and some of them assembled into 1000+ basepair sequences. Can anybody explain to me exactly how that can occur, since they (apparently) did not assemble the first time? How come they assemble now?
Thanks a lot! :)