Daniel:Notebook/GenomeMiner/2013-9-12: Difference between revisions
Jump to navigation
Jump to search
>Djacobse (Created page with "=Mock HL155 (Started 9/9/2013)= Back to Calendar ==Bowtie on Perfect Mock Sequence Data with Perfect...") |
>Djacobse |
||
Line 14: | Line 14: | ||
'''bowtie2''' | '''bowtie2''' | ||
bowtie2 -x bbuild_v4s1 -U v4s1_mockseq_perfseqbqual.fq --phred33 -S mockseq_v4s1_perfseqbqual.sam | bowtie2 -x bbuild_v4s1 -U v4s1_mockseq_perfseqbqual.fq --phred33 -S mockseq_v4s1_perfseqbqual.sam | ||
'''Alignment Results''' | |||
2959000 reads; of these: | |||
2959000 (100.00%) were unpaired; of these: | |||
0 (0.00%) aligned 0 times | |||
2953000 (99.80%) aligned exactly 1 time | |||
6000 (0.20%) aligned >1 times | |||
100.00% overall alignment rate | |||
So it looks like even with indisputable base quality and sequences, about 0.2% of sequences are so close that bowtie2 thinks they can be several sequences. |
Revision as of 18:23, 12 September 2013
Mock HL155 (Started 9/9/2013)
Bowtie on Perfect Mock Sequence Data with Perfect Base Qualities
Based on the results from yesterday I'm seeing what kind of alignment bowtie2 does on perfect sequences with perfect base qualities. Yesterday's sequences were direct copies of the .fa file, but they used base qualities randomly picked from the original sequencing data. These base qualities are all I (an extremely high quality score for phred33).
Generated using MATLAB (MockHL155_Master, Switch 3)
scp to Genome Miner
scp V4S1_mockseq_perfseqbqual.fq djacobse@132.239.135.41:/media/LTS_15T/DEJ_LTS/SeqStore/130628_HL155/mockseq/
bowtie2
bowtie2 -x bbuild_v4s1 -U v4s1_mockseq_perfseqbqual.fq --phred33 -S mockseq_v4s1_perfseqbqual.sam
Alignment Results
2959000 reads; of these: 2959000 (100.00%) were unpaired; of these: 0 (0.00%) aligned 0 times 2953000 (99.80%) aligned exactly 1 time 6000 (0.20%) aligned >1 times 100.00% overall alignment rate
So it looks like even with indisputable base quality and sequences, about 0.2% of sequences are so close that bowtie2 thinks they can be several sequences.