View Single Post
Old 05-09-2014, 11:35 AM   #2
Location: SE MN

Join Date: Oct 2013
Posts: 44

By its nature, SAM has the alignment in it (SAM=Sequence Alignment Map).

First: Did you sort your SAM file? (will need samtools for this).

Checking the manual again, it looks like it should be:

velveth $FOLDERNAME $KMER -reference ref.fasta -shortPaired -sam illumina_align.sam
Can you post the first few lines of your sam file by using

head illumina_align.sam
Since you already have to do the alignment of your reads to the reference, have you tried looking at the alignments? Are there dips in coverage or re-arrangements, things like that?
ctseto is offline   Reply With Quote