Read metadata from a pacbio bam file
Hi,
I have got a BAM file from PacBio data and would like to read some of the metadata of reads in R, e.g. ip, pw, rq. See here for pacbio BAM format: http://pacbiofileformats.readthedocs.io/en/3.0/BAM.html
I have read the BAM into R with Rsamtools::scanBam('file.bam') and I have got a BAM file in R now but I don't know how to extract the ip, pw, rq info.
I have no problem extract sequencing info, mapq, flag and usual BAM tag. It's the metadata that I have a problem with, the additional columns in the BAM file.
I hope someone can help.
Thanks.
C.
• 4,213 views
•
link
written
by
biomagician
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Extract unique reads from SAM file
written by niu.shengyong •
-
Rsamtools usage with URL
written by Friederike
-
Allow ambiguous sequence in bwa index or bwa mem to capture barcode information
written by yesitsjess •
-
Rsamtools and BAM file discordance
written by biomagician
-
How to upload my own BAM file to UCSC genome browsers?
written by sustcyucw •
-
Can I open a BAM file without a header in Rsamtools?
written by biomagician
-
HISAT2 (v2.0.0-beta) MAPQ query for allele specific expression analysis
written by prasundutta87
-
Whole exomes: Rsamtools, RAM issue
written by bioinfo8
-
How to compute MAPQ distribution with Rsamtools
written by Nicola Casiraghi
-
How to read .bam file in Rsamtools R package?
written by M K
You may be better off using SMRTanalysis v.3.x for this purpose.
I've started using SMRTLink in a browser, do you know if it's the same thing? http://www.pacb.com/wp-content/uploads/SMRT-Link-User-Guide-v4.0.0.pdf
Yes. If you want to work with PacBio specific data like what you describe above this would be the best option. It is called SMRTlink now instead of SMRTanalysis (old name).