Can GSNAP generate SAM/BAM files?

[efoss]

I am running gsnap with this command:

gsnap --gunzip -d Drosophila5 -D /home/efoss/gene_databases/GMAP_GSNAP_db/ D09NJACXX_s8_1_illumina12index_8_SL7777.fastq.gz D09NJACXX_s8_2_illumina12index_8_SL7777.fastq.gz > D09NJACXX_s8_illumina12index_8_SL7777.sam

The output clearly is not in standard sam format. Is it possible to generate standard sam/bam files with gsnap, and if so, how?

Thank you.

Eric

[cedance]

Are specific columns missing? or certain columns are 0? what's the problem you see?

[efoss]

Hi cedance,

I see output like this:


>CNAANATATATTGCACCATCCGCCGATCCGCAGGCAATTTTTAGGCCATT 0 unpaired <#0<#28899;<>><<?=<?>?@<==49=?;===?>1:;9:>>??;>??< HWI-ST156:39709NJACXX:8:1101:1790:1988

<ATCACAAACGAACACTTTAAAATCAGCTGTGATTTTGCTCGGGCTTCCTT 0 unpaired F?@@DFDDDDF;1C;BBDAEGCBFHJI<CGEIJEIEBF?D)6D@?FDHB= HWI-ST156:39709NJACXX:8:1101:1790:1988

>TNGANCTGGACCAGGTAATGGGATATACTGTTGTCCGGCGCCTGTTCAGC 0 unpaired 3#0:#2288>>?@@<@=??=<7:>9>>9@?@??>>9=88=?8?<1=???? HWI-ST156:39709NJACXX:8:1101:1980:1992

<AAAAACACTCACACCAAACATCAGCAGGATGCTGCAGCCCCCGGATGAGA 0 unpaired F<@<D+ADDDFF3??EDG?:CF?F<@;F811899?@C<@FG>AA'B;==E HWI-ST156:39709NJACXX:8:1101:1980:1992

>CTGANACTCGGGATCTGCACGGCGATAGCGGGCACCCAGGCTGGAGATGA 0 unpaired @?@F#2=BCDAA?FHIGCHIGIIIBEEGEHGGBACHEDFBCCCBABC>>> HWI-ST156:39709NJACXX:8:1101:1775:1997

This doesn't match the mandatory fields for a sam file:

http://samtools.sourceforge.net/SAM1.pdf

Thank you for the help.

Eric

[cedance]

I remember vaguely but there is a --format=sam parameter to specifically output the results in SAM format. You should check out the parameter list and add it to your command line options accordingly.

[efoss]

I remember vaguely but there is a --format=sam parameter to specifically output the results in SAM format. You should check out the parameter list and add it to your command line options accordingly.

Thanks so much, cedance. I see that that is one of my critical mistakes. I thought I had read through the options carefully, but I guess I hadn't read them carefully enough. I very much appreciate the help.

Best wishes,

Eric