BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,516,689	T→C	100%	T382A (ACG→GCG)	yfeA ←	putative diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,516,689	0	T	C	100.0%	430.8 / ‑6.7	124	T382A (ACG→GCG)	yfeA	putative diguanylate cyclase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (65/58); total (65/59)

Only 100 of 124 total aligned reads displayed.

AGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCT  >  NC_000913/2516592‑2516784

aGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCAc                                                                                               <  1:930487/100‑1 (MQ=255)

   cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                            >  1:1685974/1‑100 (MQ=255)

    gCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCa                                                                                           >  1:1853526/1‑100 (MQ=255)

     cccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCaa                                                                                          >  2:1142920/1‑100 (MQ=255)

      ccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCaaa                                                                                         <  1:4109585/100‑1 (MQ=255)

      ccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCaaa                                                                                         >  1:960624/1‑100 (MQ=255)

      cattCCGTTTCCGGCACGCTCTGTACTAACAGCAGTTCACTACCCGGTAATTGATGCAACTTTTCGTTTTCCTGCATCAACGGCAGGAGCGCACGGaaaa                                                                                         >  1:3647969/3‑100 (MQ=255)

       cTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAAt                                                                                        <  1:766671/100‑1 (MQ=255)

       cTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAAt                                                                                        >  2:353351/1‑100 (MQ=255)

       cTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAAt                                                                                        <  1:4354658/100‑1 (MQ=255)

         tCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCACTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTg                                                                                      >  1:661228/1‑100 (MQ=255)

           cGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAg                                                                                    >  1:3097543/1‑100 (MQ=255)

            gTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGc                                                                                   <  2:2359049/100‑1 (MQ=255)

            gTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGc                                                                                   >  1:3924710/1‑100 (MQ=255)

            gTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGGTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGa                                                                                     >  1:1642410/1‑98 (MQ=255)

             tttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCg                                                                                  >  2:4311341/1‑100 (MQ=255)

             tttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCg                                                                                  <  1:312335/100‑1 (MQ=255)

               tCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGa                                                                                >  2:375918/1‑100 (MQ=255)

                 nggCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                              >  1:900173/2‑100 (MQ=255)

                 cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                              >  2:3653758/1‑100 (MQ=255)

                 cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                              >  1:3963180/1‑100 (MQ=255)

                   gCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACa                                                                            >  1:116609/1‑100 (MQ=255)

                      cGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTg                                                                         <  2:357257/100‑1 (MQ=255)

                      cGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTg                                                                         >  1:1925334/1‑100 (MQ=255)

                       gCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGa                                                                        <  2:532617/100‑1 (MQ=255)

                        cTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGaa                                                                       >  1:2961386/1‑100 (MQ=255)

                         tCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAc                                                                      <  2:3793888/100‑1 (MQ=255)

                         tCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAc                                                                      <  1:1415668/100‑1 (MQ=255)

                          cAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcg                                                                     >  2:1390857/1‑100 (MQ=255)

                           aGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgc                                                                    <  2:1853533/100‑1 (MQ=255)

                           aGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgc                                                                    >  1:3475055/1‑100 (MQ=255)

                            gCACTAACCGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGGAGCAGCGTACGGCAAATTGAGCGGCTACAGTGAAcgcg                                                                   <  1:2016506/100‑1 (MQ=255)

                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGaa                                                                       >  1:3961163/1‑96 (MQ=255)

                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                   <  2:4204059/100‑1 (MQ=255)

                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                   >  2:690915/1‑100 (MQ=255)

                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                   >  2:2567257/1‑100 (MQ=255)

                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAACGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                   <  1:2461761/100‑1 (MQ=255)

                              aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                 <  2:2263876/100‑1 (MQ=255)

                              aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                 <  1:2798629/100‑1 (MQ=255)

                              aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                 <  2:1720242/100‑1 (MQ=255)

                              aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                 <  1:575138/100‑1 (MQ=255)

                              aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                 <  2:4060557/100‑1 (MQ=255)

                               cTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                 >  2:3259287/1‑99 (MQ=255)

                                tAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcacc                                                               <  2:2999650/100‑3 (MQ=255)

                                 aaCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatca                                                              >  2:336973/1‑100 (MQ=255)

                                    agcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTa                                                           <  1:1328160/100‑1 (MQ=255)

                                     gtagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTaa                                                          <  1:2666904/98‑1 (MQ=255)

                                     gcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTaa                                                          >  2:1955332/1‑100 (MQ=255)

                                        gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGcc                                                       >  1:2948499/1‑100 (MQ=255)

                                        gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGcc                                                       >  2:4273599/1‑100 (MQ=255)

                                        gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGcc                                                       <  2:1940016/100‑1 (MQ=255)

                                        gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                        >  2:852645/1‑99 (MQ=255)

                                        gTTCACTACCCAGCAATTGCTACAACTTTTCGTTTTCCTGAATCTACGGCAGCAGCGCACGGCACATCGCGCGGATACAGTGAACGTCCATCATTAAGcc                                                       >  1:4175111/1‑100 (MQ=255)

                                         ttCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCg                                                      <  2:4395103/100‑1 (MQ=255)

                                          tCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGt                                                     <  2:2723648/100‑1 (MQ=255)

                                          tCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGt                                                     <  1:375557/100‑1 (MQ=255)

                                           cACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTa                                                    <  1:3624369/100‑1 (MQ=255)

                                            aCTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTaa                                                   >  2:470785/1‑100 (MQ=255)

                                            aCTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTaa                                                   >  2:456020/1‑100 (MQ=255)

                                            aCTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTaa                                                   >  2:3376694/1‑100 (MQ=255)

                                            aCTACCCGGCAATTGATACAACGTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTaa                                                   <  1:1028368/100‑1 (MQ=255)

                                             cTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAAt                                                  >  2:3251221/1‑100 (MQ=255)

                                               aCCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATGATTAAGCCGTAATga                                                >  2:4348964/1‑100 (MQ=255)

                                               aCCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATga                                                <  1:3691653/100‑1 (MQ=255)

                                               aCCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATga                                                <  2:2105826/100‑1 (MQ=255)

                                                cccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgac                                               >  2:4330598/1‑100 (MQ=255)

                                                 ccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacg                                              <  1:3802552/100‑1 (MQ=255)

                                                   ggCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACTGTGAACGCGCATCATTAAGCCGTAATgacgac                                            <  2:4185891/100‑1 (MQ=255)

                                                   ggCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacgac                                            <  2:2564790/100‑1 (MQ=255)

                                                   ggCAAGTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacgac                                            <  2:1735354/100‑1 (MQ=255)

                                                    gCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCCCACGGCAAATTGAGCGGATACAGTGAACGCGCATCCTTAAGCCGTAATgacgac                                            >  1:2847676/1‑99 (MQ=255)

                                                      aaTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCa                                         <  2:4126768/100‑1 (MQ=255)

                                                          gggaCAACGTTGCGGTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATaaa                                     <  1:3450624/97‑1 (MQ=255)

                                                          gATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATaaa                                     <  2:1339605/100‑1 (MQ=255)

                                                           aTACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                    <  1:381024/100‑1 (MQ=255)

                                                            tACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACt                                   >  2:4297902/1‑100 (MQ=255)

                                                             aCAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCACACGTCAAATTGAGCGGATACAGTGAACGCGCACCATTAAGCCGTAATGAAGACTCATaaa                                     >  1:2781353/1‑97 (MQ=255)

                                                              cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCa                                 >  2:30342/1‑100 (MQ=255)

                                                                aCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAg                               <  2:324849/100‑1 (MQ=255)

                                                                aCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAg                               >  2:3090091/1‑100 (MQ=255)

                                                                 cTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGa                              <  2:238350/100‑1 (MQ=255)

                                                                  ttttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGAt                             <  2:3475078/100‑1 (MQ=255)

                                                                  ttttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATCGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGAt                             >  2:3270213/1‑100 (MQ=255)

                                                                   tttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGaag                            >  1:238350/1‑98 (MQ=255)

                                                                   tatCGTTTTCCTGAATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGAtt                            >  2:2974531/3‑100 (MQ=255)

                                                                    ttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTa                           <  2:557525/100‑1 (MQ=255)

                                                                    ttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTa                           <  1:2669280/100‑1 (MQ=255)

                                                                       gTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCa                        >  1:499713/1‑100 (MQ=255)

                                                                         tttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATCACGACTCATAAACTCAAGATTATCAAt                      >  1:3709450/1‑100 (MQ=255)

                                                                              tGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgca                 <  2:3464094/100‑1 (MQ=255)

                                                                                cATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcagg               >  1:3235219/1‑100 (MQ=255)

                                                                                 aTCAACGGCAGCAGCGCACGGCAAGTTGAGCGGATACAGTGAACGCGCAGCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc              <  1:2346222/100‑1 (MQ=255)

                                                                                  tCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggca             <  2:2719203/100‑1 (MQ=255)

                                                                                    aaCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagc           >  1:1695514/1‑100 (MQ=255)

                                                                                    aaCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCAGTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggaa             >  1:2965310/1‑96 (MQ=255)

                                                                                    aaCGGCAGCAGCGCAAGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACACAAGATTATCAATGCGCAGgcagc           >  1:1081143/1‑100 (MQ=255)

                                                                                        gcagcaGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaaa       >  1:2000008/1‑100 (MQ=255)

                                                                                         cagcagCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaaaa      >  1:2393724/1‑100 (MQ=255)

                                                                                             agcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActct  <  1:1955343/100‑1 (MQ=255)

                                                                                             agcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActct  <  2:1888024/100‑1 (MQ=255)

AGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCT  >  NC_000913/2516592‑2516784

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