Take the 2-minute tour ×
Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. It's 100% free, no registration required.

I'm taking an online bioinformatics class and I am stuck on a problem. The resources I have found don't help that much. Here is the problem:

"Assemble the error-free reads from a circular genome on the dataset below. You do not need to consider the reverse complement of the reads.

Dataset: CGTTCT CTAG GACGT TAGAC

As the reads do not have the same length, convert the reads to 4-mers using the approach proposed by Idury and Waterman."

I believe this means I need to convert this into sets of 4 characters each. Do I simply break them up by taking the first 4 characters and last 4 characters like this?

CGTTCT: CGTT --> TTCT
CTAG: CTAG --> CTAG
GACGT: GACG --> ACGT
TAGAC: TAGA --> AGAC

And end up with a new dataset of: "CGTT TTCT CTAG GACG ACGT TAGA AGAC"?

Or is there more to it than that? Thanks for any help.

share|improve this question
1  
I would guess that it would be all 4mers. In this case that would mean that CGTTCT would also yield GTTC. But as I say - this is a guess. –  Alan Boyd Feb 21 at 6:05
1  
It does in fact mean all possible 4mers. I don't know how much material you get from the course, but they refer to this paper: citeseerx.ist.psu.edu/viewdoc/… –  skymninge Feb 21 at 7:12
    
To be sure I'm understanding correctly, does that mean the new dataset should be: "CGTT GTTC TTCT CTAG GACG ACGT TAGA AGAC"? –  Brett Feb 21 at 21:19

1 Answer 1

Actually the new data set will be,

for k = 4 (k-mer value)

CGTTCT -> CGTT GTTC TTCT

or

enter image description here

Source: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2874646/

share|improve this answer

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.