You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
While extracting the unmapped reads from a nvBowtie bam file, the unmapped R2 is always replaced by a reverse-compliment version of R1.
This is very similar to an issue raised in 2015 https://groups.google.com/g/nvbio-users/c/is28EEvm2QE
# R1 read from bam file
CAAATGTATCTCTCTCTCTCACACACAGTATCCAGATAACTGATTACTGGAATGTGTGATAGAATAATACTACTGCAGCCACGAATGGTATCTATTTGAAAAGTCTTCCTTGAATAGAAGTCTAATGCCGTCTACAGGATGTAGTAGATG
# R2 read from bam file, which is reverse_complement of R1
CATCTACTACATCCTGTAGACGGCATTAGACTTCTATTCAAGGAAGACTTTTCAAATAGATACCATTCGTGGCTGCAGTAGTATTATTCTATCACACATTCCAGTAATCAGTTATCTGGATACTGTGTGTGAGAGAGAGAGATACATTTG
# R1 read from test_1.fq
CAAATGTATCTCTCTCTCTCACACACAGTATCCAGATAACTGATTACTGGAATGTGTGATAGAATAATACTACTGCAGCCACGAATGGTATCTATTTGAAAAGTCTTCCTTGAATAGAAGTCTAATGCCGTCTACAGGATGTAGTAGATG
# R2 read from test_2.fq
AACTTAGAGTTCACTCTGTACAGATAGATAGATACAAGTTACCACAGAGATCATACTACATCTACTACATACTGTAAACGGCATTAGACTTCTATTCAAGGACGACGATTCAAATAGATACCATTCGTGGCTGCAGAACTATAATTCTAA
The text was updated successfully, but these errors were encountered:
While extracting the unmapped reads from a nvBowtie bam file, the unmapped R2 is always replaced by a reverse-compliment version of R1.
This is very similar to an issue raised in 2015 https://groups.google.com/g/nvbio-users/c/is28EEvm2QE
Here is an example to reproduce the error.
A tar file with minimal reference file and fq files is available at https://sid.erda.dk/share_redirect/hMPcixT4XT
This results in
The text was updated successfully, but these errors were encountered: