Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Counting mapped nucleotides in a bam file

    Hi !

    Is there a way to go over a sorted bam file and count how many different matches there are for a specific nucleotide - for example for a reference nucleotide T in position x - there are :12 T , 3A, 4C 0G..?\

    Thanks ahead!

  • #2
    if it's just for some nucleotides and you need no autmoation:
    IGV shows that information when you visualize it there's a coverage track. y moving the cursor over the relevant bar you'll get that information

    Comment


    • #3
      if it's just for some nucleotides and you need no autmoation:
      IGV shows that information when you visualize it there's a coverage track. y moving the cursor over the relevant bar you'll get that information

      Comment


      • #4
        Thanks for your answer!

        I saw the anomaly in the IGV for couple of genes and now I would like to check it in about 500 genes... so I have to find a way to do it automatically ...

        Comment


        • #5
          What about traditional SNP calling (for example using GATK). Is there any reason not to do it?

          Comment


          • #6
            As far as I know snp calling won't always count a snp like one if there is not enough coverage or not sufficient quality score...
            In addition, when I used a snp caller I got the consensus nucleotide and since I would like to see how much of each nucleotide covers specific position I am not sure if it fits...

            Do you know how to overcome these issues using snp callers?
            Thanks

            Comment


            • #7
              you can set the criteria quite low and when using GATK (don't know for the others) it outputs the allele balance between the reference allele and the substituted allele.

              It might miss some snps though if evidence is very low but I guess you do search variants with biological meaning and I guess he wouldn't miss those.

              However, you would only get information of balance of 2 nucleotides. I don't think triallelic SNPs are supported yet.

              Can't think of any other solution but maybe there's a fuction in bamtools/bedtools/samtools/picard that you might use

              Comment


              • #8
                Thanks for your help!

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Exploring the Dynamics of the Tumor Microenvironment
                  by seqadmin




                  The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
                  07-08-2024, 03:19 PM
                • seqadmin
                  Exploring Human Diversity Through Large-Scale Omics
                  by seqadmin


                  In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
                  06-25-2024, 06:43 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, Yesterday, 05:49 AM
                0 responses
                15 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 07-15-2024, 06:53 AM
                0 responses
                26 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 07-10-2024, 07:30 AM
                0 responses
                37 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 07-03-2024, 09:45 AM
                0 responses
                204 views
                0 likes
                Last Post seqadmin  
                Working...
                X