Seqlogo generation and interpretation
1
0
Entering edit mode
23 months ago
jaqx008 ▴ 110

seqlogo -F pdf -C 13 -c -S -B 2 -f male.13.fa -o seqlogo13.pdf

hello everyone. below (link1) is a seqlogo I generated using the command above. I am trying to interprete the image with regards to reads that begin with either a G or an A on the 5' end. I dont know what to deduce from this image since the first and second lines dont begin with either G or A. Can someone help? Also, I want to make a seqlogo figure like bellow (link2) but the command keep complaining about uneven length when i set the length to be anything above 13. Does anyone know how I can get the image bellow since I have the same seq data to make it? A recommendation of how I can interprete it would also be helpful.

for first figure I only provided a fasta file

>1-5183
TAGACTCCGACAC
>2-2196
TCGCAGACACCGA
>3-1786
TGACGTCGTAGAC
>4-1406
TGTCGTCGCAGAC
>5-1178
TAGACTCCGATAC
>6-1127
TAGACTCCGACAC

in attempt o make second figure, I supplied seqlogo with a fasta file containing all reads up to 20, 21, 22, 23 ... 30. hope this helps

Screen-Shot-2019-12-19-at-8-03-39-AM

seqlogo 5' ID polymerase • 407 views
ADD COMMENT
0
Entering edit mode

show us the nucleotide matrix you fed seqlogo

ADD REPLY
0
Entering edit mode

for first figure I only provided a fasta file

>1-5183
TAGACTCCGACAC
>2-2196
TCGCAGACACCGA
>3-1786
TGACGTCGTAGAC
>4-1406
TGTCGTCGCAGAC
>5-1178
TAGACTCCGATAC
>6-1127
TAGACTCCGACAC

in attempt o make second figure, I supplied seqlogo with a fasta file containing all reads up to 20, 21, 22, 23 ... 30. hope this helps

ADD REPLY
0
Entering edit mode
23 months ago

This is what I got using seqlogo from weblogo.2.8.2 your fasta file. Could you be just looking at the wrong output file from a previous run?

enter image description here

ADD COMMENT
0
Entering edit mode

what I pasted was actually a small portion of the fasta file. However, I see now that seq logo gives that were as weblogo gives close to what I am looking for. Thanks alot

ADD REPLY
0
Entering edit mode

in this case seqlogo is just the CLI version of weblogo. they should be identical given the same parameters and input.

ADD REPLY

Login before adding your answer.

Traffic: 2126 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6