[Biojava-l] Parsing a BLAST file

David Waring dwaring@u.washington.edu
Fri, 2 Nov 2001 13:19:04 -0800


This is a multi-part message in MIME format.

------=_NextPart_000_0000_01C163A0.F1D463F0
Content-Type: text/plain;
	charset="us-ascii"
Content-Transfer-Encoding: 7bit

I have just been doing this type of thing. I like Keith's ssbind classes. I
have attached a demo. Just run it as 'java BlastTest blast.out
blastInput.fasta'.

This parses the blast file and builds SequenceDBSearchResults into a list.
It is a little bit more complicated than that really. But this complication
gives very great functionality. The SearchResultBuilder must have two things
that you might not expect, a SequenceDB with all the query sequences that
blast was called with, and a SequenceDBInstallation which contains a
SequenceDB with the same name as that found in the blast output file, in the
demo this is 'genome'. With these things in place you can get both the
subject, and query sequences of any hit from the SequenceDBSearchResult. I
included a little sample of how to do this below since it is not in the
demo.

But, you say, this is a blast against some foreign database, How can I have
a sequencDB with all this data. The truth is you do not really need it. You
just need an empty SequenceDB with the correct name inside your
SequenceDBInstallation. But then of course you can not get the subject
sequences from the search result.

If anyone is interested I have been working on a tool that will run blast
against a local blast database, and give results back as
SequenceDBSearchResults. It makes a system call to your local blast program
to do this. Because of this it requires a little environment specific
configuration. I have it running on Unix, but it will probably run on
Windows too. It is just about finished. Because of the system specific
issue, I have not planned on putting this into biojava, but if anyone is
interested in it, let me know.

David

example of getting subject and query sequences from a SequenceDBSearchResult
( this is typed in from scratch so I don't guarantee it will complile :)

for (int j = 0; i < resultList.size(); i++){
    SequenceDBSearchResult r = (SequenceDBSearchResult)resultList.get(i);
    Sequence querySeq = (Sequence)r.getQuerySequence();
    SequenceDB subjectDB = r.getSequenceDB();

    hits = r.getHits();
    hitsI = hits.iterator();
    hit = (SeqSimilaritySearchHit)hitsI.next();
    String hid = hit.getSequenceID();
    Sequence subjectSeq = subjectDB.getSequence(hid);
}
---------------------------------------------------------------

> -----Original Message-----

> From: biojava-l-admin@biojava.org [mailto:biojava-l-admin@biojava.org]On
> Behalf Of Susan Glass

>
> I have just downloaded the Biojava package and made a BLAST
> output file parser according to the tutorial on
> http://biojava.org/tutorials/blastlikeParsingCookBook/index.html
>
> That means I have written my own XML content handler and put the
> parsing logic into startElement, using the
> BlastLikeDataSetCollection DTD file as a guide to the tag names.
>
> This works fine, but I believe I should be able to use objects in
> the org.biojava.bio.search package to avoid writing my own
> content handler.  Is this true? If so, could someone please point
> me to some example code (or other help) that starts with a BLAST
> output file, and ends with hit objects?  Unfortunately I'm having
> trouble choosing the correct classes to use with only the
> javadocs as a guide.  I found a mail list posting from August
> that deals with
s( 
> http://biojava.org/pipermail/biojava-l/2001-August/001480.html ), 


------=_NextPart_000_0000_01C163A0.F1D463F0
Content-Type: application/x-tar;
	name="blastTest.tar"
Content-Transfer-Encoding: quoted-printable
Content-Disposition: attachment;
	filename="blastTest.tar"

BlastTest.java=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00000744 =00023547 =
=00000017 =0000000011604 07370602146 014477 =
=000=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00ustar=0000dwaring=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00users=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00000000 =00000000 =
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00/* BlastTest
*  a demo of how to use the BlastLikeSaxParser and ssbind tools to parse =
blast results
*/

import java.io.*;
import java.util.*;

import org.biojava.utils.*;
import org.biojava.bio.seq.*;
import org.biojava.bio.seq.io.*;
import org.biojava.bio.seq.db.*;
import org.biojava.bio.symbol.*;
import org.biojava.bio.search.*;
import org.biojava.bio.program.ssbind.*;
import org.biojava.bio.program.sax.*;
import org.xml.sax.*;



public class BlastTest{
    public static void main(String[] args) throws Exception {
        List resultList;
        org.xml.sax.XMLReader parser;
        SeqSimilarityAdapter cHandler;
        BlastLikeSearchBuilder sBuilder;
        FileInputStream fStream;
        BufferedReader br;
        SequenceIterator si;
        SequenceDB qDb;
        SimpleSequenceDBInstallation sDbInst;
       =20
        if(args.length < 2){
            System.out.println("I need a blast result file name, and a =
fasta filename");
            System.exit(0);
        }
       String blastFileName =3D args[0];
       String fastaFileName =3D args[1];


        //
        // set up querydb and subjectdbinstallation  Is this needed?
       =20
        // query database is needed, I must contain all the query =
sequences that were used in the
        // when blast was called.
       =20
        qDb =3D new HashSequenceDB();
        br =3D new BufferedReader(new	InputStreamReader(new =
FileInputStream(fastaFileName)));
        si =3D SeqIOTools.readFastaDNA(br);
        while (si.hasNext()){
            qDb.addSequence(si.nextSequence());
        }
       =20
        // a subject database installation must exist. It does not have =
to be populated, but if
        // it is you can get the hit sequence from the SearchResult. The =
SequenceDBInstallation
        // must contain a database with the same name as will be found =
in the blast result file
        // in the 'database' field
        sDbInst =3D new SimpleSequenceDBInstallation();
        sDbInst.addSequenceDB(new HashSequenceDB("genome"),new =
HashSet());
       =20
       =20
        // now set up the blast parser things
        parser =3D (XMLReader) new BlastLikeSAXParser();
        ((BlastLikeSAXParser)parser).setModeLazy(); // if you have a =
newer version of blast than 2.01 you need this=20

        // the BlastLikeSearchBuilder will put SeqSimilaritySearch =
results in here
        resultList =3D new Vector();
        sBuilder =3D new BlastLikeSearchBuilder(resultList);
       =20
        // the next 2 lines allows the builder to build SearhResults =
that have access
        // to the query sequences and the database sequences. THIS NOT =
OPTIONAL
        sBuilder.setQuerySeqHolder(qDb);=20
        sBuilder.setSubjectDBInstallation(sDbInst);
       =20
        // adapter converts the datafrom the parser to something the =
builder can handle contentHandler->SearchContentHandler
        cHandler =3D new SeqSimilarityAdapter();
        parser.setContentHandler(cHandler);         // parser will pump =
its data into cHandler
        cHandler.setSearchContentHandler(sBuilder); // handler will pump =
its data into sBuilder
               =20
       =20
        // now read the blast file and parse it
        fStream =3D new FileInputStream(blastFileName);
        parser.parse(new InputSource(fStream));
       =20
        // the results are now int the resultList
        int rCount =3D resultList.size();       =20
        System.out.println("You have " + rCount + " results");
       =20
       =20
        // now lets look at the results
        for (int i =3D 0; i < rCount;i++){
            SeqSimilaritySearchResult r =3D =
(SeqSimilaritySearchResult)resultList.get(i);
            lookAtResult(r);
        }
      =20
       =20
    }
    public static void lookAtResult(SeqSimilaritySearchResult r){
        List hits =3D r.getHits();
        Iterator hIter =3D hits.iterator();
       =20
        while(hIter.hasNext()){
            lookAtHit((SeqSimilaritySearchHit)hIter.next());           =20
        }
    }
   =20
    public static void lookAtHit(SeqSimilaritySearchHit hit){=20
        List sHits =3D hit.getSubHits();
        System.out.println("Subject Id " + hit.getSequenceID());
        System.out.println("start " + hit.getSubjectStart() + " end " + =
hit.getSubjectEnd() + " subhits " + sHits.size());

        Iterator i =3D sHits.iterator();
        while(i.hasNext()){
            lookAtSubHit((SeqSimilaritySearchSubHit)i.next());
            System.out.println("");
        }
    }
    public static void lookAtSubHit(SeqSimilaritySearchSubHit sHit){
        System.out.println("Score " + sHit.getScore());
        Alignment a =3D sHit.getAlignment();
        List l =3D a.getLabels();
        Iterator aIter =3D l.iterator();
        SymbolList sl;
       =20
        while (aIter.hasNext()){
            sl =3D a.symbolListForLabel(aIter.next());
            System.out.println("seq " + sl.seqString());
        }
    }
           =20
=20

}=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
BlastTest.class=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00000664 =00023547 =00000017 =
=0000000010232 07370602152 014655 =
=000=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00ustar=0000dwaring=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00users=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00000000 =00000000 =
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=CA=FE=BA=BE=00=03=00=
-=00=F2=0A=
=00N=00_	=00`=00a=08=00b=0A=
=00c=00d=0A=
=00`=00e=07=00f=0A=
=00=06=00_=07=00g=07=00h=07=00i=0A=
=00=0A=
=00j=0A=
=00	=00k=0A=
=00=08=00l=0A=
=00m=00n=0B=00o=00p=0B=00q=00r=0B=00o=00s=07=00t=0A=
=00=12=00_=08=00u=0A=
=00=06=00j=07=00v=0A=
=00=16=00_=0A=
=00=12=00w=07=00x=0A=
=00=19=00_=0A=
=00=19=00y=07=00z=0A=
=00=1C=00_=07=00{=0A=
=00=1E=00|=0A=
=00=1E=00}=0A=
=00=1E=00~=07=00=7F=0A=
=00"=00_=0B=00=80=00=81=0A=
=00"=00=82=07=00=83=0A=
=00&=00k=0B=00=80=00=84=0B=00=85=00=86=07=00=87=0A=
=00*=00_=08=00=88=0A=
=00*=00=89=0A=
=00*=00=8A=08=00=8B=0A=
=00*=00=8C=0B=00=85=00=8D=07=00=8E=0A=
=00M=00=8F=0B=002=00=90=0B=00=85=00=91=0B=00=92=00=93=07=00=94=0A=
=00M=00=95=0B=00=92=00s=0B=007=00=96=08=00=97=0B=007=00=98=08=00=99=0B=00=
7=00=9A=08=00=9B=0B=007=00=9C=08=00=9D=07=00=9E=0A=
=00M=00=9F=08=00=A0=08=00=A1=0B=00B=00=A2=0A=
=00*=00=A3=0B=00B=00=A4=0B=00=A5=00=A6=0B=00=A5=00=A7=08=00=A8=0B=00=A9=00=
=AA=07=00=AB=07=00=AC=01=00=06<init>=01=00=03()V=01=00=04Code=01=00=0FLin=
eNumberTable=01=00=04main=01=00=16([Ljava/lang/String;)V=01=00=0A=
Exceptions=07=00=AD=01=00=0ClookAtResult=01=005(Lorg/biojava/bio/search/S=
eqSimilaritySearchResult;)V=01=00	=
lookAtHit=01=002(Lorg/biojava/bio/search/SeqSimilaritySearchHit;)V=01=00=0C=
lookAtSubHit=01=005(Lorg/biojava/bio/search/SeqSimilaritySearchSubHit;)V=01=
=00=0A=
SourceFile=01=00=0EBlastTest.java=0C=00O=00P=07=00=AE=0C=00=AF=00=B0=01=00=
5I need a blast result file name, and a fasta =
filename=07=00=B1=0C=00=B2=00=B3=0C=00=B4=00=B5=01=00%org/biojava/bio/seq=
/db/HashSequenceDB=01=00=16java/io/BufferedReader=01=00=19java/io/InputSt=
reamReader=01=00=17java/io/FileInputStream=0C=00O=00=B3=0C=00O=00=B6=0C=00=
O=00=B7=07=00=B8=0C=00=B9=00=BA=07=00=BB=0C=00=BC=00=BD=07=00=BE=0C=00=BF=
=00=C0=0C=00=C1=00=C2=01=003org/biojava/bio/seq/db/SimpleSequenceDBInstal=
lation=01=00=06genome=01=00=11java/util/HashSet=0C=00=C3=00=C4=01=00.org/=
biojava/bio/program/sax/BlastLikeSAXParser=0C=00=C5=00P=01=00=10java/util=
/Vector=01=005org/biojava/bio/program/ssbind/BlastLikeSearchBuilder=0C=00=
O=00=C6=0C=00=C7=00=C8=0C=00=C9=00=CA=01=003org/biojava/bio/program/ssbin=
d/SeqSimilarityAdapter=07=00=CB=0C=00=CC=00=CD=0C=00=CE=00=CF=01=00=17org=
/xml/sax/InputSource=0C=00=D0=00=D1=07=00=D2=0C=00=D3=00=D4=01=00=16java/=
lang/StringBuffer=01=00	You have =0C=00=D5=00=D6=0C=00=D5=00=D7=01=00=08 =
results=0C=00=D8=00=D9=0C=00=DA=00=DB=01=000org/biojava/bio/search/SeqSim=
ilaritySearchResult=0C=00W=00X=0C=00=DC=00=DD=0C=00=DE=00=DF=07=00=E0=0C=00=
=E1=00=E2=01=00-org/biojava/bio/search/SeqSimilaritySearchHit=0C=00Y=00Z=0C=
=00=E3=00=DD=01=00=0BSubject Id =0C=00=E4=00=D9=01=00=06start =
=0C=00=E5=00=D4=01=00=05 end =0C=00=E6=00=D4=01=00	 subhits =
=01=000org/biojava/bio/search/SeqSimilaritySearchSubHit=0C=00[=00\=01=00=00=
=01=00=06Score =
=0C=00=E7=00=E8=0C=00=D5=00=E9=0C=00=EA=00=EB=07=00=EC=0C=00=ED=00=DD=0C=00=
=EE=00=EF=01=00=04seq =07=00=F0=0C=00=F1=00=D9=01=00	=
BlastTest=01=00=10java/lang/Object=01=00=13java/lang/Exception=01=00=10ja=
va/lang/System=01=00=03out=01=00=15Ljava/io/PrintStream;=01=00=13java/io/=
PrintStream=01=00=07println=01=00=15(Ljava/lang/String;)V=01=00=04exit=01=
=00=04(I)V=01=00=18(Ljava/io/InputStream;)V=01=00=13(Ljava/io/Reader;)V=01=
=00!org/biojava/bio/seq/io/SeqIOTools=01=00=0CreadFastaDNA=01=00@(Ljava/i=
o/BufferedReader;)Lorg/biojava/bio/seq/SequenceIterator;=01=00$org/biojav=
a/bio/seq/SequenceIterator=01=00=0CnextSequence=01=00 =
()Lorg/biojava/bio/seq/Sequence;=01=00%org/biojava/bio/seq/db/SequenceDBL=
ite=01=00=0BaddSequence=01=00!(Lorg/biojava/bio/seq/Sequence;)V=01=00=07h=
asNext=01=00=03()Z=01=00=0DaddSequenceDB=01=005(Lorg/biojava/bio/seq/db/S=
equenceDB;Ljava/util/Set;)V=01=00=0BsetModeLazy=01=00=13(Ljava/util/List;=
)V=01=00=11setQuerySeqHolder=01=00&(Lorg/biojava/bio/seq/db/SequenceDB;)V=
=01=00=18setSubjectDBInstallation=01=002(Lorg/biojava/bio/seq/db/Sequence=
DBInstallation;)V=01=00=15org/xml/sax/XMLReader=01=00=11setContentHandler=
=01=00=1F(Lorg/xml/sax/ContentHandler;)V=01=00=17setSearchContentHandler=01=
=008(Lorg/biojava/bio/program/search/SearchContentHandler;)V=01=00=05pars=
e=01=00=1C(Lorg/xml/sax/InputSource;)V=01=00=0Ejava/util/List=01=00=04siz=
e=01=00=03()I=01=00=06append=01=00,(Ljava/lang/String;)Ljava/lang/StringB=
uffer;=01=00=1B(I)Ljava/lang/StringBuffer;=01=00=08toString=01=00=14()Lja=
va/lang/String;=01=00=03get=01=00=15(I)Ljava/lang/Object;=01=00=07getHits=
=01=00=12()Ljava/util/List;=01=00=08iterator=01=00=16()Ljava/util/Iterato=
r;=01=00=12java/util/Iterator=01=00=04next=01=00=14()Ljava/lang/Object;=01=
=00=0A=
getSubHits=01=00=0DgetSequenceID=01=00=0FgetSubjectStart=01=00=0DgetSubje=
ctEnd=01=00=08getScore=01=00=03()D=01=00=1B(D)Ljava/lang/StringBuffer;=01=
=00=0CgetAlignment=01=00$()Lorg/biojava/bio/symbol/Alignment;=01=00 =
org/biojava/bio/symbol/Alignment=01=00	=
getLabels=01=00=12symbolListForLabel=01=007(Ljava/lang/Object;)Lorg/bioja=
va/bio/symbol/SymbolList;=01=00!org/biojava/bio/symbol/SymbolList=01=00	=
seqString=00!=00M=00N=00=00=00=00=00=05=00=01=00O=00P=00=01=00Q=00=00=00=1D=
=00=01=00=01=00=00=00=05*=B7=00=01=B1=00=00=00=01=00R=00=00=00=06=00=01=00=
=00=00=14=00	=
=00S=00T=00=02=00Q=00=00=01=B6=00=07=00=0F=00=00=01&*=BE=05=A2=00=0F=B2=00=
=02=12=03=B6=00=04=03=B8=00=05*=032:=0A=
*=042:=0B=BB=00=06Y=B7=00=07:=08=BB=00=08Y=BB=00	Y=BB=00=0A=
Y=19=0B=B7=00=0B=B7=00=0C=B7=00=0D:=06=19=06=B8=00=0E:=07=A7=00=11=19=08=19=
=07=B9=00=0F=01=00=B9=00=10=02=00=19=07=B9=00=11=01=00=9A=FF=EB=BB=00=12Y=
=B7=00=13:	=19	=
=BB=00=06Y=12=14=B7=00=15=BB=00=16Y=B7=00=17=B6=00=18=BB=00=19Y=B7=00=1AM=
,=C0=00=19=B6=00=1B=BB=00=1CY=B7=00=1DL=BB=00=1EY+=B7=00=1F:=04=19=04=19=08=
=B6=00 =19=04=19	=B6=00!=BB=00"Y=B7=00#N,-=B9=00$=02=00-=19=04=B6=00%=BB=00=0A=
Y=19=0A=
=B7=00=0B:=05,=BB=00&Y=19=05=B7=00'=B9=00(=02=00+=B9=00)=01=006=0C=B2=00=02=
=BB=00*Y=B7=00+=12,=B6=00-=15=0C=B6=00.=12/=B6=00-=B6=000=B6=00=04=036=0D=
=A7=00=18+=15=0D=B9=001=02=00=C0=002:=0E=19=0E=B8=003=84=0D=01=15=0D=15=0C=
=A1=FF=E7=B1=00=00=00=01=00R=00=00=00~=00=1F=00=00=00 =
=00=06=00!=00=0E=00"=00=12=00$=00=17=00%=00=1C=00.=00%=00/=00>=000=00E=00=
1=00H=002=00V=001=00`=009=00i=00:=00~=00>=00=86=00?=00=8D=00B=00=95=00C=00=
=9F=00G=00=A6=00H=00=AD=00K=00=B5=00L=00=BC=00M=00=C2=00Q=00=CD=00R=00=DC=
=00U=00=E4=00V=01=03=00Z=01	=
=00[=01=16=00\=01=1B=00Z=01%=00`=00U=00=00=00=04=00=01=00V=00	=
=00W=00X=00=01=00Q=00=00=00S=00=01=00=03=00=00=00'*=B9=004=01=00L+=B9=005=
=01=00M=A7=00=0F,=B9=006=01=00=C0=007=B8=008,=B9=009=01=00=9A=FF=EE=B1=00=
=00=00=01=00R=00=00=00=1A=00=06=00=00=00b=00=07=00c=00=0E=00e=00=11=00f=00=
=1D=00e=00&=00h=00	=
=00Y=00Z=00=01=00Q=00=00=00=BF=00=03=00=03=00=00=00=87*=B9=00:=01=00L=B2=00=
=02=BB=00*Y=B7=00+=12;=B6=00-*=B9=00<=01=00=B6=00-=B6=000=B6=00=04=B2=00=02=
=BB=00*Y=B7=00+=12=3D=B6=00-*=B9=00>=01=00=B6=00.=12?=B6=00-*=B9=00@=01=00=
=B6=00.=12A=B6=00-+=B9=00)=01=00=B6=00.=B6=000=B6=00=04+=B9=005=01=00M=A7=
=00=17,=B9=006=01=00=C0=00B=B8=00C=B2=00=02=12D=B6=00=04,=B9=009=01=00=9A=
=FF=E6=B1=00=00=00=01=00R=00=00=00&=00	=
=00=00=00k=00=07=00l=00%=00m=00_=00o=00f=00p=00i=00q=00u=00r=00}=00p=00=86=
=00t=00	=
=00[=00\=00=01=00Q=00=00=00=A5=00=04=00=05=00=00=00m=B2=00=02=BB=00*Y=B7=00=
+=12E=B6=00-*=B9=00F=01=00=B6=00G=B6=000=B6=00=04*=B9=00H=01=00L+=B9=00I=01=
=00M,=B9=005=01=00N=A7=000+-=B9=006=01=00=B9=00J=02=00:=04=B2=00=02=BB=00=
*Y=B7=00+=12K=B6=00-=19=04=B9=00L=01=00=B6=00-=B6=000=B6=00=04-=B9=009=01=
=00=9A=FF=CD=B1=00=00=00=01=00R=00=00=00&=00	=
=00=00=00v=00=1E=00w=00%=00x=00,=00y=003=00|=006=00}=00D=00~=00c=00|=00l=00=
=80=00=01=00]=00=00=00=02=00^=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00blast.out=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00000664 =00023547 =00000017 =
=0000000114221 07370577366 013604 =
=000=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00ustar=0000dwaring=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00users=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00000000 =00000000 =
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00BLASTN 2.2.1 =
[Aug-1-2001]=0A=
=0A=
=0A=
Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. =
Schaffer, =0A=
Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997), =0A=
"Gapped BLAST and PSI-BLAST: a new generation of protein database search=0A=
programs",  Nucleic Acids Res. 25:3389-3402.=0A=
=0A=
Query=3D UnnamedSeq1=0A=
         (3650 letters)=0A=
=0A=
Database: genome=0A=
           3518 sequences; 2,805,113,105 total letters=0A=
=0A=
Searching.........done=0A=
=0A=
                                                                   Score =
    E=0A=
Sequences producing significant alignments:                        =
(bits)  Value=0A=
=0A=
ref|NT_007819.4|Hs7_7976 Homo sapiens chromosome 7 working draft...  =
6593  0.0=0A=
ref|NT_025821.3|Hs10_25977 Homo sapiens chromosome 10 working dr...    =
46  0.14=0A=
ref|NT_011493.4|Hs20_11650 Homo sapiens chromosome 20 working dr...    =
44  0.54=0A=
ref|NT_010795.4|Hs17_10952 Homo sapiens chromosome 17 working dr...    =
44  0.54=0A=
ref|NT_010651.4|Hs17_10808 Homo sapiens chromosome 17 working dr...    =
44  0.54=0A=
ref|NT_010604.4|Hs16_10761 Homo sapiens chromosome 16 working dr...    =
44  0.54=0A=
ref|NT_026437.2|Hs14_26604 Homo sapiens chromosome 14 working dr...    =
44  0.54=0A=
ref|NT_024517.4|Hs13_24673 Homo sapiens chromosome 13 working dr...    =
44  0.54=0A=
ref|NT_009459.4|Hs12_9616 Homo sapiens chromosome 12 working dra...    =
44  0.54=0A=
ref|NT_022840.4|Hs4_22996 Homo sapiens chromosome 4 working draf...    =
44  0.54=0A=
ref|NT_005392.4|Hs2_5549 Homo sapiens chromosome 2 working draft...    =
44  0.54=0A=
ref|NT_004641.4|Hs1_4798 Homo sapiens chromosome 1 working draft...    =
44  0.54=0A=
ref|NT_004470.4|Hs1_4627 Homo sapiens chromosome 1 working draft...    =
44  0.54=0A=
=0A=
>ref|NT_007819.4|Hs7_7976 Homo sapiens chromosome 7 working draft =
sequence segment=0A=
          Length =3D 7744418=0A=
=0A=
 Score =3D 6593 bits (3326), Expect =3D 0.0=0A=
 Identities =3D 3340/3347 (99%)=0A=
 Strand =3D Plus / Minus=0A=
=0A=
                                                                         =
  =0A=
Query: 1       =
tttgagggtcttgcaatctacttggaagatacaaggacatgggatttgtaagatgggcca 60=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1096201 =
tttgagggtcttgcaatctacttggaagatacaaggacatgggatttgtaagatgggcca 1096142=0A=
=0A=
                                                                         =
  =0A=
Query: 61      =
actttcacttccttttatgtacctctatattgtttgcatttttcacaaggaagtatcact 120=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1096141 =
actttcacttccttttatgtacctctatattgtttgcatttttcacaaggaagtatcact 1096082=0A=
=0A=
                                                                         =
  =0A=
Query: 121     =
ttgtgcagtgaaggtgcacatgagatgcaatgggcagacggggcttcatggagagggtga 180=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1096081 =
ttgtgcagtgaaggtgcacatgagatgcaatgggcagacggggcttcatggagagggtga 1096022=0A=
=0A=
                                                                         =
  =0A=
Query: 181     =
gccttgagttcaaccttgaagaagaaactgatgtcaaactctctcctaaggcatcttttt 240=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1096021 =
gccttgagttcaaccttgaagaagaaactgatgtcaaactctctcctaaggcatcttttt 1095962=0A=
=0A=
                                                                         =
  =0A=
Query: 241     =
aacattccacaccagaagtgagctgtggaatgcaagaaagactaattacaatgggaaaat 300=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095961 =
aacattccacaccagaagtgagctgtggaatgcaagaaagactaattacaatgggaaaat 1095902=0A=
=0A=
                                                                         =
  =0A=
Query: 301     =
ggatgctgaataaaagcagagaagcccctatgcatcttgggcagtattaagcaacaaatc 360=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095901 =
ggatgctgaataaaagcagagaagcccctatgcatcttgggcagtattaagcaacaaatc 1095842=0A=
=0A=
                                                                         =
  =0A=
Query: 361     =
tcaagatcacatggcagggttttgactcgctcaacattacattactgaatacaaatgcct 420=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095841 =
tcaagatcacatggcagggttttgactcgctcaacattacattactgaatacaaatgcct 1095782=0A=
=0A=
                                                                         =
  =0A=
Query: 421     =
cacatctccttcaccggagtgcaaagataccaagaggtttggcttttgtctgaccgtcac 480=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095781 =
cacatctccttcaccggagtgcaaagataccaagaggtttggcttttgtctgaccgtcac 1095722=0A=
=0A=
                                                                         =
  =0A=
Query: 481     =
tgttggatgaaaaattaagcttgacatttctcaaccctgacttcttactgagaaattgtc 540=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095721 =
tgttggatgaaaaattaagcttgacatttctcaaccctgacttcttactgagaaattgtc 1095662=0A=
=0A=
                                                                         =
  =0A=
Query: 541     =
accattgtctgaacagttactgctgtcttaaagcaaagttaaccaccacatatattgtct 600=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095661 =
accattgtctgaacagttactgctgtcttaaagcaaagttaaccaccacatatattgtct 1095602=0A=
=0A=
                                                                         =
  =0A=
Query: 601     =
gaaacatatttgaatataattcattttacaaaacaaacacattccaggaaagctagtgat 660=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095601 =
gaaacatatttgaatataattcattttacaaaacaaacacattccaggaaagctagtgat 1095542=0A=
=0A=
                                                                         =
  =0A=
Query: 661     =
aaaatgaatctgtctaagtcagctcccactcacctcactttctcacttagcacattgatt 720=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095541 =
aaaatgaatctgtctaagtcagctcccactcacctcactttctcacttagcacattgatt 1095482=0A=
=0A=
                                                                         =
  =0A=
Query: 721     =
tagaaattcactctggtttgggagaacggacttgagaagatggaaaggtgaggagagaga 780=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095481 =
tagaaattcactctggtttgggagaacggacttgagaagatggaaaggtgaggagagaga 1095422=0A=
=0A=
                                                                         =
  =0A=
Query: 781     =
actatttcttactatcctcttctggacaaataataaattggccatttcaggaacatatta 840=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095421 =
actatttcttactatcctcttctggacaaataataaattggccatttcaggaacatatta 1095362=0A=
=0A=
                                                                         =
  =0A=
Query: 841     =
tcaaacagaaaatgttaattttctgattatgcacatgacaagttagctgcaatgaaaaaa 900=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095361 =
tcaaacagaaaatgttaattttctgattatgcacatgacaagttagctgcaatgaaaaaa 1095302=0A=
=0A=
                                                                         =
  =0A=
Query: 901     =
ttgactgaaacctaattaagaaattaggccataagttgttttacttgaggttaaaaaagc 960=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095301 =
ttgactgaaacctaattaagaaattaggccataagttgttttacttgaggttaaaaaagc 1095242=0A=
=0A=
                                                                         =
  =0A=
Query: 961     =
atgtgaaaaggaaatattttcttcctattacatctagctttttaaatggagttcatctca 1020=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095241 =
atgtgaaaaggaaatattttcttcctattacatctagctttttaaatggagttcatctca 1095182=0A=
=0A=
                                                                         =
  =0A=
Query: 1021    =
tcacatctcattctgtgtatgaatcataaatcttgcctactcctgctagaggaaggatct 1080=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095181 =
tcacatctcattctgtgtatgaatcataaatcttgcctactcctgctagaggaaggatct 1095122=0A=
=0A=
                                                                         =
  =0A=
Query: 1081    =
ccctctggtacagctctcctaccaggctaggcccagacatggctccctcctagggaaact 1140=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095121 =
ccctctggtacagctctcctaccaggctaggcccagacatggctccctcctagggaaact 1095062=0A=
=0A=
                                                                         =
  =0A=
Query: 1141    =
gcctccccgacacccctccctgacacactgcccctcccacacctgtctcaattttccctc 1200=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095061 =
gcctccccgacacccctccctgacacactgcccctcccacacctgtctcaattttccctc 1095002=0A=
=0A=
                                                                         =
  =0A=
Query: 1201    =
tccttcagacccaatccagatgtcactttctccataaaaccttatgcaataaaaccttat 1260=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1095001 =
tccttcagacccaatccagatgtcactttctccataaaaccttatgcaataaaaccttat 1094942=0A=
=0A=
                                                                         =
  =0A=
Query: 1261    =
gtatggttttcctaaacaagcctttgaagtgtgggctttggcccagggagacttgcttct 1320=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094941 =
gtatggttttcctaaacaagcctttgaagtgtgggctttggcccagggagacttgcttct 1094882=0A=
=0A=
                                                                         =
  =0A=
Query: 1321    =
gagttttttcattcacatggccacatcgtccctgggggcacaccgggtacaggcttcatt 1380=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094881 =
gagttttttcattcacatggccacatcgtccctgggggcacaccgggtacaggcttcatt 1094822=0A=
=0A=
                                                                         =
  =0A=
Query: 1381    =
gcagagactttccaagtgagttttctcccaagtgagataaccttttatctgtatctttct 1440=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094821 =
gcagagactttccaagtgagttttctcccaagtgagataaccttttatctgtatctttct 1094762=0A=
=0A=
                                                                         =
  =0A=
Query: 1441    =
taaaacaattcccattttatgcccccaaattatggtaaagtataggtgtcttaccttccc 1500=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094761 =
taaaacaattcccattttatgcccccaaattatggtaaagtataggtgtcttaccttccc 1094702=0A=
=0A=
                                                                         =
  =0A=
Query: 1501    =
tgatagagtcaaccccttgaggttgggaaccacttcctcacttgtgaatctctgacattg 1560=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094701 =
tgatagagtcaaccccttgaggttgggaaccacttcctcacttgtgaatctctgacattg 1094642=0A=
=0A=
                                                                         =
  =0A=
Query: 1561    =
cagcagagttgcacatacagtaggtgctcagtaatgaagcagttattgaacaagtgcatg 1620=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094641 =
cagcagagttgcacatacagtaggtgctcagtaatgaagcagttattgaacaagtgcatg 1094582=0A=
=0A=
                                                                         =
  =0A=
Query: 1621    =
gacagtacgagtcatgactctggcctctggtgaggctggcagtgcccacagaaggttcag 1680=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094581 =
gacagtacgagtcatgactctggcctctggtgaggctggcagtgcccacagaaggttcag 1094522=0A=
=0A=
                                                                         =
  =0A=
Query: 1681    =
tcaccttctccatcccttacatgaataccaacttgcctggctttagcatgagtgtttcca 1740=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094521 =
tcaccttctccatcccttacatgaataccaacttgcctggctttagcatgagtgtttcca 1094462=0A=
=0A=
                                                                         =
  =0A=
Query: 1741    =
taacttattacaactatatgttcaaatgcatcttaagtgagtcatccctgtcctccatat 1800=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094461 =
taacttattacaactatatgttcaaatgcatcttaagtgagtcatccctgtcctccatat 1094402=0A=
=0A=
                                                                         =
  =0A=
Query: 1801    =
tccctattctatgtcatttgtggatttaatcagaagcaagggtttatgtgatgcactgag 1860=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094401 =
tccctattctatgtcatttgtggatttaatcagaagcaagggtttatgtgatgcactgag 1094342=0A=
=0A=
                                                                         =
  =0A=
Query: 1861    =
aacctctaacttccttcttcatcaggaggaggctcagaacccagaggctactcagcattc 1920=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094341 =
aacctctaacttccttcttcatcaggaggaggctcagaacccagaggctactcagcattc 1094282=0A=
=0A=
                                                                         =
  =0A=
Query: 1921    =
agggcacgtattcccacagagaatacaagctggggggaaaactcacttggagactctctg 1980=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094281 =
agggcacgtattcccacagagaatacaagctggggggaaaactcacttggagactctctg 1094222=0A=
=0A=
                                                                         =
  =0A=
Query: 1981    =
cattgaagcatgtgcctggtgtgcccccagggacaatgcggccacgtggatgaaaacact 2040=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094221 =
cattgaagcatgtgcctggtgtgcccccagggacaatgcggccacgtggatgaaaacact 1094162=0A=
=0A=
                                                                         =
  =0A=
Query: 2041    =
cagaagcaaatctccctgggccaaagcccacacttcaaaggcttgtttaggaaaaccata 2100=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094161 =
cagaagcaaatctccctgggccaaagcccacacttcaaaggcttgtttaggaaaaccata 1094102=0A=
=0A=
                                                                         =
  =0A=
Query: 2101    =
cactcaaccttctttgtgtattgtgaaaataaacatttacagaacattggaatgtgtgaa 2160=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094101 =
cactcaaccttctttgtgtattgtgaaaataaacatttacagaacattggaatgtgtgaa 1094042=0A=
=0A=
                                                                         =
  =0A=
Query: 2161    =
aaatgtccttcgttgtgtgtttatggaaaacatttggcaagtatttgctataaatgcttc 2220=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1094041 =
aaatgtccttcgttgtgtgtttatggaaaacatttggcaagtatttgctataaatgcttc 1093982=0A=
=0A=
                                                                         =
  =0A=
Query: 2221    =
ctgtaaatgctgtgagtgagaatcagcaggaaaatgtatgagatgagagggagccaactc 2280=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093981 =
ctgtaaatgctgtgagtgagaatcagcaggaaaatgtatgagatgagagggagccaactc 1093922=0A=
=0A=
                                                                         =
  =0A=
Query: 2281    =
ccacagcaaaagaatgtccccttcccttccacacaccccacgcgctgcctccactctgca 2340=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093921 =
ccacagcaaaagaatgtccccttcccttccacacaccccacgcgctgcctccactctgca 1093862=0A=
=0A=
                                                                         =
  =0A=
Query: 2341    =
ctcctccgccgcccccctccccagataaggcccaggggatccgggagtgtgagggtgggg 2400=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093861 =
ctcctccgccgcccccctccccagataaggcccaggggatccgggagtgtgagggtgggg 1093802=0A=
=0A=
                                                                         =
  =0A=
Query: 2401    =
aacagagaagggacaatctgtgaaggggcctccaggctcagcattccatgttggagtgca 2460=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093801 =
aacagagaagggacaatctgtgaaggggcctccaggctcagcattccatgttggagtgca 1093742=0A=
=0A=
                                                                         =
  =0A=
Query: 2461    =
gaaaaggaaactgactcgtgtttcctaggtgatgacgccccctaataggatttccttccc 2520=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093741 =
gaaaaggaaactgactcgtgtttcctaggtgatgacgccccctaataggatttccttccc 1093682=0A=
=0A=
                                                                         =
  =0A=
Query: 2521    =
cttcattcttattggttaacagttctgattccaagtcccccattccagctggaaacattg 2580=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093681 =
cttcattcttattggttaacagttctgattccaagtcccccattccagctggaaacattg 1093622=0A=
=0A=
                                                                         =
  =0A=
Query: 2581    =
cctcaacacccacccggtcttccccacctgaagactggcctgcctctcagttaactcttg 2640=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093621 =
cctcaacacccacccggtcttccccacctgaagactggcctgcctctcagttaactcttg 1093562=0A=
=0A=
                                                                         =
  =0A=
Query: 2641    =
gacgctttcaggtgaaaagtaaaaagtggaaagcaagagactcaaatatatttcttcaga 2700=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093561 =
gacgctttcaggtgaaaagtaaaaagtggaaagcaagagactcaaatatatttcttcaga 1093502=0A=
=0A=
                                                                         =
  =0A=
Query: 2701    =
acaagatgcatgggacatttctcttaattcctcaacaatggacaaagacctctgatgata 2760=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093501 =
acaagatgcatgggacatttctcttaattcctcaacaatggacaaagacctctgatgata 1093442=0A=
=0A=
                                                                         =
  =0A=
Query: 2761    =
tgaactgaaaagggaagagagtgttgcaggccacagcaatgaacgctcgtggtttttacc 2820=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093441 =
tgaactgaaaagggaagagagtgttgcaggccacagcaatgaacgctcgtggtttttacc 1093382=0A=
=0A=
                                                                         =
  =0A=
Query: 2821    =
aggaagcaaccttgaccgcctccagattgccctctgggatgaaccaggtcctcccctgcc 2880=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093381 =
aggaagcaaccttgaccgcctccagattgccctctgggatgaaccaggtcctcccctgcc 1093322=0A=
=0A=
                                                                         =
  =0A=
Query: 2881    =
acatgaggtatgtagtgggtactgcagagaggctggctgctcttggttcccaagcagcct 2940=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093321 =
acatgaggtatgtagtgggtactgcagagaggctggctgctcttggttcccaagcagcct 1093262=0A=
=0A=
                                                                         =
  =0A=
Query: 2941    =
ccgtgataaacagcctcagtgttgaaatatcagggagattaattatctcgcagatgttac 3000=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093261 =
ccgtgataaacagcctcagtgttgaaatatcagggagattaattatctcgcagatgttac 1093202=0A=
=0A=
                                                                         =
  =0A=
Query: 3001    =
atgccaaaaaggctacnnnnnnncttcatttcagaattctttccttttttaagaaataca 3060=0A=
               ||||||||||||||||       =
|||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093201 =
atgccaaaaaggctactttttttcttcatttcagaattctttccttttttaagaaataca 1093142=0A=
=0A=
                                                                         =
  =0A=
Query: 3061    =
aggacacaaatccaatgtgcttaattaagcctgccaaggcgaagaggcccagctgcctgg 3120=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093141 =
aggacacaaatccaatgtgcttaattaagcctgccaaggcgaagaggcccagctgcctgg 1093082=0A=
=0A=
                                                                         =
  =0A=
Query: 3121    =
ctgtatttgtatctgcacaggtgcccagggaacaccagcctcactctcagcttcatctcc 3180=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093081 =
ctgtatttgtatctgcacaggtgcccagggaacaccagcctcactctcagcttcatctcc 1093022=0A=
=0A=
                                                                         =
  =0A=
Query: 3181    =
aggtcccctacctctgcccagaggtcctcacatctgcacctgactggagctaatgacaaa 3240=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1093021 =
aggtcccctacctctgcccagaggtcctcacatctgcacctgactggagctaatgacaaa 1092962=0A=
=0A=
                                                                         =
  =0A=
Query: 3241    =
catccttacacagcaacttttacgctgatttgcaaaatgatttattttaagtacgtttct 3300=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1092961 =
catccttacacagcaacttttacgctgatttgcaaaatgatttattttaagtacgtttct 1092902=0A=
=0A=
                                                              =0A=
Query: 3301    aaaaacataaacgttaagtcaaaaggtctcagcatttctgccactta 3347=0A=
               |||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1092901 aaaaacataaacgttaagtcaaaaggtctcagcatttctgccactta 1092855=0A=
=0A=
=0A=
 Score =3D  444 bits (224), Expect =3D e-121=0A=
 Identities =3D 224/224 (100%)=0A=
 Strand =3D Plus / Minus=0A=
=0A=
                                                                         =
  =0A=
Query: 3427    =
agtattgataaaccataatctaaaaatgatgtcatttcacacttcccaaagcaatggagg 3486=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1092775 =
agtattgataaaccataatctaaaaatgatgtcatttcacacttcccaaagcaatggagg 1092716=0A=
=0A=
                                                                         =
  =0A=
Query: 3487    =
agttccctttaaaatcatttcaccggctcacctagatttatcatgggaaagaaaaatcca 3546=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1092715 =
agttccctttaaaatcatttcaccggctcacctagatttatcatgggaaagaaaaatcca 1092656=0A=
=0A=
                                                                         =
  =0A=
Query: 3547    =
cagtctcctctgatggaagaagaaaggggctatggggccccccagactttcccaatcagc 3606=0A=
               =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1092655 =
cagtctcctctgatggaagaagaaaggggctatggggccccccagactttcccaatcagc 1092596=0A=
=0A=
                                                           =0A=
Query: 3607    cccatgctgtggctctccccagccatggaccatgacagccactc 3650=0A=
               ||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 1092595 cccatgctgtggctctccccagccatggaccatgacagccactc 1092552=0A=
=0A=
=0A=
 Score =3D  194 bits (98), Expect =3D 2e-46=0A=
 Identities =3D 134/146 (91%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                                                         =
  =0A=
Query: 1957    =
gaaaactcacttggagactctctgcattgaagcatgtgcctggtgtgcccccagggacaa 2016=0A=
               ||||||||||||||| | |||||||| |||||| ||| || =
||||||||||||||||| |=0A=
Sbjct: 1094797 =
gaaaactcacttggaaagtctctgcaatgaagcctgtacccggtgtgcccccagggacga 1094856=0A=
=0A=
                                                                         =
  =0A=
Query: 2017    =
tgcggccacgtggatgaaaacactcagaagcaaatctccctgggccaaagcccacacttc 2076=0A=
               || ||||| ||| ||||||| |||||||||||| =
||||||||||||||||||||||||||=0A=
Sbjct: 1094857 =
tgtggccatgtgaatgaaaaaactcagaagcaagtctccctgggccaaagcccacacttc 1094916=0A=
=0A=
                                         =0A=
Query: 2077    aaaggcttgtttaggaaaaccataca 2102=0A=
               ||||||||||||||||||||||||||=0A=
Sbjct: 1094917 aaaggcttgtttaggaaaaccataca 1094942=0A=
=0A=
=0A=
 Score =3D  194 bits (98), Expect =3D 2e-46=0A=
 Identities =3D 134/146 (91%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                                                         =
  =0A=
Query: 1260    =
tgtatggttttcctaaacaagcctttgaagtgtgggctttggcccagggagacttgcttc 1319=0A=
               |||||||||||||||||||||||||||||||||||||||||||||||||||| =
|||||||=0A=
Sbjct: 1094100 =
tgtatggttttcctaaacaagcctttgaagtgtgggctttggcccagggagatttgcttc 1094159=0A=
=0A=
                                                                         =
  =0A=
Query: 1320    =
tgagttttttcattcacatggccacatcgtccctgggggcacaccgggtacaggcttcat 1379=0A=
               ||||| ||||||| ||| ||||| ||| ||||||||||||||||| || ||| =
|||||| =0A=
Sbjct: 1094160 =
tgagtgttttcatccacgtggccgcattgtccctgggggcacaccaggcacatgcttcaa 1094219=0A=
=0A=
                                         =0A=
Query: 1380    tgcagagactttccaagtgagttttc 1405=0A=
               |||||||| | |||||||||||||||=0A=
Sbjct: 1094220 tgcagagagtctccaagtgagttttc 1094245=0A=
=0A=
=0A=
>ref|NT_025821.3|Hs10_25977 Homo sapiens chromosome 10 working draft =
sequence segment=0A=
          Length =3D 1368624=0A=
=0A=
 Score =3D 46.1 bits (23), Expect =3D 0.14=0A=
 Identities =3D 23/23 (100%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                     =0A=
Query: 964    tgaaaaggaaatattttcttcct 986=0A=
              |||||||||||||||||||||||=0A=
Sbjct: 519120 tgaaaaggaaatattttcttcct 519142=0A=
=0A=
=0A=
>ref|NT_011493.4|Hs20_11650 Homo sapiens chromosome 20 working draft =
sequence segment=0A=
          Length =3D 1489287=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.54=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Minus=0A=
=0A=
                                    =0A=
Query: 1573   acatacagtaggtgctcagtaa 1594=0A=
              ||||||||||||||||||||||=0A=
Sbjct: 280588 acatacagtaggtgctcagtaa 280567=0A=
=0A=
=0A=
>ref|NT_010795.4|Hs17_10952 Homo sapiens chromosome 17 working draft =
sequence segment=0A=
          Length =3D 1138774=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.54=0A=
 Identities =3D 25/26 (96%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                        =0A=
Query: 1425   ttatctgtatctttcttaaaacaatt 1450=0A=
              |||||| |||||||||||||||||||=0A=
Sbjct: 270394 ttatctctatctttcttaaaacaatt 270419=0A=
=0A=
=0A=
>ref|NT_010651.4|Hs17_10808 Homo sapiens chromosome 17 working draft =
sequence segment=0A=
          Length =3D 3054331=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.54=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                     =0A=
Query: 1573    acatacagtaggtgctcagtaa 1594=0A=
               ||||||||||||||||||||||=0A=
Sbjct: 2160182 acatacagtaggtgctcagtaa 2160203=0A=
=0A=
=0A=
>ref|NT_010604.4|Hs16_10761 Homo sapiens chromosome 16 working draft =
sequence segment=0A=
          Length =3D 2471742=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.54=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                     =0A=
Query: 847     agaaaatgttaattttctgatt 868=0A=
               ||||||||||||||||||||||=0A=
Sbjct: 1935902 agaaaatgttaattttctgatt 1935923=0A=
=0A=
=0A=
>ref|NT_026437.2|Hs14_26604 Homo sapiens chromosome 14 working draft =
sequence segment=0A=
          Length =3D 12185353=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.54=0A=
 Identities =3D 25/26 (96%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                         =0A=
Query: 276     gaaagactaattacaatgggaaaatg 301=0A=
               ||||||| ||||||||||||||||||=0A=
Sbjct: 3982427 gaaagacgaattacaatgggaaaatg 3982452=0A=
=0A=
=0A=
>ref|NT_024517.4|Hs13_24673 Homo sapiens chromosome 13 working draft =
sequence segment=0A=
          Length =3D 888756=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.54=0A=
 Identities =3D 25/26 (96%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                        =0A=
Query: 1573   acatacagtaggtgctcagtaatgaa 1598=0A=
              |||||||| |||||||||||||||||=0A=
Sbjct: 716216 acatacagaaggtgctcagtaatgaa 716241=0A=
=0A=
=0A=
>ref|NT_009459.4|Hs12_9616 Homo sapiens chromosome 12 working draft =
sequence segment=0A=
          Length =3D 1482321=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.54=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Minus=0A=
=0A=
                                     =0A=
Query: 1574    catacagtaggtgctcagtaat 1595=0A=
               ||||||||||||||||||||||=0A=
Sbjct: 1416921 catacagtaggtgctcagtaat 1416900=0A=
=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.54=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Minus=0A=
=0A=
                                     =0A=
Query: 1574    catacagtaggtgctcagtaat 1595=0A=
               ||||||||||||||||||||||=0A=
Sbjct: 1416957 catacagtaggtgctcagtaat 1416936=0A=
=0A=
=0A=
>ref|NT_022840.4|Hs4_22996 Homo sapiens chromosome 4 working draft =
sequence segment=0A=
          Length =3D 480938=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.54=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Minus=0A=
=0A=
                                    =0A=
Query: 302    gatgctgaataaaagcagagaa 323=0A=
              ||||||||||||||||||||||=0A=
Sbjct: 311390 gatgctgaataaaagcagagaa 311369=0A=
=0A=
=0A=
>ref|NT_005392.4|Hs2_5549 Homo sapiens chromosome 2 working draft =
sequence segment=0A=
          Length =3D 856227=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.54=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                    =0A=
Query: 837    attatcaaacagaaaatgttaa 858=0A=
              ||||||||||||||||||||||=0A=
Sbjct: 139872 attatcaaacagaaaatgttaa 139893=0A=
=0A=
=0A=
>ref|NT_004641.4|Hs1_4798 Homo sapiens chromosome 1 working draft =
sequence segment=0A=
          Length =3D 2011841=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.54=0A=
 Identities =3D 25/26 (96%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                         =0A=
Query: 1436    tttcttaaaacaattcccattttatg 1461=0A=
               |||||||||||||| |||||||||||=0A=
Sbjct: 1943104 tttcttaaaacaatacccattttatg 1943129=0A=
=0A=
=0A=
>ref|NT_004470.4|Hs1_4627 Homo sapiens chromosome 1 working draft =
sequence segment=0A=
          Length =3D 1587415=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.54=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                     =0A=
Query: 1779    gagtcatccctgtcctccatat 1800=0A=
               ||||||||||||||||||||||=0A=
Sbjct: 1221697 gagtcatccctgtcctccatat 1221718=0A=
=0A=
=0A=
  Database: genome=0A=
    Posted date:  Aug 29, 2001 12:46 PM=0A=
  Number of letters in database: 2,805,113,105=0A=
  Number of sequences in database:  3518=0A=
  =0A=
Lambda     K      H=0A=
    1.37    0.711     1.31 =0A=
=0A=
Gapped=0A=
Lambda     K      H=0A=
    1.37    0.711     1.31 =0A=
=0A=
=0A=
Matrix: blastn matrix:1 -3=0A=
Gap Penalties: Existence: 5, Extension: 2=0A=
Number of Hits to DB: 8,146,679=0A=
Number of Sequences: 3518=0A=
Number of extensions: 8146679=0A=
Number of successful extensions: 420=0A=
Number of sequences better than  1.0: 13=0A=
length of query: 3650=0A=
length of database: 2,805,113,105=0A=
effective HSP length: 22=0A=
effective length of query: 3628=0A=
effective length of database: 2,805,035,709=0A=
effective search space: 10176669552252=0A=
effective search space used: 10176669552252=0A=
T: 0=0A=
A: 40=0A=
X1: 6 (11.9 bits)=0A=
X2: 15 (29.7 bits)=0A=
S1: 12 (24.3 bits)=0A=
S2: 22 (44.1 bits)=0A=
BLASTN 2.2.1 [Aug-1-2001]=0A=
=0A=
=0A=
Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. =
Schaffer, =0A=
Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997), =0A=
"Gapped BLAST and PSI-BLAST: a new generation of protein database search=0A=
programs",  Nucleic Acids Res. 25:3389-3402.=0A=
=0A=
Query=3D UnnamedSeq2=0A=
         (1860 letters)=0A=
=0A=
Database: genome=0A=
           3518 sequences; 2,805,113,105 total letters=0A=
=0A=
Searching.........done=0A=
=0A=
                                                                   Score =
    E=0A=
Sequences producing significant alignments:                        =
(bits)  Value=0A=
=0A=
ref|NT_006014.5|Hs3_6171 Homo sapiens chromosome 3 working draft...   =
694  0.0=0A=
ref|NT_010705.4|Hs17_10862 Homo sapiens chromosome 17 working dr...    =
48  0.017=0A=
ref|NT_011651.4|HsX_11808 Homo sapiens chromosome X working draf...    =
44  0.27=0A=
ref|NT_024524.4|Hs13_24680 Homo sapiens chromosome 13 working dr...    =
44  0.27=0A=
ref|NT_022776.4|Hs4_22932 Homo sapiens chromosome 4 working draf...    =
44  0.27=0A=
ref|NT_005187.5|Hs2_5344 Homo sapiens chromosome 2 working draft...    =
44  0.27=0A=
=0A=
>ref|NT_006014.5|Hs3_6171 Homo sapiens chromosome 3 working draft =
sequence segment=0A=
          Length =3D 1084762=0A=
=0A=
 Score =3D  694 bits (350), Expect =3D 0.0=0A=
 Identities =3D 350/350 (100%)=0A=
 Strand =3D Plus / Minus=0A=
=0A=
                                                                         =
 =0A=
Query: 434    =
cattttcctatatctgttttggttcaaattagttcccaaaaactgtctccctcctagtca 493=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 843794 =
cattttcctatatctgttttggttcaaattagttcccaaaaactgtctccctcctagtca 843735=0A=
=0A=
                                                                         =
 =0A=
Query: 494    =
aagatgttaaaaaacagcatataatgacttacatttaaactaaaattaaggtgggaaaat 553=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 843734 =
aagatgttaaaaaacagcatataatgacttacatttaaactaaaattaaggtgggaaaat 843675=0A=
=0A=
                                                                         =
 =0A=
Query: 554    =
ataattttaaataagtttaaacgagtcaaagacaaatgttctaaatttgaaaaaagttag 613=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 843674 =
ataattttaaataagtttaaacgagtcaaagacaaatgttctaaatttgaaaaaagttag 843615=0A=
=0A=
                                                                         =
 =0A=
Query: 614    =
acaataaagctaaattaacctaacatttctgaatgcttacctaaaccacgttaattctat 673=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 843614 =
acaataaagctaaattaacctaacatttctgaatgcttacctaaaccacgttaattctat 843555=0A=
=0A=
                                                                         =
 =0A=
Query: 674    =
acgaaacagccttacttacaaaatagttacccataaatacaagtaataagaggatagatt 733=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 843554 =
acgaaacagccttacttacaaaatagttacccataaatacaagtaataagaggatagatt 843495=0A=
=0A=
                                                                =0A=
Query: 734    acaacctaaatatttttaaatggtaacgtacaaccccattgatctagatg 783=0A=
              ||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 843494 acaacctaaatatttttaaatggtaacgtacaaccccattgatctagatg 843445=0A=
=0A=
=0A=
 Score =3D  668 bits (337), Expect =3D 0.0=0A=
 Identities =3D 337/337 (100%)=0A=
 Strand =3D Plus / Minus=0A=
=0A=
                                                                         =
 =0A=
Query: 48     =
tgttcttcttaggcaatacaatttaacactaaaattaagatgttaccactaaatacagtg 107=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 844180 =
tgttcttcttaggcaatacaatttaacactaaaattaagatgttaccactaaatacagtg 844121=0A=
=0A=
                                                                         =
 =0A=
Query: 108    =
tctctcatgtttaaaatcattcaagtctttcttgcaggaagaatttaagcttcatcagac 167=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 844120 =
tctctcatgtttaaaatcattcaagtctttcttgcaggaagaatttaagcttcatcagac 844061=0A=
=0A=
                                                                         =
 =0A=
Query: 168    =
actcatgggatgaggactgaggtatgaaatctacatgttttctgaccagaacacatccaa 227=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 844060 =
actcatgggatgaggactgaggtatgaaatctacatgttttctgaccagaacacatccaa 844001=0A=
=0A=
                                                                         =
 =0A=
Query: 228    =
catggcctaaagattatcttactcttctaaattaagatcaaatagatttagagtaactta 287=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 844000 =
catggcctaaagattatcttactcttctaaattaagatcaaatagatttagagtaactta 843941=0A=
=0A=
                                                                         =
 =0A=
Query: 288    =
ctgtagcagggttcttagatttgtattttatacagacgaatcaagaaaccattctactaa 347=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 843940 =
ctgtagcagggttcttagatttgtattttatacagacgaatcaagaaaccattctactaa 843881=0A=
=0A=
                                                   =0A=
Query: 348    tataaaatatcagtagttatgaaaattattttcaacc 384=0A=
              |||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 843880 tataaaatatcagtagttatgaaaattattttcaacc 843844=0A=
=0A=
=0A=
 Score =3D  601 bits (303), Expect =3D e-169=0A=
 Identities =3D 317/324 (97%)=0A=
 Strand =3D Plus / Minus=0A=
=0A=
                                                                         =
 =0A=
Query: 1537   =
tgtactagttaatagaagtttcaactggagcaacaatacttatgaatgcatacatttaaa 1596=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 842691 =
tgtactagttaatagaagtttcaactggagcaacaatacttatgaatgcatacatttaaa 842632=0A=
=0A=
                                                                         =
 =0A=
Query: 1597   =
ggtcagatttatttgctatgtaccaagtaaaaccaaaaaagtcacattcatatagcatct 1656=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 842631 =
ggtcagatttatttgctatgtaccaagtaaaaccaaaaaagtcacattcatatagcatct 842572=0A=
=0A=
                                                                         =
 =0A=
Query: 1657   =
cattaattgtaaacacaatacgaaaaatgtttaaatactagttaaaatacagtgtgtatg 1716=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 842571 =
cattaattgtaaacacaatacgaaaaatgtttaaatactagttaaaatacagtgtgtatg 842512=0A=
=0A=
                                                                         =
 =0A=
Query: 1717   =
atactagctatnnnnnnnttaattccctacaaatttaaaatacaatctatttacaactaa 1776=0A=
              |||||||||||       =
||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 842511 =
atactagctataaaaaaattaattccctacaaatttaaaatacaatctatttacaactaa 842452=0A=
=0A=
                                                                         =
 =0A=
Query: 1777   =
atgctaccactttctccactcttcttggtaccaagtaacccctaaaattgagaggtggct 1836=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 842451 =
atgctaccactttctccactcttcttggtaccaagtaacccctaaaattgagaggtggct 842392=0A=
=0A=
                                      =0A=
Query: 1837   tagcaaactgcgtcaacatgccac 1860=0A=
              ||||||||||||||||||||||||=0A=
Sbjct: 842391 tagcaaactgcgtcaacatgccac 842368=0A=
=0A=
=0A=
 Score =3D  428 bits (216), Expect =3D e-117=0A=
 Identities =3D 232/240 (96%)=0A=
 Strand =3D Plus / Minus=0A=
=0A=
                                                                         =
 =0A=
Query: 989    =
atctaatccttagtattttgaaagnnnnnnnntacccaacatcacactgttgttttaatt 1048=0A=
              ||||||||||||||||||||||||        =
||||||||||||||||||||||||||||=0A=
Sbjct: 843239 =
atctaatccttagtattttgaaagaaaaaaaatacccaacatcacactgttgttttaatt 843180=0A=
=0A=
                                                                         =
 =0A=
Query: 1049   =
taaagaccaatctctgtcattcagcatttatgaaatactagttttactaaacattatatg 1108=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 843179 =
taaagaccaatctctgtcattcagcatttatgaaatactagttttactaaacattatatg 843120=0A=
=0A=
                                                                         =
 =0A=
Query: 1109   =
gattttctaaactcagttctccacctaaatttttcattggaactccaagcttcagaacaa 1168=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 843119 =
gattttctaaactcagttctccacctaaatttttcattggaactccaagcttcagaacaa 843060=0A=
=0A=
                                                                         =
 =0A=
Query: 1169   =
tgaatgagaacatctgttgcatttaagtaataaaatcttgcagattgccaatccaaatgg 1228=0A=
              =
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||=0A=
Sbjct: 843059 =
tgaatgagaacatctgttgcatttaagtaataaaatcttgcagattgccaatccaaatgg 843000=0A=
=0A=
=0A=
>ref|NT_010705.4|Hs17_10862 Homo sapiens chromosome 17 working draft =
sequence segment=0A=
          Length =3D 374660=0A=
=0A=
 Score =3D 48.1 bits (24), Expect =3D 0.017=0A=
 Identities =3D 24/24 (100%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                      =0A=
Query: 520    acttacatttaaactaaaattaag 543=0A=
              ||||||||||||||||||||||||=0A=
Sbjct: 130537 acttacatttaaactaaaattaag 130560=0A=
=0A=
=0A=
>ref|NT_011651.4|HsX_11808 Homo sapiens chromosome X working draft =
sequence segment=0A=
          Length =3D 2274920=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.27=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                     =0A=
Query: 63      atacaatttaacactaaaatta 84=0A=
               ||||||||||||||||||||||=0A=
Sbjct: 2173742 atacaatttaacactaaaatta 2173763=0A=
=0A=
=0A=
>ref|NT_024524.4|Hs13_24680 Homo sapiens chromosome 13 working draft =
sequence segment=0A=
          Length =3D 5607782=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.27=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                     =0A=
Query: 552     atataattttaaataagtttaa 573=0A=
               ||||||||||||||||||||||=0A=
Sbjct: 5167684 atataattttaaataagtttaa 5167705=0A=
=0A=
=0A=
>ref|NT_022776.4|Hs4_22932 Homo sapiens chromosome 4 working draft =
sequence segment=0A=
          Length =3D 1039623=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.27=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Minus=0A=
=0A=
                                    =0A=
Query: 549    aaaatataattttaaataagtt 570=0A=
              ||||||||||||||||||||||=0A=
Sbjct: 960496 aaaatataattttaaataagtt 960475=0A=
=0A=
=0A=
>ref|NT_005187.5|Hs2_5344 Homo sapiens chromosome 2 working draft =
sequence segment=0A=
          Length =3D 1557795=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.27=0A=
 Identities =3D 25/26 (96%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                         =0A=
Query: 352     aaatatcagtagttatgaaaattatt 377=0A=
               |||||||||| |||||||||||||||=0A=
Sbjct: 1187166 aaatatcagttgttatgaaaattatt 1187191=0A=
=0A=
=0A=
  Database: genome=0A=
    Posted date:  Aug 29, 2001 12:46 PM=0A=
  Number of letters in database: 2,805,113,105=0A=
  Number of sequences in database:  3518=0A=
  =0A=
Lambda     K      H=0A=
    1.37    0.711     1.31 =0A=
=0A=
Gapped=0A=
Lambda     K      H=0A=
    1.37    0.711     1.31 =0A=
=0A=
=0A=
Matrix: blastn matrix:1 -3=0A=
Gap Penalties: Existence: 5, Extension: 2=0A=
Number of Hits to DB: 3,438,566=0A=
Number of Sequences: 3518=0A=
Number of extensions: 3438566=0A=
Number of successful extensions: 184=0A=
Number of sequences better than  1.0: 6=0A=
length of query: 1860=0A=
length of database: 2,805,113,105=0A=
effective HSP length: 21=0A=
effective length of query: 1839=0A=
effective length of database: 2,805,039,227=0A=
effective search space: 5158467138453=0A=
effective search space used: 5158467138453=0A=
T: 0=0A=
A: 40=0A=
X1: 6 (11.9 bits)=0A=
X2: 15 (29.7 bits)=0A=
S1: 12 (24.3 bits)=0A=
S2: 22 (44.1 bits)=0A=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00blastInput.fasta=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00000764 =00023547 =00000017 =0000000013017 =
07370601500 015067 =
=000=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00ustar=0000dwaring=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00users=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00000000 =00000000 =
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00>UnnamedSeq1=0A=
TTTGAGGGTCTTGCAATCTACTTGGAAGATACAAGGACATGGGATTTGTA=0A=
AGATGGGCCAACTTTCACTTCCTTTTATGTACCTCTATATTGTTTGCATT=0A=
TTTCACAAGGAAGTATCACTTTGTGCAGTGAAGGTGCACATGAGATGCAA=0A=
TGGGCAGACGGGGCTTCATGGAGAGGGTGAGCCTTGAGTTCAACCTTGAA=0A=
GAAGAAACTGATGTCAAACTCTCTCCTAAGGCATCTTTTTAACATTCCAC=0A=
ACCAGAAGTGAGCTGTGGAATGCAAGAAAGACTAATTACAATGGGAAAAT=0A=
GGATGCTGAATAAAAGCAGAGAAGCCCCTATGCATCTTGGGCAGTATTAA=0A=
GCAACAAATCTCAAGATCACATGGCAGGGTTTTGACTCGCTCAACATTAC=0A=
ATTACTGAATACAAATGCCTCACATCTCCTTCACCGGAGTGCAAAGATAC=0A=
CAAGAGGTTTGGCTTTTGTCTGACCGTCACTGTTGGATGAAAAATTAAGC=0A=
TTGACATTTCTCAACCCTGACTTCTTACTGAGAAATTGTCACCATTGTCT=0A=
GAACAGTTACTGCTGTCTTAAAGCAAAGTTAACCACCACATATATTGTCT=0A=
GAAACATATTTGAATATAATTCATTTTACAAAACAAACACATTCCAGGAA=0A=
AGCTAGTGATAAAATGAATCTGTCTAAGTCAGCTCCCACTCACCTCACTT=0A=
TCTCACTTAGCACATTGATTTAGAAATTCACTCTGGTTTGGGAGAACGGA=0A=
CTTGAGAAGATGGAAAGGTGAGGAGAGAGAACTATTTCTTACTATCCTCT=0A=
TCTGGACAAATAATAAATTGGCCATTTCAGGAACATATTATCAAACAGAA=0A=
AATGTTAATTTTCTGATTATGCACATGACAAGTTAGCTGCAATGAAAAAA=0A=
TTGACTGAAACCTAATTAAGAAATTAGGCCATAAGTTGTTTTACTTGAGG=0A=
TTAAAAAAGCATGTGAAAAGGAAATATTTTCTTCCTATTACATCTAGCTT=0A=
TTTAAATGGAGTTCATCTCATCACATCTCATTCTGTGTATGAATCATAAA=0A=
TCTTGCCTACTCCTGCTAGAGGAAGGATCTCCCTCTGGTACAGCTCTCCT=0A=
ACCAGGCTAGGCCCAGACATGGCTCCCTCCTAGGGAAACTGCCTCCCCGA=0A=
CACCCCTCCCTGACACACTGCCCCTCCCACACCTGTCTCAATTTTCCCTC=0A=
TCCTTCAGACCCAATCCAGATGTCACTTTCTCCATAAAACCTTATGCAAT=0A=
AAAACCTTATGTATGGTTTTCCTAAACAAGCCTTTGAAGTGTGGGCTTTG=0A=
GCCCAGGGAGACTTGCTTCTGAGTTTTTTCATTCACATGGCCACATCGTC=0A=
CCTGGGGGCACACCGGGTACAGGCTTCATTGCAGAGACTTTCCAAGTGAG=0A=
TTTTCTCCCAAGTGAGATAACCTTTTATCTGTATCTTTCTTAAAACAATT=0A=
CCCATTTTATGCCCCCAAATTATGGTAAAGTATAGGTGTCTTACCTTCCC=0A=
TGATAGAGTCAACCCCTTGAGGTTGGGAACCACTTCCTCACTTGTGAATC=0A=
TCTGACATTGCAGCAGAGTTGCACATACAGTAGGTGCTCAGTAATGAAGC=0A=
AGTTATTGAACAAGTGCATGGACAGTACGAGTCATGACTCTGGCCTCTGG=0A=
TGAGGCTGGCAGTGCCCACAGAAGGTTCAGTCACCTTCTCCATCCCTTAC=0A=
ATGAATACCAACTTGCCTGGCTTTAGCATGAGTGTTTCCATAACTTATTA=0A=
CAACTATATGTTCAAATGCATCTTAAGTGAGTCATCCCTGTCCTCCATAT=0A=
TCCCTATTCTATGTCATTTGTGGATTTAATCAGAAGCAAGGGTTTATGTG=0A=
ATGCACTGAGAACCTCTAACTTCCTTCTTCATCAGGAGGAGGCTCAGAAC=0A=
CCAGAGGCTACTCAGCATTCAGGGCACGTATTCCCACAGAGAATACAAGC=0A=
TGGGGGGAAAACTCACTTGGAGACTCTCTGCATTGAAGCATGTGCCTGGT=0A=
GTGCCCCCAGGGACAATGCGGCCACGTGGATGAAAACACTCAGAAGCAAA=0A=
TCTCCCTGGGCCAAAGCCCACACTTCAAAGGCTTGTTTAGGAAAACCATA=0A=
CACTCAACCTTCTTTGTGTATTGTGAAAATAAACATTTACAGAACATTGG=0A=
AATGTGTGAAAAATGTCCTTCGTTGTGTGTTTATGGAAAACATTTGGCAA=0A=
GTATTTGCTATAAATGCTTCCTGTAAATGCTGTGAGTGAGAATCAGCAGG=0A=
AAAATGTATGAGATGAGAGGGAGCCAACTCCCACAGCAAAAGAATGTCCC=0A=
CTTCCCTTCCACACACCCCACGCGCTGCCTCCACTCTGCACTCCTCCGCC=0A=
GCCCCCCTCCCCAGATAAGGCCCAGGGGATCCGGGAGTGTGAGGGTGGGG=0A=
AACAGAGAAGGGACAATCTGTGAAGGGGCCTCCAGGCTCAGCATTCCATG=0A=
TTGGAGTGCAGAAAAGGAAACTGACTCGTGTTTCCTAGGTGATGACGCCC=0A=
CCTAATAGGATTTCCTTCCCCTTCATTCTTATTGGTTAACAGTTCTGATT=0A=
CCAAGTCCCCCATTCCAGCTGGAAACATTGCCTCAACACCCACCCGGTCT=0A=
TCCCCACCTGAAGACTGGCCTGCCTCTCAGTTAACTCTTGGACGCTTTCA=0A=
GGTGAAAAGTAAAAAGTGGAAAGCAAGAGACTCAAATATATTTCTTCAGA=0A=
ACAAGATGCATGGGACATTTCTCTTAATTCCTCAACAATGGACAAAGACC=0A=
TCTGATGATATGAACTGAAAAGGGAAGAGAGTGTTGCAGGCCACAGCAAT=0A=
GAACGCTCGTGGTTTTTACCAGGAAGCAACCTTGACCGCCTCCAGATTGC=0A=
CCTCTGGGATGAACCAGGTCCTCCCCTGCCACATGAGGTATGTAGTGGGT=0A=
ACTGCAGAGAGGCTGGCTGCTCTTGGTTCCCAAGCAGCCTCCGTGATAAA=0A=
CAGCCTCAGTGTTGAAATATCAGGGAGATTAATTATCTCGCAGATGTTAC=0A=
ATGCCAAAAAGGCTACTTTTTTTCTTCATTTCAGAATTCTTTCCTTTTTT=0A=
AAGAAATACAAGGACACAAATCCAATGTGCTTAATTAAGCCTGCCAAGGC=0A=
GAAGAGGCCCAGCTGCCTGGCTGTATTTGTATCTGCACAGGTGCCCAGGG=0A=
AACACCAGCCTCACTCTCAGCTTCATCTCCAGGTCCCCTACCTCTGCCCA=0A=
GAGGTCCTCACATCTGCACCTGACTGGAGCTAATGACAAACATCCTTACA=0A=
CAGCAACTTTTACGCTGATTTGCAAAATGATTTATTTTAAGTACGTTTCT=0A=
AAAAACATAAACGTTAAGTCAAAAGGTCTCAGCATTTCTGCCACTTANNN=0A=
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN=0A=
NNNNNNNNNNNNNNTATATATATATAAGTATTGATAAACCATAATCTAAA=0A=
AATGATGTCATTTCACACTTCCCAAAGCAATGGAGGAGTTCCCTTTAAAA=0A=
TCATTTCACCGGCTCACCTAGATTTATCATGGGAAAGAAAAATCCACAGT=0A=
CTCCTCTGATGGAAGAAGAAAGGGGCTATGGGGCCCCCCAGACTTTCCCA=0A=
ATCAGCCCCATGCTGTGGCTCTCCCCAGCCATGGACCATGACAGCCACTC=0A=
>UnnamedSeq2=0A=
GCATCTATTTTGCTGCACTCGCGGTTTTTTGCTAGGAGACCCCCCAGTGT=0A=
TCTTCTTAGGCAATACAATTTAACACTAAAATTAAGATGTTACCACTAAA=0A=
TACAGTGTCTCTCATGTTTAAAATCATTCAAGTCTTTCTTGCAGGAAGAA=0A=
TTTAAGCTTCATCAGACACTCATGGGATGAGGACTGAGGTATGAAATCTA=0A=
CATGTTTTCTGACCAGAACACATCCAACATGGCCTAAAGATTATCTTACT=0A=
CTTCTAAATTAAGATCAAATAGATTTAGAGTAACTTACTGTAGCAGGGTT=0A=
CTTAGATTTGTATTTTATACAGACGAATCAAGAAACCATTCTACTAATAT=0A=
AAAATATCAGTAGTTATGAAAATTATTTTCAACCAAATAAAAACAAATTT=0A=
TAAAAGAAAAGTTCAAAGGAAATGTAACAAAATCATTTTCCTATATCTGT=0A=
TTTGGTTCAAATTAGTTCCCAAAAACTGTCTCCCTCCTAGTCAAAGATGT=0A=
TAAAAAACAGCATATAATGACTTACATTTAAACTAAAATTAAGGTGGGAA=0A=
AATATAATTTTAAATAAGTTTAAACGAGTCAAAGACAAATGTTCTAAATT=0A=
TGAAAAAAGTTAGACAATAAAGCTAAATTAACCTAACATTTCTGAATGCT=0A=
TACCTAAACCACGTTAATTCTATACGAAACAGCCTTACTTACAAAATAGT=0A=
TACCCATAAATACAAGTAATAAGAGGATAGATTACAACCTAAATATTTTT=0A=
AAATGGTAACGTACAACCCCATTGATCTAGATGNNNNNNNNNNNNNNNNN=0A=
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN=0A=
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN=0A=
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN=0A=
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNATCTAATCCTTA=0A=
GTATTTTGAAAGAAAAAAAATACCCAACATCACACTGTTGTTTTAATTTA=0A=
AAGACCAATCTCTGTCATTCAGCATTTATGAAATACTAGTTTTACTAAAC=0A=
ATTATATGGATTTTCTAAACTCAGTTCTCCACCTAAATTTTTCATTGGAA=0A=
CTCCAAGCTTCAGAACAATGAATGAGAACATCTGTTGCATTTAAGTAATA=0A=
AAATCTTGCAGATTGCCAATCCAAATGGNNNNNNNNNNNNNNNNNNNNNN=0A=
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN=0A=
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN=0A=
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN=0A=
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN=0A=
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN=0A=
NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNTGTACTAGTTAATA=0A=
GAAGTTTCAACTGGAGCAACAATACTTATGAATGCATACATTTAAAGGTC=0A=
AGATTTATTTGCTATGTACCAAGTAAAACCAAAAAAGTCACATTCATATA=0A=
GCATCTCATTAATTGTAAACACAATACGAAAAATGTTTAAATACTAGTTA=0A=
AAATACAGTGTGTATGATACTAGCTATAAAAAAATTAATTCCCTACAAAT=0A=
TTAAAATACAATCTATTTACAACTAAATGCTACCACTTTCTCCACTCTTC=0A=
TTGGTACCAAGTAACCCCTAAAATTGAGAGGTGGCTTAGCAAACTGCGTC=0A=
AACATGCCAC=0A=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00 =3D =
2274920=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.27=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                     =0A=
Query: 63      atacaatttaacactaaaatta 84=0A=
               ||||||||||||||||||||||=0A=
Sbjct: 2173742 atacaatttaacactaaaatta 2173763=0A=
=0A=
=0A=
>ref|NT_024524.4|Hs13_24680 Homo sapiens chromosome 13 working draft =
sequence segment=0A=
          Length =3D 5607782=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.27=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                     =0A=
Query: 552     atataattttaaataagtttaa 573=0A=
               ||||||||||||||||||||||=0A=
Sbjct: 5167684 atataattttaaataagtttaa 5167705=0A=
=0A=
=0A=
>ref|NT_022776.4|Hs4_22932 Homo sapiens chromosome 4 working draft =
sequence segment=0A=
          Length =3D 1039623=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.27=0A=
 Identities =3D 22/22 (100%)=0A=
 Strand =3D Plus / Minus=0A=
=0A=
                                    =0A=
Query: 549    aaaatataattttaaataagtt 570=0A=
              ||||||||||||||||||||||=0A=
Sbjct: 960496 aaaatataattttaaataagtt 960475=0A=
=0A=
=0A=
>ref|NT_005187.5|Hs2_5344 Homo sapiens chromosome 2 working draft =
sequence segment=0A=
          Length =3D 1557795=0A=
=0A=
 Score =3D 44.1 bits (22), Expect =3D 0.27=0A=
 Identities =3D 25/26 (96%)=0A=
 Strand =3D Plus / Plus=0A=
=0A=
                                         =0A=
Query: 352     aaatatcagtagttatgaaaattatt 377=0A=
               |||||||||| |||||||||||||||=0A=
Sbjct: 1187166 aaatatcagttgttatgaaaattatt 1187191=0A=
=0A=
=0A=
  Database: genome=0A=
    Posted date:  Aug 29, 2001 12:46 PM=0A=
  Number of letters in database: 2,805,113,105=0A=
  Number of sequences in database:  3518=0A=
  =0A=
Lambda     K      H=0A=
    1.37    0.711     1.31 =0A=
=0A=
Gapped=0A=
Lambda     K      H=0A=
    1.37    0.711     1.31 =0A=
=0A=
=0A=
Matrix: blastn matrix:1 -3=0A=
Gap Penalties: Existence: 5, Extension: 2=0A=
Number of Hits to DB: 3,438,566=0A=
Number of Sequences: 3518=0A=
Number of extensions: 3438566=0A=
Number of successful extensions: 184=0A=
Number of sequences better than  1.0: 6=0A=
length of query: 1860=0A=
length of database: 2,805,113,105=0A=
effective HSP length: 21=0A=
effective length of query: 1839=0A=
effective length of database: 2,805,039,227=0A=
effective search space: 5158467138453=0A=
effective search space used: 5158467138453=0A=
T: 0=0A=
A: 40=0A=
X1: 6 (11.9 bits)=0A=
X2: 15 (29.7 bits)=0A=
S1: 12 (24.3 bits)=0A=
S2: 22 (44.1 bits)=0A=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00blastInput.fasta=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00000764 =00023547 =00000017 =0000000013017 =
07370601500 015067 =
=000=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00ustar=0000dwaring=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00users=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00000000 =00000000 =
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=
=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00=00
------=_NextPart_000_0000_01C163A0.F1D463F0--