Should I use -a or -g when demultiplexing ONT reads with dual barcodes?
Nobody has claimed this yet.
- Dominant language
- Python
- Stars
- 587
- Forks
- 144
- Avg merge
- 2h 37m
- Merged PRs (30d)
- 1
Description
cutadapt 4.9
I have 16S amplicon reads that were sequenced with ONT that I am trying to demultiplex. Each sample was PCR barcoded with a 13 base barcode on both ends, so I expect a read to start with a barcode and end with its reverse complement. I put together a fasta file of all my pairs, some are listed below.
>HL001_FW
ATCCGGTCGGAGA...TCTCCGACCGGAT
>HL002_FW
CTGAGGTGATCAG...CTGATCACCTCAG
>HL003_FW
AGTGTCCTGCTAG...CTAGCAGGACACT
>HL004_FW
ATAAGCAATTCGA...TCGAATTGCTTAT
The problem I run into is whether to use the -a or -g flag. Looking through the documentation I see it used almost interchangeably for linked adapters, but I get different outputs depending on which I use and I'm not sure which is correct. I used the below commands, for reference
cutadapt -e 1 -a file:barcodes_for_cutadapt.fasta -o trimmed-{name}.fastq.gz reads.fastq.gz
cutadapt -e 1 -g file:barcodes_for_cutadapt.fasta -o trimmed-{name}.fastq.gz reads.fastq.gz
Contributor guide
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start with the cutadapt documentation for the -a and -g options and linked adapters, then reproduce both command lines using the supplied barcode FASTA structure. Compare the demultiplexing results and determine what guidance the documentation should provide about choosing between the two options.
Written by the indexing model from the issue text.
Assessment
- Domain
- bioinformatics
- Issue type
- Documentation
- Difficulty
- 3/5
- Estimated time
- 1-2 days
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 35/100