Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Wrong mapping with BBmap default parameters

    I performed mapping for my RNA-seq data (150 bp reads) with BBMap, just the default parameters (the only special parameter set was low memory):
    Code:
    bbmap.sh ref=gencode.v29.transcripts.fa path=index_bb_lomem in=sampleA_1.fq.gz out=sampleA.bam overwrite=f usemodulo
    When I looked in the resulting sam file, I noticed that for many reads I had CIGAR string looking like this:
    2=4X6=11I1X7=1X3=1X14=1X2=1X1=1X15=1X3=2X9=1X9=2X21=2X23=1X5=
    This how the alignment looked (upper line is the read, lower is suggested hit):
    CAGAAGTCAAGGAGTAATTTAGGAGTAATTTTGACTTTCAAGTCTTATTACTTAATAAATACATTTCATAAAGCTACAGCTGCCATAGATAGTGATTCCTCTGATGGATCTGGGCAAAGTCAATTGAAAACCTTCTGGAAAGGAGTCACC
    ..TGGA......************G.......C...C..............T..G.G...............G...CT.........G.........GT.....................CA.......................T.....

    It is clearly incorrect mapping: this read should be just assigned to unmapped. From 1414 reads assigned to my gene of interest (attached), 949 had editing distance of 21 and more (34 reads had distance of 40 and more, up to 357!), i.e. obviously didn't belong to that gene. It means that a lot of the counts I get with the default settings are totally wrong. I am wondering if this is normal result for the mapper and what parameters I need to adjust to avoid this issue - `minratio`, `minind`, or something else?
    Attached Files

  • #2
    Bit late, but you could do this

    nmtag=f Write NM tags.

    and then filter by NM = number of mismatches later.

    This might be easier:

    editfilter=-1 Ban alignments with more than this many edits.

    OR also

    idfilter=0 Independant of minid; sets exact minimum identity
    allowed for alignments to be printed. Range 0 to 1.
    (set to 0.95 ?)

    subfilter=-1 Ban alignments with more than this many substitutions.
    (set to 5 ?)

    Comment


    • #3
      Colin, thanks a lot for your input!
      (Unfortunately for now I had to switch to STAR because it's hard to get help for BBmap questions while STAR's author is very responsive and helpful)

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Best Practices for Single-Cell Sequencing Analysis
        by seqadmin



        While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
        06-06-2024, 07:15 AM
      • seqadmin
        Latest Developments in Precision Medicine
        by seqadmin



        Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

        Somatic Genomics
        “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
        05-24-2024, 01:16 PM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, Yesterday, 06:58 AM
      0 responses
      13 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 06-06-2024, 08:18 AM
      0 responses
      20 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 06-06-2024, 08:04 AM
      0 responses
      18 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 06-03-2024, 06:55 AM
      0 responses
      13 views
      0 likes
      Last Post seqadmin  
      Working...
      X