I have to find matching of short reads of a .fq file onto human genome. I am using BowTie for that.
Here is one read of my input file:
And BowTie outputs:
And
Can anybody inform me: what does this output means ? I mean: specially + 10 50100 and 1.
Here is one read of my input file:
Code:
@10_71499258_71499890_0_1_0_0_1:0:0_2:0:0_0/1 TATTAGATCGTGTGATTATATTTGACAGGTCTTAATTGACGCGCTGTTCAGCCCTTTGAGTTCGGTTGAGTTTTGGGTTGGAGAATTTTCTTCCACAAGG + 2222222222222222222222222222222222222222222222222222222222222222222222222222222222222222222222222222
Code:
10_71499258_71499890_0_1_0_0_1:0:0_2:0:0_0/1 + 10 50100 TATTAGATCGTGTGATTATATTTGACAGGTCTTAATTGACGCGCTGTTCAGCCCTTTGAGTTCGGTTGAGTTTTGGGTTGGAGAATTTTCTTCCACAAGG 2222222222222222222222222222222222222222222222222222222222222222222222222222222222222222222222222222 1
Code:
10_71499258_71499890_0_1_0_0_1:0:0_2:0:0_0/1 + 18 4514 TATTAGATCGTGTGATTATATTTGACAGGTCTTAATTGACGCGCTGTTCAGCCCTTTGAGTTCGGTTGAGTTTTGGGTTGGAGAATTTTCTTCCACAAGG 2222222222222222222222222222222222222222222222222222222222222222222222222222222222222222222222222222 1
Comment