Importing read pairs from BAM with specific R1 and R2 start/end positions for pair
Hi there,
Is it possible to set up a ScanBamParam() to import read pairs where the R1 and R2 reads start and end at specific positions? Note that this is not a range, but rather a conditional to only import the read pair if the R1 has a certain position AND its R2 pair has a certain position.
If not with ScanBamParam(), is there another way to filter for this? Help would be very much appreciated :)
Kind regards, Emily
• 689 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
Emily Wasson •
•
written
by
Emily Wasson •
1 answer
using samjdk https://jvarkit.readthedocs.io/en/latest/SamJdk/ assuming the MC (mate cigar) attribute is available.
java -jar dist/jvarkit.jar samjdk -e 'return record.getReadPairedFlag() && !record.getReadUnmappedFlag() && record.getContig().equals("chr1") && record.getReadPairedFlag() && !record.getMateUnmappedFlag() && record.getContig().equals(record.getMateReferenceName()) && record.getStart()==12345 && SAMUtils.getMateAlignmentEnd(record)==56789;' in.bam
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
Looking for RNA-seq aligner that can identify not properly mapped paired
written by romain.lannes
-
Matching R1 and R2 in paired end reads from Illumina
written by Priyanka •
-
Extract the overlapping paired-end reads
written by ATRX
-
Identifying mismatches between read mates at position in genome
written by wpearman1996 •
-
How to filter both reads of paired end reads by AS when only 1 of the 2 reads has an AS below 50
written by jelle •
-
Paired-end reads that span a position in an assembly
written by saladi •
-
isolate reads from a paired-ended fastq files that fall within a certain bp range
written by nibhelim •
-
merging multiple pair end read files into one
written by Bioinfonext
-
How to trim both pair end reads and remove low quality reads
written by Bioinfonext
-
unpaired R1 and R2 sequences
written by kamiljaron
