analyzing samples with no Barcodes and no Linker primer squences

Skip to first unread message


Sep 20, 2016, 4:11:14 PM9/20/16
to Qiime 1 Forum
Hello all! 
I have been asked to analyze a data set that has already been processed. It is now 4 fastq files with No barcodes and No linker primer sequences. 

So I ran the file -b and got no errors

I was thinking that I would try to split them: -f 3BC_Sky_1Ab.fastq -o 3BC_Sky_1Ab -c fastq_to_fastaqual -b

And then wanted to split the libraries so that I don't need to go through and manually add the headers to the fasta file. But now I am getting a mapping file error with the command -m map1.txt -f 3BC_Sky_1Ab.fna -b 0 -o split_lib_out1 -q 3BC_Sky_1Ab.qual

Error message: 
ValueError: Errors were found with mapping file, please run to identify problems.
dyrdahlyoung@if-plp-haustorium2:~/ECBMetagenome/ECBMetagenomics/3BC_Sky_1Ab$  qiime > -m map1.txt -f 3BC_Sky_1Ab.fna -b 0 --barcode_type 'not_barcoded' -c  -o split_lib_out4 -q 3BC_Sky_1Ab.qual -l 10
Traceback (most recent call last):
  File "/usr/local/bin/", line 411, in <module>
  File "/usr/local/bin/", line 408, in main
  File "/usr/lib/python2.7/dist-packages/qiime/", line 1289, in preprocess
    barcode_type, added_demultiplex_field)
  File "/usr/lib/python2.7/dist-packages/qiime/", line 310, in check_map
    'identify problems.')
ValueError: Errors were found with mapping file, please run to identify problems.


I also tried to run the fastq files with the command:

dyrdahlyoung@if-plp-haustorium2:~/ECBMetagenome/ECBMetagenomics$  qiime > -i 3BC_Sky_1Ab.fastq -m 3BC_Sky_1Ab/map1.txt --barcode_type 'not-barcoded' -s 1 -o split_lib_fastq --sample_ids I3BC

Traceback (most recent call last):
  File "/usr/local/bin/", line 365, in <module>
  File "/usr/local/bin/", line 344, in main
    for fasta_header, sequence, quality, seq_id in seq_generator:
  File "/usr/lib/python2.7/dist-packages/qiime/", line 239, in process_fastq_single_end_read_file_no_barcode
  File "/usr/lib/python2.7/dist-packages/qiime/", line 317, in process_fastq_single_end_read_file
    parse_fastq(fastq_read_f, strict=False, phred_offset=phred_offset)):
  File "/usr/lib/python2.7/dist-packages/skbio/parse/sequences/", line 174, in parse_fastq
skbio.parse.sequences._exception.FastqParseError: Failed qual conversion for seq id: 1R3YL:00017:00042. This may be because you passed an incorrect value for phred_offset.

Thanks for your thoughts!! I have attached my -tf

Reply all
Reply to author
0 new messages