This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Sequence Logo For Different Alleles Or Generated From Sam/Bam

I have a SAM/BAM file and instead of displaying it in IGV or another browser, I would like to make a sequence logo out of it. My favourite software WebLogo does not support SAM/BAM file as an input (instead it supports CLUSTALW, FASTA, plain flatfile, MSF, NBRF, PIR, NEXUS and PHYLIP)

I would like to avoid parsing SAM file by myself and adding '-' when needed. Do you have any idea if there is already some software that could do such a thing? Thanks a lot.

The second question is about different alleles. I would like to visualize sequences (reads) from each allele with another color. I have found many coloring schemes, but none dividing sequences into groups together with corresponding coloring. Thanks a lot.

sam
• 7,893 views • link
updated by Ram Ram • written by Biomonika (Noolean) Biomonika (Noolean)

3 answers

I quickly wrote a java tool sam4weblogo: https://github.com/lindenb/jvarkit#sam4weblogo using the picard library.

enter image description here

$ java -jar dist/sam4weblogo.jar I=path/to/samtools-0.1.18/examples/sorted.bam R=seq1:80-110 2> /dev/null 

>B7_593:4:106:316:452/1
TGTTG--------------------------
>B7_593:4:106:316:452a/1
TGTTG--------------------------
>B7_593:4:106:316:452b/1
TGTTG--------------------------
>B7_589:8:113:968:19/2
TGGGG--------------------------
>B7_589:8:113:968:19a/2
TGGGG--------------------------
>B7_589:8:113:968:19b/2
TGGGG--------------------------
>EAS54_65:3:321:311:983/1
TGTGGG-------------------------
>EAS54_65:3:321:311:983a/1
TGTGGG-------------------------
>EAS54_65:3:321:311:983b/1
TGTGGG-------------------------
>B7_591:6:155:12:674/2
TGTGGGGG-----------------------
>B7_591:6:155:12:674a/2
TGTGGGGG-----------------------
>B7_591:6:155:12:674b/2
TGTGGGGG-----------------------
>EAS219_FC30151:7:51:1429:1043/2
TGTGGGGGGCGCCG-----------------
>EAS219_FC30151:7:51:1429:1043a/2
TGTGGGGGGCGCCG-----------------
>EAS219_FC30151:7:51:1429:1043b/2
TGTGGGGGGCGCCG-----------------
>B7_591:5:42:540:501/1
TGTGGGGGCCGCAGTG---------------
>EAS192_3:5:223:142:410/1
TGGGGGGGGCGCAGT----------------
>B7_591:5:42:540:501a/1
TGTGGGGGCCGCAGTG---------------
>EAS192_3:5:223:142:410a/1
TGGGGGGGGCGCAGT----------------
>B7_591:5:42:540:501b/1
TGTGGGGGCCGCAGTG---------------
>EAS192_3:5:223:142:410b/1
TGGGGGGGGCGCAGT----------------
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Pierre, you are definitely Tolkien of bioinformatics! Thanks.

• 0 views • link
written by Biomonika (Noolean) Biomonika (Noolean)

Please help me in running sam4weblogo

I downloaded jvarkit-master.zip and unzip how to use it

• 0 views • link
updated by Ram Ram • written by Tark • Tark

you need to compile it; Read the manual please.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

I am having problems installing berkeleydb.jar, could you please give me short advice how to install it? (I am ubuntu user)

• 0 views • link
written by Biomonika (Noolean) Biomonika (Noolean)

you don't need bdb, if you have filled the path to the picard jars. Just type:

ant sam4weblogo

should work.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

There are many tools to convert a SAM file to FASTA:

You should be able to use one of these to convert your alignment to FASTA, and then make use of WebLogo.

• 22 views • link
written by Giovanni M Dall'Olio Giovanni M Dall'Olio

Thank you! I will check how they deal with insertions tomorrow, I am quite curious about that.

• 0 views • link
written by Biomonika (Noolean) Biomonika (Noolean)

Almost there:) sam2fasta does what I want except dealing with insertions (I would like to gap reference when insertions are introduced). SamToFastq does not seem to do the alignment.

• 162 views • link
written by Biomonika (Noolean) Biomonika (Noolean)

To solve this another way, without custom parsing (scripts do all the work for you):

  1. Grab the FASTA files for your genomic build-of-interest (e.g., grab and extract all the gzipped FASTA files for hg19 here).
  2. Convert your SAM to BED
  3. Convert the resulting BED to FASTA, incorporating the whole-genome FASTA file you made from Step 1
  4. Run the FASTA file through WebLogo to generate the sequence logo
• 0 views • link
written by Alex Reynolds Alex Reynolds

Thank you, interesting. Unfortunately I have no reference for my species.

• 0 views • link
written by Biomonika (Noolean) Biomonika (Noolean)

Log in to answer this question.

More posts like this
  • How To Sort Sam File
    written by Chen Sun Chen Sun
    <p>I want to process the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file created by <a href='http://bio-bwa.sourceforge.net/'>BWA</a>, if it is sorted by the sequence coordinate on reference, it will be …
  • Sam To Bam - Loss Of Data Or Just Great Compression?
    written by Clare Clare
    <p>I use picard tools to take a <a href='http://bio-bwa.sourceforge.net/'>BWA</a> alignment <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file and convert it into a sorted <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. Normally this works …
  • Samtools View Of A Sam File Entire Chromosome
    written by dfernan dfernan
    <p>Hi, </p> <p>I have (what I think) is a valid <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file.</p> <p>When I do:</p> <pre><code>samtools view -S F1i1F.C57.2.Aligned.out.sam chr19_C57 </code></pre> <p>I get the …
  • Basic Paired-End Sam Questions
    written by dlv • dlv
    <p>I can't seem to find answers to 2 very basic questions.</p> <p>For a paired-end <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a>, is there a separate <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> line for mate1's …
  • Error In Converting Sam To Bam By Samtools
    written by Raghav Raghav
    <p>hello every one,</p> <p>after mapping,i encountered an error during converting my <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> out file into <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> out put file.</p> <p>[cdac@nbri <a href='http://samtools.sourceforge.net/'>samtools</a>-0.1.19]$ <a …
  • Sample Sam File
    written by Steffi Steffi
    <p>Hi,</p> <p>I have a <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file which contain alignments of about 30 million reads. I just would like to randomly sample e.g. 1 million …
  • Filtering A Sam File For Quality Scores
    written by Varun Gupta Varun Gupta
    <p>Hi I have a <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file. The 5th column of the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file tells us the quality score. In my <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> alignment …
  • In Bowtie2'S Sam Output, What Does The Yt Attribute Mean?
    written by Eric Fournier Eric Fournier
    <p>In Bowtie2's <a href='http://samtools.sourceforge.net/SAM1.pdf'>SAM</a> output, what does the YT flag mean? I can't seem to find a mention of it in Bowtie2's documentation, and the …
  • Sam To Bam Conversion Reference Fasta File
    written by Atom Smasher • Atom Smasher
    <p>Hello,</p> <p>I had some data in "eland_export.txt" format and I want to convert it to the <a href='http://samtools.sourceforge.net/SAM1.pdf'>SAM</a> format first and then later to the …
  • Soap Result Converted To Bam Not Working In Samtool
    written by Sarwar • Sarwar
    <p>I converted the output of soap2 which is map file to <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file using map2sam.pl and then converted this <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> to <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> …