Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Mark
    started a topic BBMAP paired-end alignment problem

    BBMAP paired-end alignment problem

    Hi All

    Trying to use bbmap to align some Illumina PE reads. My command line looks like this

    Code:
    bbmap.sh overwrite=t ref=A.fasta in=R1_001_val_1.100.fq in2=R2_001_val_2.100.fq minid=0.25 mappedonly=true out=out.sam  outu=unaligned_reads.fa maxindel=100000
    When I do this out.sam contains alignments of only read1 though the log indicates that read2 reads were aligned

    Code:
    Read 1 data:            pct reads       num reads       pct bases          num bases
    
    mapped:                 100.0000%             100       100.0000%              29779
    unambiguous:            100.0000%             100       100.0000%              29779
    ambiguous:                0.0000%               0         0.0000%                  0
    low-Q discards:           0.0000%               0         0.0000%                  0
    
    perfect best site:        0.0000%               0         0.0000%                  0
    semiperfect site:        35.0000%              35        35.1758%              10475
    rescued:                  0.0000%               0
    
    Match Rate:                   NA               NA        88.8176%              26449
    Error Rate:              65.0000%              65         4.4763%               1333
    Sub Rate:                65.0000%              65         0.3627%                108
    Del Rate:                 0.0000%               0         0.0000%                  0
    Ins Rate:                49.0000%              49         4.1136%               1225
    N Rate:                 100.0000%             100         6.7061%               1997
    
    
    Read 2 data:            pct reads       num reads       pct bases          num bases
    
    mapped:                 100.0000%             100       100.0000%              25636
    unambiguous:            100.0000%             100       100.0000%              25636
    ambiguous:                0.0000%               0         0.0000%                  0
    low-Q discards:           0.0000%               0         0.0000%                  0
    
    perfect best site:        0.0000%               0         0.0000%                  0
    semiperfect site:        26.0000%              26        25.3745%               6505
    rescued:                  2.0000%               2
    
    Match Rate:                   NA               NA        86.8232%              22258
    Error Rate:              74.0000%              74         4.9852%               1278
    Sub Rate:                74.0000%              74         0.8894%                228
    Del Rate:                 0.0000%               0         0.0000%                  0
    Ins Rate:                42.0000%              42         4.0958%               1050
    N Rate:                 100.0000%             100         8.1916%               2100

    If I change the command line to

    Code:
    bbmap.sh overwrite=t ref=A.fasta in=R2_001_val_2.100.fq minid=0.25 mappedonly=true out=out.sam  outu=unaligned_reads.fa maxindel=100000
    Read2 reads are indeed aligned.

    Any thoughts what I'm doing wrong not to get paired end alignments with the first command line? I am using BBMap version 35.92.

    Thanks

    Mark

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

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 04-11-2024, 12:08 PM
0 responses
35 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 10:19 PM
0 responses
38 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 09:21 AM
0 responses
33 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-04-2024, 09:00 AM
0 responses
54 views
0 likes
Last Post seqadmin  
Working...
X