SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
[Velvet,assembly] core dumped occured by runnning velvet matador3000 De novo discovery 0 12-17-2011 07:31 AM
Large K-mer Velvet NGS_user De novo discovery 3 06-10-2011 07:59 PM
Why using k-mer? papori De novo discovery 12 03-07-2011 03:30 AM
remove false positives from less than 17 mer hpgala Illumina/Solexa 0 02-26-2011 12:30 PM
Optimal k-mer and N50? AronaldJ De novo discovery 1 12-28-2010 09:03 AM

Reply
 
Thread Tools
Old 10-03-2011, 10:43 PM   #1
ramadatta.88
Junior Member
 
Location: singapore

Join Date: Sep 2011
Posts: 7
Default Query on K-mer using in velvet

Hi, I am newbie to NGS. TO assemble a genome i started using velvet. During which I came across a term called 'K-mer' which means the overlapping segment in the reads. I would like to know that if the K-mer length assigned is suppose 21, is it that the overlapping is done exactly for 21 nucleotides at the ends or is it a minimum requirement for the reads to match? Could i kindly know the difference..Thanks in advance..
ramadatta.88 is offline   Reply With Quote
Old 10-04-2011, 12:54 AM   #2
boetsie
Senior Member
 
Location: NL, Leiden

Join Date: Feb 2010
Posts: 245
Default

Hi!

A read is hashed into your defined k-mers. And only ends are considered of the k-mers, simply because you are stripping the read.

Say you have a 10bp read, and a k-mer of 5bp, the hashing of the read into kmers will be;

read: ACGATGATAGTA
kmer: ACGAT
kmer: .CGATG
kmer: ..GATGA
kmer: ...ATGAT
kmer: ....GATAG
kmer: .....ATAGT
kmer: ......TAGTA

The program than finds k-1 overlap (thus here: 5-1=4 bp overlap) at the ends of the k-mers.

E.g:
ACGAT
.CGATG


So you find overlap at the ends of the k-mers, but within the reads.

Have a look at this poster from Velvet for maybe detailed information;

http://www.ebi.ac.uk/~zerbino/velvet/velvet_poster.pdf

Regards,
Boetsie
boetsie is offline   Reply With Quote
Old 10-04-2011, 07:23 PM   #3
ramadatta.88
Junior Member
 
Location: singapore

Join Date: Sep 2011
Posts: 7
Default

thank you..My doubt is cleared..
ramadatta.88 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 05:14 PM.


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