Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bam file of RNA-seq

    Hi, forks

    When I import my RNA-seq in BAM format to the CLC genomic workbench, it automatically turns out to be two files with different file name. I thougth these two files are the same. But when I did reads-mapping of the two files separately, I found that total reads of two files are different, also the percentage of mapped reads varied with the same parameter settings.

    Very confused with these two files. Could anyone help to explain it? Thanks!

  • #2
    Are there different read groups, perhaps two different samples in the BAM? How was it produced (tool, command example, etc.)? Maybe you could show examples of SAM lines from each of the two files ...

    Comment


    • #3
      The bam file is generated by sequencing company. And they are indicated as one sample.
      Could you tell me how to show SAM lines from the two files?
      Thanks!

      Comment


      • #4
        If you're comfortable working on the command line (linux/unix, Mac OS X terminal?) and have samtools installed, use 'samtools view' ...
        If not, can you get helpful details about the aligned reads by hovering over features in CLC?

        Comment


        • #5
          I get samtools installed and it is a huge file.
          which command can be used to see the file page by page?
          I tried less | samtools view myfile.bam.
          It didnot work.

          Comment


          • #6
            Maybe CLC Bio partitions the files into paired and non-paired files. Could this be waht happened? (It does this with FASTQ files, after trimming and discarding bad reads)

            samtools view myfile.bam |less

            should work for you.
            Last edited by kopi-o; 04-12-2012, 01:05 PM.

            Comment


            • #7
              Thank you, kopi-o!

              Both of them are pointing to paired files. One is *** [no read group](paired), the other is ***(paired).

              Comment


              • #8
                OK, so maybe your BAM file contains some alignments for reads tagged with a read group (@RG) and some without a read group. It could be a merged file of alignments from two different instrument runs, for instance. Could that be it?

                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, 07-10-2024, 07:30 AM
                0 responses
                30 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 07-03-2024, 09:45 AM
                0 responses
                202 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 07-03-2024, 08:54 AM
                0 responses
                212 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 07-02-2024, 03:00 PM
                0 responses
                194 views
                0 likes
                Last Post seqadmin  
                Working...
                X