PacificBiosciences / pbbioconda

PacBio Secondary Analysis Tools on Bioconda. Contains list of PacBio packages available via conda.
BSD 3-Clause Clear License
243 stars 44 forks source link

iso-seq reads and pbfusion #688

Closed mariachiara-github closed 2 months ago

mariachiara-github commented 2 months ago

Hi! Pbfusion uses aligned Iso-Seq reads. Does this mean that I can use flnc.bam reads (generated from the Iso-Seq pipeline), aligne them to hg38 using the pbmm2 tool (pbmm2 index --preset ISOSEQand pbmm2 align --preset ISOSEQ), and then run pbfusion on the aligned flnc reads?

Or is it mandatory to cluster the flnc.bam reads first with the isoseq cluster2tool, then align the clustered reads, and lastly perform the pbfusion?

Thank you for the help! :)

zeeev commented 2 months ago

Hi @meryyr,

I'd recommend that you start with the FLCN reads. The tool also works on the clustered reads, however, we find better performance using the reads. The distribution of reads allows us to capture variability in mapping.

mariachiara-github commented 2 months ago

Hi @meryyr,

I'd recommend that you start with the FLCN reads. The tool also works on the clustered reads, however, we find better performance using the reads. The distribution of reads allows us to capture variability in mapping.

Thank you!!