In the process of data analysis, it is always necessary to customize scripts to accomplish specific tasks.
This program processes paired-end bam file sorted by read name, calculates the fragment length for a pair of reads aligned to the same chromosome, and visualizes it.
- python3
- numpy
- matplotlib
- HTSeq
- logging
- argparse
python Fragment_length.py -h