SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
AVA-difference betweeen "computation" reads and "global align" reads CCBIO 454 Pyrosequencing 0 01-17-2014 04:13 AM
htseq-count: many reads are "no_feature". zzhao2 Bioinformatics 5 01-16-2014 09:23 AM
"not_aligment" from the output of HTseq-count flyingoyster Bioinformatics 5 09-03-2013 01:21 AM
HTseq count - all reads are "alignment_not_unique" dmsoanes Bioinformatics 2 07-12-2012 06:07 AM
SEQanswers second "publication": "How to map billions of short reads onto genomes" ECO Literature Watch 0 06-30-2009 12:49 AM

Reply
 
Thread Tools
Old 01-27-2014, 12:37 PM   #1
zzhao2
Member
 
Location: Dallas

Join Date: Aug 2012
Posts: 21
Default htseq-count: reads with letter "N"

Hi,
I've just found that a read starting with letter "N" in my sam file was marked as "no_feature" by htseq-count. I hacked that read by changing the "N" to the reference letter and also modifying the affected fields such as FLAG and MD, and got htseq-count worked with that read. My question is: does htseq-count discard any reads with "N"s in anywhere, or just those starting with "N"s? It looks like that in that read all other letters of the sequence were same with the reference, so I would rather keep that read.
I raised this question because by using FastQC I've found that many of my reads have "N"s in their tails (3'?). Will htseq-count output better counts if I trim those tails?

Thanks,
Sylvia

Last edited by zzhao2; 01-27-2014 at 12:41 PM.
zzhao2 is offline   Reply With Quote
Old 01-27-2014, 01:01 PM   #2
zzhao2
Member
 
Location: Dallas

Join Date: Aug 2012
Posts: 21
Default

Please ignore my initial post. The real reason why my read is not mapped is not the "N", it's because in the FLAG field it's marked as "0x10 SEQ being reverse complemented". When I change this single bit to 0, the read got counted. Does any body have any idea?
zzhao2 is offline   Reply With Quote
Reply

Tags
htseq-count

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 07:47 PM.


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