BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,735,733	(G)_6→7	intergenic (+67/‑32)	pheL → / → pheA	PheL/chorismate mutase‑P

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,735,727	1	.	G	94.9%	668.2 / 26.3	217	intergenic (+61/‑38)	pheL/pheA	PheL/chorismate mutase‑P
Reads supporting (aligned to +/- strand): ref base . (9/2); new base G (96/110); total (105/112)
Fisher's exact test for biased strand distribution p-value = 2.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

TTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTG  >  NC_000913/2735635‑2735815

 tttCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑g                                                                                         <  2:5007836/93‑1 (MQ=255)

 tttCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑g                                                                                         <  2:5915761/93‑1 (MQ=255)

 tttCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑g                                                                                         <  2:2880783/93‑1 (MQ=255)

  ttcttcGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑gg                                                                                        >  1:707007/1‑93 (MQ=255)

   tcttcGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑ggg                                                                                       >  1:4995002/1‑93 (MQ=255)

    cttcGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑gggg                                                                                      >  2:3320656/1‑93 (MQ=255)

     ttcGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑ggggg                                                                                     <  1:349695/93‑1 (MQ=255)

     ttcGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑ggggg                                                                                     <  1:5960987/93‑1 (MQ=255)

      tcGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑gggggg                                                                                    <  1:185538/93‑1 (MQ=255)

       cccATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑ggggggg                                                                                   <  1:4932004/91‑2 (MQ=255)

        gCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑gggggggc                                                                                  >  1:5530233/1‑91 (MQ=255)

        gCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑gggggggc                                                                                  >  1:1306963/1‑91 (MQ=255)

         cATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑gggggggcc                                                                                 >  2:4588952/1‑90 (MQ=255)

          aTTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑gggggggcct                                                                                >  2:4431919/1‑89 (MQ=255)

            tCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑GGGGGGGCcttt                                                                              >  1:4706677/1‑93 (MQ=255)

            tCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑GGGGGGGCcttt                                                                              >  2:562920/1‑93 (MQ=255)

             cTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑GGGGGGGCctttt                                                                             >  1:3162911/1‑93 (MQ=255)

               tttttACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCtttttt                                                                            >  1:4333011/1‑93 (MQ=255)

                ttttACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTa                                                                           <  2:1077102/93‑1 (MQ=255)

                  ttACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTAtt                                                                         <  2:429235/93‑1 (MQ=255)

                  ttACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTAtt                                                                         <  2:2740206/93‑1 (MQ=255)

                   tACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTg                                                                        >  2:534003/1‑93 (MQ=255)

                   tACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTg                                                                        >  2:3759440/1‑93 (MQ=255)

                     ccTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGAt                                                                      >  1:4439816/1‑93 (MQ=255)

                     ccTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGAt                                                                      <  2:446588/93‑1 (MQ=255)

                      cTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATa                                                                     >  2:1591683/1‑93 (MQ=255)

                      cTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATa                                                                     <  2:6057726/93‑1 (MQ=255)

                      cTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATa                                                                     <  1:5733854/93‑1 (MQ=255)

                       ttCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTTATaa                                                                    >  2:2267937/1‑93 (MQ=255)

                       ttCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATaa                                                                    >  2:2692506/1‑93 (MQ=255)

                       ttCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATaa                                                                    >  2:2994308/1‑93 (MQ=255)

                       ttCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATaa                                                                    >  1:1438718/1‑93 (MQ=255)

                       ttCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATaa                                                                    <  1:1771856/93‑1 (MQ=255)

                        tCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATaaa                                                                   >  2:5487570/1‑92 (MQ=255)

                        tCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAAc                                                                   <  1:4719939/93‑1 (MQ=255)

                           ccTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACaaa                                                                <  2:5237350/93‑1 (MQ=255)

                           ccTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACaaa                                                                <  1:2995258/93‑1 (MQ=255)

                           ccTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACaaa                                                                >  2:3668294/1‑93 (MQ=255)

                            cTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACaaaa                                                               >  1:2063898/1‑93 (MQ=255)

                            cTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACaaaa                                                               >  2:2373732/1‑93 (MQ=255)

                            cTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACaaaa                                                               >  2:4148956/1‑93 (MQ=255)

                             tGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACaaaaa                                                              >  2:1101373/1‑93 (MQ=255)

                              gAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAg                                                             >  2:18540/1‑93 (MQ=255)

                               aaTGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAgg                                                            <  2:428721/93‑1 (MQ=255)

                               aaTGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAgg                                                            <  1:6224053/93‑1 (MQ=255)

                                 tGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCa                                                          <  2:5802936/93‑1 (MQ=255)

                                  gggAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCaa                                                         <  1:4920933/93‑1 (MQ=255)

                                  gggAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCaa                                                         <  2:1961735/93‑1 (MQ=255)

                                  gggAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCaa                                                         >  2:1690789/1‑93 (MQ=255)

                                  gggAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCaa                                                         <  1:6051597/93‑1 (MQ=255)

                                    gAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAaca                                                       <  1:5119935/93‑1 (MQ=255)

                                     aGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAacac                                                      >  1:4812706/1‑93 (MQ=255)

                                     aGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAacac                                                      <  2:2156457/93‑1 (MQ=255)

                                        cGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTAt                                                   <  2:6033303/93‑1 (MQ=255)

                                         gTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATg                                                  >  2:621700/1‑93 (MQ=255)

                                          tttCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATGG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGa                                                 <  1:4442008/93‑1 (MQ=255)

                                          tttCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGa                                                 <  2:7882/93‑1 (MQ=255)

                                          tttCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGa                                                 >  2:3486784/1‑93 (MQ=255)

                                          tatcgtcgTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGa                                                 >  2:3976108/3‑93 (MQ=255)

                                            tcgtcgTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACa                                               >  1:315521/1‑93 (MQ=255)

                                              gtcgtGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATc                                             >  2:5360900/1‑93 (MQ=255)

                                               tcgtGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCg                                            >  2:479533/1‑93 (MQ=255)

                                               tcgtGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCg                                            <  2:4651718/93‑1 (MQ=255)

                                               tcgtGTGAAACAGAACGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATc                                             <  1:854214/93‑1 (MQ=255)

                                                 gtgtGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGa                                          >  1:6198271/1‑93 (MQ=255)

                                                 gtgtGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGa                                          <  1:1508718/93‑1 (MQ=255)

                                                  tgtgAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGaa                                         <  1:3087169/93‑1 (MQ=255)

                                                    tgAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGaaaa                                       <  2:4141839/93‑1 (MQ=255)

                                                     gAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAAc                                      <  1:4028441/93‑1 (MQ=255)

                                                      aaaCAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAAcc                                     >  1:2419356/1‑93 (MQ=255)

                                                         cAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGt                                  >  1:3435669/1‑93 (MQ=255)

                                                          aGAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGtt                                 <  2:2956919/93‑1 (MQ=255)

                                                           gAATGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTa                                >  1:2201550/1‑93 (MQ=255)

                                                            aaTGCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTAc                               <  1:5859218/93‑1 (MQ=255)

                                                               gCGAAGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTgg                            <  2:2620472/93‑1 (MQ=255)

                                                                  aaGACGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGgcgc                         <  1:4575679/93‑1 (MQ=255)

                                                                     aCGAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGc                      <  2:5415611/93‑1 (MQ=255)

                                                                       gAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCga                    <  1:4660544/93‑1 (MQ=255)

                                                                       gAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCga                    >  2:4740698/1‑93 (MQ=255)

                                                                       gAACAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCga                    >  2:1732220/1‑93 (MQ=255)

                                                                         aCAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCgaga                  <  1:544881/93‑1 (MQ=255)

                                                                          cAATAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCgagag                 <  1:3812436/93‑1 (MQ=255)

                                                                           aaTAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCgagaga                >  1:2882402/1‑93 (MQ=255)

                                                                             tAAGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGaaa              <  1:2254186/93‑1 (MQ=255)

                                                                               aGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAAt            >  2:3244800/1‑93 (MQ=255)

                                                                               aGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAAt            <  1:1611416/93‑1 (MQ=255)

                                                                               aGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAAt            >  2:3447760/1‑93 (MQ=255)

                                                                               aGGCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAAt            <  2:3607760/93‑1 (MQ=255)

                                                                                 gCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCa          >  2:5592776/1‑93 (MQ=255)

                                                                                 gCCTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCa          >  2:1517814/1‑93 (MQ=255)

                                                                                  ccTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAg         <  2:1393708/93‑1 (MQ=255)

                                                                                   cTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgc        <  1:4796458/93‑1 (MQ=255)

                                                                                   cTCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgc        >  1:366606/1‑93 (MQ=255)

                                                                                    tCCCAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgcg       >  2:409865/1‑93 (MQ=255)

                                                                                     cccAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgcgc      <  1:1553845/93‑1 (MQ=255)

                                                                                      ccAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgcgcg     >  2:2553386/1‑93 (MQ=255)

                                                                                       cAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgcgcgc    >  1:3908209/1‑93 (MQ=255)

                                                                                       cAAATCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgcgcgc    <  1:2741813/93‑1 (MQ=255)

                                                                                        aaaTCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCt   <  1:1724111/93‑1 (MQ=255)

                                                                                        aaaTCG‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCt   >  2:849064/1‑93 (MQ=255)

TTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑‑GGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTG  >  NC_000913/2735635‑2735815

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`