Thanks for replies.
I took a look at SAM format specification. It was more clear to me than bowtie2 manual.
I am now testing bits 0x1, 0x2 and 0x8 from flags and if they are on I am determining that current match is R2 and previous match is R1 (of the same pair).
Hopefully it works correctly now.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
As long as we're talking best practices, I would definitely pipe the output through samtools and save a BAM rather than a SAM file. There's no reason a SAM file should ever exist on a hard drive.
Leave a comment:
-
Originally posted by GenoMax View PostIt appears that sampie omitted "-S filename_to_hold_alignments.sam" option on the command line to write the alignments to a file.
-S <hit>
File to write SAM alignments to. By default, alignments are written to the "standard out" or "stdout" filehandle (i.e. the console).
Leave a comment:
-
It appears that sampie omitted "-S filename_to_hold_alignments.sam" option on the command line to write the alignments to a file.
Leave a comment:
-
Bowtie2 should be outputting BAM/SAM files. Read up on that format. If you have further questions after reading about SAM files -- especially the bitwise flag field -- then ask a question.
Leave a comment:
-
Bowtie2 output with paired end reads
Hi,
I am trying to use bowtie2 to align paired end Illumina read files.
I am using following command:
./bowtie2 -a --threads 8 --n-ceil 42 --very-sensitive-local --quiet --no-sq --no-hd --np 0 --ma 1 --mp 0,0 --rdg 1,1 --rfg 1,1 --score-min L,0,0.95 -x rdp/bt2_index/rdp.db -q -1 sR1.fastq -2 sR2.fastq > bt2out.txt
The index is made from RDP II database and currently sR1.fastq and sR2.fastq contain only one read:
sR1.fastq:
@1_M01363:14:000000000-A5K8J:1:1101:8617:1211 1:N:0:16
TACGGAGGATCCGAGCGTTATCCGGA...
sR2.fastq:
@2_M01363:14:000000000-A5K8J:1:1101:8617:1211 2:N:0:16
TTCACCGTTGCCGGCGTACTCCCCAGGTGGA
The bt2out.txt file shows mostly results with sR1.fastq read header.
1_M01363:14:000000000-A5K8J:1:1101:8617:1211....
1_M01363:14:000000000-A5K8J:1:1101:8617:1211....
1_M01363:14:000000000-A5K8J:1:1101:8617:1211....
.
.
.
I am not sure how to read bowtie2 results. How can I obtain information how each matching read aligned as a pair?
ThanksTags: None
Latest Articles
Collapse
-
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...-
Channel: Articles
07-08-2024, 03:19 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
Yesterday, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
160 views
0 likes
|
Last Post
by seqadmin
07-19-2024, 07:20 AM
|
||
Started by seqadmin, 07-16-2024, 05:49 AM
|
0 responses
127 views
0 likes
|
Last Post
by seqadmin
07-16-2024, 05:49 AM
|
Leave a comment: