BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,271,477	G→A	G251D (GGC→GAC)	pdeN →	putative c‑di‑GMP phosphodiesterase PdeN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,271,477	0	G	A	100.0%	1384.5 / NA	328	G251D (GGC→GAC)	pdeN	putative c‑di‑GMP phosphodiesterase PdeN
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (166/162); total (166/162)

Only 100 of 328 total aligned reads displayed.

GGGCATTCCTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGGCCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGCGTATCAACC  >  NC_000913/2271379‑2271573

gggCATTCCTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGAcc                                                                                                >  2:4232709/1‑101 (MQ=255)

 ggCATTCCTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACct                                                                                               <  2:4232164/101‑1 (MQ=255)

 ggCATTCCTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACct                                                                                               <  2:4232240/101‑1 (MQ=255)

  gCATTCCTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACctg                                                                                              <  2:4232357/101‑1 (MQ=255)

   cATTCCTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACctgc                                                                                             >  1:4232713/1‑101 (MQ=255)

    aTTCCTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACctgct                                                                                            <  2:4232281/101‑1 (MQ=255)

    aTTCCTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACctgct                                                                                            >  2:4232715/1‑101 (MQ=255)

      tCCTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACctgctct                                                                                          <  1:4232251/101‑1 (MQ=255)

      tCCTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACctgctct                                                                                          >  1:4232721/1‑101 (MQ=255)

      tCCTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACctgctct                                                                                          <  1:4232380/101‑1 (MQ=255)

       ccTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACctgctctg                                                                                         <  1:4232145/101‑1 (MQ=255)

        ctctGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGAATGAGTGGAACTGTCGTTGACctgctctg                                                                                         >  2:4232724/1‑100 (MQ=255)

         tctGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACctgctctgct                                                                                       <  1:4232135/101‑1 (MQ=255)

          ctGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTa                                                                                      <  2:4232215/101‑1 (MQ=255)

           tGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTAt                                                                                     >  1:4232728/1‑101 (MQ=255)

           tGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTAt                                                                                     >  2:4232730/1‑101 (MQ=255)

            gACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTAtt                                                                                    <  2:4232423/101‑1 (MQ=255)

            gACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTAtt                                                                                    <  1:4232199/101‑1 (MQ=255)

            gACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTAtt                                                                                    >  1:4232734/1‑101 (MQ=255)

             aCCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTa                                                                                   <  1:4231837/101‑1 (MQ=255)

               cGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTAcc                                                                                 >  1:4232739/1‑101 (MQ=255)

                gTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCt                                                                                >  2:4232740/1‑101 (MQ=255)

                gTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCt                                                                                >  1:4232743/1‑101 (MQ=255)

                  tCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCtga                                                                              >  2:4232749/1‑101 (MQ=255)

                    ggCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCtgatg                                                                            <  2:4231869/101‑1 (MQ=255)

                    ggCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCtgatg                                                                            >  1:4232751/1‑101 (MQ=255)

                     gCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCtgatga                                                                           <  1:4231929/101‑1 (MQ=255)

                      cTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAg                                                                          >  2:4232757/1‑101 (MQ=255)

                      cTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAg                                                                          >  1:4232756/1‑101 (MQ=255)

                       tttATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGc                                                                         >  2:4232760/1‑101 (MQ=255)

                       tttATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGc                                                                         >  2:4232764/1‑101 (MQ=255)

                         tATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGcgt                                                                       <  1:4232157/101‑1 (MQ=255)

                         tATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGcgt                                                                       <  1:4232183/101‑1 (MQ=255)

                          aTGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGcgta                                                                      >  2:4232767/1‑101 (MQ=255)

                           tGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGcgtac                                                                     <  1:4232466/101‑1 (MQ=255)

                           tGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGcgtac                                                                     <  2:4232311/101‑1 (MQ=255)

                            gCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGcgtacg                                                                    >  2:4232771/1‑101 (MQ=255)

                              aGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGcgtacgta                                                                  >  1:4232775/1‑101 (MQ=255)

                               gATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTAt                                                                 >  2:4232776/1‑101 (MQ=255)

                                aTGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATg                                                                <  2:4232480/101‑1 (MQ=255)

                                aTGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATg                                                                <  2:4232302/101‑1 (MQ=255)

                                aTGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATg                                                                >  1:4232779/1‑101 (MQ=255)

                                 tGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATgc                                                               <  1:4232244/101‑1 (MQ=255)

                                 tGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGAATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATgc                                                               >  2:4232784/1‑101 (MQ=255)

                                   aCTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATgcgc                                                             <  1:4232270/101‑1 (MQ=255)

                                    cTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGcc                                                            >  2:4232789/1‑101 (MQ=255)

                                     tGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGccc                                                           >  1:4232790/1‑101 (MQ=255)

                                        aCATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCgg                                                        >  1:4232796/1‑101 (MQ=255)

                                        aCATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCgg                                                        >  2:4232795/1‑101 (MQ=255)

                                          aTGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCa                                                      <  2:4232198/101‑1 (MQ=255)

                                          aTGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCa                                                      >  1:4232800/1‑101 (MQ=255)

                                            ggAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCaga                                                    >  2:4232803/1‑101 (MQ=255)

                                               aCGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGaa                                                 <  1:4232045/101‑1 (MQ=255)

                                                cGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGaaa                                                >  2:4232807/1‑101 (MQ=255)

                                                 gATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCag                                                     <  1:4232324/95‑1 (MQ=255)

                                                 gATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAAt                                               >  1:4232811/1‑101 (MQ=255)

                                                  aTGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATc                                              >  2:4232815/1‑101 (MQ=255)

                                                   tgtgtgGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCa                                             <  1:4232313/101‑1 (MQ=255)

                                                       tgGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGAc                                         >  1:4232817/1‑101 (MQ=255)

                                                             gCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCAt                                   <  1:4232361/101‑1 (MQ=255)

                                                             gCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCAt                                   >  2:4232822/1‑101 (MQ=255)

                                                             gCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGACCCGGCAGAGAAATCATGACCGCCAt                                   >  1:4232821/1‑101 (MQ=255)

                                                              cATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATc                                  <  2:4232094/101‑1 (MQ=255)

                                                               aTTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCa                                 >  1:4232825/1‑101 (MQ=255)

                                                               aTTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCa                                 <  2:4232267/101‑1 (MQ=255)

                                                                tttttACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCaa                                <  1:4232689/101‑1 (MQ=255)

                                                                  tttACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGc                                      <  1:4232217/93‑1 (MQ=255)

                                                                   ttACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAgcg                             >  2:4232828/1‑101 (MQ=255)

                                                                    tACTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAgcgc                            <  2:4232291/101‑1 (MQ=255)

                                                                     aCTGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGcgcg                           <  2:4232494/101‑1 (MQ=255)

                                                                       tGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGaa                         >  1:4232835/1‑101 (MQ=255)

                                                                       tGGGCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGaa                         >  2:4232836/1‑101 (MQ=255)

                                                                        gggCGGCATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAAc                        >  2:4232839/1‑101 (MQ=255)

                                                                         ggcggcATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACa                       >  2:4232843/1‑101 (MQ=255)

                                                                           cggcATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAAt                     <  2:4232395/101‑1 (MQ=255)

                                                                            ggcATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAgcg                             <  1:4232445/92‑1 (MQ=255)

                                                                             gcATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAAttt                   <  2:4232576/101‑1 (MQ=255)

                                                                             gcATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAAttt                   >  1:4232846/1‑101 (MQ=255)

                                                                             gcATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAAttt                   <  2:4232338/101‑1 (MQ=255)

                                                                             gcATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAAttt                   <  2:4232404/101‑1 (MQ=255)

                                                                              cATGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAAtttt                  <  2:4232041/101‑1 (MQ=255)

                                                                               aTGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTa                 <  1:4232416/101‑1 (MQ=255)

                                                                                tGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTAc                >  1:4232850/1‑101 (MQ=255)

                                                                                tGAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGAATGCGCCCCGGCa                                                      >  2:4232852/1‑63 (MQ=255)

                                                                                 gAGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACg               <  2:4232436/101‑1 (MQ=255)

                                                                                  aGTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGt              <  1:4232080/101‑1 (MQ=255)

                                                                                   gTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTg             <  2:4232474/101‑1 (MQ=255)

                                                                                   gTGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTg             >  2:4232858/1‑101 (MQ=255)

                                                                                    tGGAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTgg            <  1:4232352/101‑1 (MQ=255)

                                                                                     ggAATTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGc           >  1:4232860/1‑101 (MQ=255)

                                                                                     ggAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGc           <  1:4232297/101‑1 (MQ=255)

                                                                                     ggAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGc           <  2:4232367/101‑1 (MQ=255)

                                                                                      gAACTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGCg          >  2:4232862/1‑101 (MQ=255)

                                                                                       aaCTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGCGt         <  2:4232378/101‑1 (MQ=255)

                                                                                        aCTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGCGTa        <  2:4232332/101‑1 (MQ=255)

                                                                                         cTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGCGTAt       <  2:4232259/101‑1 (MQ=255)

                                                                                         cTGTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGCGTAt       >  2:4232866/1‑101 (MQ=255)

                                                                                           gTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGCGTATCa     >  1:4232870/1‑101 (MQ=255)

                                                                                           gTCGTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGCGTATCa     <  2:4232522/101‑1 (MQ=255)

                                                                                              gTTGACCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGCGTATCAAcc  >  1:4232874/1‑101 (MQ=255)

GGGCATTCCTCTGACCGTTCGGCTTTATGCAGATGACTGGACATGGAACGATGTGTGGTACGCATTTTTACTGGGCGGCATGAGTGGAACTGTCGTTGGCCTGCTCTGCTATTACCTGATGAGCGTACGTATGCGCCCCGGCAGAGAAATCATGACCGCCATCAAGCGCGAACAATTTTACGTGGCGTATCAACC  >  NC_000913/2271379‑2271573

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