Can Mlst (Multi Locus Sequence Typing) Profiles Be Built From Raw Reads
2
1
Entering edit mode
11.2 years ago
Lee Katz ★ 3.2k

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 • 2.4k views
ADD COMMENT
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 REPLY
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 REPLY
1
Entering edit mode
11.2 years ago

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 COMMENT
0
Entering edit mode

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

ADD REPLY
1
Entering edit mode
11.2 years ago
Lee Katz ★ 3.2k

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 COMMENT

Login before adding your answer.

Traffic: 1872 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6