Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Wildcard Nucleotide Bases

    I have some pileup data from samtools and it is using wildcards to represent a combination of bases I'm unsure of.

    I know the 5 standard bases (ACGTU) and the well known pyrimidine and purine wildcards (YR), but don't know these other wildcards "SMWK". Anyone know what they definitively mean?

    Here's some sample pileup output:
    Code:
    NC_000913.2     82      T       W       24      24      60      19     .AAAA.....,,.,...,.     RBBBBBgBhgeBfahdffc
    NC_000913.2     3622    A       M       23      23      60      34      .,,C,C.CC,,C,,,.,,,,...g,,.g,,....      ]egBhBaBBhdBhcahdccfBhfBT^[Bb`hhcf
    NC_000913.2     6609    G       S       45      45      60      50      .$.,,.C,.,.,C.C,,CC.,..,C,.,.,.,,.,,.,,,.,,,c.c,,^~,^~, a_gh`BhYhZfBBBffBJhXgBdB`gfWYhcBYheBebdfaecBfBhfBB
    NC_000913.2     7073    C       Y       21      21      60      40      ,,,.,..,.,,.,,T..,.,.,....,.,,.t.g,ttt.g        ghhhgeQhVhg^ehBee\gBhBhW^hWhVBhBhBBBBBhB
    NC_000913.2     1258    G       R       26      26      60      36      ,.,A.,,ACC$,,,C....,,,,,,.,,,,...aat.   h\hB]ghBBBfhfB`cgc_a_a^afhffhdafBBBh
    NC_000913.2     3622    A       M       23      23      60      34      .,,C,C.CC,,C,,,.,,,,...g,,.g,,....      ]egBhBaBBhdBhcahdccfBhfBT^[Bb`hhcf
    NC_000913.2     13748   T       K       31      31      60      48      .,..,.G,GG,.,...,.,..,.,,......,,,.......,.g,ggg        dgeRhSBhBBhBhfBfhB]abBdOadba\ghWB_ghadchhghBfBBB

  • #2
    IUPAC ambiguity codes...



    S = strong (g/c)
    M = aMino (a/c)
    W = weak (a/t)
    K = keto (g/t)

    Comment


    • #3
      Exactly what I was looking for, thanks very much.

      Comment


      • #4
        Can anyone tell me how to expand this ambiguity codes with biopython?
        I mean, if I call with a sequence like 'YATCG' it returns 2 sequences 'TATCG' and 'CATCG'

        Thanks a lot.

        Comment


        • #5
          I got it. Thank you anyway.

          Code:
          from Bio.Data import IUPACData
          
          #-------------------------------------------
          #	Functions
          #-------------------------------------------
          
          def expand_wildcards(seq):
          	bl = [seq]
          	for ix in range(len(seq)):
          		sw = True
          		while sw:
          			iseq = bl.pop(0)
          			sl = list(iseq)
          			v = IUPACData.ambiguous_dna_values[sl[ix]]
          			if len(v) == 1:
          				sw = False
          				bl.append(iseq)
          			else:
          				for jx in range(len(v)):
          					sl[ix] = v[jx]
          					iseq = (''.join(sl))
          					bl.append(iseq)
          
          	bl.sort()	# no need indeed
          	return bl
          Last edited by GenoMax; 11-15-2015, 01:34 PM. Reason: Added CODE tags

          Comment


          • #6
            ^^^^^
            It is correctly indented but I do not know how to leave it well. It did not aloow me to upload it as attachment.

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Choosing Between NGS and qPCR
              by seqadmin



              Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
              10-18-2024, 07:11 AM
            • seqadmin
              Non-Coding RNA Research and Technologies
              by seqadmin




              Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

              Nobel Prize for MicroRNA Discovery
              This week,...
              10-07-2024, 08:07 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 11-01-2024, 06:09 AM
            0 responses
            18 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 10-30-2024, 05:31 AM
            0 responses
            18 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 10-24-2024, 06:58 AM
            0 responses
            24 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 10-23-2024, 08:43 AM
            0 responses
            53 views
            0 likes
            Last Post seqadmin  
            Working...
            X