Extract unique reads from SAM file · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Extract unique reads from SAM file

I'm now trying to develop some packages written in R. I've conducted mapping through Rsubread and produce SAM file by the function of asSam in Rsamtools. However, I don't know how to extract the unique mapping reads in the SAM file. In the previous approach, I know people would extract the unique reads by the tag of XT:A:U, but I don't have the tag in the current version of SAM. Any suggestion? Thanks!

Simon

sam unique reads rsamtools

• 1,609 views • link

written 8.3 years ago by niu.shengyong •

Here are some thoughts: BWA mem how to know if a read is mapped uniquely?

• 1 views • link

written 8.3 years ago by cschu181

0 answers

No answers yet.

Log in to answer this question.

More posts like this

Question: how to get low quality read from sam file
written 5.4 years ago by minghuiguo448 •

Hello, l am a starter.I used hisat2 to generate a sam file and have extracted unmapping reads. l want to extract reads with a mapping …
The error generated by Rsamtools
written 6.5 years ago by wangdp123

Hi there, I use Rsubread to produce the BAM files and then use Rsamtools to manipulate the BAM files. A few errors were thrown out: …
Error of SAMSequenceDictionary
written 8.2 years ago by niu.shengyong •

Hi there, after I sort and index my sam file by igvtools, I load my sam file into igv and got the following error: Error …
SAM CIGAR parsing and position counter in R
written 8.3 years ago by niu.shengyong •

I have a SAM file like this: SRR400619.13 0 NC_000913.3 3440402 40 1S25M10S * 0 0 NTGCCGTTGGTTTCCATTTCGATGACNNNNNNNNNN %-:9986986666:5##################### HI:i:1 NH:i:1 NM:i:0 SRR400619.14 0 NC_000913.3 3442282 …
How to unique bwa mem reads?
written 8.3 years ago by niu.shengyong •

I use bwa mem with my paired-ended files and mixed reference genome (concatenate by hg19 and mm10) as the following command: bwa mem -t 4 …
Read metadata from a pacbio bam file
written 9.2 years ago by biomagician

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, …
HISAT2 (v2.0.0-beta) MAPQ query for allele specific expression analysis
written 9.3 years ago by prasundutta87

Hello, For an allele specific expression analysis, one of the ways to reduce ambiguous mapping reads is to extract or use only unique mapping (not …
Subsetting data from BAM file using R tools
written 10.5 years ago by Ilya IN •

Good day to all! This is my first post in Biostars. I have BAM file with Paired-End RNA seq reads, and I need filter this …
How to filter unique reads from sam file?
written 11.4 years ago by Medhat

How to extract unique mapped reads from sam alignment file of pair end reads The mapping was done using bowtie2. There are different posts for …
Bwa: "Xt:A:U" And Mapq Of 0 At The Same Time
written 14.4 years ago by Steffi

<p>I map RNA-Seq Data with <a href='http://bio-bwa.sourceforge.net/'>BWA</a> to the genome. The output files from <a href='http://bio-bwa.sourceforge.net/'>BWA</a> in <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a>-format contain reads that have on the …