Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • IGV shows bowtie2 mapped reads not corresponding to genomic sequence

    Hi,
    I have used Bowtie2 to map my reads to my genome. The reads are similar to ChIP seq reads. I filtered my reads for a MAPQ score > 10. When I look at the alignments using IGV, reads with a very high MAPQ score (say, 42 for example), seem to not match very well at all on the base-pair level for the genomic regions of interest I'm looking at (regulatory regions). How can this be? Doesn't a high MAPQ score from Bowtie2 mean it is certain that read maps there?

    My Bowtie command:
    bowtie2 --phred33 --fr -I 100 -X 500 -p 8 --seed 123 -q -t -x Sea_Urchin_Bowtie2_ref_index -1 clean_R1.fastq -2 clean_R2.fastq -S Sea_Urchin_sample_1_bowtie2.sam

    [IMG][/IMG]

    Any suggestions or comments or help?

    Thanks!

  • #2
    There's a discrepancy between what genome you aligned to, and what genome is in IGV, or the reference index is messed up.

    So I'd start by removing that reference genome from IGV, close IGV, to make sure it's not remembering the old index, remaking the index with samtools faidx, and then reimporting it into IGV, and trying that. That's pretty fast.

    If that doesn't work, reindex the genome with bowtie, and realign, to make 100% sure that the genome you are aligning to is consistent all the way through.

    Comment


    • #3
      Thank you, that makes sense. Even though I thought they were both using the 3.0 version of the S. Purpuratus genome, I just used whatever IGV had rather than importing my own.

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Current Approaches to Protein Sequencing
        by seqadmin


        Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
        04-04-2024, 04:25 PM
      • seqadmin
        Strategies for Sequencing Challenging Samples
        by seqadmin


        Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
        03-22-2024, 06:39 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, Yesterday, 12:08 PM
      0 responses
      11 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 04-10-2024, 10:19 PM
      0 responses
      17 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 04-10-2024, 09:21 AM
      0 responses
      14 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 04-04-2024, 09:00 AM
      0 responses
      43 views
      0 likes
      Last Post seqadmin  
      Working...
      X