SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Write the subset of reads from BAM file into new SAM/BAM file, using R tools. (http://seqanswers.com/forums/showthread.php?t=65845)

Old Pioneer 01-27-2016 04:41 AM

Write the subset of reads from BAM file into new SAM/BAM file, using R tools.
 
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!


All times are GMT -8. The time now is 04:37 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.