Originally posted by dpryan
View Post
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
-
The bowtie2 commands I specified were these:
bowtie2-build ref base_name
bowtie2 --sensitive-local -x base_name -1 fq1 -2 fq2 -S sam
I'm aware that Nterm and Cterm have particular meanings, I'm just not too worried about it at this point. Once I get past this section in my testing I'll make sure the orientation is taken into account.
Leave a comment:
-
A couple points. Firstly, what was the exact bowtie2 command that you specified? Depending on the flags you used, the unmapped reads may or may not even be included. Secondly, I hope you take the orientation of your genes into account when you start calling the files Nterm and Cterm. They'll be swapped if the gene is on the - strand.
BTW, the most straight forward solution is to just find a couple reads like this from the bwa alignment and use grep to see what bowtie2 did with them. Presuming they map the same, you'll then know if bowtie2 is doing something unexpected with the flags.
Leave a comment:
-
Unable to find flag in SAM with bowtie2 - but can with BWA
I'm currently using Bowtie2 to map my reads to a small region of a genome to find what region it's in by collecting the unmapped pairs of the reads which have mapped.
I'm currently using samtools flags to get the reads I want:
---a--- ......... -------------
|------------------------gene----------------|
-------..............----b-----
What we want is
a: retain unmapped reads whose mate is on reverse strand (-f 36)
b: retain unmapped reads (-f 4) excluding reads whose mate is on reverse strand or unmapped (-F 40)
samtools view -Sb -f 36 sam > NtermBam
samtools view -Sb -f 4 -F 40 sam > CtermBam
Previously when mapping with BWA both of these commands would work and I would get an N terminus bam and a C terminus bam. However, when I map with Bowtie2 (because the mapping itself is better), I get no reads in the N terminus bam and all the reads go to the C terminus bam.
Does anyone know why this might be? I'm pretty sure the samtools commands are right because it works with BWA. I'm thinking it has something to do with Bowtie2, maybe not recording flags correctly?Tags: None
Latest Articles
Collapse
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
-
by seqadmin
Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...-
Channel: Articles
03-22-2024, 06:39 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
31 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
33 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
28 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
53 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Leave a comment: