Answers.Body: <p>I was able to get the hg002 annotation files (<a href="https://s3-us-west-2.amazonaws.com/human-pangenomics/index.html?prefix=T2T/HG002/assemblies/annotation/" rel="nofollow noreferrer">link</a>) hosted on human pangenome S3 bucket (<a href="https://s3-us-west-2.amazonaws.com/human-pangenomics/index.html" rel="nofollow noreferrer">https://s3-us-west-2.amazonaws.com/human-pangenomics/index.html</a>).</p>
  <p>It was due to this article cited below that I was able to figure out the annotation directory.</p>
  <p>Liao, WW., Asri, M., Ebler, J. et al. A draft human pangenome reference. Nature 617, 312–324 (2023). <a href="https://doi.org/10.1038/s41586-023-05896-x" rel="nofollow noreferrer">https://doi.org/10.1038/s41586-023-05896-x</a></p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2024-04-23T11:06:10.397Z
Answers.LastActivityDate: 2024-04-23T11:06:10.397Z
Answers.UserDisplayName: Panda_1996
Answers.Body: <p>If you click through on UCSC you can get to the <a href="https://hgdownload.soe.ucsc.edu/hubs/GCA/018/852/605/GCA_018852605.2/" rel="nofollow noreferrer">UCSC source files for hg002</a>.</p>
  <pre><code>    GCA_018852605.2.2bit                          2023-12-24 09:52  732M  
        GCA_018852605.2.2bit.bpt                      2023-12-24 09:52  4.5K  
        GCA_018852605.2.agp.gz                        2023-12-14 07:12  512   
        GCA_018852605.2.chrom.sizes.txt               2023-12-14 07:12  474   
        GCA_018852605.2.chromAlias.bb                 2024-02-15 12:08   39K  
        GCA_018852605.2.chromAlias.txt                2024-02-15 12:08  806   
        GCA_018852605.2.fa.gz                         2023-12-24 09:52  878M  
        GCA_018852605.2.repeatMasker.out.gz           2023-12-22 12:36  160M  
        GCA_018852605.2.repeatMasker.version.txt      2023-12-22 09:51  1.7K  
        GCA_018852605.2.repeatModeler.2bit            2023-12-25 19:37  729M  
        GCA_018852605.2.repeatModeler.families.fa.gz  2023-12-25 14:49  205K  
        GCA_018852605.2.repeatModeler.families.stk.gz 2023-12-25 14:49  4.6M  
        GCA_018852605.2.repeatModeler.out.gz          2023-12-25 17:38  135M  
        GCA_018852605.2.repeatModeler.version.txt     2023-12-25 14:55  1.4K  
        GCA_018852605.2.rmod.log.txt                  2023-12-25 14:49  5.1K  
        GCA_018852605.2.rmsk.customLib.fa.gz          2023-12-25 14:49  205K  
        GCA_018852605.2.trans.gfidx                   2023-12-24 09:52  4.3G  
        GCA_018852605.2.untrans.gfidx                 2023-12-24 09:52  2.0G  
        GCA_018852605.2_assembly_report.txt           2023-12-23 08:58  2.9K  
        bbi/                                          2024-01-10 14:37    -   
        genes/                                        2024-01-10 14:37    -   
        groups.txt                                    2024-03-09 14:52  542   
        html/                                         2024-03-09 14:54    -   
        hub.txt                                       2024-03-09 14:52  5.3K  
        ixIxx/                                        2024-01-10 14:37    -   
        trackDb.txt                                   2024-01-10 14:37  4.4K  
  </code></pre>
  <p>In the BigBed (<code>bbi/</code>) <a href="https://hgdownload.soe.ucsc.edu/hubs/GCA/018/852/605/GCA_018852605.2/bbi/" rel="nofollow noreferrer">directory</a> you can see several annotation tracks:</p>
  <pre><code>      GCA_018852605.2_ASM1885260v2.allGaps.bb              2023-12-24 08:50   19K  
        GCA_018852605.2_ASM1885260v2.assembly.bb             2023-12-14 07:12   22K  
        GCA_018852605.2_ASM1885260v2.augustus.bb             2023-12-24 19:30  2.0M  
        GCA_018852605.2_ASM1885260v2.cpgIslandExt.bb         2023-12-24 16:56  1.0M  
        GCA_018852605.2_ASM1885260v2.cpgIslandExtUnmasked.bb 2023-12-24 16:55  1.7M  
        GCA_018852605.2_ASM1885260v2.cytoBand.bb             2023-12-14 07:12   21K  
        GCA_018852605.2_ASM1885260v2.gap.bb                  2023-12-14 07:12   24K  
        GCA_018852605.2_ASM1885260v2.gc5Base.bw              2023-12-14 07:12  1.5G  
        GCA_018852605.2_ASM1885260v2.rModel.align.bb         2023-12-25 19:30  1.1G  
        GCA_018852605.2_ASM1885260v2.rModel.bb               2023-12-25 19:23  313M  
        GCA_018852605.2_ASM1885260v2.rmsk.align.bb           2023-12-22 13:59  1.0G  
        GCA_018852605.2_ASM1885260v2.rmsk.bb                 2023-12-22 13:53  320M  
        GCA_018852605.2_ASM1885260v2.simpleRepeat.bb         2023-12-24 08:48   48M  
        GCA_018852605.2_ASM1885260v2.tandemDups.bb           2023-12-14 07:12   50M  
        GCA_018852605.2_ASM1885260v2.windowMasker.bb         2023-12-24 09:51  110M  
        GCA_018852605.2_ASM1885260v2.xenoRefGene.bb          2023-12-24 18:00  3.4M  
  </code></pre>
  <p>I don't see a specific homopolymer bed file but they do have RepeatMasker outputs which you should be able to extract this info from.</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2024-04-18T17:32:02.677Z
Answers.Text: The tandemdups one might be the OP's best bet. UCSC also has a ["segmental duplications"](https://genome.ucsc.edu/cgi-bin/hgTrackUi?hgsid=2138075028_NjmytRhLArrBbu0AcfAEXbM6GXxw&db=hg38&c=chr7&g=genomicSuperDups) track which I don't see in the list here but is probably even better.
Answers.UserDisplayName: terdon
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2024-04-19T06:06:05.067Z
Answers.Text: I was able to work with this [resource](https://s3-us-west-2.amazonaws.com/human-pangenomics/index.html?prefix=T2T/HG002/assemblies/annotation/). It contains homopolymer regions bed file in [microsatellites](https://s3-us-west-2.amazonaws.com/human-pangenomics/index.html?prefix=T2T/HG002/assemblies/annotation/microsatellites/) directory with the name v1.0.1.mnr10.bb. Though there is no explanation for the name 'mnr' but I guess it means mononucleotide repeats.
Answers.UserDisplayName: Panda_1996
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2024-04-19T16:15:14.030Z
Answers.Text: @Panda_1996 if that solved your problem, you should feel free to post an answer to your own question with this information.
Answers.UserDisplayName: Maximilian Press
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2024-04-17T18:12:54.373Z
Answers.LastActivityDate: 2024-04-17T18:12:54.373Z
Answers.UserDisplayName: Maximilian Press
Body: <p><em>This question was also asked on <a href="https://www.biostars.org/p/9592730/" rel="nofollow noreferrer">Biostars</a></em></p>
  <p>I am doing an experiment where I am trying to analyze the errors in the homopolymer regions between the polished reference <a href="https://github.com/marbl/HG002" rel="nofollow noreferrer">hg002</a> genome and hifiasm assembly done inhouse using just long ONT duplex reads.</p>
  <p>I am currently working on Haplotype aware long read error correction mechanisms, so Hg002(diploid) genome and its associated interesting regions (homopolymer, SDs and Tandem repeats) are important for this study. Assemblies only using long ONT duplex or HiFi reads tend to give a lower quality consensus in these regions.</p>
  <p>I am not able to find the .bed files for the hg002 genome. All I was able to find was the homopolymer .bed files for chm13, GRCh38, GRCH37 reference genomes <a href="https://ftp-trace.ncbi.nlm.nih.gov/ReferenceSamples/giab/release/genome-stratifications/v3.3/" rel="nofollow noreferrer">here</a>.</p>
  <p>Thanks a lot!</p>
  
Tags: <assembly><bed><human-genome><repeat-elements>
Title: Where to find the homopolymer regions bed file for Hg002 genome?
ThreadText: # Where to find the homopolymer regions bed file for Hg002 genome?
  
  (Asked by: username_0 on 2024-04-16)
  
  *This question was also asked on [Biostars](https://www.biostars.org/p/9592730/)*
  
  I am doing an experiment where I am trying to analyze the errors in the homopolymer regions between the polished reference [hg002](https://github.com/marbl/HG002) genome and hifiasm assembly done inhouse using just long ONT duplex reads.
  
  I am currently working on Haplotype aware long read error correction mechanisms, so Hg002(diploid) genome and its associated interesting regions (homopolymer, SDs and Tandem repeats) are important for this study. Assemblies only using long ONT duplex or HiFi reads tend to give a lower quality consensus in these regions.
  
  I am not able to find the .bed files for the hg002 genome. All I was able to find was the homopolymer .bed files for chm13, GRCh38, GRCH37 reference genomes [here](https://ftp-trace.ncbi.nlm.nih.gov/ReferenceSamples/giab/release/genome-stratifications/v3.3/).
  
  Thanks a lot!
  
  **Tags:** assembly, bed, human-genome, repeat-elements
  
  
  ---
  
  ## 2 answers
  
  ### Answer by username_0 on 2024-04-23. Score: 0 | Accepted answer
  
  I was able to get the hg002 annotation files ([link](https://s3-us-west-2.amazonaws.com/human-pangenomics/index.html?prefix=T2T/HG002/assemblies/annotation/)) hosted on human pangenome S3 bucket (<https://s3-us-west-2.amazonaws.com/human-pangenomics/index.html>).
  
  It was due to this article cited below that I was able to figure out the annotation directory.
  
  Liao, WW., Asri, M., Ebler, J. et al. A draft human pangenome reference. Nature 617, 312–324 (2023). <https://doi.org/10.1038/s41586-023-05896-x>
  
  
  ---
  
  ### Answer by username_1 on 2024-04-17. Score: 1
  
  If you click through on UCSC you can get to the [UCSC source files for hg002](https://hgdownload.soe.ucsc.edu/hubs/GCA/018/852/605/GCA_018852605.2/).
  
  ```
      GCA_018852605.2.2bit                          2023-12-24 09:52  732M  
        GCA_018852605.2.2bit.bpt                      2023-12-24 09:52  4.5K  
        GCA_018852605.2.agp.gz                        2023-12-14 07:12  512   
        GCA_018852605.2.chrom.sizes.txt               2023-12-14 07:12  474   
        GCA_018852605.2.chromAlias.bb                 2024-02-15 12:08   39K  
        GCA_018852605.2.chromAlias.txt                2024-02-15 12:08  806   
        GCA_018852605.2.fa.gz                         2023-12-24 09:52  878M  
        GCA_018852605.2.repeatMasker.out.gz           2023-12-22 12:36  160M  
        GCA_018852605.2.repeatMasker.version.txt      2023-12-22 09:51  1.7K  
        GCA_018852605.2.repeatModeler.2bit            2023-12-25 19:37  729M  
        GCA_018852605.2.repeatModeler.families.fa.gz  2023-12-25 14:49  205K  
        GCA_018852605.2.repeatModeler.families.stk.gz 2023-12-25 14:49  4.6M  
        GCA_018852605.2.repeatModeler.out.gz          2023-12-25 17:38  135M  
        GCA_018852605.2.repeatModeler.version.txt     2023-12-25 14:55  1.4K  
        GCA_018852605.2.rmod.log.txt                  2023-12-25 14:49  5.1K  
        GCA_018852605.2.rmsk.customLib.fa.gz          2023-12-25 14:49  205K  
        GCA_018852605.2.trans.gfidx                   2023-12-24 09:52  4.3G  
        GCA_018852605.2.untrans.gfidx                 2023-12-24 09:52  2.0G  
        GCA_018852605.2_assembly_report.txt           2023-12-23 08:58  2.9K  
        bbi/                                          2024-01-10 14:37    -   
        genes/                                        2024-01-10 14:37    -   
        groups.txt                                    2024-03-09 14:52  542   
        html/                                         2024-03-09 14:54    -   
        hub.txt                                       2024-03-09 14:52  5.3K  
        ixIxx/                                        2024-01-10 14:37    -   
        trackDb.txt                                   2024-01-10 14:37  4.4K  
  
  ```
  
  In the BigBed (`bbi/`) [directory](https://hgdownload.soe.ucsc.edu/hubs/GCA/018/852/605/GCA_018852605.2/bbi/) you can see several annotation tracks:
  
  ```
        GCA_018852605.2_ASM1885260v2.allGaps.bb              2023-12-24 08:50   19K  
        GCA_018852605.2_ASM1885260v2.assembly.bb             2023-12-14 07:12   22K  
        GCA_018852605.2_ASM1885260v2.augustus.bb             2023-12-24 19:30  2.0M  
        GCA_018852605.2_ASM1885260v2.cpgIslandExt.bb         2023-12-24 16:56  1.0M  
        GCA_018852605.2_ASM1885260v2.cpgIslandExtUnmasked.bb 2023-12-24 16:55  1.7M  
        GCA_018852605.2_ASM1885260v2.cytoBand.bb             2023-12-14 07:12   21K  
        GCA_018852605.2_ASM1885260v2.gap.bb                  2023-12-14 07:12   24K  
        GCA_018852605.2_ASM1885260v2.gc5Base.bw              2023-12-14 07:12  1.5G  
        GCA_018852605.2_ASM1885260v2.rModel.align.bb         2023-12-25 19:30  1.1G  
        GCA_018852605.2_ASM1885260v2.rModel.bb               2023-12-25 19:23  313M  
        GCA_018852605.2_ASM1885260v2.rmsk.align.bb           2023-12-22 13:59  1.0G  
        GCA_018852605.2_ASM1885260v2.rmsk.bb                 2023-12-22 13:53  320M  
        GCA_018852605.2_ASM1885260v2.simpleRepeat.bb         2023-12-24 08:48   48M  
        GCA_018852605.2_ASM1885260v2.tandemDups.bb           2023-12-14 07:12   50M  
        GCA_018852605.2_ASM1885260v2.windowMasker.bb         2023-12-24 09:51  110M  
        GCA_018852605.2_ASM1885260v2.xenoRefGene.bb          2023-12-24 18:00  3.4M  
  
  ```
  
  I don't see a specific homopolymer bed file but they do have RepeatMasker outputs which you should be able to extract this info from.
  
  > **Comments:**
  > **username_2**: The tandemdups one might be the OP's best bet. UCSC also has a ["segmental duplications"](https://genome.ucsc.edu/cgi-bin/hgTrackUi?hgsid=2138075028_NjmytRhLArrBbu0AcfAEXbM6GXxw&db=hg38&c=chr7&g=genomicSuperDups) track which I don't see in the list here but is probably even better.  
  > **username_0**: I was able to work with this [resource](https://s3-us-west-2.amazonaws.com/human-pangenomics/index.html?prefix=T2T/HG002/assemblies/annotation/). It contains homopolymer regions bed file in [microsatellites](https://s3-us-west-2.amazonaws.com/human-pangenomics/index.html?prefix=T2T/HG002/assemblies/annotation/microsatellites/) directory with the name v1.0.1.mnr10.bb. Though there is no explanation for the name 'mnr' but I guess it means mononucleotide repeats.  
  > **username_1**: @username_0 if that solved your problem, you should feel free to post an answer to your own question with this information.  
  
  ---
  
  