Changes

From Genome Analysis Wiki
Jump to: navigation, search

Examples of Read Mapping with Karma and BWA

150 bytes added, 10:06, 11 December 2009
m
no edit summary
# Some instructions for read mapping and variant calling using the
# University of Michigan tools and procedures. Please view this as source.
# -- Paul Anderson and Tom Blackwell, December 11, 2009 --
# software components of this process are:
# bwa index # -- Do this only when a new version of the genome reference # sequence is released. Takes just under two hours. # bwa aln # -- Run this on every .fastq file individually. Highly # variable timings -- takes between 10 minutes and many # hours per .fastq file. Better data runs quicker. # bwa samse # -- Converts a single .fastq / .sai pair to .sam alignment # format. Usually under 1 minute per file. # bwa sampe # -- Converts paired end .fastq / .sai pairs (four files total) # to a single .sam alignment file. 1 - 2 minutes per run. # (In general, the CEU trio chromosome 20 test data set # contains VERY small .fastq files. Chromosome 20 is just # over 2% of the entire genome.)
# Sample command lines. This is csh syntax.
29
edits

Navigation menu