text: user2442: I have used canu (correct) and smartdenovo to assembly PacBio reads. Next I am going to polish my assembly using illumina pair end reads by Pilon. However, I found there is high heterozygosity in my pair end reads when I used jellyfish to do k-mer analysis. There are two peaks in the k-mer frequencies distribution plot. The size of my genome is about 190M, and it is a diploid organism.
  
  What I am doubting is that can I use this pair end data to polish my assembly?
  If not, I have to do sequencing again. 
  
  user2085: Resequencing your sample is not going to make it any less heterozygous. I am assuming that you are working with a diploid organism? Please update your question and specify the genome size and ploidy of your organism (the type of organism will help us update our answers, too).
  The answer, yes, you can use Pilon + short reads to correct the PacBio assembly even with a highly heterozygous sample. To do so you should use the --diploid option in Pilon to let it know that there will be heterozygous sites. The way that it corrects heterozygous SNPs and heterozygous indels is different, and it seems to perform best with SNPs and sometimes does not correctly identify heterozygous indels (source).  
  Pilon itself doesn't handle phasing, but just corrects with whatever reads are in the bam file that it is provided with. The end result of providing a bam with reads from both haplotypes is that the corrected sequence will be a mosaic of the dominant allele at the heterozygous site as was observed in the short read data. This may or may not select for some alleles based on sequencing coverage bias of Illumina sequencers.
  
original_text: 2442: <p>I have used canu (correct) and smartdenovo to assembly PacBio reads. Next I am going to polish my assembly using illumina pair end reads by Pilon. However, I found there is high heterozygosity in my pair end reads when I used jellyfish to do k-mer analysis. There are two peaks in the k-mer frequencies distribution plot. The size of my genome is about 190M, and it is a diploid organism.<br>
  <a href="https://i.stack.imgur.com/y93rV.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/y93rV.png" alt="enter image description here"></a>
  What I am doubting is that can I use this pair end data to polish my assembly?</p>
  
  <p>If not, I have to do sequencing again. </p>
   2085: <p>Resequencing your sample is not going to make it any less heterozygous. I am assuming that you are working with a diploid organism? Please update your question and specify the genome size and ploidy of your organism (the type of organism will help us update our answers, too).</p>
  
  <p>The answer, yes, you can use Pilon + short reads to correct the PacBio assembly even with a highly heterozygous sample. To do so you should use the <code>--diploid</code> option in Pilon to let it know that there will be heterozygous sites. The way that it corrects heterozygous SNPs and heterozygous indels is different, and it seems to perform best with SNPs and sometimes does not correctly identify heterozygous indels (<a href="http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0112963" rel="nofollow noreferrer">source</a>).  </p>
  
  <p>Pilon itself doesn't handle phasing, but just corrects with whatever reads are in the <code>bam</code> file that it is provided with. The end result of providing a <code>bam</code> with reads from both haplotypes is that the corrected sequence will be a mosaic of the dominant allele at the heterozygous site as was observed in the short read data. This may or may not select for some alleles based on sequencing coverage bias of Illumina sequencers.</p>
  