MAINTENANCE OUTAGE: The University Wiki Service will undergo maintenance on September 26th, 2017, from 6 pm to 8 pm. During this 2 hour time period https://wikis.utexas.edu may be unavailable. Users are advised to save content locally that may be needed during this time and to otherwise save all edits as unsaved work may be lost. Please contact the UT Service Desk at 512-475-9400 for any questions.
The University Wiki Service has upgraded the Confluence Server software, from version 5.9.14 to 5.10.8. Please refer to the knowledge base article, KB0015891, for a high level summary of upgrade changes. Thank you!
Skip to end of metadata
Go to start of metadata

Using Maq 0.6.6 for colorspace alignment:

Note: Compatibility issues with Maq 0.7.1 and Maqview 0.2.3 have been experienced. So, use Maq 0.6.6 and maqview 0.2.3 together.

With colorspace reads file readsF3.csfasta and readsF3_QV.qual and reference file ref.fasta

1. solid2fastq.pl reads out : this converts the csfasta and quality file into mock based fastq file called out.fastq

2. maq fastq2bfq out.fastq out.bfq : this converts the fastq file into binary format

3. maq fasta2csfa ref.fasta > ref.csfa  : this converts the reference file into mock based reference

4. maq fasta2bfa ref.csfa > ref.csbfa : this converts the mock based reference file into binary format

5. maq fasta2bfa ref.fasta > ref.bfa  : this converts the original reference file into binary format (for use in step 7 and 8)

6. maq map -c aln.cs.map ref.csbfa out.bfq 2>aln.log : this maps the reads to the reference- the mapping is done in colorspace.

7. maq csmap2nt aln.nt.map ref.bfa aln.cs.map : this converts the color alignment to nucleotide alignment (aln.nt.map). After this step, the color reference file and color alignment will not need to be used.

8. maq assemble out.cns ref.bfa aln.nt.map : this creates a consensus sequence (out.cns) from the mapping results

Additionally, to convert any mapping result file (.map) to a readable tab delimited file, use:

maq mapview aln.nt.map > aln.nt.info

  • No labels