Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • swbarnes2
    replied
    Thanks for posting the answer on the thread.

    Leave a comment:


  • Michael.James.Clark
    replied
    Pretty straightforward answer, but I'll post it anyway in case anyone else encounters this in the future (not that there are that many people out there dealing with both Illumina and SOLiD, but here you go).

    It was indeed the indexes. When indexing the first time, I indexed normal and colorspace in the same folder, colorspace second, using default output. It seems some of the resulting indexes, therefore, overwrite. Of course the colorspace indexes don't work with Illumina data.

    Second time around, I indexed them with different names (in different folders, actually), and now things are aligning beautifully.

    Leave a comment:


  • Michael.James.Clark
    replied
    Thanks Jon!

    I think I do have those. Here's my whole (top secret) reference folder:
    Code:
    lrwxrwxrwx 1 mjclark rpm   56 Jan 18 22:02 human_g1k_v37.dict -> ../GATK/human_g1k_v37.dict
    lrwxrwxrwx 1 mjclark rpm   57 Jan 18 22:02 human_g1k_v37.fasta -> ../GATK/human_g1k_v37.fasta
    -rw-r--r-- 1 mjclark rpm 6.5K Feb 28 20:15 human_g1k_v37.fasta.amb
    -rw-r--r-- 1 mjclark rpm 6.7K Feb 28 20:15 human_g1k_v37.fasta.ann
    -rw-r--r-- 1 mjclark rpm 1.1G Feb 28 21:04 human_g1k_v37.fasta.bwt
    lrwxrwxrwx 1 mjclark rpm   61 Jan 18 22:02 human_g1k_v37.fasta.fai -> ../GATK/human_g1k_v37.fasta.fai
    -rw-r--r-- 1 mjclark rpm 6.5K Feb 28 20:14 human_g1k_v37.fasta.nt.amb
    -rw-r--r-- 1 mjclark rpm 6.7K Feb 28 20:14 human_g1k_v37.fasta.nt.ann
    -rw-r--r-- 1 mjclark rpm 740M Feb 28 20:14 human_g1k_v37.fasta.nt.pac
    -rw-r--r-- 1 mjclark rpm 740M Feb 28 20:15 human_g1k_v37.fasta.pac
    -rw-r--r-- 1 mjclark rpm 1.1G Feb 28 21:05 human_g1k_v37.fasta.rbwt
    -rw-r--r-- 1 mjclark rpm 740M Feb 28 20:15 human_g1k_v37.fasta.rpac
    -rw-r--r-- 1 mjclark rpm 370M Feb 28 21:22 human_g1k_v37.fasta.rsa
    -rw-r--r-- 1 mjclark rpm 370M Feb 28 21:13 human_g1k_v37.fasta.sa
    -rwxr--r-- 1 mjclark rpm 6.1G Oct 15 16:35 human_g1k_v37.nix
    Maybe it's that I did the colorspace indexing in the same folder. I'll try re-doing it seperate from one-another in lieu of another idea.

    Leave a comment:


  • Jon_Keats
    replied
    Looks like some of the index files are missing. This is an example of what I have in my bwa index directory.

    Code:
    -rw-r--r-- 1 jkeats domainuser 3142044949 Feb  3 16:23 hg18.fasta
    -rw-r--r-- 1 jkeats domainuser       6152 Feb  3 16:23 hg18.fasta.amb
    -rw-r--r-- 1 jkeats domainuser        946 Feb  3 16:23 hg18.fasta.ann
    -rw-r--r-- 1 jkeats domainuser 1155163564 Feb  3 16:23 hg18.fasta.bwt
    -rw-r--r-- 1 jkeats domainuser  770109014 Feb  3 16:23 hg18.fasta.pac
    -rw-r--r-- 1 jkeats domainuser 1155163564 Feb  3 16:23 hg18.fasta.rbwt
    -rw-r--r-- 1 jkeats domainuser  770109014 Feb  3 16:23 hg18.fasta.rpac
    -rw-r--r-- 1 jkeats domainuser  385054532 Feb  3 16:23 hg18.fasta.rsa
    -rw-r--r-- 1 jkeats domainuser  385054532 Feb  3 16:23 hg18.fasta.sa

    Leave a comment:


  • Michael.James.Clark
    started a topic BWA, mostly unmapped reads

    BWA, mostly unmapped reads

    Hi all,
    Hoping someone may know exactly what I did wrong off the bat here since I think we have a lot of BWA gurus here. This is my first time using BWA. Previously I've used Novoalign on the same exome-seq data to great success, aligning the majority of the reads. So I was surprised after running BWA on it that less than 1% of the data mapped and most of it was unmapped.

    The data in question are single lanes of HiSeq human exome-seq data.

    I indexed the reference genome:
    Code:
    bwa index -a bwtsw human_g1k_v37.fasta
    That created (in the same folder):
    Code:
    human_g1k_v37.fasta.amb
    human_g1k_v37.fasta.ann
    human_g1k_v37.fasta.pac
    human_g1k_v37.fasta.rpac
    (I also indexed for colorspace in the same directory since I have SOLiD data I need to align in a few days.)

    Then I ran BWA as follows:
    Code:
    $bwa aln -t 8 $ref $f1 > $out.aln_sa1.sai
    $bwa aln -t 8 $ref $f2 > $out.aln_sa2.sai
    $bwa sampe -r "$rg" $ref $out.aln_sa1.sai $out.aln_sa2.sai $f1 $f2 > $out.sam
    There were no errors while it ran except it mapped almost nothing.

    Can anyone see a glaring problems in my commands here that would lead to tons of unmapped reads? Any help appreciated!
    Last edited by Michael.James.Clark; 03-02-2011, 12:45 PM.

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