Reads mapping to mitochondrial DNA or unassigned sequences were discarded.
Sequence data were aligned to the specified National Center for Biotechnology Information human reference sequence after low-quality sequences were discarded.
OTUs with <0.005% of the total number of sequences were discarded as recommended by Navas-Molina et al. [12] for downstream analysis, resulting in 1011 OTUs.
Highly similar (≥95% identical), and incomplete sequences were discarded.
Metazoans, marine euryarchaeote group II sequences (obtained with EukA primer), and short (<100 bp) sequences were discarded.
If a sequence differed at more than one nucleotide from the highest scoring alignment, then the sequences were discarded into a separate dataset.
BLAST analysis of the contigs was performed to verify the α-gliadin identity of the contigs and short (<100 bp) and bad sequences were discarded.
The remaining cDNAs were compared to the LSUrdb, all 23S putative sequences were discarded and the remaining fraction was used to determine the functional content of the sample.
In each cluster, those sequences with more detailed taxonomic annotations were selected to represent the cluster, and other sequences were discarded.
Sequences were discarded as mapping to multiple sites when they had more than one match on the human genome differing in identity less than 2%.
Matched sequences were discarded.
