Answers.Body: <p>to understand how to interpret sanger sequencing chromatograms I'd like to give a brief introduction to how this data are generated (sorry, if you already know this)</p>
  
  <p>The target region is amplified by PCR. The following sequencing reaction is quite similar to a PCR, despite two important differences:</p>
  
  <ol>
  <li>Only one strand is used as a template, by adding just 1 primer instead of two.</li>
  <li>In addition two normal dNTPs, we have modified ddNTPS. Those are fluorescence labeled and the 3'-OH is missing, so there is no chance to add more dNTPS.</li>
  </ol>
  
  <p>The result of this reaction are single stranded DNA fragments of different length. The last base is fluorescence labeled. So we have fragments that looks like this:</p>
  
  <pre><code>A
  AT
  ATC
  ATCG
  ATCGA
  </code></pre>
  
  <p>Of course there a many copies of these fragments.</p>
  
  <p>Now these fragments are separated by capillary electrophoresis. Small fragments are detected first. We can differ the several bases by its fluorescence color.</p>
  
  <p>Let's imagine a DNA strand like this:</p>
  
  <pre><code>ATCGA
  </code></pre>
  
  <p>What happens if we had on the second strand a SNP?</p>
  
  <pre><code>ATCGA
  ATAGA
  </code></pre>
  
  <p>Now the sequence would detect two signal of equal intensity on base 3 - one for C and one for A. This have two effects in the chromatogram.</p>
  
  <ol>
  <li>Two peaks of equal heights with different colors</li>
  <li>The relative peak height is just half the height, compared to those sequences that have just one base at this position.</li>
  </ol>
  
  <p>Let's see what happen in case of a heterozygous deletion. An alignment would look like this</p>
  
  <pre><code>ATCGA
  AT-GA
  </code></pre>
  
  <p>But our capillary electrophoresis doesn't know about alignment. It only knows about fragment length. So this is what it will detect</p>
  
  <pre><code>ATCGA
  ATGA
  </code></pre>
  
  <p>With the beginning of the deletion it will detect to different signals. Just like in case of the SNP, the relativ signal intensity is now just about the half.</p>
  
  <p>With this knowledge let's go back to your initial questions:</p>
  
  <blockquote>
    <p>How can evidence of (het) indels like this be distinguished from noise
    in the chromatogram? My guess is it includes things like</p>
    
    <ul>
    <li>the overlapping peaks last until the end of the sample (and don't start at the beginning, but somewhere in the middle); and</li>
    </ul>
  </blockquote>
  
  <p>They don't start at the beginning, if your indel is not right after the primer, correct. In case of an insertion or deletion you have strands of different size. So if your target region is quite small you may have luck, that you can see the end of the overlap at the end of the sequence, because the shorter strand reaches its end.</p>
  
  <blockquote>
    <ul>
    <li>the peaks are still at the same horizontal positions; i.e., one is not shifted right or left relative to the other?</li>
    </ul>
  </blockquote>
  
  <p>Yes. And one of the best things to distinguish between noise and true heterozygous variant is a relativ drop off of the signal intensities by approximately 50% compared to homozygous samples.</p>
  
  <p>Concerning your second chromatogram: This is not an indel. As the other already said, this is a mistake made by the polymerase. Invitro polymerases doesn't like homopolymer stretches. They begin to add extra bases or "forget" one. As we do PCR this error accumulates. In my experience up to 9 time the same base is not a big problem, 10 only if you have luck, from 11 upwards your sequence are terrible. If you are interested in the sequence part behind this stretch, you have to sequence from the other direction.</p>
  
  <p>fin swimmer</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-17T06:23:43.857Z
Answers.Text: Thanks for your thoroughness! Interesting that you bring up the relative intensity of heterozygous peaks relative to a het peak at the same spot--I'm not quite sure I'll be able to use this, though, because I've only sequenced one sample, so I don't have multiple sequencing results that cover the same locus. 
  
  I'm interested in to know how you're confident that the chromatogram with the homopolymer does not show an indel—surely it's *possible* that there is an indel in a homopolymer region, right?
Answers.UserDisplayName: Randoms
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-17T06:33:58.407Z
Answers.Text: You see no significant drop-off of the average peak height and the additional base peaks are to small. That's why it's unlikely that there is a heterozygous indel. But as I wrote, the homopolymer stretch is to long to be sure. I wouldn't give any definitively answer here.
Answers.UserDisplayName: finswimmer
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-18T07:41:58.733Z
Answers.Text: Downvote without comment? Not a good style :(
Answers.UserDisplayName: finswimmer
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-17T05:22:17.720Z
Answers.LastActivityDate: 2018-12-17T05:22:17.720Z
Answers.UserDisplayName: finswimmer
Answers.Body: <p>There are specialized tools to deconvolute (separate) heterozygous InDels in a chromatogram trace using a reference genome. For instance, <a href="https://www.gear-genomics.com/indigo/" rel="nofollow noreferrer">Indigo</a> is a web application for this, the command-line application is called <a href="https://github.com/gear-genomics/tracy" rel="nofollow noreferrer">Tracy</a>.</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2019-08-01T09:27:52.080Z
Answers.LastActivityDate: 2019-09-24T11:47:37.817Z
Answers.LastEditDate: 2019-09-24T11:47:37.817Z
Answers.UserDisplayName: tora
Answers.Body: <p>Sanger sequencing produces a bit of messy signal at the start and a lot at the end of the sequence, but should have a region of clean signal in the middle. If the signal is getting mixed up in the middle, it suggests that there are sequences that are out of phase with each other. That <em>could</em> be a heterozygous INDEL, but it could also indicate that multiple template sequences have been amplified by the same primer.</p>
  
  <p>Sanger sequencing is frequently carried out in both directions (forward and reverse complement), and if similar things are happening in both directions at the same place, it makes it more likely that there is a heterozygous INDEL and not something else funny going on. It's a very good idea to have homozygous templates for all expected variants, if possible. If the secondary peaks in the signal match what is observed in the homozygous variants, that can also be used as evidence. In the absence of control sequences, it is also possible to eyeball the peaks to find self-similarity.</p>
  
  <p>Where triple peaks are seen, that could point to a region where there are two INDEL variants in close proximity, or multiple regions getting amplified at the same time, or even a single INDEL variant with more than two variant states (e.g. a homopolymer with 8 As that is sometimes 7 As, and sometimes 9 As).</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-17T06:19:31.687Z
Answers.Text: Thanks for your response! Can you expand on what you mean by "positive and negative homozygous variant templates"?
Answers.UserDisplayName: Randoms
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-17T10:38:34.303Z
Answers.Text: Sorry, I was thinking of positive and negative controls; I've changed it to clarify as "homozygous templates for all expected variants"
Answers.UserDisplayName: gringer
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-16T18:00:30.307Z
Answers.LastActivityDate: 2018-12-17T10:37:45.160Z
Answers.LastEditDate: 2018-12-17T10:37:45.160Z
Answers.UserDisplayName: gringer
Answers.Body: <p>Honestly I would not try to call heterozygous indels with Sanger sequencing.  It's just not the right technology.  </p>
  
  <p>If you were trying to confirm a known indel that didn't flank a homopolymer stretch, and you had sequence from the other side that got messy at exactly the same point, maybe you could safely conclude there was something going on there, and possibly you could eyeball the exact right sequence of mixed peaks, but I don't see that what's what you have in your examples at all.  It just looks like noise.</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-18T05:23:30.533Z
Answers.Text: Thanks for the input. By the other side do you mean the other strand—like forward and reverse Sanger sequencing—or the same strand, just on the other side of the variant?
Answers.UserDisplayName: Randoms
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-19T09:39:18.440Z
Answers.Text: For me sanger sequencing is still the gold standard for calling indels :) Compared to short-read sequencing is it quite easy to discover indels of any length.
Answers.UserDisplayName: finswimmer
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-19T17:31:51.167Z
Answers.Text: But mixed indels?  Unambiguously no matter what the surrounding sequence?  If there is an indel in that poly-A region, It would be exceedingly difficult to prove.
Answers.UserDisplayName: swbarnes2
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-17T23:17:26.713Z
Answers.LastActivityDate: 2018-12-17T23:17:26.713Z
Answers.UserDisplayName: swbarnes2
Answers.Body: <p>One important point that @gringer briefly touched on is that there is often off-target amplification when performing PCR on whole-genome DNA.</p>
  
  <p>So, the <strong>noise in the chromatograms very well could be from another part of the genome and not related to indels at all</strong>. </p>
  
  <p>The only way to ensure that this is not a problem is to clone your PCR products into a vector, transform into something like E. coli, grow up a single colony, cut the product with restriction enzymes, size-select, then sequence the correct-sized insert. Or you could just re-sequence your PCR product with Illumina or PacBio to accurately sequence the complete population of molecules -- indels and all. Check out:</p>
  
  <p><a href="https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-018-4611-3" rel="nofollow noreferrer">Hebert, Paul DN, et al. "A Sequel to Sanger: amplicon sequencing that scales." BMC genomics 19.1 (2018): 219.</a></p>
  
  <p>The safe thing to do is to not read too much into the noise of Sanger chromatograms unless it is a simple and very clean SNP.</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-17T20:06:36.150Z
Answers.Text: Before making such a big stunt with cloning, NGS, PacBio etc. one should check the  specificity of the primer, make changes to the PCR protocol, make sure that the clean up steps work, ... If you are working on region that have no homologues, off-target amplification isn't normal.
Answers.UserDisplayName: finswimmer
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-17T19:49:40.823Z
Answers.LastActivityDate: 2018-12-17T20:48:19.580Z
Answers.LastEditDate: 2018-12-17T20:48:19.580Z
Answers.UserDisplayName: conchoecia
Answers.Body: <p>If you want to call heterozygous SNPs and INDELs, you may try <a href="https://droog.gs.washington.edu/polyphred/" rel="nofollow noreferrer">polyphred</a>. I believe there are a couple of other similar tools, but I have forgotten their names.</p>
  
  <blockquote>
    <p>How can evidence of (het) indels like this be distinguished from noise in the chromatogram?</p>
  </blockquote>
  
  <p>Your two points are correct. Basically, at some position, you start to see two streaks of overlapping peaks. Once identified, the two streaks can be aligned to each other if you shift one of them. This is a strong constraint, but it is hard to tell by eye because it is not obvious how to assign peaks to streaks without running an automatic tool. Doing this is even more difficult given noisy traces like your last figure. I think you'd better run polyphred or similar.</p>
  
  <blockquote>
    <p>the putative indel seems to start in a homopolymer region—cast doubt on this?</p>
  </blockquote>
  
  <p>The INDEL mutation rate in long homopolymer runs is orders of magnitude higher than elsewhere in the genome. This has two implications: 1) you expect to see more homopolymer indels; 2) PCR polymerase is more likely to introduce errors in homopolymers––a higher fraction of homopolymer indels tend to be wrong.</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-17T06:14:54.140Z
Answers.Text: Thanks! I'll look into polyphred and other tools. To be clear, what's your take on the indel in the homopolymer region—that the chromatogram is inconclusive itself? What can I do to verify it?
Answers.UserDisplayName: Randoms
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-17T13:03:11.903Z
Answers.Text: @Randoms Overlapping peaks alone are not strong evidence. You need to check if the overlapping peaks can be separated into two sets and the two sets can be aligned with a shift. I can't do this in my head. Even if there is an indel there, the chance of seeing a PCR error is still high. Run polyphred first. PS: finswimmer is right that overlapping peaks should be twice as "low". This makes the second figure less likely to be true.
Answers.UserDisplayName: user172818
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-18T21:23:05.003Z
Answers.Text: I think you're right that running a tool is really the way to go. Thanks! I'm going to try Poly Peak Parser: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4525701/
Answers.UserDisplayName: Randoms
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-12-16T14:56:18.683Z
Answers.LastActivityDate: 2018-12-16T14:56:18.683Z
Answers.UserDisplayName: user172818
Body: <p>I have around 100 chromatograms (<code>.ab1</code> files) from Sanger sequencing a genome at loci believed to have an indel. </p>
  
  <p>I'm new to interpreting this kind of data in general, but I've read a bit on the general idea—mostly in guides like <a href="https://seqcore.brcf.med.umich.edu/sites/default/files/html/interpret.html" rel="noreferrer">this</a>. What I'm not able to figure out from the resources I've checked is <strong>how to identify insertions and deletions in a chromatogram</strong>.</p>
  
  <p>I'm most concerned with  heterozygous indels, and these seem to be simpler. (And since my samples are (supposed to be) larger indels that are more rare and more likely to be het, I think.) A heterozygous indel, according to my understanding, would like something like consecutive heterozygous SNPs—overlapping peaks, maybe not with the same height, but with the same horizontal placement—until the end of the sample. </p>
  
  <p>In other words, something like the top track in the open window in this image. </p>
  
  <p><a href="https://i.sstatic.net/BFkKb.png" rel="noreferrer"><img src="https://i.sstatic.net/BFkKb.png" alt="Indel Chromatogram"></a>
  (Source: <a href="https://www.codoncode.com/aligner/mutation-detection.htm" rel="noreferrer">CodonCode</a>.)</p>
  
  <p>(I'm trying to do this manually, and not with packages like this for several reasons, including: I want to see how this is actually done, I want to do it as thoroughly as possible—I will probably check my results with an automated package—and I was asked to do it like this.)</p>
  
  <p><em>How can evidence of (het) indels like this be distinguished from noise in the chromatogram?</em> My guess is it includes things like</p>
  
  <ul>
  <li>the overlapping peaks last until the end of the sample (and don't start at the beginning, but somewhere in the middle); and </li>
  <li>the peaks are still at the same horizontal positions; i.e., one is not shifted right or left relative to the other?</li>
  </ul>
  
  <p>Is this correct?</p>
  
  <hr>
  
  <p>With my own chromatograms, I have a lot of doubt about distinguishing noise from indel. Here's an example, visualized with 4Peaks. (Full chromatogram <a href="http://s000.tinyupload.com/?file_id=41522383324429952799" rel="noreferrer">here</a>.)</p>
  
  <p><a href="https://i.sstatic.net/Zii0P.png" rel="noreferrer"><img src="https://i.sstatic.net/Zii0P.png" alt="enter image description here"></a></p>
  
  <p>On the left, the chromatogram appears high quality. There are sharp, regularly spaced, single peaks. On the right, there appear to be two sets of peaks lined up with each other. Does this constitute evidence for an indel? Does the homopolymer region—I've noticed in many of these samples, the putative indel seems to start in a homopolymer region—cast doubt on this? </p>
  
  <hr>
  
  <p>I've also noticed that in some situations like this, there appear to be not two but <em>three</em> peaks. </p>
  
  <p><a href="https://i.sstatic.net/5efEm.png" rel="noreferrer"><img src="https://i.sstatic.net/5efEm.png" alt="enter image description here"></a></p>
  
  <p>Is one of these just noise and can be disregarded? Or does this cast doubt on the validity of the data—which is supposed to be coming from the diploid human genome?</p>
  
Tags: <phylogenetics><ngs><genome><variant-calling>
Title: Identifying Indels from Chromatograms
ThreadText: # Identifying Indels from Chromatograms
  
  (Asked by: username_0 on 2018-12-16)
  
  I have around 100 chromatograms (`.ab1` files) from Sanger sequencing a genome at loci believed to have an indel.
  
  I'm new to interpreting this kind of data in general, but I've read a bit on the general idea—mostly in guides like [this](https://seqcore.brcf.med.umich.edu/sites/default/files/html/interpret.html). What I'm not able to figure out from the resources I've checked is **how to identify insertions and deletions in a chromatogram**.
  
  I'm most concerned with heterozygous indels, and these seem to be simpler. (And since my samples are (supposed to be) larger indels that are more rare and more likely to be het, I think.) A heterozygous indel, according to my understanding, would like something like consecutive heterozygous SNPs—overlapping peaks, maybe not with the same height, but with the same horizontal placement—until the end of the sample.
  
  In other words, something like the top track in the open window in this image.
  
  [![Indel Chromatogram](https://i.sstatic.net/BFkKb.png)](https://i.sstatic.net/BFkKb.png)
  (Source: [CodonCode](https://www.codoncode.com/aligner/mutation-detection.htm).)
  
  (I'm trying to do this manually, and not with packages like this for several reasons, including: I want to see how this is actually done, I want to do it as thoroughly as possible—I will probably check my results with an automated package—and I was asked to do it like this.)
  
  *How can evidence of (het) indels like this be distinguished from noise in the chromatogram?* My guess is it includes things like
  
  * the overlapping peaks last until the end of the sample (and don't start at the beginning, but somewhere in the middle); and
  * the peaks are still at the same horizontal positions; i.e., one is not shifted right or left relative to the other?
  
  Is this correct?
  
  ---
  
  With my own chromatograms, I have a lot of doubt about distinguishing noise from indel. Here's an example, visualized with 4Peaks. (Full chromatogram [here](http://s000.tinyupload.com/?file_id=41522383324429952799).)
  
  [![enter image description here](https://i.sstatic.net/Zii0P.png)](https://i.sstatic.net/Zii0P.png)
  
  On the left, the chromatogram appears high quality. There are sharp, regularly spaced, single peaks. On the right, there appear to be two sets of peaks lined up with each other. Does this constitute evidence for an indel? Does the homopolymer region—I've noticed in many of these samples, the putative indel seems to start in a homopolymer region—cast doubt on this?
  
  ---
  
  I've also noticed that in some situations like this, there appear to be not two but *three* peaks.
  
  [![enter image description here](https://i.sstatic.net/5efEm.png)](https://i.sstatic.net/5efEm.png)
  
  Is one of these just noise and can be disregarded? Or does this cast doubt on the validity of the data—which is supposed to be coming from the diploid human genome?
  
  **Tags:** phylogenetics, ngs, genome, variant-calling
  
  
  ---
  
  ## 6 answers
  
  ### Answer by username_6 on 2018-12-16. Score: 4 | Accepted answer
  
  If you want to call heterozygous SNPs and INDELs, you may try [polyphred](https://droog.gs.washington.edu/polyphred/). I believe there are a couple of other similar tools, but I have forgotten their names.
  
  > How can evidence of (het) indels like this be distinguished from noise in the chromatogram?
  
  Your two points are correct. Basically, at some position, you start to see two streaks of overlapping peaks. Once identified, the two streaks can be aligned to each other if you shift one of them. This is a strong constraint, but it is hard to tell by eye because it is not obvious how to assign peaks to streaks without running an automatic tool. Doing this is even more difficult given noisy traces like your last figure. I think you'd better run polyphred or similar.
  
  > the putative indel seems to start in a homopolymer region—cast doubt on this?
  
  The INDEL mutation rate in long homopolymer runs is orders of magnitude higher than elsewhere in the genome. This has two implications: 1) you expect to see more homopolymer indels; 2) PCR polymerase is more likely to introduce errors in homopolymers––a higher fraction of homopolymer indels tend to be wrong.
  
  > **Comments:**
  > **username_0**: Thanks! I'll look into polyphred and other tools. To be clear, what's your take on the indel in the homopolymer region—that the chromatogram is inconclusive itself? What can I do to verify it?  
  > **username_6**: @username_0 Overlapping peaks alone are not strong evidence. You need to check if the overlapping peaks can be separated into two sets and the two sets can be aligned with a shift. I can't do this in my head. Even if there is an indel there, the chance of seeing a PCR error is still high. Run polyphred first. PS: username_1 is right that overlapping peaks should be twice as "low". This makes the second figure less likely to be true.  
  > **username_0**: I think you're right that running a tool is really the way to go. Thanks! I'm going to try Poly Peak Parser: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4525701/  
  
  ---
  
  ### Answer by username_1 on 2018-12-17. Score: 3
  
  to understand how to interpret sanger sequencing chromatograms I'd like to give a brief introduction to how this data are generated (sorry, if you already know this)
  
  The target region is amplified by PCR. The following sequencing reaction is quite similar to a PCR, despite two important differences:
  
  1. Only one strand is used as a template, by adding just 1 primer instead of two.
  2. In addition two normal dNTPs, we have modified ddNTPS. Those are fluorescence labeled and the 3'-OH is missing, so there is no chance to add more dNTPS.
  
  The result of this reaction are single stranded DNA fragments of different length. The last base is fluorescence labeled. So we have fragments that looks like this:
  
  ```
  A
  AT
  ATC
  ATCG
  ATCGA
  
  ```
  
  Of course there a many copies of these fragments.
  
  Now these fragments are separated by capillary electrophoresis. Small fragments are detected first. We can differ the several bases by its fluorescence color.
  
  Let's imagine a DNA strand like this:
  
  ```
  ATCGA
  
  ```
  
  What happens if we had on the second strand a SNP?
  
  ```
  ATCGA
  ATAGA
  
  ```
  
  Now the sequence would detect two signal of equal intensity on base 3 - one for C and one for A. This have two effects in the chromatogram.
  
  1. Two peaks of equal heights with different colors
  2. The relative peak height is just half the height, compared to those sequences that have just one base at this position.
  
  Let's see what happen in case of a heterozygous deletion. An alignment would look like this
  
  ```
  ATCGA
  AT-GA
  
  ```
  
  But our capillary electrophoresis doesn't know about alignment. It only knows about fragment length. So this is what it will detect
  
  ```
  ATCGA
  ATGA
  
  ```
  
  With the beginning of the deletion it will detect to different signals. Just like in case of the SNP, the relativ signal intensity is now just about the half.
  
  With this knowledge let's go back to your initial questions:
  
  > How can evidence of (het) indels like this be distinguished from noise
  > in the chromatogram? My guess is it includes things like
  >
  > * the overlapping peaks last until the end of the sample (and don't start at the beginning, but somewhere in the middle); and
  
  They don't start at the beginning, if your indel is not right after the primer, correct. In case of an insertion or deletion you have strands of different size. So if your target region is quite small you may have luck, that you can see the end of the overlap at the end of the sequence, because the shorter strand reaches its end.
  
  > * the peaks are still at the same horizontal positions; i.e., one is not shifted right or left relative to the other?
  
  Yes. And one of the best things to distinguish between noise and true heterozygous variant is a relativ drop off of the signal intensities by approximately 50% compared to homozygous samples.
  
  Concerning your second chromatogram: This is not an indel. As the other already said, this is a mistake made by the polymerase. Invitro polymerases doesn't like homopolymer stretches. They begin to add extra bases or "forget" one. As we do PCR this error accumulates. In my experience up to 9 time the same base is not a big problem, 10 only if you have luck, from 11 upwards your sequence are terrible. If you are interested in the sequence part behind this stretch, you have to sequence from the other direction.
  
  fin swimmer
  
  > **Comments:**
  > **username_0**: Thanks for your thoroughness! Interesting that you bring up the relative intensity of heterozygous peaks relative to a het peak at the same spot--I'm not quite sure I'll be able to use this, though, because I've only sequenced one sample, so I don't have multiple sequencing results that cover the same locus.   I'm interested in to know how you're confident that the chromatogram with the homopolymer does not show an indel—surely it's *possible* that there is an indel in a homopolymer region, right?  
  > **username_1**: You see no significant drop-off of the average peak height and the additional base peaks are to small. That's why it's unlikely that there is a heterozygous indel. But as I wrote, the homopolymer stretch is to long to be sure. I wouldn't give any definitively answer here.  
  > **username_1**: Downvote without comment? Not a good style :(  
  
  ---
  
  ### Answer by username_2 on 2019-08-01. Score: 1
  
  There are specialized tools to deconvolute (separate) heterozygous InDels in a chromatogram trace using a reference genome. For instance, [Indigo](https://www.gear-genomics.com/indigo/) is a web application for this, the command-line application is called [Tracy](https://github.com/gear-genomics/tracy).
  
  
  ---
  
  ### Answer by username_3 on 2018-12-16. Score: 1
  
  Sanger sequencing produces a bit of messy signal at the start and a lot at the end of the sequence, but should have a region of clean signal in the middle. If the signal is getting mixed up in the middle, it suggests that there are sequences that are out of phase with each other. That *could* be a heterozygous INDEL, but it could also indicate that multiple template sequences have been amplified by the same primer.
  
  Sanger sequencing is frequently carried out in both directions (forward and reverse complement), and if similar things are happening in both directions at the same place, it makes it more likely that there is a heterozygous INDEL and not something else funny going on. It's a very good idea to have homozygous templates for all expected variants, if possible. If the secondary peaks in the signal match what is observed in the homozygous variants, that can also be used as evidence. In the absence of control sequences, it is also possible to eyeball the peaks to find self-similarity.
  
  Where triple peaks are seen, that could point to a region where there are two INDEL variants in close proximity, or multiple regions getting amplified at the same time, or even a single INDEL variant with more than two variant states (e.g. a homopolymer with 8 As that is sometimes 7 As, and sometimes 9 As).
  
  > **Comments:**
  > **username_0**: Thanks for your response! Can you expand on what you mean by "positive and negative homozygous variant templates"?  
  > **username_3**: Sorry, I was thinking of positive and negative controls; I've changed it to clarify as "homozygous templates for all expected variants"  
  
  ---
  
  ### Answer by username_4 on 2018-12-17. Score: 1
  
  Honestly I would not try to call heterozygous indels with Sanger sequencing. It's just not the right technology.
  
  If you were trying to confirm a known indel that didn't flank a homopolymer stretch, and you had sequence from the other side that got messy at exactly the same point, maybe you could safely conclude there was something going on there, and possibly you could eyeball the exact right sequence of mixed peaks, but I don't see that what's what you have in your examples at all. It just looks like noise.
  
  > **Comments:**
  > **username_0**: Thanks for the input. By the other side do you mean the other strand—like forward and reverse Sanger sequencing—or the same strand, just on the other side of the variant?  
  > **username_1**: For me sanger sequencing is still the gold standard for calling indels :) Compared to short-read sequencing is it quite easy to discover indels of any length.  
  > **username_4**: But mixed indels?  Unambiguously no matter what the surrounding sequence?  If there is an indel in that poly-A region, It would be exceedingly difficult to prove.  
  
  ---
  
  ### Answer by username_5 on 2018-12-17. Score: 0
  
  One important point that @username_3 briefly touched on is that there is often off-target amplification when performing PCR on whole-genome DNA.
  
  So, the **noise in the chromatograms very well could be from another part of the genome and not related to indels at all**.
  
  The only way to ensure that this is not a problem is to clone your PCR products into a vector, transform into something like E. coli, grow up a single colony, cut the product with restriction enzymes, size-select, then sequence the correct-sized insert. Or you could just re-sequence your PCR product with Illumina or PacBio to accurately sequence the complete population of molecules -- indels and all. Check out:
  
  [Hebert, Paul DN, et al. "A Sequel to Sanger: amplicon sequencing that scales." BMC genomics 19.1 (2018): 219.](https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-018-4611-3)
  
  The safe thing to do is to not read too much into the noise of Sanger chromatograms unless it is a simple and very clean SNP.
  
  > **Comments:**
  > **username_1**: Before making such a big stunt with cloning, NGS, PacBio etc. one should check the  specificity of the primer, make changes to the PCR protocol, make sure that the clean up steps work, ... If you are working on region that have no homologues, off-target amplification isn't normal.  
  
  ---
  
  