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