According to TopHat documentation ( for --library-type option) if a strand-specific protocol was used in paired-end sequencing there should be an XS attribute tag assigned to every read. Right - I got my reads in FASTA format - where is my XS tag? If it's not there, where would it normally be? First four reads from raw FASTA file:
Code:
@A:3:1101:1441:2249#CTTGCAAC/1 CGCTGGGGCACAGTGACAGGAAACTGAGCAAAACAGCTCGTTGGCACTGATATGCCCCTCCCACCTCCCAGCTCCTGACTCCTATGGGCTGTTGGCTGCC +A:3:1101:1441:2249#CTTGCAAC/1 bbbeeeeegggggfghiiiihhhhghdghhhiiihiiidgghghiihhifggfhcegfhiggggeeeebdddbbcbbbcc]bR]b]acbccbYbc^ @A:3:1101:1593:2030#CTTGCAAC/1 NTGGAAGCAGCAGGTGGAATCCAAGAATATGCCCTTTCAGGATGGCCAAGAATTTGAACTGAGCATCTCAGTGCTGCCAGATAAGTACCAGGTAATGGTC +A:3:1101:1593:2030#CTTGCAAC/1 BP\c`acae^eegdffdeabgfdgddgff^egfghPI^cafhhhhZacfgfhghfghggbe_dghbfdgdgdgabZ^`bbbbbbabbbcbbb_
Comment