Entering edit mode
6.4 years ago
Yuyayuya
▴
250
Hello everyone, I'm doing de novo assembly from 150 bp whole genome shotgun reads. My pipeline is:
- Trimmomatic (trimming and quality filtering)
- SOAPdenovo2 (de novo assembly)
- mumMER3 (align to reference sequence).
After aligning contigs to reference sequence, I connect the contigs according to their order. My problem is how to make all contigs in the same orientation?
Thank you!
There is a range of contig reordering tools available:
However, if you're already using Mummer, then you could look into the show-tiling option as it can also fix order and orientation and can generate a pseudo-sequence when
-p filename Output a pseudo molecule of the query contigs to file
option is specified.I think one possible way could be the reverse complement of your contig.