# Problems with q2-cutadapt multiplexed paired end sequences without sequence ID lines or quality scores

**URL:** <https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583>\
**Category:** User Support\
**Created:** [August 14, 2018, 10:12pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583 "2018-08-14T22:12:17Z")\
**Posts on this page:** 16\
**Page:** 1

<div class="post-metadata">

**Author:** ![ropolifroni](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/r/cab0a1/32.png) [@ropolifroni](https://forum.qiime2.org/u/ropolifroni)\
**Post date:** [August 14, 2018, 10:12pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/1 "2018-08-14T22:12:17Z")

</div>

Hi everyone;  
I have sequences in two files (R1 and R2) that I assume as forward and reverse senses. I followed the q2-cutadapt tutorial in [Demultiplexing and Trimming Adapters from Reads with q2-cutadapt](https://forum.qiime2.org/t/demultiplexing-and-trimming-adapters-from-reads-with-q2-cutadapt/2313).  
I tried with:  
qiime tools import \

> --type MultiplexedPairedEndBarcodeInSequence   
> --input-path foward.fastq.gz   
> --input-path reverse.fastq.gz   
> --output-path multiplexed-seqs.qza  
> Usage: qiime tools import [OPTIONS]

Error: Option --input-path was specified multiple times in the command.

Then I tried with forward file only as a single end barcode in sequence and I have error again:

$ qiime tools import \

> --type MultiplexedSingleEndBarcodeInSequence   
> --input-path forward.fastq.gz   
> --output-path multiplexed-seqs.qza

There was a problem importing forward.fastq.gz:

forward.fastq.gz is not a(n) FastqGzFormat file:

Header on line 5 is not FASTQ, records may be misaligned

I checked the file and didn't find nothing. I don't realize where is the mistake.

Could you give me an idea what can I do?

Thank a lot.

---

<div class="post-metadata">

**Author:** ![ebolyen](https://forum.qiime2.org/user_avatar/forum.qiime2.org/ebolyen/32/11_2.png) [@ebolyen](https://forum.qiime2.org/u/ebolyen)\
**Post date:** [August 15, 2018, 5:26pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/2 "2018-08-15T17:26:26Z")

</div>



---

<div class="post-metadata">

**Author:** ![Mehrbod\_Estaki](https://forum.qiime2.org/user_avatar/forum.qiime2.org/mehrbod_estaki/32/4001_2.png) [@Mehrbod\_Estaki](https://forum.qiime2.org/u/Mehrbod_Estaki)\
**Post date:** [August 15, 2018, 10:03pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/3 "2018-08-15T22:03:45Z")

</div>

Hi @ropolifroni,

The first error you receive during your paired-end import is because you are providing 2 lines of `--input-path` when the script is expecting one.

See this note at the bottom of the link you posted:

# Regarding paired-end reads

- The import format for paired-end reads with the barcodes still in the sequence is `MultiplexedPairedEndBarcodeInSequence` - this format expects two files in a directory ( `forward.fastq.gz` and `reverse.fastq.gz` ).

So you need to create a directory that contains your `forward.fastq.qz`, and `reverse.fastq.gz` files and refer to this folder's path in your instead in your `--input-path` parameter.

As for your second error I'm not quite sure what is happening. Could you try running this command and post the output of it:

`zcat forward.fastq.gz | head`

Finally, keep in mind one another point from that tutorial:

- Demultiplexing currently only works if the barcodes are in the forward reads — we plan to support dual-indexing strategies in a future release of QIIME 2.

So if you have barcodes in your forward AND reverse reads (dual indexed) you will not be able to use cutadapt to demultiplex them. You would have to do this outside of qiime2 then import back here. You can do this in qiime1 or some other software of your choice.

Keep us posted!

---

<div class="post-metadata">

**Author:** ![ropolifroni](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/r/cab0a1/32.png) [@ropolifroni](https://forum.qiime2.org/u/ropolifroni)\
**Post date:** [August 16, 2018, 3:45pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/4 "2018-08-16T15:45:36Z")

</div>

> [@Mehrbod\_Estaki](#):
>
> zcat forward.fastq.gz | head

Hi Mehrbod Estaki for your advices.  
I've created a directory where I copied R1 and R2 sequences. Then I run the script that you said. I run the comand in qiime2 and the result was:

$ zcat forward.fastq.gz | head  
NAGTCTCAGTGCAGCCGCCGCGGTAATACGGAGGGTGCAAGCGTTGTTCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGCTTTTTAAGTCAGATGTGAAAGTCCACGGCTCAACCGTGGAGGTGCATTTGAAACTGGGAAGCTTGAGTATTGGAGGGAGTAGTGGAATTCCCGGTGTAGAGGTGAAATTCGTAGATATCGGGAGGAATACCAGTGGCGAAGGCGACTACCTGGCCAAATACTGACACTGAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCGGGTAGTCC  
NAGTGTGAGTGCCAGCAGCCGCGGTAATACAGAGGTCTCGAGCGTTAATCGGAGTCACTGGGCTTAAAGGGTGCGTAGGCGGCGCTGCAAGTGCCTTGTGAAATCCCATCGCTCAACGATGGAAATGCTTGGTAGACTGTAGTGCTTGAGGTTGGTAGGGGCTGTCGGAACTCTAGGTGGAGCGGTGAAATGCGTAGATATCTAGAGGAACGCCAAAGGTGAAGACAGACAGCTGGGCCAATTCTGACGCTGAGGCACGAAAGCGTGGGGAGCAAACGGGATTAGATACCCGCGTAGTCC  
NAGTGTGAGTGCCAGCAGCCGCGGTAATACAGAGGTCTCGAGCGTTAATCGGAGTCACTGGGCTTAAAGGGTGCGTAGGCGGCGCTATAAGTGCCTTGTGAAATCCCATCGCTCAACGATGGAATTGCTTGGTAGACTGTAGTGCTTGAGGTTGGTAGGGGCTGTCGGAACTCTAGGTGGAGCGGTGAAATGCGTAGATATCTAGAGGAACGCCAAAGGTGAAGACAGACAGCTGGGCCAATTCTGACGCTGAGGCACGAAAGCGTGGGGAGCAAACGGGATTAGATACCCGCGTAGTCCC  
GAGTGTGAGTGCCAGCCGCCGCGGTAATACAGAGGTCTCGAGCGTTAATCGGAGTCACTGGGCTTAAAGGGTGCGTAGGCGGCGCTACAAGTGCCTTGTGAAATCCCATCGCTCAACGATGGAATTGCTTGGTAGACTGTAGTGCTTGAGGTTGGTAGGGGCTGTCGGAACTCTAGGTGGAGCGGTGAAATGCGTAGATATCTAGAGGAACGCCAAAGGTGAAGACAGACAGCTGGGCCAATTCTGACGCTGAGGCACGAAAGCGTGGGGAGCAAACGGGATTAGATACCCGCGTAGTCC  
GAGTAGTGGTGCCAGCCGCCGCGGTAATACAGAGGGGGCTAGCGTTGTTCGGAATTACTGGGCGTAAAGCGCACGTAGGCGGATCGGAAAGTTGGGGGTGAAATCCCGGGGCTCAACCCCGGAACTGCCTTCAAAACTCCCGGTCTTGAGTTCGAGAGAGGTGAGTGGAATTCCGAGTGTAGAGGTGAAATTCGTAGATATTCGGAGGAACACCAGTGGCGAAGGCGGCTCACTGGCTCGATACTGACGCTGAGGTGCGAAAGTGTGGGGAGCAAACAGGATTAGATACCCGGGTAGTCC  
GAGTTCACGTGCCAGCAGCCGCGGTAATACGAACTGCGCAAACGTTATTCGGAATCACTGGGCTTACAGGGTGCGTAGGCGGTCCACCAGGTCAGATGTGAAATCCCACAGCTCAACTGGGGAACTGCATTTGAAACCGGCGGACTTGAGTAGACTAGAGGTACACGGAACTTCAGGTGGAGCGGTGAAATGTGTTGAGATCTGAAGGAACACCAGAGGCGAAAGCGGTGTACTGGAGTCTGTCTGACGCTGAGGCACGAAAGCTAGGGGAGCGAACGGGATTAGATACCCGTGTAGTCC  
GAGTCACTGTGCCAGCAGCCGCGGTAATACGGAGGGGGCTAGCGTTGTTCGGAATTACTGGGCGTAAAGCGCACGTAGGCGGCTATCCAAGTTGGGGGTGAAATCCCTGGGCTTAACCCAGGAACTGCCTTCAAAACTGGATGGCTAGAGTCCGAGAGAGGTGAGTGGAATTTCCAGTGTAGAGGTGAAATTCGTAGATATTGGAAAGAACACCAGTGGCGAAGGCGGCTCACTGGCTCGGTACTGACGCTGAGGTGCGACAGCGTGGGGAGCAAACAGGATTAGATACCCGCGTAGTCC  
GAGTCACTGTGCCAGCAGCCGCGGTAATACGGAGGGGGCTAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGAATGTGAAAGCCCCGGGCTCAACCTGGGAACTGCATCCAAAACTGGCAAGCTAGAGTATGGCAGAGGGTGGTGGAATTTCCTGTGTAGCGGTGAAATGCGTAGATATAGGAAGGAACACCAGTGGCGAAGGCGACCACCTGGGCTAATACTGACACTGAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCCGGTAGTCC  
GAGTCAGAGTGCCAGCAGCCGCGGTAATACGGAGGGGGCTAGCGTTGTTCGGAATTACTGGGCGTAAAGCGCACGTAGGCGGATCGGAAAGTTGGGGGTGAAATCCCGGGGCTCAACCCCGGAACTGCCTTCAAAACTCCCGGTCTTGAGTTCGAGAGAGGTGAGTGGAATTCCGAGTGTAGAGGTGAAATTCGTAGATATTCGGAGGAACACCAGTGGCGAAGGCGGCTCACTGGCTCGATACTGACGCTGAGGTGCGAAAGTGTGGGGAGCAAACAGGATTAGATACCCGCGTAGTCCC  
GAGTAGTGGTGCCAGCAGCCGCGGTAATACGAAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGAATGTGAAAGCCCCGGGCTCAACCTGGGAACTGCATCCAAAACTGGCAAGCTAGAGTATGGCAGAGGGTGGTGGAATTTCCTGTGTAGCGGTGAAATGCGTAGATATAGGAAGGAACACCAGTGGCGAAGGCGACCACCTGGGCTAATACTGACACTGAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCCCGTAGTCC

I'll be waiting your answer.  
Thanks in advance.

Ro

---

<div class="post-metadata">

**Author:** ![Mehrbod\_Estaki](https://forum.qiime2.org/user_avatar/forum.qiime2.org/mehrbod_estaki/32/4001_2.png) [@Mehrbod\_Estaki](https://forum.qiime2.org/u/Mehrbod_Estaki)\
**Post date:** [August 16, 2018, 4:07pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/5 "2018-08-16T16:07:41Z")

</div>

Hmm, this certainly doesn't look like what I expected. Could you also run the same as above for your other file, the R2, whatever its called now.

---

<div class="post-metadata">

**Author:** ![ropolifroni](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/r/cab0a1/32.png) [@ropolifroni](https://forum.qiime2.org/u/ropolifroni)\
**Post date:** [August 16, 2018, 4:20pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/6 "2018-08-16T16:20:44Z")

</div>

> [@ropolifroni](#):
>
> $ zcat forward.fastq.gz | head

Here is the result:

zcat reverse.fastq.gz | head  
zcat reverse.fastq.gz | head  
GAGTAGACGTGCCAGCCGCCGCGGTAATACGAAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCGCGTAGGTGGTTTGTTAAGTTGAATGTGAAAGCCCCGGGCTCAACCTGGGAACTGCATCCATAACTGGCAAGCTAGAGTATGGCAGAGGGTGGTGGAATTTCCTGTGTAGCGGTGAAATGCGTAGATATAGGAAGGAACACCAGTGGCGAAGGCGACCACCTGGGCTAATACTGCCACTGAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATATCCGTGTTGTCCC  
GAGTAGACGTGCCAGCAGCCGCGGTAATACGAAGGGTGCAAGCGTTACTCGGAATTACTGGGCGTAAAGCGTGCGTAGGTGGTGGTTTAAGTCCGTTGTGAAAGCCCTGGGCTCAACCTGGGAATTGCAGTGGATACTGGATCACTCGAGTGTGGTAGAGGATAGTGGAATTTCCGGTGTAGCAGTGAAATGCGTAGAGATCGGAAGGAACATCTGTGGCGAAGGCGACTATCTGGGCCAACACTGACACCGAGGCACGCAAGCGTGGGGAGCCAACAGGATTAGCAACCCGGTTAGTCAC  
GAGTAGTGGTGCCAGCCGCCGCGGTAATACGAAGGGTGCAAGCGTTAATCGGAATTACTGGGCGTAAAGCGCACGTAGGCGGATCGGAAAGTTGGGGGTGAAATCCCGGGGCTCAACCCCGGAACTGCCTTCAAAACTCCCGGTCTTGAGTTCGAGAGAGGTGAGTGGAATTCCGAGTGTAGAGGTGAAATTCGTAGATATTCGGAGGAACACCAGTGGCGAAGGCGGCTCACTGGCTCGATACTGACGCTGAGGTGCGAAAGTGTGGGGAGCAAAGAGGATTAGATACCCCGGGAGTCCC  
GAGTAGACGTGCCAGCCGCCGCGGTAATACGTAGGGTGCGAGCGTTGTCCGGAATTACTGGGCGTAAAGAGCTCGTAGGCGGTTTGTCACGTCGTCTGTGAAATCCCTCGGCTTAACCGGGGGCGTGCAGGCGATACGGGCAGACTTGAGTACTACAGGGGAGACTGGAATTCCTGGTGTAGCGGTGAAATGCGCAGATATCAGGAGGAACACCGGTGGCGAAGGTGGGTCTCTGTGTAGTAACTGACGCTGAGGAGGGAAAGCATGGGTAGCGAACGGGATTAGGTATCCCGGTGGTCCC  
AGTGCCAGCAGCCGCGGTAATACGAACTGCGCAAACGTTATTCGGAATCACTGGGCTTACAGGGTGCGTAGGCGGTCCACCAGGTCAGATGTGAAATCCCACAGCTCAACTGTGGAACTGCATTTGAAACCGGCGGACTTGAGTAGACTAGAGGTACACGGAACTTCAGGTGGAGCGGTGAAATGTGTTGAGATCTGAAGGAACACCAGAGGCGAAAGCGGTGTACTGGAGTCTGTCTGACGCTGAGGCACGAAAGCTAGGGGAGCGAACGGGATTAGATACCCCTGTAGTACAGAACGGG  
GAGTCACTGTGCCAGCAGCCGCGGTAATACGGAGGGGGCTAGCGTTGTTCGGAATTACTGGGCGTAAAGCGCACGTAGGCGGATCGGAAAGTTGGGGGTGAAATCCCGGGGCTCAACCCCGGAACTGCCTTCAAAACTCCCGGTCTTGAGTTCGAGAGAGGTGAGTGGAATTCCGAGGGTAGAGGTGAAATTCGTAGATATTCGGAGGAACACCAGTGGCGAAGGCGGCTCACTGGCTCGATACTGCCGCTGAGGTGCGAAAGTGTGGGGAGCAAACAGGATTAGAGACCCCTGCAGTCCG  
GAGTAGTGGTGCCAGCAGCCGCGGTAATACGGAGGGGGCTAGCGTTGTTCGGAATTACCGGGCGTAAAGCGCACGTAGGCGGATCGGAAAGTTGGGGGTGAAATCCCGGGGCTCAACCCCGGAACTGCCTTCAAAACTCCCGGTCTTGAGTTCGAGAGAGGTGAGTGGAATTCCGAGTGTAGAGGTGAACTTCGTAGATATTCGGAGGAACACCAGTGGCGAAGGCGGCTCACTGGCTCGATCCTGCCGCTGAGGTGCGAAAGGGTGGGGGGCAAACCGGCTCAGATACCCTCGTAGGCCC  
GAGTCACTGTGCCAGCCGCCGCGGTAATACGGAGGGAGCTAGCGTTGTTCGGAATTACTGGGCGTAAAGCGCACGTAGGCGGCTTTTCAAGTCAGGGGTGAAATCCCGGGGCTCAACCCCGGAACTGCCCTTGAAACTGGATGGCTAGAATACTGGAGAGGTGAGTGGAATTCCGAGTGTAGAGGTGAAATTCGTAGATATTCGGAAGAACACCAGTGGCGAAGGCGACTCACTGGACCGTTATTGACGCTGAGGTGCGAAAGCGTGGGGAGCAAACAGGAATAGATCCCCAGGTAGGCCC  
GAGTCAGAGTGCCAGCAGCCGCGGTAATACGGAGGGGGCTAGCGTTGTTCGGAATTACTGGGCGTAAAGCGCACGTAGGCGGCTATCCAAGTTGGGGGTGAAATCCCTGGGCTTAACCCAGGAACTGCCTTCAAAACTGGATGGCTAGAGTCCGAGAGAGGTGAGTGGAATTTCCAGTGTAGAGGTGAAATTCGTAGATATTGGAAAGAACACCAGTGGCGAAGGCGGCTCACTGGCTCGGTACTGACGCTGAGGTGCGACAGCGGGGGGAGCAAACAGGATTAGATACCCGGGTAGTCCC  
GAGTCACTGTGCCAGCCGCCGCGGTAATACGAACTGCGCAAACGTTATTCGGAATCACTGGGCTTACAGGGTGCGTAGGCGGTCCACCAGGTCAGATGTGAAATCCCACAGCTCAACTGTGGAACTGCATTTGAAACCGGCGGACTTGAGTAGACTAGAGGTACACGGAACTTCAGGTGGAGCGGTGAAATGTGTTGAGATCTGAAGGAACACCAGAGGCGAAAGCGGTGTACTGGAGTCTGTCTGACGCTGAGGCACGAAAGCTAGGGGAGCGAACGGGATTAGATCCCCCGGTAGTCCC

What do you see wrong?

---

<div class="post-metadata">

**Author:** ![Mehrbod\_Estaki](https://forum.qiime2.org/user_avatar/forum.qiime2.org/mehrbod_estaki/32/4001_2.png) [@Mehrbod\_Estaki](https://forum.qiime2.org/u/Mehrbod_Estaki)\
**Post date:** [August 16, 2018, 5:02pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/9 "2018-08-16T17:02:55Z")

</div>

HI @ropolifroni,  
The format in which your files are in are not what we usually see. For example, a typical fastq file would look something like this:

```
@M02501:63:000000000-BT3BH:1:1101:14311:1491 1:N:0: CGTACTAG+TATCCTCT
GTGCCAGCAGCCGCGGTAACACGATTAACCCAAGTCAATAGAAGCCGGCGTAAAGAGTGTTTTGGATCACCCCCTCCCCAATAAAGCTAAAACTCACCTGAGTTGTAAAAAACTCCAGTTGACACAAAATAGACTACGAAAGTGGCTTTAACATATCTGAACACACAATAGCTAAGACCCAAACTGGGATTAGAAACCCTTGTAGTCC
+
CCCCCGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGCFFGGGGGGGGGGGEFGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG&lt;FEGGGGGGGGGFGGGGGGGGGGGGGGGGGGGGGGGGFGCDGGGGGDFGGGGGDGFGGGFGFGEGGCFGGGGGFGFGGFFFFFFFGF5FGGFCEG;,5= 

```

Where there are 4 lines associated with each read. 1=header, 2=sequences, 3=filler/blank, 4=quality scores. Your reads seem to be purely sequences without quality scores in 1 line format. Are you sure there are no other files associated with these 2 files you have? Something else that could map these files to individual samples or barcodes that could correspond with these?  
Some additional questions: How many samples are meant to be represented in these reads? What sequencing platform are these from? Has there been any other types of pre-processing or quality control done to these reads prior to this point, either by your or the sequencing facility? In particular, can you locate your barcodes in these lines you posted? Anything that could help us demultiplex these...because as it is I can't see a way to do that without that info.

---

<div class="post-metadata">

**Author:** ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)\
**Post date:** [August 16, 2018, 5:03pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/10 "2018-08-16T17:03:58Z")

</div>



---

<div class="post-metadata">

**Author:** ![Mehrbod\_Estaki](https://forum.qiime2.org/user_avatar/forum.qiime2.org/mehrbod_estaki/32/4001_2.png) [@Mehrbod\_Estaki](https://forum.qiime2.org/u/Mehrbod_Estaki)\
**Post date:** [August 16, 2018, 5:27pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/13 "2018-08-16T17:27:23Z")

</div>



---

<div class="post-metadata">

**Author:** ![Mehrbod\_Estaki](https://forum.qiime2.org/user_avatar/forum.qiime2.org/mehrbod_estaki/32/4001_2.png) [@Mehrbod\_Estaki](https://forum.qiime2.org/u/Mehrbod_Estaki)\
**Post date:** [August 16, 2018, 5:27pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/14 "2018-08-16T17:27:25Z")

</div>



---

<div class="post-metadata">

**Author:** ![ropolifroni](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/r/cab0a1/32.png) [@ropolifroni](https://forum.qiime2.org/u/ropolifroni)\
**Post date:** [August 16, 2018, 5:59pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/15 "2018-08-16T17:59:49Z")

</div>

I have 17 samples that I want to analyze in 3 goups or experiments (ExpA= 4 samples; ExpB= 8 samples and ExpC= 5 samples). Now I'd like to work with ExpA samples (4 samples).  
The barcodes of these 4 samples are:  
GAGTAGAC  
GAGTAGTG  
GAGTCACT  
GAGTCAGA  
I get them from the mapping.file provided by the sequence service.  
The platform were came from the sequences is Illumina MiSeq and they are paired, so I have 2 files with the sequences from the 17 samples!!! That is why I want to cut the 4 samples from ExpA.  
SAM1-17\_S5\_L001\_R1\_001.fastq and SAM1-17\_S5\_L001\_R2\_001.fastq are 2.1 Gb each. I get this 2 files from Illumina Basespace.  
I extracted the ExpA sequences from the SAM..... R1 and SAM.....R2 files using the barcodes named before.  
AHHHH!!! Now I'm realize!!! The extracts aren't in .fastq format!! Are only sequences... Sorry...

How I can get the barcode. gz file from the SAM....R1 and SAM....R2 files, in order to can go ahead?  
I tried the Qiime1 comand:  
extract\_barcodes.py -f inseqs\_R1.fastq -r inseqs\_R2.fastq -c barcode\_paired\_end --bc1\_len 6 --bc2\_len 6 -o processed\_seqs

And the results were:

Traceback (most recent call last):  
File "/home/rosana/miniconda3/envs/qiime1/bin/extract\_barcodes.py", line 4, in   
**import** ('pkg\_resources').run\_script('qiime==1.9.1', 'extract\_barcodes.py')  
File "/home/rosana/miniconda3/envs/qiime1/lib/python2.7/site-packages/pkg\_resources/ **init**.py", line 658, in run\_script  
self.require(requires)[0].run\_script(script\_name, ns)  
File "/home/rosana/miniconda3/envs/qiime1/lib/python2.7/site-packages/pkg\_resources/ **init**.py", line 1438, in run\_script  
exec(code, namespace, namespace)  
File "/home/rosana/miniconda3/envs/qiime1/lib/python2.7/site-packages/qiime-1.9.1-py2.7.egg-info/scripts/extract\_barcodes.py", line 12, in   
from skbio.util import create\_dir  
File "/home/rosana/miniconda3/envs/qiime1/lib/python2.7/site-packages/skbio/ **init**.py", line 15, in   
import [skbio.io](http://skbio.io)  
File "/home/rosana/miniconda3/envs/qiime1/lib/python2.7/site-packages/skbio/io/ **init**.py", line 309, in   
import\_module('skbio.io.clustal')  
File "/home/rosana/miniconda3/envs/qiime1/lib/python2.7/importlib/ **init**.py", line 37, in import\_module  
**import** (name)  
File "/home/rosana/miniconda3/envs/qiime1/lib/python2.7/site-packages/skbio/io/clustal.py", line 123, in   
from skbio.alignment import Alignment  
File "/home/rosana/miniconda3/envs/qiime1/lib/python2.7/site-packages/skbio/alignment/ **init**.py", line 230, in   
from .\_alignment import Alignment, SequenceCollection, StockholmAlignment  
File "/home/rosana/miniconda3/envs/qiime1/lib/python2.7/site-packages/skbio/alignment/\_alignment.py", line 21, in   
from skbio.stats.distance import DistanceMatrix  
File "/home/rosana/miniconda3/envs/qiime1/lib/python2.7/site-packages/skbio/stats/ **init**.py", line 45, in   
from .\_subsample import subsample, subsample\_counts, isubsample  
File "/home/rosana/miniconda3/envs/qiime1/lib/python2.7/site-packages/skbio/stats/\_subsample.py", line 22, in   
from .\_\_subsample import \_subsample\_counts\_without\_replacement  
File " **init**.pxd", line 155, in init skbio.stats.\_\_subsample (skbio/stats/\_\_subsample.c:3964)  
ValueError: numpy.dtype has the wrong size, try recompiling

What can I do in order to demultiplex the 17 samples, but could work in 3 independent experiments of 4, 8 and 5 samples each?  
I'll appreciate very much your help.

THANKS A LOT

Ro

---

<div class="post-metadata">

**Author:** ![ropolifroni](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/r/cab0a1/32.png) [@ropolifroni](https://forum.qiime2.org/u/ropolifroni)\
**Post date:** [August 16, 2018, 7:56pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/16 "2018-08-16T19:56:13Z")

</div>

I resolved the issue with an external software. Now I have forward.fastq.gz, reverse.fastq.gz and barcodes.fastq.gz files in order to process my data following the Atacama soil tutorial.

Thanks a lot.

Ro

---

<div class="post-metadata">

**Author:** ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)\
**Post date:** [August 16, 2018, 8:29pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/17 "2018-08-16T20:29:53Z")

</div>

Hi @ropolifroni,  
Glad you were able to convert your data into a usable format!

Would you mind posting the steps you used to convert your data? Just in case future users run into the same problems...

Thanks!

---

<div class="post-metadata">

**Author:** ![ropolifroni](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/r/cab0a1/32.png) [@ropolifroni](https://forum.qiime2.org/u/ropolifroni)\
**Post date:** [August 16, 2018, 8:44pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/18 "2018-08-16T20:44:14Z")

</div>

Hi Nicholas,  
I used the free software FASTQProcessor ([FASTQ Processor](http://www.mrdnalab.com/16freesoftware/fastq-processor.html)) (I not sure that I can share this http 🤨  
Was pretty easy, but it ran in windows.  
But my issue didn't finish, because I get the 3 files (forward.fastq.gz, reverse.fastq.gz and barecodes.fastq.gz) for the 17 samples, and not for only 4 samples from ExpA.  
Nonetheless, I ran the next command for the 17 samples:  
qiime demux emp-paired   
--m-barcodes-file sample-metadata.tsv   
--m-barcodes-column BarcodeSequence   
--i-seqs emp-paired-end-sequences.qza   
--o-per-sample-sequences demux   
--p-rev-comp-mapping-barcodes

The terminal give me an error:  
Plugin error from demux:

No sequences were mapped to samples. Check that your barcodes are in the correct orientation (see the rev\_comp\_barcodes and/or rev\_comp\_mapping\_barcodes options).

Debug info has been saved to /tmp/qiime2-q2cli-err-hcc9svrn.log

Here I live the error file content:  
Traceback (most recent call last):  
File "/home/rosana/miniconda3/envs/qiime2-2018.6/lib/python3.5/site-packages/q2cli/commands.py", line 274, in **call**  
results = action(\*\*arguments)  
File "", line 2, in emp\_paired  
File "/home/rosana/miniconda3/envs/qiime2-2018.6/lib/python3.5/site-packages/qiime2/sdk/action.py", line 232, in bound\_callable  
output\_types, provenance)  
File "/home/rosana/miniconda3/envs/qiime2-2018.6/lib/python3.5/site-packages/qiime2/sdk/action.py", line 367, in _callable\_executor_  
output\_views = self.\_callable(\*\*view\_args)  
File "/home/rosana/miniconda3/envs/qiime2-2018.6/lib/python3.5/site-packages/q2\_demux/\_demux.py", line 378, in emp\_paired  
raise ValueError('No sequences were mapped to samples. Check that '  
ValueError: No sequences were mapped to samples. Check that your barcodes are in the correct orientation (see the rev\_comp\_barcodes and/or rev\_comp\_mapping\_barcodes options).

Could you tell me what is wrong?

Thanks a lot.

Ro

---

<div class="post-metadata">

**Author:** ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)\
**Post date:** [August 20, 2018, 3:18pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/20 "2018-08-20T15:18:27Z")

</div>

Hi @ropolifroni,  
See this part of the error message:

> [@ropolifroni](#):
>
> No sequences were mapped to samples. Check that your barcodes are in the correct orientation (see the rev\_comp\_barcodes and/or rev\_comp\_mapping\_barcodes options).

Make sure your forward and reverse barcodes are also corresponding to the correct read files.

I hope that helps!

---

<div class="post-metadata">

**Author:** ![system](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/3X/2/1/21af5fe23cb6f4579467c66a9ed94e55274ca7bd.svg) [@system](https://forum.qiime2.org/u/system)\
**Post date:** [September 20, 2018, 9:26pm UTC](https://forum.qiime2.org/t/problems-with-q2-cutadapt-multiplexed-paired-end-sequences-without-sequence-id-lines-or-quality-scores/5583/21 "2018-09-20T21:26:12Z")

</div>

This topic was automatically closed 31 days after the last reply. New replies are no longer allowed.
