QC metrics from pacbio long reads data
1
0
Entering edit mode
7 weeks ago
eli_bayat ▴ 50

Hello,

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?

Thanks

coverage QC Long_Reads smrtlink Pacbio • 249 views
ADD COMMENT
1
Entering edit mode
6 weeks ago
gconcepcion ▴ 140

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

ADD COMMENT

Login before adding your answer.

Traffic: 1622 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6