Question: ABySS assembly: got unitig instead of contig
gravatar for philousap10
4 days ago by
philousap100 wrote:


I have run ABySS (abyss-pe) command to assemble a de novo genome. I got only unitigs.fa file and not a contig.fa file. Is it normal? Thank you

commad: abyss-pe k=30 name=AA lib='MI.M03992_0274.001.BioOHT_59.AA_R2_Q20HL50 MI.M03992_0274.001.BioOHT_11.AA_R1_Q20HL50' \ MI.M03992_0274.001.BioOHT_11.AA_R1_Q20HL50='MI.M03992_0274.001.BioOHT_11.AA_R1_Q20HL50_R1.fq MI.M03992_0274.001.BioOHT_11.AA_R1_Q20HL50_R2.fq' MI.M03992_0274.001.BioOHT_59.AA_R2_Q20HL50='MI.M03992_0274.001.BioOHT_59.AA_R2_Q20HL50_R1.fq MI.M03992_0274.001.BioOHT_59.AA_R2_Q20HL50_R2.fq' \ se='MI.M03992_0274.001.BioOHT_11.AA_R1_Q20HL50_single.fq MI.M03992_0274.001.BioOHT_59.AA_R2_Q20HL50_single.fq'

assembly • 101 views
ADD COMMENTlink written 4 days ago by philousap100
gravatar for Carambakaracho
4 days ago by
Carambakaracho360 wrote:

Maybe abyss doesn't understand the library name. I'd try again with something like lib='BioOHT_59 BioOHT_11'. Other than that, your command looks fine, however I recommend using the Code (101/010 - button) formatting. The only thing is the k=30 option - is this the ideal kmer size? What read lengths do you use? What is the insert size of your PE library? Maybe unitigs is all abyss can do with your data.

ADD COMMENTlink written 4 days ago by Carambakaracho360


Thanks for your reply! I found it... That happens when something goes wrong during the assembly. Indeed, i got a error warning during the assembly.
I re-run my command fixing the mistake and it worked perfectly! Thank you!

ADD REPLYlink written 4 days ago by philousap100

I realized this had better been an answer. When a post is helpful a thumbs up is generally appreciated and will help others in the future to reevaluate an answer. In case the answer solved your problem, please mark it as accepted for pretty much the same reasons.

ADD REPLYlink written 3 days ago by Carambakaracho360
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1000 users visited in the last hour