BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,145,905	G→A	S880S (TCG→TCA)	dgcE →	putative diguanylate cyclase DgcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,145,905	0	G	A	100.0%	1222.5 / NA	291	S880S (TCG→TCA)	dgcE	putative diguanylate cyclase DgcE
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (146/145); total (146/145)

Only 100 of 295 total aligned reads displayed.

GCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCT  >  NC_000913/2145807‑2146002

gCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCAcc                                                                                                 <  1:3123353/101‑1 (MQ=255)

  tGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTcaccac                                                                                               <  2:3123423/101‑1 (MQ=255)

  tGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTcaccac                                                                                               >  1:3123809/1‑101 (MQ=255)

   gCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACg                                                                                              <  1:3123318/101‑1 (MQ=255)

     cGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGa                                                                                            <  2:3123379/101‑1 (MQ=255)

     cGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGa                                                                                            >  2:3123811/1‑101 (MQ=255)

         caTAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATccc                                                                                        <  1:3123390/101‑1 (MQ=255)

          aTAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCg                                                                                       >  1:3123816/1‑101 (MQ=255)

            aGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGa                                                                                     <  2:3123408/101‑1 (MQ=255)

            aGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGa                                                                                     >  2:3123818/1‑101 (MQ=255)

             gCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGaa                                                                                    <  1:3123330/101‑1 (MQ=255)

             gCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGaa                                                                                    <  2:3123359/101‑1 (MQ=255)

             gCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGaa                                                                                    <  1:3123382/101‑1 (MQ=255)

              cGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAg                                                                                   >  2:3123821/1‑101 (MQ=255)

                   gggCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAgcgc                                                                                >  2:3123825/1‑99 (MQ=255)

                     gcggcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAgcgc                                                                                >  1:3123827/1‑97 (MQ=255)

                      cggcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGa                                                                                            <  2:3123479/84‑1 (MQ=255)

                       ggcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAAtt                                                                          <  2:3123618/101‑1 (MQ=255)

                       ggcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAAtt                                                                          >  1:3123830/1‑101 (MQ=255)

                       ggcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAAtt                                                                          <  1:3123702/101‑1 (MQ=255)

                        gcgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAAttt                                                                         <  1:3123299/101‑1 (MQ=255)

                         cgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTtg                                                                        >  2:3123835/1‑101 (MQ=255)

                         cgATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTtg                                                                        >  1:3123833/1‑101 (MQ=255)

                          gATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTtgt                                                                       <  2:3123351/101‑1 (MQ=255)

                            tGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTgg                                                                     >  1:3123841/1‑101 (MQ=255)

                            tGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTgg                                                                     >  2:3123840/1‑101 (MQ=255)

                              tCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCt                                                                   >  1:3123845/1‑101 (MQ=255)

                                gCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGa                                                                 >  1:3123848/1‑101 (MQ=255)

                                  ttGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGAtt                                                               >  1:3123851/1‑101 (MQ=255)

                                   tgatgaACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGAttt                                                              <  2:3123397/101‑1 (MQ=255)

                                   tgatgaACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGAttt                                                              >  1:3123855/1‑101 (MQ=255)

                                     atgaACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCa                                                            >  2:3123858/1‑101 (MQ=255)

                                      tgaACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCAc                                                           >  1:3123860/1‑101 (MQ=255)

                                        aaCAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACtt                                                         <  1:3123527/101‑1 (MQ=255)

                                        aaCAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACtt                                                         <  2:3123631/101‑1 (MQ=255)

                                         aCAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTa                                                        <  2:3123449/101‑1 (MQ=255)

                                          cAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa                                                       <  1:3123574/101‑1 (MQ=255)

                                          cAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTaa                                                       >  1:3123864/1‑101 (MQ=255)

                                           aGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTaa                                                                            >  2:3123867/1‑79 (MQ=255)

                                           aGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAg                                                      <  1:3123489/101‑1 (MQ=255)

                                              ggCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGctc                                                   <  2:3123570/101‑1 (MQ=255)

                                              ggCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGctc                                                   <  2:3123233/101‑1 (MQ=255)

                                               gCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGctct                                                  >  2:3123871/1‑101 (MQ=255)

                                               gCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGctct                                                  <  2:3123285/101‑1 (MQ=255)

                                                cGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTg                                                 >  2:3123874/1‑101 (MQ=255)

                                                cGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTg                                                 >  1:3123875/1‑101 (MQ=255)

                                                cGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTg                                                 <  2:3123490/101‑1 (MQ=255)

                                                  gatgatTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGa                                               >  1:3123878/1‑101 (MQ=255)

                                                   atgatTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAg                                              <  1:3123422/101‑1 (MQ=255)

                                                    tgatTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGt                                             >  1:3123882/1‑101 (MQ=255)

                                                     gatTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGtt                                            <  1:3123475/101‑1 (MQ=255)

                                                      atTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTg                                           >  2:3123884/1‑101 (MQ=255)

                                                      atTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTg                                           <  1:3123539/101‑1 (MQ=255)

                                                       ttAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGc                                          <  1:3123615/101‑1 (MQ=255)

                                                        tAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCg                                         >  1:3123888/1‑101 (MQ=255)

                                                        tAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCg                                         <  1:3123446/101‑1 (MQ=255)

                                                          aaGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGaa                                       <  2:3123503/101‑1 (MQ=255)

                                                          aaGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGaa                                       >  2:3123890/1‑101 (MQ=255)

                                                          aaGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGaa                                       <  1:3123222/101‑1 (MQ=255)

                                                           agagAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAg                                      <  2:3123442/101‑1 (MQ=255)

                                                           agagAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAg                                      >  2:3123893/1‑101 (MQ=255)

                                                             agaATCAGTTGATGATGCTCGCCCACGGTGTCACTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTCAAGCTCTGGAGTTGCGAAGGc                                    <  1:3123342/101‑1 (MQ=255)

                                                              gaATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCg                                   >  1:3123897/1‑101 (MQ=255)

                                                               aaTCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCga                                  >  1:3123900/1‑101 (MQ=255)

                                                                aTCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTCGCGAAGGCgag                                 <  2:3123430/101‑1 (MQ=255)

                                                                 tCAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCgaga                                >  1:3123903/1‑101 (MQ=255)

                                                                  cAGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGat                               >  2:3123904/1‑101 (MQ=255)

                                                                   aGTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGatt                              <  1:3123752/101‑1 (MQ=255)

                                                                    gTTGATGATGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGatta                             <  2:3123512/101‑1 (MQ=255)

                                                                    gTTGATGATGCTCGCCCACGGGGTCGCTTCACCACGGATCCCGGa                                                                                     >  2:3123908/1‑45 (MQ=255)

                                                                      tgatgaTGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGattatt                           <  2:3123306/101‑1 (MQ=255)

                                                                      tgatgaTGCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGattatt                           <  2:3123462/101‑1 (MQ=255)

                                                                         tgatgCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtgat                        >  2:3123914/1‑101 (MQ=255)

                                                                         tgatgCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtgat                        >  1:3123912/1‑101 (MQ=255)

                                                                           atgCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtgatga                      <  1:3123399/101‑1 (MQ=255)

                                                                             gCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaac                    >  2:3123918/1‑101 (MQ=255)

                                                                             gCTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaac                    >  1:3123919/1‑101 (MQ=255)

                                                                              cTCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaaca                   <  1:3123021/101‑1 (MQ=255)

                                                                               tCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCa                                                            <  2:3123597/59‑1 (MQ=255)

                                                                               tCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACtt                                                         >  1:3123925/1‑62 (MQ=255)

                                                                               tCGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGatt                              >  2:3123924/1‑89 (MQ=255)

                                                                                cGCCCACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaacaaa                 <  1:3123648/101‑1 (MQ=255)

                                                                                   ccACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAt              >  1:3123930/1‑101 (MQ=255)

                                                                                   ccACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAt              >  2:3123931/1‑101 (MQ=255)

                                                                                    cACGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAtt             <  2:3123466/101‑1 (MQ=255)

                                                                                     aCGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAttt            >  2:3123937/1‑101 (MQ=255)

                                                                                     aCGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAttt            >  1:3123936/1‑101 (MQ=255)

                                                                                     aCGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAttt            <  1:3123428/101‑1 (MQ=255)

                                                                                      cGGTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTtc           <  2:3123553/101‑1 (MQ=255)

                                                                                       ggTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTtcg          >  2:3123942/1‑101 (MQ=255)

                                                                                       ggTGTCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTtcg          >  1:3123943/1‑101 (MQ=255)

                                                                                        gtgtCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTtcgt         <  2:3123679/101‑1 (MQ=255)

                                                                                        gtgtCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTtcgt         <  1:3123510/101‑1 (MQ=255)

                                                                                         tgtCGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTtcgtc        <  2:3123525/101‑1 (MQ=255)

                                                                                            cGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTa     >  1:3123948/1‑101 (MQ=255)

                                                                                            cGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTa     <  1:3123587/101‑1 (MQ=255)

                                                                                            cGCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTCAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTa     <  2:3123586/101‑1 (MQ=255)

                                                                                             gCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAg    >  2:3123949/1‑101 (MQ=255)

                                                                                             gCTTCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAg    <  2:3123538/101‑1 (MQ=255)

                                                                                               ttCACCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGGTTATTGATGAACAAACATTTCGTCGTAGCt  >  1:3123953/1‑101 (MQ=255)

GCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGCGGATGATTAAAGAGAATCAGTTGATGATGCTCGCCCACGGTGTCGCTTCGCCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCT  >  NC_000913/2145807‑2146002

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