Hi everybody,
I have reads from transcriptome and the reference genome, but the reference is fragmented (more than 20k contigs), Somebody knows if is possible to improve the genome assembly using the RNA-SEQ reads? I`m sure that it is not possible with IDBA, SOAP and SPADES
Thanks!
You may be able to improve local regions with RNAseq reads but chances of finding spanning reads (to bridge contigs, assuming you have paired-end transcriptome reads) would be smaller. Even if you do find them there would be no way to fill the sequence in between.
Since the reference is that fragmented I guess it can help bridging contigs, correctly orientating them, but indeed the sequence in between will not be filled in.