question-mark
Stuck on an issue?

Lightrun Answers was designed to reduce the constant googling that comes with debugging 3rd party libraries. It collects links to all the places you might be looking at while hunting down a tough bug.

And, if you’re still stuck at the end, we’re happy to hop on a call to see how we can help out.

/1 and /2 suffixes in paired-end reads

See original GitHub issue

Hi, Can umi-tools analyze paired-end reads with /1 and /2 suffixes instead of 1:N and 2:N ?

trying with v0.5.5 and get the error message below with umi-tools extract:

umi-tools extract 
--bc-pattern=XXXXXX 
--bc-pattern2 CCCCCCNNNN 
--stdin H619_RG_C000QSP_6_1_CC9HDACXX_S0_L001_R1_001.fastq.gz 
--stdout H619_RG_C000QSP_6_1_CC9HDACXX_S0_L001_R1_001_extracted.fastq.gz 
--read2-in H619_RG_C000QSP_6_1_CC9HDACXX_S0_L001_R2_001.fastq.gz 
--read2-out H619_RG_C000QSP_6_1_CC9HDACXX_S0_L001_R2_001_extracted.fastq.gz 
--filter-cell-barcode 
--whitelist=whitelist.txt

Read pairs do not match H3:1:CC9HDACXX:6:2209:1046:2076/1 != H3:1:CC9HDACXX:6:2209:1046:2076/2

Thanks,

Issue Analytics

  • State:closed
  • Created 5 years ago
  • Comments:17 (13 by maintainers)

github_iconTop GitHub Comments

1reaction
IanSudberycommented, Mar 15, 2019

I wonder if instead of having two options we should just have one, that takes the delimiter and only uses it if one is provided. In fact I reckon that we already have a default, which is space.

On Fri, 15 Mar 2019, 4:16 pm Tom Smith, notifications@github.com wrote:

Great. Let me know when you’ve run the alignment and deduplication and confirmed this all works OK. Then I’ll merge this into the master and close the issue.

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/CGATOxford/UMI-tools/issues/325#issuecomment-473349042, or mute the thread https://github.com/notifications/unsubscribe-auth/AFJFjo81ivtg5ebDxY-BGKhS2q1nLZ36ks5vW8dOgaJpZM4btnh9 .

0reactions
TomSmithCGATcommented, Sep 14, 2020

This is now resolved in the master branch with option --ignore-read-pair-suffixes and will be available in the next release shortly

Read more comments on GitHub >

github_iconTop Results From Across the Web

adding suffix /1 & /2 to PE data - abyss input data - Biostars
Hi Abyssers, I have run trimmomatric on my PE data which generates R1 & R2 files that do not have any suffixes. Now...
Read more >
Read Names for Paired-End Reads - Support Illumina
By a common convention, read names can include suffixes (such as /1 or /2), which indicate the end of a pair the read...
Read more >
seqtools - McDonald Lab
Validates order of paired-end reads in given fastq files. ... Adds \1 and \2 suffix (tags) to the first (left) and second (right)...
Read more >
User guide — Cutadapt 4.2 documentation - Read the Docs
This is useful, for example, if you work with merged overlapping paired-end reads. Add the $ character to the end of an adapter...
Read more >
Manual - Bowtie
Reads files containing paired-end reads will sometimes name the reads according to whether they are the #1 or #2 mates by appending a...
Read more >

github_iconTop Related Medium Post

No results found

github_iconTop Related StackOverflow Question

No results found

github_iconTroubleshoot Live Code

Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start Free

github_iconTop Related Reddit Thread

No results found

github_iconTop Related Hackernoon Post

No results found

github_iconTop Related Tweet

No results found

github_iconTop Related Dev.to Post

No results found

github_iconTop Related Hashnode Post

No results found