Separate BAM reads by the number and type of mismatchs on the read
I'm interesting in analyzing pair-end SLAM-seq data in a custom manner.
I want to split an aligned RNA-seq into 2 files, every read with more than one mismatch of a given type (every T>C mutation) into one sam/bam, and the others into another sam/bam.
I'm sure I could do it with pysamtools, but was wondering if there is an easier way by directly parsing the MD tag. Could probably get partially there with a quick filter on the NM tag, but I also want specific types of mismatches.
• 674 views
•
link
written
by
Anna-Leigh •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to discard alignments from SAM/BAM file that have mismatches
written by Mo •
-
What's the absolute easiest way of visualizing stranded RNA-seq data on custom genomes?
written by Trivas •
-
What are the columns from pysam fetch?
written by Eric •
-
Adding information tags with types, positions and counts of matches and mismatches in every read in…
written by lechu •
-
Program/script to count how many sequences have mismatches on each position of my primers, and also…
written by gastonviarengo •
-
Filtering bam files generated by bwa mem
written by kulvait
-
Deeptools sample scaling
written by Constantine
-
questions regarding Oncotator and Annovar
written by mangfu100
-
Gsnap: How To Output Only Alignments With Minimal Number Of Mismatches?
written by Herpson •
-
Next Generation Read Sequence Assembly Visualization Software
written by Neeraj •
