BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	78,533	A→C	M305L (ATG→CTG)	setA →	sugar exporter SetA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	78,533	0	A	C	98.5%	183.6 / NA	68	M305L (ATG→CTG)	setA	sugar exporter SetA
Reads supporting (aligned to +/- strand): ref base A (1/0); new base C (30/37); total (31/37)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGACT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTCTTTAATAGCCGTATGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAAC  >  NC_000913/78288‑78777

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)

                                                                                                                                                                                                                                          ttatttcaTAGCCGTCTGGTG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGGTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGGTATTCAGTGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGGCTACTGGGTAATTGCGGTTATTTCTGTTATCGCat            <  1:913607/243‑1 (MQ=255)

                                                                                                                                                                                                                                           tatttcaTAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCatt           <  2:194276/244‑1 (MQ=255)

                                                                                                                                                                                                                                                caTAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTAttttt     >  1:898630/2‑251 (MQ=255)

                                                                                                                                                                                                                                                   aGCCGTCTGGCG‑TTGATGACGCTGCACCTTTTTAACGCTGTATTTATCGGCATTGATGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGCAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCCTTATTTTTAAc  >  1:857693/1‑251 (MQ=255)

                                                                                                                                                                                                                                                   aGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGATATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGTTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGAATTATTTTTAAc  >  2:802764/1‑251 (MQ=255)

CCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGACT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTCTTTAATAGCCGTATGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAAC  >  NC_000913/78288‑78777

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: `‑`