Skip to content

Workflow for collapsing of standard binary alignment mapping file generated by MSK-ACCESS

License

Notifications You must be signed in to change notification settings

msk-access/bam_collapsing

Repository files navigation

description
Collapsing of reads a aligned to the genome having same unique molecular identifier.

BAM Collapsing

PyPI version Build Status Updates Python 3

The bam collapsing workflow for MSK-ACCESS assay takes in a standard processed bam from standard processing workflow available here.

Workflow

Features

Given a uncollapsed Unique Molecular Index (UMI) tagged BAM file:

  1. Group Reads using UMI and generate a BAM file
  2. Call Duplex Consensus Reads from the BAM and generate an unmapped BAM file
  3. Convert unmapped BAM to FASTQ and do standard BAM alignment.
  4. Merge the unmapped BAM with the aligned BAM and this generated the collapsed BAM.
  5. Generate Duplex and Simplex BAM files from the collapsed BAM.

Credits

  • CMO cfDNA Informatics Team
  • This package was created with Cookiecutter and the audreyr/cookiecutter-pypackage project template.