Guardat en:
| Autor principal: | |
|---|---|
| Format: | Recurso digital |
| Idioma: | |
| Publicat: |
Zenodo
2022
|
| Accés en línia: | https://doi.org/10.5281/zenodo.16333929 |
| Etiquetes: |
Afegir etiqueta
Sense etiquetes, Sigues el primer a etiquetar aquest registre!
|
Taula de continguts:
- <h1>seq-format-conversion</h1> <p>Workflows for converting between sequence data formats</p> <h3>cram-to-bam :</h3> <p>This script should convert a CRAM to SAM to BAM and outputs a BAM, BAM Index, and validation report. The reason this approach was chosen instead of converting CRAM to BAM directly using Samtools is because Samtools 1.3 produces incorrect bins due to an old version of htslib included in the package. Samtools versions 1.4 & 1.5 have an NM issue that causes them to not validate with Picard.</p> <h4>Requirements/expectations</h4> <ul> <li>Cram file</li> </ul> <h4>Outputs</h4> <ul> <li>Bam file and index</li> <li>Validation report</li> </ul> <h3>paired-fastq-to-unmapped-bam :</h3> <p>This WDL converts paired FASTQ to uBAM and adds read group information</p> <h4>Requirements/expectations</h4> <ul> <li>Pair-end sequencing data in FASTQ format (one file per orientation)</li> <li>The following metada descriptors per sample: <ul> <li>readgroup</li> <li>sample_name</li> <li>library_name</li> <li>platform_unit</li> <li>run_date</li> <li>platform_name</li> <li>sequecing_center</li> </ul> </li> </ul> <h4>Outputs</h4> <ul> <li>Unmapped BAM</li> </ul> <h3>bam-to-unmapped-bams :</h3> <p>This WDL converts BAM to unmapped BAMs</p> <h4>Requirements/expectations</h4> <ul> <li>BAM file</li> </ul> <h4>Outputs</h4> <ul> <li>Sorted Unmapped BAMs</li> <li>Text file listing the unmapped file paths (FOFN)</li> </ul> <h3>interleaved-fastq-to-paired-fastq :</h3> <p>This WDL takes in a single interleaved(R1+R2) FASTQ file and separates it into separate R1 and R2 FASTQ (i.e. paired FASTQ) files. Paired FASTQ files are the input format for the tool that generates unmapped BAMs (the format used in most GATK processing and analysis tools).</p> <h4>Requirements/expectations</h4> <ul> <li>Interleaved Fastq file</li> </ul> <h4>Outputs</h4> <ul> <li>Separate R1 and R2 FASTQ files (i.e. paired FASTQ)</li> </ul> <h3>interleaved-fastq-to-paired-fastq :</h3> <p>This WDL takes in a single interleaved(R1+R2) FASTQ file and separates it into separate R1 and R2 FASTQ (i.e. paired FASTQ) files. Paired FASTQ files are the input format for the tool that generates unmapped BAMs (the format used in most GATK processing and analysis tools).</p> <h4>Requirements/expectations</h4> <ul> <li>Interleaved Fastq file</li> </ul> <h4>Outputs</h4> <ul> <li>Separate R1 and R2 FASTQ files (i.e. paired FASTQ)</li> </ul> <h3>Software version requirements :</h3> <ul> <li>GATK4 or later</li> <li>Samtools 1.3.1</li> <li>Picard 2.8.3</li> <li>Cromwell version support <ul> <li>Successfully tested on v50</li> </ul> </li> </ul> <h3>Important Notes :</h3> <ul> <li>Runtime parameters are optimized for Broad's Google Cloud Platform implementation.</li> <li>The provided JSON is a ready to use example JSON template of the workflow. Users are responsible for reviewing the <a href="https://gatk.broadinstitute.org/hc/en-us/categories/360002310591">GATK Tool and Tutorial Documentations</a> to properly set the reference and resource variables.</li> <li>For help running workflows on the Google Cloud Platform or locally please view the following tutorial <a href="https://gatk.broadinstitute.org/hc/en-us/articles/360035530952">(How to) Execute Workflows from the gatk-workflows Git Organization</a>.</li> <li>Relevant reference and resources bundles can be accessed in <a href="https://gatk.broadinstitute.org/hc/en-us/articles/360035890811">Resource Bundle</a>.</li> </ul> <h3>Contact Us :</h3> <ul> <li>The following material is provided by the Data Science Platforum group at the Broad Institute. Please direct any questions or concerns to one of our forum sites : <a href="https://gatk.broadinstitute.org/hc/en-us/community/topics">GATK</a> or <a href="https://support.terra.bio/hc/en-us/community/topics/360000500432">Terra</a>.</li> </ul> <h3>LICENSING :</h3> <p>Copyright Broad Institute, 2019 | BSD-3 This script is released under the WDL open source code license (BSD-3) (full license text at https://github.com/openwdl/wdl/blob/master/LICENSE). Note however that the programs it calls may be subject to different licenses. Users are responsible for checking that they are authorized to run all programs before running this script.</p>