Quantcast
Channel: Post Feed
Viewing all articles
Browse latest Browse all 3231

Filter Paired-End Sam File For Xt:A:U

$
0
0

Dear all,

I have a sam file (BWA output, paired-end reads). I would like to retain only reads which are "properly paired". This I would do by:

samtools view -f 0x002 file.sam > file_filtered.sam

Additionally I would like to retain only those pairs of reads where both reads have the XT:A:U tag. It is important to me that after the filterting step I still have the pairs together (so read1, read2, read1, read2, ...).

Any ideas how to do so?

Thanks for any help! Stefanie


Viewing all articles
Browse latest Browse all 3231