Ylaine/2009-7-15

From ZhangLabWiki
Jump to navigation Jump to search

Comparing Maq/Maq and Bwa/Sam

Outputs from the two aligning/consensus calling methods often have different read numbers associated with SNPs at the same location, suggesting that the mapping algorithms produce different results. This is confirmed by direct examination of the output 'pileup' files.

  • Run Sam on Maq output. Created new directory "MaqSamTarget"
/Users/kunzhang/Downloads/samtools/misc/maq2sam-long ../MaqTarget/all.map > all.sam
/Users/kunzhang/Downloads/samtools/samtools import /Users/kunzhang/WorkSpace/Exome/targets/ccdsExon_unique_miRNA_200bp.fa.combined.fai all.sam all.bam