SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Can GSNAP generate SAM/BAM files? (http://seqanswers.com/forums/showthread.php?t=14773)

efoss 10-15-2011 06:24 AM

Can GSNAP generate SAM/BAM files?
 
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 10-15-2011 07:22 AM

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

efoss 10-15-2011 07:46 AM

Hi cedance,

I see output like this:


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

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

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

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

>CTGANACTCGGGATCTGCACGGCGATAGCGGGCACCCAGGCTGGAGATGA 0 unpaired @?@F#2=BCDAA?FHIGCHIGIIIBEEGEHGGBACHEDFBCCCBABC>>> HWI-ST156:397:D09NJACXX: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 10-15-2011 08:02 AM

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 10-16-2011 08:11 PM

Quote:

Originally Posted by cedance (Post 54040)
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


All times are GMT -8. The time now is 09:30 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.