pip install lapa

lapa --alignment {rep1.bam},{rep2.bam},{rep3.bam} \
	--fasta {fasta} \
	--annotation {gtf} \
	--chrom_sizes {chrom_sizes} \
	--output_dir {output}

$ lapa --help

Usage: lapa [OPTIONS]

  CLI interface for lapa polyA cluster calling.

Options:
  --alignment TEXT                Single or multiple bam file paths are
                                  separated with a comma.Alternatively, CSV
                                  file with columns of sample, dataset, path
                                  where the sample columns contains the name
                                  of the sample, the dataset is the group of
                                  samples replicates of each other, and path
                                  is the path of bam file.  [required]
  --fasta TEXT                    Genome reference (GENCODE or ENSEMBL fasta)
                                  [required]
  --annotation TEXT               Standart genome annotation (GENCODE or
                                  ENSEMBL gtf). GENCODE gtf file do not
                                  contains annotation for `five_prime_utr` and
                                  `three_prime_utr` so need to be corrected
                                  with `gencode_utr_fix` (see https://github.c
                                  om/MuhammedHasan/gencode_utr_fix.git).
                                  [required]
  --chrom_sizes TEXT              Chrom sizes files (can be generated with
                                  `faidx fasta -i chromsizes > chrom_sizes`)
                                  [required]
  --output_dir TEXT               Output directory of LAPA. See
                                  lapa.readthedocs.io/en/latest/output.html)
                                  for the details of the directory structure
                                  and file format.  [required]
  ...

sample,dataset,path
ENCFF772LYG,myoblast,ENCFF772LYG.bam
ENCFF421MIL,myoblast,ENCFF421MIL.bam
ENCFF699KOR,myotube,ENCFF699KOR.bam
ENCFF731HHB,myotube,ENCFF731HHB.bam

lapa --alignment samples.csv \
	--fasta {fasta} \
	--annotation {gtf} \
	--chrom_sizes {chrom_sizes} \
	--output_dir {output}
 ...

lapa_tss --alignment samples.csv \
	--fasta {fasta} \
	--annotation {gtf} \
	--chrom_sizes {chrom_sizes} \
	--output_dir {output}

$ lapa_tss --help

Usage: lapa_tss [OPTIONS]

  CLI interface for lapa tss cluster calling.

Options:
  --alignment TEXT                Single or multiple bam file paths are
                                  separated with a comma.Alternatively, CSV
                                  file with columns of sample, dataset, path
                                  where the sample columns contains the name
                                  of the sample, the dataset is the group of
                                  samples replicates of each other, and path
                                  is the path of bam file.  [required]
  --fasta TEXT                    Genome reference (GENCODE or ENSEMBL fasta)
                                  [required]
  --annotation TEXT               Standart genome annotation (GENCODE or
                                  ENSEMBL gtf). GENCODE gtf file do not
                                  contains annotation for `five_prime_utr` and
                                  `three_prime_utr` so need to be corrected
                                  with `gencode_utr_fix` (see https://github.c
                                  om/MuhammedHasan/gencode_utr_fix.git)
                                  [required]
  --chrom_sizes TEXT              Chrom sizes files (can be generated
                                  with)`faidx fasta -i chromsizes >
                                  chrom_sizes`)  [required]
  --output_dir TEXT               Output directory of LAPA. See
                                  lapa.readthedocs.io/en/latest/output.html)
                                  for the details of the directory structure
                                  and file format.  [required]

@article{celik2022analysis,
  title={Analysis of alternative polyadenylation from long-read or short-read RNA-seq with LAPA},
  author={Celik, Muhammed Hasan and Mortazavi, Ali},
  journal={bioRxiv},
  year={2022},
  publisher={Cold Spring Harbor Laboratory}
}