Dear All,
I am trying to map strand-specific pair-end reads using Tophat v1.1.4, but in the result BAM file 'ISIZE' or 'TLEN' field is identical to ZERO all the time. Take some hits for example:
A8054YABXX:8:1:17161:19551#ATCACGAT 163 chr1 395709 255 100M = 395803 0 GCGCCAGTCGCAGAGGAAGATGGCGTCGCGAAGGCTGGGGTCGTTTCTGCGCGGGCTGCCGGCGCTGCTACTGCTGATGCTGTTCTTGCTTGGGCCCTGG GGGGGGGGGGGGGGGGBGGFDGGGGDFFFFDFGGGGGGGGAFFDFEEEEEGGGFGGGGDGGG?E?EFDAAEC=CAEBBDCEEAE=ABA?########### NM:i:0 NH:i:1 XS:A:+
A8054YABXX:8:4:14320:96976#ATCACGAT 163 chr1 395709 255 100M = 395803 0 GCGCCAGTCGCAGAGGAAGATGGCGTCGCGAAGGCTGGGGTCGTTTCTGCGCGGGCTGCCGGCGCTGCTACTGCTGATGCTGTTCTTGCTTGGGCCCCGG GGGGGGGEGGGGGDGGGFGGGGGFGAFFFFFFDFFFEGGGBEEAEFFFFFGGGEGG?GDG?DDDBBA?=CEBCAECB:A::??C5?B############# NM:i:1 NH:i:1 XS:A:+
The Run Command Line is:
I do not know what's wrong. And any help would be appreciated. Thanks ahead!
I am trying to map strand-specific pair-end reads using Tophat v1.1.4, but in the result BAM file 'ISIZE' or 'TLEN' field is identical to ZERO all the time. Take some hits for example:
A8054YABXX:8:1:17161:19551#ATCACGAT 163 chr1 395709 255 100M = 395803 0 GCGCCAGTCGCAGAGGAAGATGGCGTCGCGAAGGCTGGGGTCGTTTCTGCGCGGGCTGCCGGCGCTGCTACTGCTGATGCTGTTCTTGCTTGGGCCCTGG GGGGGGGGGGGGGGGGBGGFDGGGGDFFFFDFGGGGGGGGAFFDFEEEEEGGGFGGGGDGGG?E?EFDAAEC=CAEBBDCEEAE=ABA?########### NM:i:0 NH:i:1 XS:A:+
A8054YABXX:8:4:14320:96976#ATCACGAT 163 chr1 395709 255 100M = 395803 0 GCGCCAGTCGCAGAGGAAGATGGCGTCGCGAAGGCTGGGGTCGTTTCTGCGCGGGCTGCCGGCGCTGCTACTGCTGATGCTGTTCTTGCTTGGGCCCCGG GGGGGGGEGGGGGDGGGFGGGGGFGAFFFFFFDFFFEGGGBEEAEFFFFFGGGEGG?GDG?DDDBBA?=CEBCAECB:A::??C5?B############# NM:i:1 NH:i:1 XS:A:+
The Run Command Line is:
Code:
tophat-1.1.4.Linux_x86_64/tophat --mate-inner-dist 100 --min-anchor-length 8 --splice-mismatches 0 --max-multihits 100 --output-dir tophat-114_out --num-threads 6 --solexa1.3-quals --coverage-search --microexon-search --library-type fr-firststrand --segment-mismatches 2 --segment-length 25 bowtie_index L1_1.fq L1_2.fq
Comment