QC metrics from pacbio long reads data
Entering edit mode
19 months ago
eli_bayat ▴ 90


I have two samples that were sequenced using pacbio technology. These samples are different in terms of their shearing methods. I also have SMRTLink installed through command line interface.

I want to run some primary analysis to compare these samples, for example coverage bias, alignment summary metrics, and any thing else that can help with understanding differences between these samples.

I have unaligned bam file, pbi and xml files for each sample. I looked into pacbio guid but I couldn't find anything that helps me with these metrics.

Is there any tool I can use within SMRTLink that can do these analysis?


coverage QC Long_Reads smrtlink Pacbio • 762 views
Entering edit mode
19 months ago
gconcepcion ▴ 360

You're looking for the Mapping Analysis application - cromwell.workflows.pb_align_ccs on page 103 of the SMRTLink 10.1 docs you linked.

Running pbcromwell run pb_align_ccs ... will allow you to run the GUI mapping analysis on the command line.

or pbcromwell run pb_align_subreads... if you have subread data.

You can also just use pbmm2 that's included w/ the SMRTlink install to do the mapping manually if you are so inclined. https://github.com/PacificBiosciences/pbmm2


Login before adding your answer.

Traffic: 1073 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6