View Single Post
Old 02-01-2011, 04:41 PM   #1
khunny7
Junior Member
 
Location: UT

Join Date: Sep 2010
Posts: 2
Default Different starting position in Genbank blast result?

I am new to this field and this probably a very ignorant question.
However, when I did a blast search with following sequence, "TGTCTTTGGACATGTAAGAATTGGAGGAAAATAAATGTGGATTTGGGAAACTTTGAGG" blast returned me this following result. As seen in the result, matched sequences are in same chromosome yet seemed to have different indices. Thus thinking those are repeated regions, I was trying to locate those. However, I could find only one not even two. My assumption is that those numbers are not really indices. Can anyone help me to understand this problem? Thank you.


>ref|NT_167247.1| Homo sapiens chromosome 6 genomic contig, GRCh37.p2 reference
assembly alternate locus group ALT_REF_LOCI_5
Length=4833398


Features in this part of subject sequence:
large proline-rich protein BAT2

Score = 108 bits (58), Expect = 7e-22
Identities = 58/58 (100%), Gaps = 0/58 (0%)
Strand=Plus/Plus

Query 1 TGTCTTTGGACATGTAAGAATTGGAGGAAAATAAATGTGGATTTGGGAAACTTTGAGG 58
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct 2970579 TGTCTTTGGACATGTAAGAATTGGAGGAAAATAAATGTGGATTTGGGAAACTTTGAGG 2970636


>ref|NT_167245.1| Homo sapiens chromosome 6 genomic contig, GRCh37.p2 reference
assembly alternate locus group ALT_REF_LOCI_3
Length=4610396


Features in this part of subject sequence:
large proline-rich protein BAT2

Score = 108 bits (58), Expect = 7e-22
Identities = 58/58 (100%), Gaps = 0/58 (0%)
Strand=Plus/Plus

Query 1 TGTCTTTGGACATGTAAGAATTGGAGGAAAATAAATGTGGATTTGGGAAACTTTGAGG 58
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct 2876461 TGTCTTTGGACATGTAAGAATTGGAGGAAAATAAATGTGGATTTGGGAAACTTTGAGG 2876518


>ref|NT_113891.2| Homo sapiens chromosome 6 genomic contig, GRCh37.p2 reference
assembly alternate locus group ALT_REF_LOCI_2
Length=4795371


Features in this part of subject sequence:
large proline-rich protein BAT2

Score = 102 bits (55), Expect = 3e-20
Identities = 57/58 (99%), Gaps = 0/58 (0%)
Strand=Plus/Plus

Query 1 TGTCTTTGGACATGTAAGAATTGGAGGAAAATAAATGTGGATTTGGGAAACTTTGAGG 58
|||||||||||| |||||||||||||||||||||||||||||||||||||||||||||
Sbjct 3100500 TGTCTTTGGACACGTAAGAATTGGAGGAAAATAAATGTGGATTTGGGAAACTTTGAGG 3100557
khunny7 is offline   Reply With Quote