BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,989,898	A→G	Q378R (CAG→CGG)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,989,898	0	A	G	97.0%	514.2 / 2.8	168	Q378R (CAG→CGG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base A (3/2); new base G (81/82); total (84/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

Only 100 of 172 total aligned reads displayed.

GCGCAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCG  >  NC_000913/3989808‑3989987

gcgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCggg                                                                                         <  1:2801801/93‑1 (MQ=255)

 cgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGCAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGt                                                                                        <  1:2314667/93‑1 (MQ=255)

  gcAGTGTTTAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTCACTGGAAGATTGCTCGGGTg                                                                                       >  2:3291469/1‑93 (MQ=255)

  gcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTg                                                                                       >  2:2434429/1‑93 (MQ=255)

    aGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCg                                                                                     >  1:48664/1‑93 (MQ=255)

    aGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCg                                                                                     >  1:2571920/1‑93 (MQ=255)

       gtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTgg                                                                                  >  2:247434/1‑92 (MQ=255)

       gtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGa                                                                                  >  1:4696010/1‑93 (MQ=255)

        ttGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAg                                                                                 <  2:1317306/93‑1 (MQ=255)

         tGAGCAAGTTAGTGAGCGAGTGGGGTTGGGCCGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAgg                                                                                <  1:4498888/93‑1 (MQ=255)

          gAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGgc                                                                               <  2:2323894/93‑1 (MQ=255)

            gCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGgcgc                                                                             >  1:3397636/1‑93 (MQ=255)

              cAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGGGGCGaac                                                                           >  1:3982755/1‑93 (MQ=255)

              cAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGcac                                                                           >  1:3982756/1‑93 (MQ=255)

               aGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGcaca                                                                          <  2:185428/93‑1 (MQ=255)

                  tAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACg                                                                       <  2:1441999/93‑1 (MQ=255)

                    gTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAg                                                                     <  1:4397325/93‑1 (MQ=255)

                    gTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAg                                                                     <  1:4187002/93‑1 (MQ=255)

                     tGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGt                                                                    >  1:2401770/1‑93 (MQ=255)

                        gCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGc                                                                 <  2:3293372/93‑1 (MQ=255)

                        gCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGc                                                                 <  2:635372/93‑1 (MQ=255)

                        gCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGc                                                                 >  1:3055602/1‑93 (MQ=255)

                         ggAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCt                                                                <  1:534300/92‑1 (MQ=255)

                         cGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCt                                                                <  2:4193681/93‑1 (MQ=255)

                         cGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCt                                                                <  2:1217713/93‑1 (MQ=255)

                         cGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCt                                                                >  2:1624503/1‑93 (MQ=255)

                          gAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTc                                                               >  1:2704367/1‑93 (MQ=255)

                            gTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGa                                                             <  2:2893247/93‑1 (MQ=255)

                            gTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGa                                                             <  2:3274513/93‑1 (MQ=255)

                               gggTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcgc                                                          <  1:1177547/93‑1 (MQ=255)

                                    gggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatga                                                     >  1:4134214/1‑93 (MQ=255)

                                    gggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatga                                                     <  1:1958936/93‑1 (MQ=255)

                                    gggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatga                                                     >  2:2423172/1‑93 (MQ=255)

                                     ggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatgat                                                    <  1:1700037/93‑1 (MQ=255)

                                     ggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatgat                                                    <  2:2997425/93‑1 (MQ=255)

                                     ggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatgat                                                    <  1:1754556/93‑1 (MQ=255)

                                     ggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatgat                                                    <  1:1940261/93‑1 (MQ=255)

                                       aCGAAGCTTGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGc                                                  >  2:4346229/1‑93 (MQ=255)

                                        cGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCa                                                 >  2:4684973/1‑93 (MQ=255)

                                         gAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCag                                                >  1:3924918/1‑93 (MQ=255)

                                         gAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCag                                                <  2:3410428/93‑1 (MQ=255)

                                         gAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGAGGATGCag                                                >  2:1280077/1‑93 (MQ=255)

                                          aaGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCaga                                               <  1:136748/93‑1 (MQ=255)

                                          aaGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCaga                                               >  2:2633956/1‑93 (MQ=255)

                                           aGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCagag                                              >  2:1005348/1‑93 (MQ=255)

                                           aGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGAGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCagag                                              >  1:4029253/1‑93 (MQ=255)

                                            gCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGc                                             >  1:4457815/1‑93 (MQ=255)

                                             cTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCt                                            <  2:2404933/93‑1 (MQ=255)

                                              tCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTa                                           >  1:4551652/1‑93 (MQ=255)

                                              tCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTa                                           <  1:1413614/93‑1 (MQ=255)

                                               cGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTAc                                          <  1:2741492/93‑1 (MQ=255)

                                               cGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTAc                                          >  1:1655106/1‑93 (MQ=255)

                                                 tCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCg                                        >  2:4615458/1‑93 (MQ=255)

                                                 tCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCg                                        >  1:4056915/1‑93 (MQ=255)

                                                   tGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTa                                      >  1:3880142/1‑93 (MQ=255)

                                                   tGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTa                                      >  1:3827678/1‑93 (MQ=255)

                                                   tGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTa                                      <  2:2948951/93‑1 (MQ=255)

                                                   tGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTa                                      <  1:2630688/93‑1 (MQ=255)

                                                   tGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTa                                      >  2:4024062/1‑93 (MQ=255)

                                                    ggCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTaa                                     >  1:2250222/1‑93 (MQ=255)

                                                      cAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATc                                   <  1:4289105/93‑1 (MQ=255)

                                                       aaTGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCt                                  >  2:3333601/1‑93 (MQ=255)

                                                       aaTGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCt                                  <  1:65022/93‑1 (MQ=255)

                                                        aTGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTg                                 >  2:3074954/1‑93 (MQ=255)

                                                        aTGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTg                                 >  2:674006/1‑93 (MQ=255)

                                                         tGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGa                                >  1:2137697/1‑93 (MQ=255)

                                                          gCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGAt                               >  2:861555/1‑93 (MQ=255)

                                                           cTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATc                              <  1:981462/93‑1 (MQ=255)

                                                           cTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATc                              >  2:690404/1‑93 (MQ=255)

                                                           cTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATc                              <  1:3710574/93‑1 (MQ=255)

                                                            tCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATcc                             >  1:2042482/1‑93 (MQ=255)

                                                            tCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATcc                             >  1:1080149/1‑93 (MQ=255)

                                                              gATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGc                           >  2:199671/1‑93 (MQ=255)

                                                              gATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGc                           <  2:1696377/93‑1 (MQ=255)

                                                               aTAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCt                          <  2:1248358/93‑1 (MQ=255)

                                                               aTAACAGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCt                          >  2:487601/1‑93 (MQ=255)

                                                                tAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCtt                         <  2:4739490/93‑1 (MQ=255)

                                                                 aaCCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCAGATCCGCttt                        >  1:2882307/1‑93 (MQ=255)

                                                                 aaCCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGACCTACCGTAATCTGATCCGCttt                        <  1:831003/93‑1 (MQ=255)

                                                                  aCCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTg                       >  2:3666571/1‑93 (MQ=255)

                                                                  aCCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTg                       >  1:3515112/1‑93 (MQ=255)

                                                                  aCCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGCTGATGCAGAGCTACCGTAATCTGATCCGCtgtg                       >  1:3800768/1‑93 (MQ=255)

                                                                     gcgcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgcgc                    >  2:4718256/1‑93 (MQ=255)

                                                                       gcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGt                  <  2:47241/93‑1 (MQ=255)

                                                                       gcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGt                  >  1:3096678/1‑93 (MQ=255)

                                                                       gcTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGt                  >  2:4020303/1‑93 (MQ=255)

                                                                        cTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTc                 <  1:2020123/93‑1 (MQ=255)

                                                                        cTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTc                 >  2:2792898/1‑93 (MQ=255)

                                                                          aaCTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGc               >  2:1971374/1‑93 (MQ=255)

                                                                          aaCTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGc               >  2:1971375/1‑93 (MQ=255)

                                                                           aCTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCa              >  2:3689209/1‑93 (MQ=255)

                                                                             tGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAAt            >  1:2316124/1‑93 (MQ=255)

                                                                                aaGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAc         <  1:4679747/93‑1 (MQ=255)

                                                                                 aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        <  2:2079933/93‑1 (MQ=255)

                                                                                 aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        <  2:2254606/93‑1 (MQ=255)

                                                                                   aTTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCtt      >  1:4395124/1‑93 (MQ=255)

                                                                                      gATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCTCTTTGCGCGTCGCAATAACCTTAGc   >  2:2564596/1‑93 (MQ=255)

                                                                                      gATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGc   <  2:820439/93‑1 (MQ=255)

                                                                                       aTCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCg  >  2:3803473/1‑93 (MQ=255)

                                                                                       aTCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCg  >  1:1048911/1‑93 (MQ=255)

GCGCAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCG  >  NC_000913/3989808‑3989987

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