analyzing samples with no Barcodes and no Linker primer squences

31 views
Skip to first unread message

Roody_UF

unread,
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 validate_map.py file -b and got no errors

I was thinking that I would try to split them: 

convert_fastaqual_fastq.py -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 split_libraries.py command

split_libraries.py -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 validate_mapping_file.py to identify problems.
dyrdahlyoung@if-plp-haustorium2:~/ECBMetagenome/ECBMetagenomics/3BC_Sky_1Ab$  qiime > split_libraries.py -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/split_libraries.py", line 411, in <module>
    main()
  File "/usr/local/bin/split_libraries.py", line 408, in main
    truncate_ambi_bases=opts.truncate_ambi_bases)
  File "/usr/lib/python2.7/dist-packages/qiime/split_libraries.py", line 1289, in preprocess
    barcode_type, added_demultiplex_field)
  File "/usr/lib/python2.7/dist-packages/qiime/split_libraries.py", line 310, in check_map
    'identify problems.')
ValueError: Errors were found with mapping file, please run validate_mapping_file.py to identify problems.

~~~~~

I also tried to run the fastq files with the split_libraries_fastq.py command:

dyrdahlyoung@if-plp-haustorium2:~/ECBMetagenome/ECBMetagenomics$  qiime > split_libraries_fastq.py -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/split_libraries_fastq.py", line 365, in <module>
    main()
  File "/usr/local/bin/split_libraries_fastq.py", line 344, in main
    for fasta_header, sequence, quality, seq_id in seq_generator:
  File "/usr/lib/python2.7/dist-packages/qiime/split_libraries_fastq.py", line 239, in process_fastq_single_end_read_file_no_barcode
    phred_offset=phred_offset):
  File "/usr/lib/python2.7/dist-packages/qiime/split_libraries_fastq.py", 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/fastq.py", line 174, in parse_fastq
    seqid)
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 print_config_all.py -tf

-Roo
print_qiime_config.py.txt
Reply all
Reply to author
Forward
0 new messages