Skip to content

Generate one PDF with different plots including read length, read % ID PMD scores, Damage plot and Eveness of coverage along the refence genome

Notifications You must be signed in to change notification settings

BenjaminGuinet/DamagePlots_from_BAM

Repository files navigation

BAM TO DAMAGE PLOTS

Description

This script converts a bam file into a PDF file with multiple plots including reads as well as damage informations.

The plots were hilghy inspired from the ameta outputs (Pochon et al 2024).

You will need Amber, Picard, samtools, PDFtools as well as multiple python packages (see header of the .py script) to be installed.

Usage Example

python3 Make_multiple_plotV2.py -b Yersinia_pestis_control.bam 

Main outputfiles :

  • NoDup_sorted_filterBAM_quality_filtred_Yersinia_pestis_control_dna_damage.csv : CSV file output with many important informations.
  • NoDup_sorted_filterBAM_quality_filtred_Yersinia_pestis_control_dna_damage.pdf : PDF file output with the plots.

Other files correspond to amber and pydamage outputs.


Output PDF example :

View the PDF

About

Generate one PDF with different plots including read length, read % ID PMD scores, Damage plot and Eveness of coverage along the refence genome

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages