Did you align the second set of fastq sequences? All I see is
bwa align ref fq1 > 1.sai
but not
bwa align ref fq2 > 2.sai
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
.sams are human readable. So start by looking at those.
You did index the genome too, right?
Leave a comment:
-
newbie BWA alignment question
another newbie question - I am trying to use BWA to align paired end reads to human ref genome.
here is the sequence of what i did -
1) alignment
bwa align ref fq1 > 1.sai
2) make sam
bwi sampe ref 1.sai 2.sai fq1 fq2 > .sam
3) convert to bam
samtools view -S -b .sam -o .bam
4) sort bam
samtools sort .bam .sorted
5) index bam
samtools index .sorted
6) check stats
samtools flagstat .bam > .stat
everything finished running without any error, but when i check the .stat file, this is what i see
0 + 0 duplicates
0 + 0 mapped (0.00%:nan%)
61091524 + 0 paired in sequencing
30545762 + 0 read1
30545762 + 0 read2
0 + 0 properly paired (0.00%:nan%)
0 + 0 with itself and mate mapped
0 + 0 singletons (0.00%:nan%)
0 + 0 with mate mapped to a different chr
0 + 0 with mate mapped to a different chr (mapQ>=5)
in every sample, there is very few to no aligned reads. can anyone comment on what could possibly be the issue here and if there is a way to check the output after each stage in this process. thanks.Tags: None
Latest Articles
Collapse
-
by seqadmin
Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.
Long-Read Sequencing
Long-read sequencing has...-
Channel: Articles
Today, 01:49 PM -
-
by seqadmin
The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben MartÃnez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...-
Channel: Articles
11-06-2024, 07:24 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Today, 09:29 AM
|
0 responses
13 views
0 likes
|
Last Post
by seqadmin
Today, 09:29 AM
|
||
Started by seqadmin, Today, 09:06 AM
|
0 responses
11 views
0 likes
|
Last Post
by seqadmin
Today, 09:06 AM
|
||
Started by seqadmin, Today, 08:03 AM
|
0 responses
11 views
0 likes
|
Last Post
by seqadmin
Today, 08:03 AM
|
||
Started by seqadmin, 11-22-2024, 07:36 AM
|
0 responses
61 views
0 likes
|
Last Post
by seqadmin
11-22-2024, 07:36 AM
|
Leave a comment: