Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • cedance
    replied
    you can use "samtools view" to output a *.bam file on the go.
    however using R + IRanges seems to be the best option for me. Else you'll have to find modules on interval trees (if any) and implement in perl. But why reinvent the wheel if the purpose here is to get the job done.

    Leave a comment:


  • Dario1984
    replied
    Do you use R ? It is easy to do. You need the packages Rsamtools and GenomicRanges installed.

    Rsamtools has a function called readBamGappedAlignments. Use it on each file to read the BAM file into memory.

    Then, use the findOverlaps function in the package GenomicRanges to find which reads in one sample overlap with reads in the other.

    You don't need to create unnecessary SAM files on your hard drive this way.

    Leave a comment:


  • fabrice
    replied
    Yes. I want to compare how many mapped reads are overlap in two method.

    Thanks

    Originally posted by upendra_35 View Post
    Is your aim is to compare the % of mapped reads? Then you should convert this to sam file. There is a related query at this thread: http://seqanswers.com/forums/showthr...hlight=upendra.

    Good luck!

    Leave a comment:


  • upendra_35
    replied
    Is your aim is to compare the % of mapped reads? Then you should convert this to sam file. There is a related query at this thread: http://seqanswers.com/forums/showthr...hlight=upendra.

    Good luck!

    Leave a comment:


  • fabrice
    started a topic Compare two bam files.

    Compare two bam files.

    I have used two different map software to map my pair-end reads.

    I want to how many reads are reported as mapped by these two map software and how many reads only reported in one map software.

    Is there some easy ways to do this?

    Thank you

Latest Articles

Collapse

  • seqadmin
    Exploring the Dynamics of the Tumor Microenvironment
    by seqadmin




    The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
    07-08-2024, 03:19 PM
  • seqadmin
    Exploring Human Diversity Through Large-Scale Omics
    by seqadmin


    In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
    06-25-2024, 06:43 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 05:49 AM
0 responses
12 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-15-2024, 06:53 AM
0 responses
23 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-10-2024, 07:30 AM
0 responses
36 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-03-2024, 09:45 AM
0 responses
204 views
0 likes
Last Post seqadmin  
Working...
X