Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • lbthrice
    replied
    whoops, should have checked that...thanks!

    Leave a comment:


  • gringer
    replied
    Yeah, mpileup doesn't generate lines for regions it has no reads on. If you need to work out whole-region coverage, then you should be looking at the column that gives the base number. For covered regions, this will increase by 1 per line (assuming no inserts), but will jump multiple bases when there are no reads.

    Leave a comment:


  • swbarnes2
    replied
    The simplest answer is that you just don't have any reads across part of your region.

    Try getting the .sam file for the same region +- 500 bases, or try looking at the whole .bam in IGV. You can zoom in on your region.

    Leave a comment:


  • lbthrice
    started a topic mpileup: specified region truncated.

    mpileup: specified region truncated.

    Hello gentlepeople,

    I am using the samtools 'mpileup' program.
    I wish to generate a pileup file for a specific region.
    I have several separate .bam input files that I am generating the pileups for.
    Unfortunately, mpileup truncates the region and I only get data for the first part of the region.
    I have ran the same command on different .bam files and get different results.

    When I command:

    samtools mpileup -f 'hg19.fa' -r chr6:27114408-27115845 'experiment01.bam' > out_PILEUP2.pile
    [mpileup] 1 samples in 1 input files
    <mpileup> Set max per-file depth to 8000

    I get a file with 746 lines.
    I am expecting 1438 lines because that is the size of the interval I have specified with the '-r' flag
    If I use mpileup with the same region but with a different .bam file I get a different size file.

    samtools mpileup -f 'hg19.fa' -r chr6:27114408-27115845 'experiment02.bam' > out_PILEUP3.pile
    [mpileup] 1 samples in 1 input files
    <mpileup> Set max per-file depth to 8000

    I get a file with 623 lines.
    Again, I am expecting 1438 lines because that is the size of the interval I have specified with the '-r' flag.

    The first line in both files is for chr6:27114408, as expected.
    BUT, the files are truncated at different positions.
    With experiment01.bam the final line of the pileup file is chr6:27115577.
    With experiment02.bam the final line of the pileup file is chr6:27115344.

    This must have something to do with the .bam file.
    Can someone tell me which parameter needs to be adjusted?
    I have consulted the manual but I could not identify anything that sounded applicable.

    Thanks for you time,
    Lionel (Lee) Brooks 3rd
    Dartmouth Genetics Grad Student

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-19-2024, 07:20 AM
0 responses
25 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-16-2024, 05:49 AM
0 responses
41 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-15-2024, 06:53 AM
0 responses
46 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-10-2024, 07:30 AM
0 responses
42 views
0 likes
Last Post seqadmin  
Working...
X