User Tools

Site Tools


archive:bioinformatic_tools:bwa

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
archive:bioinformatic_tools:bwa [2011/05/30 06:21]
karplus [Using BWA] Mentioning gzip capability.
archive:bioinformatic_tools:bwa [2011/06/08 18:22]
svohr
Line 88: Line 88:
  
 {{:​bioinformatic_tools:​run2_bc08_seqprep_histogram_r2.png|}} {{:​bioinformatic_tools:​run2_bc08_seqprep_histogram_r2.png|}}
 +
 +===== After Quake =====
 +This process was repeat on the data after Quake correction. In addition to the paired reads and merge reads, Quake separates the reads whose pair could not be successfully corrected. ​
 +
 +{{:​bioinformatic_tools:​run1_quake_histogram.png|}}
 +{{:​bioinformatic_tools:​run2_bc07_quake_histogram.png|}}
 +{{:​bioinformatic_tools:​run2_bc08_quake_histogram.png|}}
 +
 +==== Mean Lengths ====
 +
 +Run 1
 +| Template | 151.7 |
 +| Merged ​  | 113.3 |
 +| Single ​  | 65.7  |
 +
 +Run 2 bc07
 +| Template | 250.3 |
 +| Merged ​  | 142.5 |
 +| Single ​  ​| ​ 88.8 |
 +
 +Run 2 bc08
 +| Template |  94.5 |
 +| Merged ​  | 104.2 |
 +| Single ​  ​| ​ 61.3 |
 +
 +
 +===== After Assembly =====
 +This process was repeated using the [[bioinformatic_tools:​soapdenovo|SOAPdenovo]] contigs from ''​assemblies/​slug/​SOAPdenovo-assembly2/​k47_w_454_contigs''​ instead of the 454 reads as the reference. The shapes of these distributions follow the same patterns as the ones found using the 454 reads, although there are more mapped pairs because of the higher coverage of the contigs and we see some longer templates sizes than before. For the run1 and run2 barcode 7 pairs, the results involve some self-reference,​ since the previous estimates were used to build the contigs to which we mapped the pairs. However, the pairs for barcode 8 were not used in the assembly because of their negative mean
 +insert size (most pairs overlap but not in a way that SeqPrep could merge). When these were mapped to the contigs, we saw the same pattern as before but with more reads that mapped successfully.
 +
 +{{:​bioinformatic_tools:​run1_assembly_histogram.png|}}
 +
 +{{:​bioinformatic_tools:​run2_bc07_assembly_histogram.png|}}
 +
 +{{:​bioinformatic_tools:​run2_bc08_assembly_histogram.png|}}
 +
  
archive/bioinformatic_tools/bwa.txt · Last modified: 2015/09/04 09:06 by 68.180.228.52