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	100.0%	1241.0 / NA	393	Q378R (CAG→CGG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (179/214); total (179/214)

Only 100 of 407 total aligned reads displayed.

GCGCAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCG  >  NC_000913/3989808‑3989987

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

  gcAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGTTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTg                                                                                       >  1:4780201/1‑93 (MQ=255)

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

   cAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGc                                                                                      >  2:3053822/1‑93 (MQ=255)

     gtgtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGTTCGTCTGGCAATGCTCGATAACCGCGCTNANTGGAAGATTGATCGGGTGCGt                                                                                    >  1:4038772/1‑93 (MQ=255)

     gtgtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGt                                                                                    >  2:7562214/1‑93 (MQ=255)

      tgtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTg                                                                                   >  1:5460707/1‑93 (MQ=255)

      tgtTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTg                                                                                   >  2:5227998/1‑93 (MQ=255)

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

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

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

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

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

                gTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACaa                                                                         <  2:1439412/93‑1 (MQ=255)

                gTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGGGCGTGAGGCGcacac                                                                         >  1:1881588/1‑92 (MQ=255)

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

                 ttAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCAGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAAc                                                                        >  2:4955020/1‑93 (MQ=255)

                    gTGAGCGAGGGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAg                                                                     <  2:5895015/93‑1 (MQ=255)

                      gAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTCCGTGAGGCGCACAACGAGtt                                                                   <  1:1424409/93‑1 (MQ=255)

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

                            gTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGa                                                             >  2:2829363/1‑93 (MQ=255)

                            gTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGa                                                             >  1:6112419/1‑93 (MQ=255)

                             tGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAc                                                            >  2:4942740/1‑93 (MQ=255)

                             tGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAc                                                            <  1:7828449/93‑1 (MQ=255)

                             tGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAc                                                            >  2:6617519/1‑93 (MQ=255)

                             tGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAc                                                            <  2:7175328/93‑1 (MQ=255)

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

                               gggTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGAcgc                                                          >  1:3889271/1‑93 (MQ=255)

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

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

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

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

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

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

                                    gggACGAAGTTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGCTTGATCGGGTGCGTGAGGCGCACAACGAGGTGCTCGACGCgaggg                                                     >  2:563271/1‑90 (MQ=255)

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

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

                                     ggACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCgatgat                                                    >  1:4529427/1‑93 (MQ=255)

                                       aCGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGc                                                  <  1:1057869/93‑1 (MQ=255)

                                        cGAAGTTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCa                                                 <  1:791104/93‑1 (MQ=255)

                                        cGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCa                                                 <  2:6917114/93‑1 (MQ=255)

                                        cGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCa                                                 <  1:6364358/93‑1 (MQ=255)

                                        cGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCa                                                 <  2:3864916/93‑1 (MQ=255)

                                          aaGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCGCAACGAGTTGCTCGACGCGATGATGCaga                                               <  1:94808/93‑1 (MQ=255)

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

                                            gCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGc                                             <  1:3568037/93‑1 (MQ=255)

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

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

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

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

                                                 tCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGa                                                                                  >  1:3067767/1‑51 (MQ=255)

                                                     gCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAAt                                    >  2:6084696/1‑93 (MQ=255)

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

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

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

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

                                                         tGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCCCAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGa                                >  1:2521416/1‑93 (MQ=255)

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

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

                                                         tGCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGa                                <  2:6155368/93‑1 (MQ=255)

                                                          gCTCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGAt                               <  2:3528971/93‑1 (MQ=255)

                                                           gtCGATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATc                              <  1:6668234/92‑1 (MQ=255)

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

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

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

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

                                                              gATAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTCCCGTAATCTGATCCGc                           >  1:6022131/1‑93 (MQ=255)

                                                                tAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCtt                         >  2:5034086/1‑93 (MQ=255)

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

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

                                                                tAACCGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCGCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCtt                         <  2:5488883/93‑1 (MQ=255)

                                                                   ccGCGCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgc                      <  2:480376/93‑1 (MQ=255)

                                                                    cgcgCTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCACAGCTACCGTAATCTGATCCGCTTTgcg                     <  2:1763014/93‑1 (MQ=255)

                                                                     gcgcTAACTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTgcgc                    <  2:7361755/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                          aaCTGGAAGATTGATCGGGTGCGGGAGCCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGc               <  1:582397/93‑1 (MQ=255)

                                                                            cTGGAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCaa             >  1:3205194/1‑93 (MQ=255)

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

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

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

                                                                               gAAGATTGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATaa          <  2:3605540/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                                   aTTGATCGGATGCGTGAGGCGCACAACGAGATGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGATTTGCGCGTCGCAATAACCtt      >  1:5926488/1‑93 (MQ=255)

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

                                                                                     tGATCGGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAg    >  2:4287837/1‑93 (MQ=255)

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

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

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

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

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

GCGCAGTGTTGAGCCAGTTAGTGAGCGAGTGGGGTTGGGACGAAGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCG  >  NC_000913/3989808‑3989987

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