SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Minia: ultra-low memory contigs assembly rayanc De novo discovery 15 04-01-2013 05:31 AM
Strange Result in Deseq Output cburke11 Bioinformatics 2 12-04-2012 08:35 AM
SAMTools index stats strange result? tatinhawk Bioinformatics 3 09-01-2011 08:29 AM
Strange sequencing result of Hi-seq 2000 angelbo Sample Prep / Library Generation 0 08-09-2011 11:46 PM
strange result of Double-stranded cDNA Synthesis wayland Illumina/Solexa 1 11-29-2010 08:23 AM

Reply
 
Thread Tools
Old 01-29-2013, 06:47 AM   #1
yzzhang
Member
 
Location: florida

Join Date: Jan 2013
Posts: 67
Default minia assembly strange result

Dear all,
I am a newbie in de novo assembly. I just used software minia 1.4683 to assembly my data (plant genome about 380 Mb,Hiseq 2000, 101 bp reads, totally 139 million paired end reads). At first, I trimmed the low quality reads and adapter contaminated reads , after trimming, the average reads quality is above 20 revealed by fastqc.
But the assembly result from minia is strange, I tried Kmer 31,35, 39, 43, from the assembly results, most contigs started with ploy A (about twenty bp), and with the Kmer grow bigger,the GC content grows down . Does anyone face things like that?
Thank you!
yun
yzzhang is offline   Reply With Quote
Old 02-02-2013, 02:25 PM   #2
rchikhi
Member
 
Location: France

Join Date: Jan 2013
Posts: 13
Default

I communicated via private messages with yun.

I checked a human assembly made with Minia, and indeed, the first contigs outputted often start with polyA. This is not an artifact of yun's data. It happens when the polyA kmer is present in the data, and this is NOT an assembly error. This has to do with how Minia it chooses a starting kmer, and should not impact assembly correctness.

To verify this assertion, I just computed full-length contigs alignments (in my case, to the human genome) to make sure that the ployA start of the contigs was correct. Visual inspection of the first dozen of contigs shows that it is.

So, to sum up, that polyA start is normal and the contigs are likely correct.

The slight variation of GC content should not be specific to Minia, and is likely a consequence of lower coverage due to higher kmer length.

Last edited by rchikhi; 02-02-2013 at 02:31 PM.
rchikhi is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 08:32 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO