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	99.6%	734.9 / ‑5.2	235	Q378R (CAG→CGG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (108/126); total (109/126)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

Only 100 of 252 total aligned reads displayed.

GCGCAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCG  >  NC_000913/3989808‑3989987

gcgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTTGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAgg                                                                                         >  2:665924/1‑93 (MQ=255)

gcgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCggg                                                                                         >  1:1018204/1‑93 (MQ=255)

gcgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCggg                                                                                         >  1:2746504/1‑93 (MQ=255)

gcgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCggg                                                                                         >  2:2485449/1‑93 (MQ=255)

 cgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCgggg                                                                                        >  2:2965723/1‑92 (MQ=255)

 cgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGt                                                                                        <  2:288193/93‑1 (MQ=255)

 cgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGt                                                                                        <  2:3336862/93‑1 (MQ=255)

 cgcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGt                                                                                        <  2:1627126/93‑1 (MQ=255)

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

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

   cAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGc                                                                                      <  2:53777/93‑1 (MQ=255)

   cAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGc                                                                                      <  2:326536/93‑1 (MQ=255)

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

     gtgtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATGACCGCGCTAACTGGAAGATTGATCGGGGGCgg                                                                                    >  2:1212349/1‑92 (MQ=255)

     gtgtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCgg                                                                                    >  1:4347107/1‑92 (MQ=255)

      tgtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTg                                                                                   <  1:4048471/93‑1 (MQ=255)

       gtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGa                                                                                  <  1:3458950/93‑1 (MQ=255)

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

          gAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGgc                                                                               <  1:1777896/93‑1 (MQ=255)

          gAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGgc                                                                               <  2:3553699/93‑1 (MQ=255)

              cAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGcac                                                                           <  2:3880863/93‑1 (MQ=255)

               aGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGcaca                                                                          >  2:2151071/1‑93 (MQ=255)

                gTTAGTGAGCGAGTGGGTTTGGGACGACGCTCGTCTGGCAATGCTCGATAACCGCGCTAATTGTAAGAGTGATCGGGTGCGTGAGGCGCACaa                                                                         <  2:1396897/93‑1 (MQ=255)

                 ttAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAAc                                                                        >  2:2544103/1‑93 (MQ=255)

                   aGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGa                                                                      <  2:1091535/93‑1 (MQ=255)

                     tGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCCCAACGAGt                                                                    <  2:4549753/93‑1 (MQ=255)

                     tGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGt                                                                    <  2:2106734/93‑1 (MQ=255)

                       aGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCCCAACGAGTTg                                                                  >  2:2943753/1‑93 (MQ=255)

                       aGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTg                                                                  >  2:603649/1‑93 (MQ=255)

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

                           aGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCg                                                              >  2:2674818/1‑93 (MQ=255)

                           aGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCg                                                              <  1:451353/93‑1 (MQ=255)

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

                              ggggTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcg                                                           <  2:3646333/93‑1 (MQ=255)

                              ggggTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcg                                                           >  2:2079409/1‑93 (MQ=255)

                               gggTTTGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcgc                                                          <  1:2480794/93‑1 (MQ=255)

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

                                ggTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcgcg                                                         >  1:167032/1‑93 (MQ=255)

                                ggTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcgcg                                                         <  1:3269550/93‑1 (MQ=255)

                                 gTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCga                                                        <  2:3996665/93‑1 (MQ=255)

                                 gTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCga                                                        >  1:4489279/1‑93 (MQ=255)

                                 gTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCga                                                        <  2:3396321/93‑1 (MQ=255)

                                  ttGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgat                                                       >  2:85703/1‑93 (MQ=255)

                                  ttGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgat                                                       <  1:1199094/93‑1 (MQ=255)

                                  ttGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgat                                                       >  1:2505788/1‑93 (MQ=255)

                                  ttGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgat                                                       >  1:649087/1‑93 (MQ=255)

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

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

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

                                      gACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGatgatg                                                   <  2:4224035/93‑1 (MQ=255)

                                      gACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGatgatg                                                   >  1:3302955/1‑93 (MQ=255)

                                      gACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGatgatg                                                   >  2:1714789/1‑93 (MQ=255)

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

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

                                           aGCTCGTCTGGCAATGCGCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCAATGATGCagag                                              <  1:808177/93‑1 (MQ=255)

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

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

                                             cTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCt                                            >  1:362060/1‑93 (MQ=255)

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

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

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

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

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

                                                         tGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATATGa                                >  1:194741/1‑93 (MQ=255)

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

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

                                                            tCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATcc                             <  2:706818/93‑1 (MQ=255)

                                                             cGATAACCTCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCg                            >  2:4577767/1‑93 (MQ=255)

                                                             cGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCg                            >  1:853736/1‑93 (MQ=255)

                                                               aTAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCt                          >  1:4434472/1‑93 (MQ=255)

                                                                 aaCCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCttt                        <  2:911936/93‑1 (MQ=255)

                                                                  aCCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTg                       <  1:143052/93‑1 (MQ=255)

                                                                    cgcgCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGGAGAGCTACCGTAATCTGATCCGCTTTgcg                     >  2:737213/1‑93 (MQ=255)

                                                                     ccgcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgcgc                    <  1:2693639/92‑1 (MQ=255)

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

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

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

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

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

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

                                                                            cTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCaa             <  1:319037/93‑1 (MQ=255)

                                                                             tGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAAt            <  1:2104878/93‑1 (MQ=255)

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

                                                                              ggAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATa           <  2:4684187/93‑1 (MQ=255)

                                                                               gAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATaa          >  2:3762958/1‑93 (MQ=255)

                                                                               gAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATaa          >  2:817212/1‑93 (MQ=255)

                                                                                 aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        >  1:4260360/1‑93 (MQ=255)

                                                                                 aGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAAcc        >  2:1600953/1‑93 (MQ=255)

                                                                                  gATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCt       <  2:2434890/93‑1 (MQ=255)

                                                                                  gATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCt       <  1:149580/93‑1 (MQ=255)

                                                                                  gATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCt       <  1:1380404/93‑1 (MQ=255)

                                                                                  gATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCt       <  1:1282289/93‑1 (MQ=255)

                                                                                     tGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAg    <  1:2357472/93‑1 (MQ=255)

                                                                                     tGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAg    <  2:988424/93‑1 (MQ=255)

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

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

                                                                                      gATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGc   >  2:1960492/1‑93 (MQ=255)

                                                                                      gATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGc   >  2:1289271/1‑93 (MQ=255)

                                                                                      gATCGGGTGCGGGAGGCGCACAAAGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGc   <  1:1878530/93‑1 (MQ=255)

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

GCGCAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCG  >  NC_000913/3989808‑3989987

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