Changes between Version 24 and Version 25 of SOPs/qc_shortReads
- Timestamp:
- 04/28/16 14:56:49 (9 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
SOPs/qc_shortReads
v24 v25 220 220 Many functions \\ 221 221 222 = Analyzing potential species composition of reads in a fastq file =222 = Analyzing potential species composition/contamination of reads in a fastq file = 223 223 224 224 * To get a quick preview of what genomes/collections (or vectors or other contaminants) of sequences the reads in a fastq file can map to, one can use the FastQ Screen tool. … … 226 226 * FastQ Screen allows you to screen your fastq file against a set of libraries, which can be set up to contain the genomes of interest, potential contaminating genomes, vectors, sequencing adaptors, ribosomal RNAs, or other contaminants commonly seen in sequencing experiments. This allows you to see if composition of the reads matches with what you expect and, where the contamination might have come from, if any. 227 227 * FastQ Screen uses bowtie or bowtie2 for mapping. The libraries you wish to screen against need to be bowtie-indexed. 228 * Paths to the bowtie indexed libraries need to be specified in the fastq_screen.conf file, which is called either from the same directory where the fastq_screen program is (by default) or from a manually specified location using the -conf /path/to/.conf file optio in.228 * Paths to the bowtie indexed libraries need to be specified in the fastq_screen.conf file, which is called either from the same directory where the fastq_screen program is (by default) or from a manually specified location using the -conf /path/to/.conf file option 229 229 * Sample commands are: 230 230 {{{ 231 231 232 "bsub fastq_screen [OPTION] file.fq" 232 Usage: fastq_screen [OPTIONS] file.fq 233 eg. fastq_screen --aligner bowtie2 myFastQ.txt 233 234 234 235 Commonly used options are: … … 237 238 [--outdir] Specify a directory in which to save output files. 238 239 [--illumina1_3] Assume that the quality values are in encoded in Illumina v1.3 format. Defaults to Sanger format if this flag is not specified. 239 [--conf] Manually specify a location for the configuration file to be used for this run. If not specified then the file will be taken from the same directory as the fastq_screen program. 240 [--conf] Manually specify a location for the configuration file to be used for this run. If not specified then the file will be taken from the same directory as the fastq_screen program. Note: the config file is already setup on our internal server 240 241 241 242 }}}