Skip site navigation (1)Skip section navigation (2)

FreeBSD Manual Pages

  
 
  

home | help
samtools-fqidx(1)	      Bioinformatics tools	       samtools-fqidx(1)

NAME
     samtools fqidx - indexes or queries regions from a fastq file

SYNOPSIS
     samtools fqidx ref.fastq [region1 [...]]

DESCRIPTION
     Index  reference  sequence  in the FASTQ format or extract subsequence from
     indexed reference sequence. If no region is specified, fqidx will index the
     file and create <ref.fastq>.fai on the disk. If regions are specified,  the
     subsequences will be retrieved and printed to stdout in the FASTQ format.

     The  input and output can be files compressed in the BGZF format. When out-
     put is compressed, the default compression level is 4.

     The sequences in the input file should all have different names.	If  they
     do  not,  indexing  will  emit  a warning about duplicate sequences and re-
     trieval will only produce subsequences from the first sequence with the du-
     plicated name.

     samtools fqidx should only be used on fastq files with a  small  number  of
     entries.	Trying to use it on a file containing millions of short sequenc-
     ing reads will produce an index that is almost as big as the original file,
     and searches using the index will be very slow and use a lot of memory.

OPTIONS
     -o, --output FILE
	     Write FASTQ to file rather than to stdout. If FILE ends  with  .gz,
	     .bgz or .bgzf then it will be BGZF compressed.

     -n, --length INT
	     Length  for  FASTQ  sequence line wrapping.  If zero, this means do
	     not line wrap.  Defaults to the line length in the input file.

     -c, --continue
	     Continue working if a non-existent region is requested.

     -r, --region-file FILE
	     Read regions from a file. Format is chr:from-to, one per line.

     -i, --reverse-complement
	     Output the sequence as the reverse complement.  When this option is
	     used, "/rc" will be appended to the sequence names.  To  turn  this
	     off or change the string appended, use the --mark-strand option.

     --mark-strand TYPE
	     Append strand indicator to sequence name.	TYPE can be one of:

	     rc     Append  '/rc'  when writing the reverse complement.  This is
		    the default.

	     no     Do not append anything.

	     sign   Append '(+)' for forward strand or '(-)' for reverse comple-
		    ment.  This matches the output of "bedtools getfasta -s".

	     custom,<pos>,<neg>
		    Append string <pos> to names when writing the forward strand
		    and <neg> when writing the reverse strand.	Spaces are  pre-
		    served,  so  it  is  possible to move the indicator into the
		    comment part of the description line by including a  leading
		    space in the strings <pos> and <neg>.

     --fai-idx FILE
	     Read/Write to specified index file.

     --gzi-idx FILE
	     Read/Write  to  specified	compressed  file  index  (used	with .gz
	     files).

     -h, --help
	     Print help message and exit.

     --output-fmt-option OPT=VAL
	     Set the output format options, level=0..9 for compression	level  0
	     to 9.

     --write-index
	     Create index for the output sequence data along with the output, in
	     same  path  as  <output name>.fai, <outputname>.gzi. This option is
	     valid only for file output.

     -@, --threads N
	     Set the number of extra threads for operations on compressed files.

AUTHOR
     Written by Heng Li,  with	modifications  by  Andrew  Whitwham  and  Robert
     Davies, all from the Sanger Institute.

SEE ALSO
     samtools(1), samtools-faidx(1), samtools-fasta(1), samtools-fastq(1)

     Samtools website: <http://www.htslib.org/>

samtools-1.23.1 		  18 March 2026 	       samtools-fqidx(1)

Want to link to this manual page? Use this URL:
<https://man.freebsd.org/cgi/man.cgi?query=samtools-fqidx&sektion=1&manpath=FreeBSD+Ports+15.1.quarterly>

home | help