<?xml version="1.0" encoding="UTF-8"?><rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom"
	xmlns:sy="http://purl.org/rss/1.0/modules/syndication/"
		>
<channel>
	<title>Comments for Heliconius Homepage</title>
	<atom:link href="http://www.heliconius.org/comments/feed/" rel="self" type="application/rss+xml" />
	<link>http://www.heliconius.org</link>
	<description>Home of the passion vine butterflies</description>
	<lastBuildDate>Fri, 10 May 2013 09:22:42 +0000</lastBuildDate>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
	<generator>http://wordpress.org/?v=3.5.1</generator>
	<item>
		<title>Comment on Morning sun by Chris Jiggins</title>
		<link>http://www.heliconius.org/2012/morning-sun/comment-page-1/#comment-546</link>
		<dc:creator>Chris Jiggins</dc:creator>
		<pubDate>Fri, 10 May 2013 09:22:42 +0000</pubDate>
		<guid isPermaLink="false">http://www.heliconius.org/?p=524#comment-546</guid>
		<description><![CDATA[Sorry, but I should call this Heliconius erato demophoon, although even that name may now be invalid - see the recent paper by Hill et al (doi:10.1111/mec.12290)]]></description>
		<content:encoded><![CDATA[<p>Sorry, but I should call this Heliconius erato demophoon, although even that name may now be invalid &#8211; see the recent paper by Hill et al (doi:10.1111/mec.12290)</p>
]]></content:encoded>
	</item>
	<item>
		<title>Comment on Consortium Meeting by Jordan</title>
		<link>http://www.heliconius.org/2013/consortium-meeting/comment-page-1/#comment-540</link>
		<dc:creator>Jordan</dc:creator>
		<pubDate>Tue, 07 May 2013 07:35:53 +0000</pubDate>
		<guid isPermaLink="false">http://www.heliconius.org/?p=802#comment-540</guid>
		<description><![CDATA[Hi there just wanted to give you a brief heads up and let you know 
a few of the images aren&#039;t loading correctly. I&#039;m not sure 
why but I think its a linking issue. I&#039;ve tried it in two different browsers and both show the same outcome.]]></description>
		<content:encoded><![CDATA[<p>Hi there just wanted to give you a brief heads up and let you know<br />
a few of the images aren&#8217;t loading correctly. I&#8217;m not sure<br />
why but I think its a linking issue. I&#8217;ve tried it in two different browsers and both show the same outcome.</p>
]]></content:encoded>
	</item>
	<item>
		<title>Comment on Consortium Meeting by jim mallet</title>
		<link>http://www.heliconius.org/2013/consortium-meeting/comment-page-1/#comment-539</link>
		<dc:creator>jim mallet</dc:creator>
		<pubDate>Tue, 07 May 2013 03:01:40 +0000</pubDate>
		<guid isPermaLink="false">http://www.heliconius.org/?p=802#comment-539</guid>
		<description><![CDATA[Heliconius biology and genomics meeting, Cmabridge MA, 12-13 July 2013

Attendance of this meeting is free, except for possible charges for any group lunch/coffee consumption. However, travel and accommodation cost must be covered personally; some help for transatlantic travel may be available from Chris Jiggins&#039;s BBSRC meeting grant. A block of rooms will be reserved at the nearby Sheraton Commander Hotel from the evening of the 11th to the morning of the 14th. Harvard has discount with the Sheraton Commander of $199.00/per night. For more information please visit http://www.starwoodhotels.com/sheraton/index.html. If you would like us to add you to the reservation list please include a note to that effect when submitting your registration confirmation. If you do not included that you would like to be added to the reservation list we will assume that you will be making other accommodations for your stay and a room will not be reserved for you.

Some other lodging possibilities include the very conveniently located Irving House at Harvard, http://irvinghouse.com/about/story/ . This offers a range of rooming options with nightly costs upwards from $135.00 a night. If you are looking for something more cost efficient please consider searching for a room in a local shared apartment or house on Airbnb, https://www.airbnb.com/. This is also a good resource if you are traveling with a group that would like to “rent” the entirety of a private residence for the duration of your stay.

For the first day (morning of 12th) at least, there will be a limited number of talk slots available over the course of the meeting. If you would like to present a talk at the meeting please include a working title for your talk and a brief description of the talk with your registration confirmation. Please note talk slots will be filled on a first-come first-serve basis, with some preference given to active or recent work from each group internationally. After the talks we will arrange free periods for discussion of future directions in genomics and other international collaborative ventures.

Please send all attendance confirmations, including if you would like to have a room reserved and your talk information to Nikki Hughes, nhughes@fas.harvard.edu, no later than May 30, 2013. If you have any question or concerns, do not hesitate to direct your inquiries to Nikki at the email address above or at 617-495-0508.

Look forward to seeing you all in July.

Regards,

Nikki Hughes and Jim Mallet]]></description>
		<content:encoded><![CDATA[<p>Heliconius biology and genomics meeting, Cmabridge MA, 12-13 July 2013</p>
<p>Attendance of this meeting is free, except for possible charges for any group lunch/coffee consumption. However, travel and accommodation cost must be covered personally; some help for transatlantic travel may be available from Chris Jiggins&#8217;s BBSRC meeting grant. A block of rooms will be reserved at the nearby Sheraton Commander Hotel from the evening of the 11th to the morning of the 14th. Harvard has discount with the Sheraton Commander of $199.00/per night. For more information please visit <a href="http://www.starwoodhotels.com/sheraton/index.html" rel="nofollow">http://www.starwoodhotels.com/sheraton/index.html</a>. If you would like us to add you to the reservation list please include a note to that effect when submitting your registration confirmation. If you do not included that you would like to be added to the reservation list we will assume that you will be making other accommodations for your stay and a room will not be reserved for you.</p>
<p>Some other lodging possibilities include the very conveniently located Irving House at Harvard, <a href="http://irvinghouse.com/about/story/" rel="nofollow">http://irvinghouse.com/about/story/</a> . This offers a range of rooming options with nightly costs upwards from $135.00 a night. If you are looking for something more cost efficient please consider searching for a room in a local shared apartment or house on Airbnb, <a href="https://www.airbnb.com/" rel="nofollow">https://www.airbnb.com/</a>. This is also a good resource if you are traveling with a group that would like to “rent” the entirety of a private residence for the duration of your stay.</p>
<p>For the first day (morning of 12th) at least, there will be a limited number of talk slots available over the course of the meeting. If you would like to present a talk at the meeting please include a working title for your talk and a brief description of the talk with your registration confirmation. Please note talk slots will be filled on a first-come first-serve basis, with some preference given to active or recent work from each group internationally. After the talks we will arrange free periods for discussion of future directions in genomics and other international collaborative ventures.</p>
<p>Please send all attendance confirmations, including if you would like to have a room reserved and your talk information to Nikki Hughes, <a href="mailto:nhughes@fas.harvard.edu">nhughes@fas.harvard.edu</a>, no later than May 30, 2013. If you have any question or concerns, do not hesitate to direct your inquiries to Nikki at the email address above or at 617-495-0508.</p>
<p>Look forward to seeing you all in July.</p>
<p>Regards,</p>
<p>Nikki Hughes and Jim Mallet</p>
]]></content:encoded>
	</item>
	<item>
		<title>Comment on Consortium Meeting by jim mallet</title>
		<link>http://www.heliconius.org/2013/consortium-meeting/comment-page-1/#comment-515</link>
		<dc:creator>jim mallet</dc:creator>
		<pubDate>Sat, 20 Apr 2013 01:59:03 +0000</pubDate>
		<guid isPermaLink="false">http://www.heliconius.org/?p=802#comment-515</guid>
		<description><![CDATA[You should contact me if you have not done so already for a place at this meeting. The meeting is at Harvard University.]]></description>
		<content:encoded><![CDATA[<p>You should contact me if you have not done so already for a place at this meeting. The meeting is at Harvard University.</p>
]]></content:encoded>
	</item>
	<item>
		<title>Comment on Mimicry among Heliconius biologists by Nicola Nadeau</title>
		<link>http://www.heliconius.org/2013/mimicry-among-heliconius-biologists/comment-page-1/#comment-337</link>
		<dc:creator>Nicola Nadeau</dc:creator>
		<pubDate>Mon, 18 Feb 2013 16:49:09 +0000</pubDate>
		<guid isPermaLink="false">http://www.heliconius.org/?p=785#comment-337</guid>
		<description><![CDATA[There also seems to be some mimicry of names occurring, perhaps best exemplified by Simon Martin who is mimetic in both first and last names. Perhaps the Heliconius community has just got large enough that these things are expected by chance, someone should do a statistical analysis...]]></description>
		<content:encoded><![CDATA[<p>There also seems to be some mimicry of names occurring, perhaps best exemplified by Simon Martin who is mimetic in both first and last names. Perhaps the Heliconius community has just got large enough that these things are expected by chance, someone should do a statistical analysis&#8230;</p>
]]></content:encoded>
	</item>
	<item>
		<title>Comment on Aligning Heliconius short read sequences by John Davey</title>
		<link>http://www.heliconius.org/2013/aligning-heliconius-short-read-sequences/comment-page-1/#comment-308</link>
		<dc:creator>John Davey</dc:creator>
		<pubDate>Mon, 28 Jan 2013 12:40:05 +0000</pubDate>
		<guid isPermaLink="false">http://www.heliconius.org/?p=700#comment-308</guid>
		<description><![CDATA[Thanks Nicola - Platypus may be worth a try, but the only data I can find on it suggests it has large false positive rates compared to GATK (http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0049110). Testing SNP callers is a whole other post!]]></description>
		<content:encoded><![CDATA[<p>Thanks Nicola &#8211; Platypus may be worth a try, but the only data I can find on it suggests it has large false positive rates compared to GATK (<a href="http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0049110" rel="nofollow">http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0049110</a>). Testing SNP callers is a whole other post!</p>
]]></content:encoded>
	</item>
	<item>
		<title>Comment on Aligning Heliconius short read sequences by Nicola Nadeau</title>
		<link>http://www.heliconius.org/2013/aligning-heliconius-short-read-sequences/comment-page-1/#comment-307</link>
		<dc:creator>Nicola Nadeau</dc:creator>
		<pubDate>Mon, 28 Jan 2013 12:26:45 +0000</pubDate>
		<guid isPermaLink="false">http://www.heliconius.org/?p=700#comment-307</guid>
		<description><![CDATA[The Stampy developers have put out their own SNP caller, Platypus (http://www.well.ox.ac.uk/platypus). Ludovik Duvaux, a postdoc with Roger Butlin, recommended it to us and has been using it on his Stampy aligned data from aphids. I haven&#039;t tried it yet because I&#039;ve been wanting to stick with tried and tested methods, but it might deal with some of Stampy&#039;s peculiarities better than GATK. Just though I would mention it because Megan brought up the issue of how well Stampy and GATK work together.]]></description>
		<content:encoded><![CDATA[<p>The Stampy developers have put out their own SNP caller, Platypus (<a href="http://www.well.ox.ac.uk/platypus" rel="nofollow">http://www.well.ox.ac.uk/platypus</a>). Ludovik Duvaux, a postdoc with Roger Butlin, recommended it to us and has been using it on his Stampy aligned data from aphids. I haven&#8217;t tried it yet because I&#8217;ve been wanting to stick with tried and tested methods, but it might deal with some of Stampy&#8217;s peculiarities better than GATK. Just though I would mention it because Megan brought up the issue of how well Stampy and GATK work together.</p>
]]></content:encoded>
	</item>
	<item>
		<title>Comment on Aligning Heliconius short read sequences by Megan</title>
		<link>http://www.heliconius.org/2013/aligning-heliconius-short-read-sequences/comment-page-1/#comment-302</link>
		<dc:creator>Megan</dc:creator>
		<pubDate>Wed, 23 Jan 2013 23:58:43 +0000</pubDate>
		<guid isPermaLink="false">http://www.heliconius.org/?p=700#comment-302</guid>
		<description><![CDATA[Here are the BWA parameters I am using for 100 bp paired end reads:
n=8 (maximum edit distance)
o=2 (maximum number of gap opens)
e=3 (maximum number of gap extensions)
l=35 (seed length)
k=2 (maximum edit distance in the seed)

We still have a ton of alignment problems, but this is the best I could come up with.  It allows more mismatches, a little less so in the seed (that could probably stand to be relaxed).  This also increases the gap parameters, but BWA still performs poorly in that respect.  Based on what I have read, BWA really is not designed for this level of divergence and we see a lot of coverage gaps across the reference.  For now I think we are stuck with this until we have an erato genome.]]></description>
		<content:encoded><![CDATA[<p>Here are the BWA parameters I am using for 100 bp paired end reads:<br />
n=8 (maximum edit distance)<br />
o=2 (maximum number of gap opens)<br />
e=3 (maximum number of gap extensions)<br />
l=35 (seed length)<br />
k=2 (maximum edit distance in the seed)</p>
<p>We still have a ton of alignment problems, but this is the best I could come up with.  It allows more mismatches, a little less so in the seed (that could probably stand to be relaxed).  This also increases the gap parameters, but BWA still performs poorly in that respect.  Based on what I have read, BWA really is not designed for this level of divergence and we see a lot of coverage gaps across the reference.  For now I think we are stuck with this until we have an erato genome.</p>
]]></content:encoded>
	</item>
	<item>
		<title>Comment on Aligning Heliconius short read sequences by John Davey</title>
		<link>http://www.heliconius.org/2013/aligning-heliconius-short-read-sequences/comment-page-1/#comment-300</link>
		<dc:creator>John Davey</dc:creator>
		<pubDate>Wed, 23 Jan 2013 19:13:36 +0000</pubDate>
		<guid isPermaLink="false">http://www.heliconius.org/?p=700#comment-300</guid>
		<description><![CDATA[Thanks Megan - the problems with BAC mapping are important I think. I suspect several aligners, certainly Stampy included, make assumptions about how likely a read is to be placed on a genome, and if only a partial genome is used as a reference, those assumptions will almost certainly be wrong.

What parameters are you using for BWA that work well in erato?]]></description>
		<content:encoded><![CDATA[<p>Thanks Megan &#8211; the problems with BAC mapping are important I think. I suspect several aligners, certainly Stampy included, make assumptions about how likely a read is to be placed on a genome, and if only a partial genome is used as a reference, those assumptions will almost certainly be wrong.</p>
<p>What parameters are you using for BWA that work well in erato?</p>
]]></content:encoded>
	</item>
	<item>
		<title>Comment on Aligning Heliconius short read sequences by Megan</title>
		<link>http://www.heliconius.org/2013/aligning-heliconius-short-read-sequences/comment-page-1/#comment-299</link>
		<dc:creator>Megan</dc:creator>
		<pubDate>Wed, 23 Jan 2013 18:57:58 +0000</pubDate>
		<guid isPermaLink="false">http://www.heliconius.org/?p=700#comment-299</guid>
		<description><![CDATA[In the special case of mapping whole genome erato data to erato BACs (0.5% of the genome) I am still using BWA, with relaxed mapping parameters, because I think the sensitive aligners don&#039;t handle this case well.  I did an &quot;eyeball&quot; analysis of Stampy and found that it there were too many cases where it mapped junk with high quality (even when accounting for its inflated mapping quality scores).  There were also genomic regions where Stampy was clearly superior, picking up what appeared to be true indels.  There were also regions where I couldn&#039;t tell if it was a true haplotype that BWA had missed, or a misalignment from another portion of the genome. 

I think including GATK in the assessment is useful if that is the SNP caller used downstream (it is for me).  It is important how well the two work together.

As for further testing of aligners, I think this is only useful if we can generate some sort of &quot;golden set&quot; with known, or at least high confidence, answers.  

I think simulations could be useful to give us a known answer and allow us to address some of the known issues.  Although there are always caveats and biases.

We are attempting to improve alignments with de novo assembly and if it works, that may be useful to get us at least closer to a golden set.  Of course you are then adding assembler biases into the mix.]]></description>
		<content:encoded><![CDATA[<p>In the special case of mapping whole genome erato data to erato BACs (0.5% of the genome) I am still using BWA, with relaxed mapping parameters, because I think the sensitive aligners don&#8217;t handle this case well.  I did an &#8220;eyeball&#8221; analysis of Stampy and found that it there were too many cases where it mapped junk with high quality (even when accounting for its inflated mapping quality scores).  There were also genomic regions where Stampy was clearly superior, picking up what appeared to be true indels.  There were also regions where I couldn&#8217;t tell if it was a true haplotype that BWA had missed, or a misalignment from another portion of the genome. </p>
<p>I think including GATK in the assessment is useful if that is the SNP caller used downstream (it is for me).  It is important how well the two work together.</p>
<p>As for further testing of aligners, I think this is only useful if we can generate some sort of &#8220;golden set&#8221; with known, or at least high confidence, answers.  </p>
<p>I think simulations could be useful to give us a known answer and allow us to address some of the known issues.  Although there are always caveats and biases.</p>
<p>We are attempting to improve alignments with de novo assembly and if it works, that may be useful to get us at least closer to a golden set.  Of course you are then adding assembler biases into the mix.</p>
]]></content:encoded>
	</item>
</channel>
</rss>
