Changes

*  Binary reference genome in nucleotide space (see [[#Input file requirement|Input file requirement]]})  

*  Binary reference genome in nucleotide space (see [[#Input file requirement|Input file requirement]]})  

*  Binary reference genome and word index in color space (see \ref{sec:2})  

*  Binary reference genome and word index in color space (see [[#Build Binary Reference Genome and Word Index|Build Binary Reference Genome and Word Index]])  

*  Color space reads in valid color space FASTQ format (see [[#Input file requirement|Input file requirement]] for file specification)  

*  Color space reads in valid color space FASTQ format (see [[#Input file requirement|Input file requirement]] for file specification)  

*  Color space reads are longer than minimum length requirement. (see [[#Minimum read length requirement|Minimum read length requirement]])  

*  Color space reads are longer than minimum length requirement. (see [[#Minimum read length requirement|Minimum read length requirement]])  

*&nbsp; Specify color space parameter when starting KARMA (see \ref{sec:3})<br>  

*&nbsp; Specify color space parameter when starting KARMA (see [[#Map Color Space Reads|Map Color Space Reads]])<br>  

Please note the hardwares requirment for KARMA are:

Please note the hardwares requirment for KARMA are:

*30G disk space  

*30G disk space  

<br> We listed a complete example reviewing the whole procedure from building word index to mapping color space reads in \ref{sec:5}.

<br> We listed a complete example reviewing the whole procedure from building word index to mapping color space reads in [[#A Complete Example|A Complete Example]].

= Build Binary Reference Genome and Word Index<br> =

= Build Binary Reference Genome and Word Index<br> =

&nbsp; First, we need to build binary reference genome (option: --createReference)<br> &nbsp; \footnote{To let KARMA map nucleotide space reads, you need to use ``--createIndex''to create the word index file.}''<br>

&nbsp; First, we need to build binary reference genome (option: --createReference)<br> &nbsp;  

\footnote{To let KARMA map nucleotide space reads, you need to use ``--createIndex''to create the word index file.}''<br>

&nbsp; in nucleotide space. Assume NCBI36.fa is a FASTA file contains sequences of all chromosomes.<br> &nbsp; The command to invoke is:<br>

&nbsp; in nucleotide space. Assume NCBI36.fa is a FASTA file contains sequences of all chromosomes.<br> &nbsp; The command to invoke is:<br>

   karma --colorSpace --createReference --createIndex --reference NCBI36CS.fa

   karma --colorSpace --createReference --createIndex --reference NCBI36CS.fa

&nbsp; An important parameter is the size of words for indexing.<br> &nbsp; We recommand 15 (default value) for human reference genome.<br> &nbsp; Specifiy ``--wordSize N`` if you like to use $N$ as word size.<br> &nbsp; Typically you will observe performance change (see [[#Choose an appropriate size for word index|Choose an appropriate size for word index]]  for more discussion).<br> &nbsp;<br> &nbsp;<br> &nbsp; Note, multiple chromosomes are supported.<br> &nbsp; In current version, KARMA can take one FASTA file which contains sequences of all chromosomes.<br>

&nbsp; An important parameter is the size of words for indexing.<br> &nbsp; We recommand 15 (default value) for human reference genome.<br> &nbsp; Specifiy ``--wordSize N`` if you like to use ''N'' as word size.<br> &nbsp; Typically you will observe performance change (see [[#Choose an appropriate size for word index|Choose an appropriate size for word index]]  for more discussion).<br> &nbsp;<br> &nbsp;<br> &nbsp; Note, multiple chromosomes are supported.<br> &nbsp; In current version, KARMA can take one FASTA file which contains sequences of all chromosomes.<br>

<br>

<br>

= Map Color Space Reads =

=Map Color Space Reads=

&nbsp; KARMA takes valid color space FASTQ files inputs.<br> &nbsp; We usually use suffix .csfastq to distinguish it from nucleotide space reads.<br> &nbsp; For single end color space read, we can invoke command:<br>

&nbsp; KARMA takes valid color space FASTQ files inputs.<br> &nbsp; We usually use suffix .csfastq to distinguish it from nucleotide space reads.<br> &nbsp; For single end color space read, we can invoke command:<br>

Size for word index is sensitive to mapping performance. A small size of word index will increase the number of calculation cycles for a single read and duplications of a single word. On the other side, a big size will require much larger memory. Please also keep in mind that appropriate size is related to your hardware architecture. For practically purpose, we found size of 15 is optimal.

Size for word index is sensitive to mapping performance. A small size of word index will increase the number of calculation cycles for a single read and duplications of a single word. On the other side, a big size will require much larger memory. Please also keep in mind that appropriate size is related to your hardware architecture. For practically purpose, we found size of 15 is optimal.

= A Complete Example =

=A Complete Example=

<br>

<br>