SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Site Reliability Engineer, Disk & Compute Farm Complete Genomics Careers Industry Jobs! 0 10-01-2010 11:34 AM

Reply
 
Thread Tools
Old 11-20-2009, 07:01 AM   #1
Layla
Member
 
Location: London

Join Date: Sep 2008
Posts: 58
Default Reliability

I am comparing the output from Bowtie and maq and using the command below I get a total of 147-e6 alignments for 34-e6 reads. After rigorous filtering, I have half the numbers of reads mapped to chr1 as I did using maq. One of the issues is this:
HWI-EAS261:6:61:525:1199#0/2 - gi|89161185|ref|NC_000001.9|NC_000001 78 TAACCCTAACCCTAACCCTAACCCTAACCCAACCCTAACCCTAAC
HWI-EAS261:6:61:525:1199#0/2 + gi|89161205|ref|NC_000003.10|NC_000003 199384671 GTTAGGGTTAGGGTTGGGTTAGGGTTAGGGTTAGGGTTAGGGTTA

Maq flags the second alignment with an MAQ score<10 and is thus easy to eliminate maintaining the 1st alignment.
Is it possible to know which alignment is correct from Bowtie?

./bowtie -t h_sapiens_asm -1 sanger1.fq -2 sanger2.fq -n 2 -l 28 -e 70 -I 0 -X 250 -a > output.txt --un unmapped_reads.fq --al reads_at_least_1_alignment.txt
#Reported 73,823,285 paired-end alignments to 1 output stream(s)

Layla
Layla is offline   Reply With Quote
Old 11-20-2009, 07:25 AM   #2
Ben Langmead
Senior Member
 
Location: Baltimore, MD

Join Date: Sep 2008
Posts: 200
Default

Hi Layla,

If filtering of repetitive alignments is what you're after, take a look at the -m and --strata options. Together, these options allow you to define what hits count as "repetitive", then filter them out. The --max option additionally diverts the repetitively-aligned reads to a separate output file.

Also, I can't see how you're running Maq, but you might be in a situation where you'd actually like to run Bowtie in unpaired mode separately on your two mate inputs. I say this because the Maq output you show seems to have aligned the mates to separate chromosomes. Perhaps I'm misinterpreting.

Thanks,
Ben
Ben Langmead is offline   Reply With Quote
Old 11-20-2009, 08:37 AM   #3
Layla
Member
 
Location: London

Join Date: Sep 2008
Posts: 58
Default

Hi Ben,

Thanks for your reply. I've got myself a little muddled with the output. I used the -m1 -k1 option but I was worried I might be missing reads that had 1MM or 2MM as these would go unreported. Hence I asked everything to be reported and will clean the data myself. I will re-align the unmapped reads using unpaired option. The above command for 34-e6 reads (17-e6 reads per file) took 6hours on a 32GB machine, is this reasonable?

I might get back to you once I have looked at the output some more.

Thanks,
L
Layla is offline   Reply With Quote
Old 11-20-2009, 09:20 AM   #4
Layla
Member
 
Location: London

Join Date: Sep 2008
Posts: 58
Default

Quick question,

How does Bowtie deal with paired end reads where both ends align to different chromosomes? I believe maq gives various flags to identify such reads.

Thanks,
L
Layla is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 11:43 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO