Uniquely mapped reads on Bowtie for colorspace RNA-seq
0
0
Entering edit mode
2.9 years ago
IsiI • 0

Hi! I mapped RNA-seq of Synechococcus sp PCC7002 and used this code:

bowtie -p 2 -C -S SynechococcusPCC7002 SRR308187.fastq > $SRR308187.sam

while it was mapping it gave me this information:

SRR308187

-reads processed: 26949795

-reads with at least one reported alignment: 13593071 (50.44%)

-reads that failed to align: 13356724 (49.56%)

Reported 13593071 alignments

No other file apart from SRR308187.sam were created. How can I get just the uniquely mapped reads?

I tried these codes:

samtools view -c -q 1 SRR308187.bam #(it gave me the number of reported alignments)

samtools view -bq 1 SRR308187.bam #(couldn't read it)

grep 'XS:' testfile.sam \
  | wc -l 
#(suppoused to give the number of alignments with more than 1 best align, but don't know if I have to just substract this number to the amount of reported alignments)

please help!

uniquely bowtie RNA-seq samtools • 574 views
ADD COMMENT

Login before adding your answer.

Traffic: 3389 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6