This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Hybrid Assembly +Gapcloser With 454 And Illumina Data

Hi everyone, I am fairly new in genome assembly. I have some 454 reads as well as illumina reads for the same baacterial strain. When i am doing assembly separtely of 454 data and illumina data. But there is some gap in genome after assembly. I want to minimize the gap by insilico method. Can anyone suggest some idea how to do it. If i do hybrid assembly how much it will be helpful and which software will be useful.

Or else i will go for separate assembly and try to fill the gap by some gap filler tool. Which will be ideal way to handle such problem.

Thank you.

454 illumina
• 6,070 views • link
updated by rtliu rtliu • written by HG HG

1 answer

You can have a try with Minimus2 to merge contigs, or use other assembler like MaSuRCA for a combined assembly.

Minimus2 is a modified version of the minimus pipeline designed for merging one or two sequence sets (S1,S2). It uses a nucmer based overlap detector which is much faster than the Smith-Waterman hash-overlap program used by minimus.

MaSuRCA is whole genome assembly software. It combines the efficiency of the de Bruijn graph and Overlap-Layout-Consensus (OLC) approaches. MaSuRCA can assemble data sets containing only short reads from Illumina sequencing or a mixture of short reads and long reads (Sanger, 454).

• 0 views • link
written by rtliu rtliu

I did 454 assembly with Newler which follow OLC and with my illumina data set i used Spades which follow de Bruijn grap. Now if i merge all the contig with Minimus2 will it be same as the algotithnm MaSuRCA follow? or result will be different whats your opinion. I did not try with any hybrid assembly approach because either they use OLC or de Bruijn . Could you please help me out

• 0 views • link
written by HG HG

Dear rtliu , Thanks for your suggestion. Now i can run with illumina pair end reads. Could you please let me know how i can include my 454 reads also in config file and make a hybrid assembly. My config file looks like this and i have only .ace file from 454 reads. please help me out. 

PATHS
JELLYFISH_PATH=/home/hiren/Desktop/MaSuRCA-2.0.3.1/bin/
SR_PATH=/home/hiren/Desktop/MaSuRCA-2.0.3.1/bin/
CA_PATH=/home/hiren/Desktop/MaSuRCA-2.0.3.1/CA/Linux-amd64/bin
END
DATA
PE= pe 250 20 /home/hiren/Desktop/MaSuRCA-2.0.3.1/data/ResetH45_S5_L001_R1_001.fastq /home/hiren/Desktop/spades/MaSuRCA-2.0.3.1/data/ResetH45_S5_L001_R2_001.fastq
END
PARAMETERS
GRAPH_KMER_SIZE=auto
USE_LINKING_MATES=1
JF_SIZE=1800000000
DO_HOMOPOLYMER_TRIM=0
NUM_THREADS=12
END
• 0 views • link
written by HG HG

I already request for .sff file . But could you please tell me how i put it inside command line after sffToCA. ??

• 0 views • link
written by HG HG

Dear rtliu, Now i got .sff file for all the strain please let me know how can i feed the assembler this two different file fastq & .sff ???

Thank you..

• 0 views • link
written by HG HG

using above sffToCA script and .sff file, generating an output file, say, 454reads.frg then add "OTHER" line to your config file (config reference: sr_config_example.txt)

"DATA

PE= pe 180 20 /FULL_PATH/frag_1.fastq /FULL_PATH/frag_2.fastq

OTHER=/FULL_PATH/454reads.frg

END"

then run the 'runSRCA.pl' script against your config file to generate a script called 'assemble.sh' finally run the 'assemble.sh' (e.g. >bash assemble.sh)

• 0 views • link
written by rtliu rtliu

Log in to answer this question.

More posts like this