EnglishFrenchSpanish

Ad


OnWorks favicon

fastq-mcf - Online in the Cloud

Run fastq-mcf in OnWorks free hosting provider over Ubuntu Online, Fedora Online, Windows online emulator or MAC OS online emulator

This is the command fastq-mcf that can be run in the OnWorks free hosting provider using one of our multiple free online workstations such as Ubuntu Online, Fedora Online, Windows online emulator or MAC OS online emulator

PROGRAM:

NAME


fastq-mcf - ea-utils: detect levels of adapter presence, compute likelihoods and locations
of the adapters

SYNOPSIS


fastq-mcf [options] <adapters.fa> <reads.fq> [mates1.fq ...]

DESCRIPTION


Version: 1.04.676

Detects levels of adapter presence, computes likelihoods and locations (start, end) of the
adapters. Removes the adapter sequences from the fastq file(s).

Stats go to stderr, unless -o is specified.

Specify -0 to turn off all default settings

If you specify multiple 'paired-end' inputs, then a -o option is required for each. IE:
-o read1.clip.q -o read2.clip.fq

OPTIONS


-h This help

-o FIL Output file (stats to stdout)

-s N.N Log scale for adapter minimum-length-match (2.2)

-t N % occurance threshold before adapter clipping (0.25)

-m N Minimum clip length, overrides scaled auto (1)

-p N Maximum adapter difference percentage (10)

-l N Minimum remaining sequence length (19)

-L N Maximum remaining sequence length (none)

-D N Remove duplicate reads : Read_1 has an identical N bases (0)

-k N sKew percentage-less-than causing cycle removal (2)

-x N 'N' (Bad read) percentage causing cycle removal (20)

-q N quality threshold causing base removal (10)

-w N window-size for quality trimming (1)

-H remove >95% homopolymer reads (no)

-X remove low complexity reads (no)

-0 Set all default parameters to zero/do nothing

-U|u Force disable/enable Illumina PF filtering (auto)

-P N Phred-scale (auto)

-R Don't remove N's from the fronts/ends of reads

-n Don't clip, just output what would be done

-C N Number of reads to use for subsampling (300k)

-S Save all discarded reads to '.skip' files

-d Output lots of random debugging stuff

Quality adjustment options:
--cycle-adjust
CYC,AMT Adjust cycle CYC (negative = offset from end) by amount AMT

--phred-adjust
SCORE,AMT Adjust score SCORE by amount AMT

--phred-adjust-max
SCORE Adjust scores > SCORE to SCOTE

Filtering options*:
--[mate-]qual-mean
NUM Minimum mean quality score

--[mate-]qual-gt
NUM,THR At least NUM quals > THR

--[mate-]max-ns
NUM Maxmium N-calls in a read (can be a %)

--[mate-]min-len
NUM Minimum remaining length (same as -l)

--homopolymer-pct
PCT Homopolymer filter percent (95)

--lowcomplex-pct
PCT Complexity filter percent (95)

If mate- prefix is used, then applies to second non-barcode read only

Adapter files are 'fasta' formatted:

Specify n/a to turn off adapter clipping, and just use filters

Increasing the scale makes recognition-lengths longer, a scale of 100 will force
full-length recognition of adapters.

Adapter sequences with _5p in their label will match 'end's, and sequences with _3p in
their label will match 'start's, otherwise the 'end' is auto-determined.

Skew is when one cycle is poor, 'skewed' toward a particular base. If any nucleotide is
less than the skew percentage, then the whole cycle is removed. Disable for methyl-seq,
etc.

Set the skew (-k) or N-pct (-x) to 0 to turn it off (should be done for miRNA, amplicon
and other low-complexity situations!)

Duplicate read filtering is appropriate for assembly tasks, and never when read length <
expected coverage. -D 50 will use 4.5GB RAM on 100m DNA reads - be careful. Great for RNA
assembly.

*Quality filters are evaluated after clipping/trimming

Homopolymer filtering is a subset of low-complexity, but will not be separately tracked
unless both are turned on.

Use fastq-mcf online using onworks.net services


Free Servers & Workstations

Download Windows & Linux apps

  • 1
    ExpressLuke GSI
    ExpressLuke GSI
    This SourceForge download page was to
    grant users to download my source built
    GSIs, based upon phhusson's great
    work. I build both Android Pie and
    Android 1...
    Download ExpressLuke GSI
  • 2
    PyQt
    PyQt
    PyQt is the Python bindings for
    Digia's Qt cross-platform
    application development framework. It
    supports Python v2 and v3 and Qt v4 and
    Qt v5. PyQt is avail...
    Download PyQt
  • 3
    garlic3
    garlic3
    Gambas3* coding tools and examples for
    use in SMB or other organizations that
    make internal programs. E.g. from
    database fields to class sourcecode with
    DbInfo...
    Download garlic3
  • 4
    Sardi
    Sardi
    Sardi is a complete restyling and
    optimisation of svg code. 6 choices for
    your applications and 10 kind of folders
    to use in your file manager. The sardi
    icons...
    Download Sardi
  • 5
    LMMS Digital Audio Workstation
    LMMS Digital Audio Workstation
    LMMS is a free cross-platform software
    which allows you to produce music with
    your computer. If you like this project
    consider getting involved in the project
    h...
    Download LMMS Digital Audio Workstation
  • 6
    FreeRTOS Real Time Kernel (RTOS)
    FreeRTOS Real Time Kernel (RTOS)
    FreeRTOS is a market-leading real-time
    operating system (RTOS) for
    microcontrollers and small
    microprocessors. Distributed freely
    under the MIT open source lice...
    Download FreeRTOS Real Time Kernel (RTOS)
  • 7
    Avogadro
    Avogadro
    Avogadro is an advanced molecular
    editor designed for cross-platform use
    in computational chemistry, molecular
    modeling, bioinformatics, materials
    science and ...
    Download Avogadro
  • More »

Linux commands

Ad