Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • how to sort out mate no mapped reads

    Hi,
    I am now doing some analysis about chromosome structure Viaration. Please help me about following questions.
    1. Is there any method to sort out reads which mate is no mapped by IGV dispaly. I have dreamed that pair end all mapped reads are not displayed, so I could see what I wanted.
    2.If IGV can't do it, is there any methods to sort out single mapped reads from SAM or BAM files. As to the sinle mapped reads,some single mapped reads may scatter around the genome, some may concentrated on one loci. the later situation may reflect a real SV. So is there any method to find concentrated single mapped reads.

  • #2
    with samtools you can perform the analysis :

    if you have paired-end reads : "samtools view -f 8 -F 4 file.bam" will give you mapped reads with unmapped mate.

    if you have unpaired + paired reads: "samtools view -F1 file.bam" will give you unpaired mapped reads

    reads concentration, as you say, means deep coverage. To check the coverage along the genome you could use genomeCoverageBed:

    genomeCoverageBed -ibam sortedBamFile.bam -g genome.txt > coverage.txt

    Paired reads with a bad insert size or mis-oriented also could be evidence of structural variation.

    please see this image:



    When you map your reads to some reference you will see similar patterns in regions with structural variation.

    Comment


    • #3
      fastq indicate "table not found"

      Who can tell me why err like this after I enter "fastq-dump -X 5 -Z SRR616966"
      2015-04-13T15:23:47 fastq-dump.2.1.7 err: table not found while opening manager within database module - failed to open 'SRR616966'

      Comment

      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, 07-19-2024, 07:20 AM
      0 responses
      119 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 07-16-2024, 05:49 AM
      0 responses
      105 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 07-15-2024, 06:53 AM
      0 responses
      103 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 07-10-2024, 07:30 AM
      0 responses
      43 views
      0 likes
      Last Post seqadmin  
      Working...
      X