Hello every one.
I am new in bioinformatic, and now I try solve the next problem. I have sorted BAM file with paired-end RNA-Seq data (several tens of thousands reads). My work tool is R. I need: 1) Filter this file by "isize" field: extract the reads with isize not equal to NA and with absolute value of isize lower than threshold (in this case, 642 bp); 2) Write this subset into SAM file, because this type is required for package which I will use.
I tried use R package Rsamtools for this, function filterBam(). Remove reads with isize == NA is possible by means scanBamFlag:
flag <- scanBamFlag(isPaired = TRUE, hasUnmappedMate = FALSE).
But I have no idea about filtration by threshold (abs(bamFile$isize) <= 642). How I can do it in R? I do not believe that effective way for this does not exist. Unfortunately, my experience in R and programming for now is inadequate for solve this problem quickly. Can someone to help me? Thanks in advance!
I am new in bioinformatic, and now I try solve the next problem. I have sorted BAM file with paired-end RNA-Seq data (several tens of thousands reads). My work tool is R. I need: 1) Filter this file by "isize" field: extract the reads with isize not equal to NA and with absolute value of isize lower than threshold (in this case, 642 bp); 2) Write this subset into SAM file, because this type is required for package which I will use.
I tried use R package Rsamtools for this, function filterBam(). Remove reads with isize == NA is possible by means scanBamFlag:
flag <- scanBamFlag(isPaired = TRUE, hasUnmappedMate = FALSE).
But I have no idea about filtration by threshold (abs(bamFile$isize) <= 642). How I can do it in R? I do not believe that effective way for this does not exist. Unfortunately, my experience in R and programming for now is inadequate for solve this problem quickly. Can someone to help me? Thanks in advance!