Pick De Novo OTUs Qiime
2
1
Entering edit mode
5.3 years ago
Ruth M. ▴ 10

Hi Everyone,

I'm trying to pick_de_novo_otus.py on 50 samples (.fa files) that I have trimmed and aligned (not using QIIME). It seemed straightforward, but I can't get the input format of the data right. I thought it could be comma delineated, but it says that the file doesn't exist 'home/Data/KK.01.fa,KK.02.fa,KK.03B.fa,KK.04.fa,KK.05.fa' (command below)

pick_de_novo_otus.py -i home/Data/KK.01.fa,KK.02.fa,KK.03B.fa,KK.04.fa,KK.05.fa -o home/Data/OTU&

I also tried using -i *.fa instead, but that gave the same result. When processing a large number of samples through the OTU picking, what is the proper way to input your files?

Thank you!!!

qiime OTU • 1.9k views
ADD COMMENT
0
Entering edit mode

Thanks, I'm unable to run split_libraries_fastq, so I have alternatively trimmed my sequences and merged the paired reads so that I now have individual .fa files for each sample rather than the .fna file. I ended up writing out pick_de_novo_otus.py script for each one in a shell script and running it using bash, and it seems to be working.

ADD REPLY
0
Entering edit mode
5.3 years ago
Picasa ▴ 610

You need to run split_libraries_fastq before.

http://qiime.org/scripts/split_libraries_fastq.html

ADD COMMENT
0
Entering edit mode
2.9 years ago
ankit hinsu ▴ 10

Hi,

I think that by this time, you might have received solution to your problem. Nevertheless, here is the solution.

  • QIIME1 requires input file in proper way (called QIIME-compatible format). Their tutorial "454 tutorial for de novo OTU picking" has a demultiplexing step. This step does three things: denmultiplexing, quality-filtering and combine all the sequences from all samples in QIIME-compatible way.

What you require is alternate to this step.

The answer is "add_qiime_labels.py" script. This script asks for a folder with all you input files and a mapping file having information about "which file belongs to which sample" [same mapping file used in QIIME with one extra column mentioning name of the file]. On running this script, it will combine all the sequences in QIIME-compatible format and produce a single fasta file (which you can use for OTU picking).

Hope it helps!!

ADD COMMENT

Login before adding your answer.

Traffic: 3086 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6