BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	78,524	A→C	N302H (AAT→CAT)	setA →	sugar exporter SetA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	78,524	0	A	C	93.9%	201.2 / ‑3.9	67	N302H (AAT→CAT)	setA	sugar exporter SetA
Reads supporting (aligned to +/- strand): ref base A (1/0); major base C (24/38); minor base T (4/0); total (29/38)
Fisher's exact test for biased strand distribution p-value = 2.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGACT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTCTTTAATAGCCGTATGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTC  >  NC_000913/78280‑78756

aTTTGTCGCCTCGACGTTAATGTGGTCCTTCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGcc                                                                                                                                                                                                                                     <  1:724043/251‑1 (MQ=255)

aTTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGcc                                                                                                                                                                                                                                     <  1:826818/251‑1 (MQ=255)

     gcGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTTGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGtctg                                                                                                                                                                                                                                <  1:26311/250‑4 (MQ=255)

      cGCCTCGACGTTAATGTGGACCTGCAAAACCATGTACATTATTGATCTGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGGTGGGGACCGACGCTGGGCT‑GGAAATACCAGCAATGGTTCTGGGTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAAAGATGGTCATAGCAGTGGCGGCAGGGG‑GACTGTTTTAAA‑CCGG‑GTTGTTTTTATTTTATAGCCGtcttg                                                                                                                                                                                                                               >  1:365023/1‑247 (MQ=255)

        ccgcggcGTTAATGTGGACCTTCAACACCATGTACATTATTGATATGCCGTTGGGGATCAGGAGCGAGGTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGGCGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCg                                                                                                                                                                                                                             <  1:115374/245‑1 (MQ=255)

           cGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑tt                                                                                                                                                                                                                          <  2:895513/250‑1 (MQ=255)

                     gtgtgcctgcactcccATGTACATTATTGATTTGCCGTTGTGTATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGGCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGc                                                                                                                                                                                                               <  2:805934/237‑1 (MQ=255)

                        ggcctgcaaccccattgtcaatattGAAATGGCGGTGTGGGTCAGTAGCGAGGTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAAc                                                                                                                                                                                                            <  1:855559/230‑1 (MQ=255)

                             gCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCTGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACtttt                                                                                                                                                                                                        <  1:895058/250‑1 (MQ=255)

                             gCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGAGGGGGACGGCAGCTGGTCTCGGAAA‑ACCAGCAATGATTCTCGCCAGCTACTATGTCAAAAGATTTGGTACGCGGCGAATGGGGGGCATAGACGTGGCGGGAGGAGGTA‑TATTTTACAACCGG‑‑TTGATTTTATTTCCTAGCCCTCTGGCG‑Tttatgtctctgtgacttttt                                                                                                                                                                                                       >  1:947285/1‑236 (MQ=255)

                                acacCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAAc                                                                                                                                                                                                    <  1:346017/251‑1 (MQ=255)

                                     aTGTACATTATTGATATGCCGTGTTGGATCAGTAGCGAGTGAGGATTGCCAGACAAACTGGCTGGTTTCCTGATGGGGACGGCAGCTTGTCT‑GGGAATACCAGCAATGATTCTGTCTGGCTACTAGGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGGAATTTCGGCGGAAG‑TACTGTTTTGCA‑CCGG‑ATTGATGTTAGTGCATAGCCGGCTGACG‑TTGATGACGCTGCAATTTTTTAACGCTGt                                                                                                                                                                                               <  2:503253/251‑1 (MQ=255)

                                             tagtGATATGCCGTTGTGGATCAGTAGCGAGTTAGGAGTTCCAGACAAATTGTCGGGTTTCCTGATGGGGAGGGCAGCTGGTCT‑GGAAATACCAGTAATGATTCTGGCGTGCTACGATGTCATACGTTATGGTAAGCGGCGAATGATTGTCAGAGTAGTGGCTGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGAGTGTATTTCAGAGCCGGCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTTTATTTATc                                                                                                                                                                                        <  2:222350/247‑1 (MQ=255)

                                                 gATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCAt                                                                                                                                                                                   <  1:691682/251‑1 (MQ=255)

                                                    atGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCAttg                                                                                                                                                                                 <  1:445877/250‑1 (MQ=255)

                                                           tgtgGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑GACTGTTTTTCA‑CCGG‑ATTGATTTTATTTTATAGCCGCCTGGCG‑TTGATGACGGTGCCACTTTTTAACGCTGTATTTTTCGGGATTTTTGCgggg                                                                                                                                                                         >  1:617770/1‑250 (MQ=255)

                                                            gtgGATCATTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTGCCTGATGGGGACGGCAGCTTGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATTGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTa                                                                                                                                                                        <  2:852372/251‑1 (MQ=255)

                                                              ggATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATATCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCACTGTTACGGGTAt                                                                                                                                                                       <  2:427912/250‑1 (MQ=255)

                                                                     tAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGGG‑TACTGTTTTTCA‑CCGG‑GTTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTTACGCTGTATTTATCGGCATTGTTGCGGGGATTGGGATgg                                                                                                                                                               >  1:898372/1‑250 (MQ=255)

                                                                          aGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTAAA‑CCGGGTTTGATTTT‑TTTCATAGCCGTCTGGCG‑TTGATGGCGCTTCAACTTTTTTACGCTGTTTTTTTCGGCATTGTTGGGGGGAtttgggagctttgg                                                                                                                                                          >  1:185990/1‑239 (MQ=255)

                                                                               ggATTGCCAGACAAACTGGCGGGTTTCCTTATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGGAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑GACTTTTTTACA‑CCGG‑ATTGATTTTATTTTATAGCCGTCTGGCGGTTG‑TGGCGCTGCAACTTTTTAACGCTGTATTTTTCGGCATTGTTGCGGGGATTTGGATGCTATGGTTTc                                                                                                                                                      >  1:898703/1‑250 (MQ=255)

                                                                                gATTGCCAGACAAACTGGCGGGTTCCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAg                                                                                                                                                    <  1:206560/251‑1 (MQ=255)

                                                                                            aaCTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCt                                                                                                                                        <  2:898372/251‑1 (MQ=255)

                                                                                                  cGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAaga                                                                                                                                  <  1:567028/251‑1 (MQ=255)

                                                                                                       ttCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCTTTGTTGCGGGTATTGGGGTGCTATGGTTTCAGGATTTAATGCCTGGAAGGGCggg                                                                                                                             >  2:411814/1‑251 (MQ=255)

                                                                                                           tGATGGGGAAGGCAGCTGGTCT‑GAAAATACCAGCAATGATTATGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGTG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTAATAACTGTCTGGAG‑TTGGTGACGGTGCAAATTTTTAACGCTGTTTTTATCGGCATTGTTGAGGGGATTGGGATGCTATGGTTTCAGGAGTTCATGCCTGGAAGGGGGGGGGc                                                                                                                          >  2:913607/1‑250 (MQ=255)

                                                                                                            gATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAAAGTTGTATTTATCGGATTTGTTGCGGGGATTGGGATGGTATGGTTTCAGGATTTAATGCCTGGAAGAGCggggga                                                                                                                          >  1:698149/1‑248 (MQ=255)

                                                                                                               ggggatggccggtgaact‑gtaAATACCTGAAATGATTGTGGCTGGCGTCTAGGTCAAACGTTATGGTAAGTGGCGAATGATGGTTATAGCAGTGGCGGCAGGAG‑TACCGTTTTACA‑CCGG‑ATTGAGTGTATTTCATAGCCGTCTGGCG‑TTGATGACGCTTCAACTTTTTAACGCTGTATGTATCGGCATTGTTGCGGGGATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCtt                                                                                                                     <  2:575355/236‑2 (MQ=255)

                                                                                                                gggACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGACTGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGATGCGGGTATTGGGTGGCTTTGGTTTCAGGTGTTAATGCCTggtaggaggggggaagctac                                                                                                                    >  2:474392/1‑233 (MQ=255)

                                                                                                                gggACGGCAGCTGGGCT‑GGAAATACCAGCAATGATGCTGGCTGGCTACTATGTCAAACGTTATTGTAAGCGGCGAAGGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTac                                                                                                                    <  2:117932/251‑1 (MQ=255)

                                                                                                                     ggCAGCTGGTAT‑GGAAATACCAGCAATGATTCTGGCTGGCTAATATGTCAACCGTTATGGTAAGCGGCGAATGATGGTAATAGCCGTGGCGGCAGGGG‑TGCTGTTTTCCA‑CAGG‑ATTGCTTTTGTTTCATAGCAGTCTGGCG‑TTGGTGACGCAGCACCTTTTTAATGCTTTATTTATTGGCATTGGTGAGGGTATTGGGCTGATATGGTTTCATGGTTTAATGACTGGAAGAGCGGGagcagcaaaaaact                                                                                                               >  2:874838/1‑243 (MQ=255)

                                                                                                                        aGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGGCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATTGTTTCAGGATTTAATGCCTGGAAGGGCGGGGGCAGCTACCCCCTTAt                                                                                                            >  2:304552/1‑251 (MQ=255)

                                                                                                                          ctggtct‑ggaAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATGGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCgg                                                                                                                              <  1:316024/231‑1 (MQ=255)

                                                                                                                          ctggtct‑ggaAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGAAGGAG‑TACTATTTTACC‑CCGG‑CTTGATTTTATTTCATAGCCGTATGGCG‑TTGATTACGCTGCTAATTTTTAACGCTGTATTTATCGGCATTGTTGGGGGTATTGGGATGATTTGGATTTAGGTATTAATGCATGGAAGAGAgggagtgagaagaggggagtgt                                                                                                          >  2:316024/1‑232 (MQ=255)

                                                                                                                           tggtct‑ggaAATACCAGCAATGATTCTAGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑AACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTAGCG‑TTGATGACGCTGCAACTTTTTTACGCAGTATTTATAGGCATTGTTGCGGGTGTTGGGATGCTATGGTTTCAGGATTTAATACCTGGAAGAGCGGGGGCAGcataccccttctttt                                                                                                         >  2:80177/5‑237 (MQ=255)

                                                                                                                              tct‑ggaAATACCAGCAATGATTCTGGCTGGCTACTATGTCAACCGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACAGTTTTACA‑CCGG‑GTTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTTTTTATCGGCCATGTTGCGGGGATTGGGAAGCTATGGGTTCAGGGTTTAAATCCTGGAAGGGCGGGGGGAGCTacaccctttattact                                                                                                       >  2:911051/2‑237 (MQ=255)

                                                                                                                                  ggaAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAAc                                                                                                    >  1:265185/1‑250 (MQ=255)

                                                                                                                                      aTACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAAAGCTGTATTTATCGGCATTGTTGCGGGTATTGGGTTGTTATTGTTTCAGGATTTATTGCCAGGAAGAGCGGGGGCAGGTACCACCTTATTTACTAACAGttt                                                                                               >  2:572793/1‑249 (MQ=255)

                                                                                                                                          cAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGGTGCTATGGTTTCAGGATTTAATGCCTGGCAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCt                                                                                           >  1:210756/1‑251 (MQ=255)

                                                                                                                                          cAGCAATGATTCGGGCTGGCTACTATGTCAAACGGTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCt                                                                                           <  1:298967/251‑1 (MQ=255)

                                                                                                                                           aGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTa                                                                                          <  1:758944/251‑1 (MQ=255)

                                                                                                                                             cAATGAGTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTAcc                                                                                        <  2:898703/251‑1 (MQ=255)

                                                                                                                                              aaTGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCg                                                                                       <  2:400660/251‑1 (MQ=255)

                                                                                                                                                tGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCACCTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGCTGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAACTACCACCTTATTTACTAACAGTATTTCTACCggg                                                                                     >  1:347714/1‑251 (MQ=255)

                                                                                                                                                  aTTCTGGCTGGCTACTCTGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑GACTGGTTTACA‑CCGG‑GTTGATTTTATTTTATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTTTTTATCGGCATTGTTGCGGGTATTGGCGTGCTTTGGTTTCAGGATTTAATTCCTGGAAGAGCGGGGGCAGGTACCACCTTATTTACTAAAAGTATTTATACCggggg                                                                                   >  2:188578/1‑250 (MQ=255)

                                                                                                                                                  aTTCTGGCTGGCTACTATGTCAAACGTTATGGAAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCCTTGTTGCGGGAATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTAACACCTTTTTTACTAACAGTATTTCTACCGGGGt                                                                                   >  1:949142/1‑251 (MQ=255)

                                                                                                                                                       ggctggctACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTaa                                                                                 <  1:109688/248‑1 (MQ=255)

                                                                                                                                                       ggctggctACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTc                                                                              >  1:787733/1‑251 (MQ=255)

                                                                                                                                                                         aCGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCTGGGGAATTTTGGCTGGCGTTATTCAgg                                                            >  2:659957/1‑251 (MQ=255)

                                                                                                                                                                              aTGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCa                                                       <  2:117086/251‑1 (MQ=255)

                                                                                                                                                                               tGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCCTTGTTGCGGGTATTGGGATGCTATGGTTTCAGGTTTTCATGCCTGGAAGAGCGGGGGCCGCTACCTCCTTATTTTCTAACAGTATTCCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCa                                                       >  2:348915/1‑250 (MQ=255)

                                                                                                                                                                                     gcggcgAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTTTTTATCGGCCTTTTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAAAAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTccgg                                                            >  2:891862/1‑236 (MQ=255)

                                                                                                                                                                                     gcggcgAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAgg                                                            <  1:891862/239‑1 (MQ=255)

                                                                                                                                                                                      cggcgAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGcac                                                >  1:656708/1‑250 (MQ=255)

                                                                                                                                                                                       ggcgAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTTCA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACaa                                              <  1:71345/251‑1 (MQ=255)

                                                                                                                                                                                            atgatgGTCATAGCAGGTGCGGCATGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTGTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGTGCAGCGTCCACCTTATTTACTAACAGTATTTCTACCTGGGTAATTCTGGCTGGCGGTATTCAGGGAGCAATGGCACAAATTTg                                         <  2:626798/251‑1 (MQ=255)

                                                                                                                                                                                                   tCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCAc                                  >  1:964614/1‑251 (MQ=255)

                                                                                                                                                                                                       aGCAGTGGCGGTAGGAG‑TACTGTGTTCCA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTAGCGGCAGTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACttt                               <  2:617770/250‑1 (MQ=255)

                                                                                                                                                                                                       aGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACttt                               <  1:523889/250‑1 (MQ=255)

                                                                                                                                                                                                           gTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGt                          <  1:587185/251‑1 (MQ=255)

                                                                                                                                                                                                             ggcggcAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTc                         <  1:411814/250‑1 (MQ=255)

                                                                                                                                                                                                                             gTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAAGTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCg          <  1:401280/250‑1 (MQ=255)

                                                                                                                                                                                                                                ttACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGtt       <  1:304552/250‑1 (MQ=255)

                                                                                                                                                                                                                                  aca‑cCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTAt     <  1:230120/250‑1 (MQ=255)

                                                                                                                                                                                                                                  aca‑cCGG‑AATGATTTTATTTCATAACCGGCTGGCG‑TTGATGACGCTTCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACAGGGGTAATTCTGGCTGGCGTTATTCATGGAGCAATTGCACACAGTTGGGGGCACTTTGCTGTATACTGGGTAATTGCGGTTAtt    >  2:716131/1‑251 (MQ=255)

                                                                                                                                                                                                                                   ca‑cCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCATTTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTAttt   >  2:370498/1‑251 (MQ=255)

                                                                                                                                                                                                                                    a‑cCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTGTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTTCCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGTGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGTGTAATTGCGGTTATTTc  <  2:538369/251‑1 (MQ=255)

ATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGACT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTCTTTAATAGCCGTATGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTC  >  NC_000913/78280‑78756

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`