Answers.Body: <p>Genome-In-A-Bottle (GIAB; version 3.3.2) contains 3.21M SNPs on auto+X chromosomes and 0.51M INDELs in 2.58Gb confident regions. The ins:del ratio is 0.92. On the CHM1-CHM13 pacbio assembly (European ancestry), there are 3.57M auto+X SNPs and 0.58M INDELs in 2.71Gb confident regions [<a href="https://rdcu.be/232R" rel="nofollow noreferrer">reference</a>] with an ins:del ratio 0.99. These give you an idea of relative short INDELs.</p>
  
  <p>For indels longer than 50bp, I would highly recommend to read <a href="https://www.nature.com/nature/journal/v517/n7536/full/nature13907.html" rel="nofollow noreferrer">this paper</a> by Chaisson et al. This call set is constructed from the whole-genome assembly of the CHM1 genome. It is far more comprehensive and probably more accurate than all the other call sets. The following is Table 1 from the paper:</p>
  
  <p><a href="https://i.sstatic.net/om410.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/om410.png" alt="CHM1-SVs"></a></p>
  
  <p>I have to say that the high ins-to-del ratio is surprising to me, but this is what the data tells us. I have done a similar analysis and reached a similar ratio.</p>
  
  <p>PS: These numbers are essentially from one haplotype. They will be higher for a diploid sample.</p>
  
  <p>EDIT: some numbers in the original post are wrong.</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-07-19T03:28:59.020Z
Answers.Text: "GIAB ignores hard regions that tend to harbor more indels" -- that's interesting. Do you have a source for that or could you explain more?
Answers.UserDisplayName: Randoms
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-07-20T01:35:16.817Z
Answers.Text: @Randoms Sorry. Some of my numbers are wrong. Please see the update. As to your question, GIAB excludes long STRs, while half of indels come from STRs. You can see that "(0.58-0.51)/(2.71-2.58)=0.54" indels per kb outside the GIAB confident region is much higher than "0.51/2.58=0.20" inside its confident region.
Answers.UserDisplayName: user172818
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-07-20T04:04:00.540Z
Answers.Text: Ok thanks. Where are those numbers in that coming from?
Answers.UserDisplayName: Randoms
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-07-20T15:56:58.377Z
Answers.Text: @Randoms These VCFs are public. You can count by yourself.
Answers.UserDisplayName: user172818
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-07-23T04:32:28.867Z
Answers.Text: Could you explain where the numbers for the CHM1-CHM13 dataset are coming from? The Syndip paper you linked refers to "0.38 million 2–50bp INDELs". The `rep2.37.broad.hc.raw.vcf.gz` VCF, by my count, has 936587 indels—493900 are within [2, 50]. The `full.37m.vcf.gz` VCF has 1145454 indels—556256 in [2, 50].
Answers.UserDisplayName: Randoms
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-07-23T04:52:27.670Z
Answers.Text: Intersecting the VCFs with the `full.37m.bed.gz`, which I guess I should have been doing, yields lower estimates—though still not what quite you're quoting, or what's in the paper strangely.
Answers.UserDisplayName: Randoms
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2018-07-24T01:36:39.603Z
Answers.Text: @Randoms There is a lot to counting: the version of data and manuscript, and the way you take intersection and count. Comments are not the best place for discussions.
Answers.UserDisplayName: user172818
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2017-06-20T02:17:29.127Z
Answers.LastActivityDate: 2018-07-20T01:28:43.777Z
Answers.LastEditDate: 2018-07-20T01:28:43.777Z
Answers.UserDisplayName: user172818
Answers.Body: <p>One of the <a href="https://www.nature.com/nature/journal/v526/n7571/full/nature15394.html" rel="noreferrer">2015 papers from the 1000 genomes project</a> has a nice figure (figure 1) showing the size distribution of medium to large sized insertions and deletions:</p>
  
  <p><a href="https://i.sstatic.net/bFMCX.jpg" rel="noreferrer"><img src="https://i.sstatic.net/bFMCX.jpg" alt="enter image description here"></a></p>
  
  <p>From <a href="https://www.nature.com/nature/journal/v526/n7571/full/nature15393.html" rel="noreferrer">another 2015 1000 genomes paper</a>, one can see that the absolute number of smaller indels is much larger, though an exact size range isn't given (as far as I saw). If you really want to know that, just download <a href="ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/" rel="noreferrer">the most recent 1000 genomes VCF file(s)</a> and compute the exact median size and/or distribution that you want.</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-06-19T21:43:05.247Z
Answers.LastActivityDate: 2017-06-19T21:43:05.247Z
Answers.UserDisplayName: Devon Ryan
Body: <p>My understanding is that indels are from 1bp to 10Kb, and a healthy genome has ~400K-500K Indels. Surely most of these are small. </p>
  
  <p>What is the distribution of insertion sizes in a healthy human genome? What is the distribution of deletion sizes? What is an average ratio of insertions:deletions? </p>
  
  <p>I have not found a definitive reference.</p>
  
Tags: <variation><indel><literature-search>
Title: What is the distribution of indel sizes in a healthy human genome? of insertion:deletion ratios?
ThreadText: # What is the distribution of indel sizes in a healthy human genome? of insertion:deletion ratios?
  
  (Asked by: username_0 on 2017-06-19)
  
  My understanding is that indels are from 1bp to 10Kb, and a healthy genome has ~400K-500K Indels. Surely most of these are small.
  
  What is the distribution of insertion sizes in a healthy human genome? What is the distribution of deletion sizes? What is an average ratio of insertions:deletions?
  
  I have not found a definitive reference.
  
  **Tags:** variation, indel, literature-search
  
  
  ---
  
  ## 2 answers
  
  ### Answer by username_3 on 2017-06-19. Score: 6 | Accepted answer
  
  One of the [2015 papers from the 1000 genomes project](https://www.nature.com/nature/journal/v526/n7571/full/nature15394.html) has a nice figure (figure 1) showing the size distribution of medium to large sized insertions and deletions:
  
  [![enter image description here](https://i.sstatic.net/bFMCX.jpg)](https://i.sstatic.net/bFMCX.jpg)
  
  From [another 2015 1000 genomes paper](https://www.nature.com/nature/journal/v526/n7571/full/nature15393.html), one can see that the absolute number of smaller indels is much larger, though an exact size range isn't given (as far as I saw). If you really want to know that, just download [the most recent 1000 genomes VCF file(s)](ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/) and compute the exact median size and/or distribution that you want.
  
  
  ---
  
  ### Answer by username_1 on 2017-06-20. Score: 5
  
  Genome-In-A-Bottle (GIAB; version 3.3.2) contains 3.21M SNPs on auto+X chromosomes and 0.51M INDELs in 2.58Gb confident regions. The ins:del ratio is 0.92. On the CHM1-CHM13 pacbio assembly (European ancestry), there are 3.57M auto+X SNPs and 0.58M INDELs in 2.71Gb confident regions [[reference](https://rdcu.be/232R)] with an ins:del ratio 0.99. These give you an idea of relative short INDELs.
  
  For indels longer than 50bp, I would highly recommend to read [this paper](https://www.nature.com/nature/journal/v517/n7536/full/nature13907.html) by Chaisson et al. This call set is constructed from the whole-genome assembly of the CHM1 genome. It is far more comprehensive and probably more accurate than all the other call sets. The following is Table 1 from the paper:
  
  [![CHM1-SVs](https://i.sstatic.net/om410.png)](https://i.sstatic.net/om410.png)
  
  I have to say that the high ins-to-del ratio is surprising to me, but this is what the data tells us. I have done a similar analysis and reached a similar ratio.
  
  PS: These numbers are essentially from one haplotype. They will be higher for a diploid sample.
  
  EDIT: some numbers in the original post are wrong.
  
  > **Comments:**
  > **username_2**: "GIAB ignores hard regions that tend to harbor more indels" -- that's interesting. Do you have a source for that or could you explain more?  
  > **username_1**: @username_2 Sorry. Some of my numbers are wrong. Please see the update. As to your question, GIAB excludes long STRs, while half of indels come from STRs. You can see that "(0.58-0.51)/(2.71-2.58)=0.54" indels per kb outside the GIAB confident region is much higher than "0.51/2.58=0.20" inside its confident region.  
  > **username_2**: Ok thanks. Where are those numbers in that coming from?  
  > **username_1**: @username_2 These VCFs are public. You can count by yourself.  
  > **username_2**: Could you explain where the numbers for the CHM1-CHM13 dataset are coming from? The Syndip paper you linked refers to "0.38 million 2–50bp INDELs". The `rep2.37.broad.hc.raw.vcf.gz` VCF, by my count, has 936587 indels—493900 are within [2, 50]. The `full.37m.vcf.gz` VCF has 1145454 indels—556256 in [2, 50].  
  > **username_2**: Intersecting the VCFs with the `full.37m.bed.gz`, which I guess I should have been doing, yields lower estimates—though still not what quite you're quoting, or what's in the paper strangely.  
  > **username_1**: @username_2 There is a lot to counting: the version of data and manuscript, and the way you take intersection and count. Comments are not the best place for discussions.  
  
  ---
  
  