text: user14146: Hi I'm new to bioinformatics and have just completed my first run on the MinION (long read sequencing Oxford Nanopore Technologies). I was hoping someone could direct me towards R packages, workflow, tutorials or guides that will help me identify species that are present in my sample mainly for fungi as my amplicons had ITS1 and ITS2 regions.
  I have been watching Youtube tutorials on RStudio for over a week and I feel like I'm no closer than I was when I started. I have both fastq and fast5 files and I tried to run the fastq files using EPI2ME, however, the database is lacking in fungi. I have tried using NCBI's BLAST but I can only run <5kb reads and most of my sequencing reads are >6kb. Also, I've tried using Galaxy to convert my fastq files to fasta but I'm not sure how to run my long reads for species identification using Galaxy or if that is even possible.
  Any help would be greatly appreciated. Thanks.
  
  user92: There are a number of useful metagenomics tools out there, including Centrifuge and Kraken2, as @gringer has mentioned.
  Also, FWIW, the WIMP tool from ONT, a Centrifuge based metageonomic classifier, does include fungus sequences derived from RefSeq, though it may not include the species you are looking for.
  If you want to use NCBI BLAST, you could split your long reads into pieces that are <5kb (there are a number of tools out there for FASTQ manipulation - my current favorite tool for this is seqkit, where you could do something like this):
  seqkit sliding -s 4900 -W 4900 test1.fastq
  The above will split your FASTQ reads into chunks of 4.9kb. You can adjust the -s parameter if you instead want overlapping chunks.
  From a GUI standpoint, you can use Galaxy for the above by making a workflow with the following steps:
  
  Trim Sequences (enter reasonable numbers here, to trim each of your reads to <5kb)
  FASTQ to FASTA converter
  NCBI BLAST+
  
  Another possibility is QIIME2, which can be used with the UNITE fungal training set. See recent reference (looking at the ITS1/ITS2 regions), UNITE, and QIIME2.
  Long term, you may want to look into using the command line, as it is much more powerful and customizable, but if BLAST gives you reasonable results for your reads, the above process in Galaxy may work for you.
  
original_text: 14146: <p>Hi I'm new to bioinformatics and have just completed my first run on the MinION (long read sequencing Oxford Nanopore Technologies). I was hoping someone could direct me towards R packages, workflow, tutorials or guides that will help me identify species that are present in my sample mainly for fungi as my amplicons had ITS1 and ITS2 regions.</p>
  <p>I have been watching Youtube tutorials on RStudio for over a week and I feel like I'm no closer than I was when I started. I have both fastq and fast5 files and I tried to run the fastq files using EPI2ME, however, the database is lacking in fungi. I have tried using NCBI's BLAST but I can only run &lt;5kb reads and most of my sequencing reads are &gt;6kb. Also, I've tried using Galaxy to convert my fastq files to fasta but I'm not sure how to run my long reads for species identification using Galaxy or if that is even possible.</p>
  <p>Any help would be greatly appreciated. Thanks.</p>
   92: <p>There are a number of useful metagenomics tools out there, including Centrifuge and Kraken2, as @gringer has mentioned.</p>
  <p>Also, FWIW, the <a href="https://epi2me.nanoporetech.com/workflow-docs-1" rel="nofollow noreferrer">WIMP tool from ONT</a>, a Centrifuge based metageonomic classifier, does include fungus sequences derived from RefSeq, though it may not include the species you are looking for.</p>
  <p>If you want to use NCBI BLAST, you could split your long reads into pieces that are &lt;5kb (there are a number of tools out there for FASTQ manipulation - my current favorite tool for this is <a href="https://bioinf.shenwei.me/seqkit/" rel="nofollow noreferrer">seqkit</a>, where you could do something like this):</p>
  <p><code>seqkit sliding -s 4900 -W 4900 test1.fastq</code></p>
  <p>The above will split your FASTQ reads into chunks of 4.9kb. You can adjust the <code>-s</code> parameter if you instead want overlapping chunks.</p>
  <p>From a GUI standpoint, you can use Galaxy for the above by making a workflow with the following steps:</p>
  <ol>
  <li>Trim Sequences (enter reasonable numbers here, to trim each of your reads to &lt;5kb)</li>
  <li>FASTQ to FASTA converter</li>
  <li>NCBI BLAST+</li>
  </ol>
  <p>Another possibility is QIIME2, which can be used with the UNITE fungal training set. See <a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7032056/" rel="nofollow noreferrer">recent reference</a> (looking at the ITS1/ITS2 regions), <a href="https://unite.ut.ee/repository.php" rel="nofollow noreferrer">UNITE</a>, and <a href="https://docs.qiime2.org/2021.11/" rel="nofollow noreferrer">QIIME2</a>.</p>
  <p>Long term, you may want to look into using the command line, as it is much more powerful and customizable, but if BLAST gives you reasonable results for your reads, the above process in Galaxy may work for you.</p>
  