SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Read counts from SAM file mapped to de novo assembled transcripts using HTSeq-count alan_sm RNA Sequencing 2 06-12-2015 09:54 PM
Sorted BAM read count >2x total FASTQ count pkMyt1 Bioinformatics 8 05-13-2015 10:43 AM
Tools that report read count AND read names that map to genomic features. foolishbrat Bioinformatics 1 02-05-2014 01:21 AM
fastq.gz stats READ-COUNT BASE-COUNT jgibbons1 Bioinformatics 9 10-30-2013 06:24 AM
multiBamCov or htseq-count to count read per feature ? NicoBxl Bioinformatics 1 07-03-2012 03:05 AM

Reply
 
Thread Tools
Old 05-13-2015, 04:00 AM   #1
vinumanikandan
Junior Member
 
Location: cochin

Join Date: Feb 2011
Posts: 4
Default views Bowtie2 : aligned result read count not matching with sam Read count

bowtie2 -f -N 1 -p 8 -x REF.fa -1 sample01_1.fasta -2 sample01_2.fasta -S sample01.sam --al sample01_al.fasta --al-conc sample01_alcon.fasta


9260783 reads; of these:
9260783 (100.00%) were paired; of these:
9260311 (99.99%) aligned concordantly 0 times
472 (0.01%) aligned concordantly exactly 1 time
0 (0.00%) aligned concordantly >1 times
----
9260311 pairs aligned concordantly 0 times; of these:
43 (0.00%) aligned discordantly 1 time
----
9260268 pairs aligned 0 times concordantly or discordantly; of these:
18520536 mates make up the pairs; of these:
18520180 (100.00%) aligned 0 times
356 (0.00%) aligned exactly 1 time
0 (0.00%) aligned >1 times
0.01% overall alignment rate


Result:

Number of reads in sample01_alcon.1.fasta: 472
Number of reads in sample01_alcon.2.fasta: 472
Number of reads in sample01_al.fasta:0


To find number of reads matching a Particular Reference from sam file I run following :

perl -ne ' if (/^\@SQ/) { @F = split(/\t|:/, $_); print $F[2]."\n" } ' sample01.sam > agi_list.txt
perl -ne ' chomp($_); print $_."\t".`grep -c "\t$_" sample01.sam ` ' agi_list.txt >1a.counts


I have Couple of questions:

1) So In theory I should get (2*aligned concordantly)+(2*discordantly) =(2*472)+(2*43)=1030 . Right? but I got another number:1742 which is like (2*472)+(2*43)+(2*356)=1742.

2) I could not find the 43 aligned discordantly sequences and the 356 sequences in my alinged file why is it so and how I could find it?

3) What I need to do is calculate RPKM and also retrieve those number of sequences


What I need to do is calculate RPKM and also retrieve those number of sequences

Last edited by vinumanikandan; 05-13-2015 at 04:05 AM.
vinumanikandan is offline   Reply With Quote
Reply

Tags
bowtie 2, bowtie alignment stats, rpkm, sam file

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 10:10 PM.


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