Biostar Beta. Not for public use.
Can Mlst (Multi Locus Sequence Typing) Profiles Be Built From Raw Reads
1
Entering edit mode
20 months ago
Lee Katz ♦ 2.9k
Atlanta, GA

Hi, I am wondering if anyone has used short read data to make MLST profiles. I imagine it would be more time-efficient and maybe even more accurate to map instead of making a de novo assembly. I've only seen a few hints at it in a blog post from Pathogenomics. These include using ICORN and SRST. I think that the other methods on the post use de novo assembly, which I am trying to avoid. ICORN seems like it could be rather involved, but I'm not sure yet, and SRST reportedly doesn't work with newer versions of Samtools.

So has anyone actually hammered out a set of scripts to map reads against an MLST database and determine a profile?

One difficulty that I can imagine is getting the right sequence on the ends of each MLST fragment. Maybe there are more difficulties.

mapping illumina • 1.4k views
ADD COMMENTlink
0
Entering edit mode

I found that SRST has been updated since that blog post and am trying it out now, but any other choices would be helpful. Thank you.

ADD REPLYlink
0
Entering edit mode

Apart from the totally unhelpful answer below, this remains unanswered. I am very hopeful that someone has the right solution.

I have worked with SRST a little because it is ideally the right answer, but the software has a few bugs. So far I have altered the main executable in the following ways:

  1. Point to the correct path for samtools.pl
  2. Append an empty string to the array "t" around 429 if its length is less than 5

However, it is choking on samtools.pl pileup2fq. I'll see if I can fix it and will write down how to fix the software when/if I can. I believe that I will have change samtools pileup to samtools mpileup somewhere, but I am not sure yet.

ADD REPLYlink
1
Entering edit mode
16 months ago
Rockville, MD

As Nick posted on his blog -- Roll Your Own (based on Anthony Underwood's suggestion. Nick even posted a python script of his own to do this.

ADD COMMENTlink
0
Entering edit mode

That would be smart to make my own and help myself. Thanks!

ADD REPLYlink
1
Entering edit mode
20 months ago
Lee Katz ♦ 2.9k
Atlanta, GA

I fixed SRST with a couple of changes which I listed on their site, https://sourceforge.net/p/srst/discussion/general/thread/dca3f134/

First, I checked an array to make sure it didn't read beyond the last element.

Second, I fixed a subroutine that expected an old version of Samtools.

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1