How could I include the unmapped reads of a BAM file for further analysis?
0
1
Entering edit mode
L ▴ 50
@dd482574
Last seen 8 months ago
Germany

Dear all,

My goal is to extract from all reads (mapped and unmapped) a part of a sequence that aligns to a specific area in the genome.

I can do this by using this excellent command from Genomic Alignments.

please find the file in here https://1drv.ms/u/s!ArOzUuixE-mggfxdOGO9Hvpp90-qTA?e=ExC5hI

stack <- stackStringsFromBam("alignment.sorted.bam", 
                             param=GRanges("Reference_barcodes:54-78"))

stack

One sequence in the bam file is considered unmapped, but I see it maps well in the designated [54-78]. How can I retrieve the unmapped sequences from the bam file?

GenomicAlignments Biostrings Rsamtools • 923 views
ADD COMMENT

Login before adding your answer.

Traffic: 554 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