BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,870,765	+G	coding (476/1764 nt)	dgoA ←	2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,870,764	1	.	G	100.0%	863.9 / NA	272	A159A (GCA→GCC)	dgoA	2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (120/152); total (120/152)

Only 100 of 285 total aligned reads displayed.

GCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCT‑GCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATG  >  NC_000913/3870672‑3870851

 cgGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCT‑g                                                                                        <  1:969932/93‑1 (MQ=255)

  ggTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCT‑gg                                                                                       >  2:2478206/1‑92 (MQ=255)

        ctctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAc                                                                                  <  1:567537/93‑1 (MQ=255)

         tctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcg                                                                                 >  2:3354234/1‑93 (MQ=255)

          ttACGGATGGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgc                                                                                <  2:1691782/92‑1 (MQ=255)

           tACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgcc                                                                               <  2:2261200/93‑1 (MQ=255)

           tACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgcc                                                                               >  1:2553288/1‑93 (MQ=255)

              ggATTGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCTCCAGGTTTTCTGGCGTCACgccgcc                                                                            >  1:546183/1‑93 (MQ=255)

              ggATTGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACgccgcc                                                                            >  2:1198576/1‑93 (MQ=255)

                 ttGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAAc                                                                         >  1:2786094/1‑93 (MQ=255)

                  tGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACg                                                                        <  2:2881563/93‑1 (MQ=255)

                    cccGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                                      >  1:1409098/1‑93 (MQ=255)

                     ccGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCa                                                                     <  1:1561219/93‑1 (MQ=255)

                       ggCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCaaa                                                                   <  1:4489870/93‑1 (MQ=255)

                        gcgcGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACACCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCTCCAACGGCAAAg                                                                  <  1:4088057/93‑1 (MQ=255)

                        gcgcGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAg                                                                  <  1:1720016/93‑1 (MQ=255)

                        gcgcGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAg                                                                  >  1:736102/1‑93 (MQ=255)

                        gcgcGATAGAGAGCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAg                                                                  <  2:3374530/93‑1 (MQ=255)

                            gATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTg                                                              <  1:837659/93‑1 (MQ=255)

                            gATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTg                                                              >  2:4367010/1‑93 (MQ=255)

                             aTAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                                             <  2:1330187/93‑1 (MQ=255)

                             aTAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                                             <  1:1578752/93‑1 (MQ=255)

                             aTAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                                             <  2:2378309/93‑1 (MQ=255)

                             aTAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                                             >  1:768866/1‑93 (MQ=255)

                              tAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCg                                                            >  1:1358990/1‑93 (MQ=255)

                               agagATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGa                                                           <  1:302169/93‑1 (MQ=255)

                                 agaTCTCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATg                                                         <  1:510136/93‑1 (MQ=255)

                                 agaTCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATg                                                         >  2:4540928/1‑93 (MQ=255)

                                 agaTCGCT‑CCGAAGCCCCGCCCCTGCACAACCGGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATg                                                         <  2:1620918/93‑1 (MQ=255)

                                  gaTCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                                        <  1:4200284/93‑1 (MQ=255)

                                  gaTCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                                        <  1:1176731/93‑1 (MQ=255)

                                   aTCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTc                                                       >  1:2102516/1‑93 (MQ=255)

                                    tCGCTGCCTAAG‑CGCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCCACGGCAAAGACTGCGATGttc                                                      >  2:4061900/1‑91 (MQ=255)

                                    tCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                      >  1:4403888/1‑93 (MQ=255)

                                    tCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                      <  2:4428654/93‑1 (MQ=255)

                                     cGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCg                                                     <  2:3711289/93‑1 (MQ=255)

                                      gCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTTCgt                                                    >  2:540023/1‑92 (MQ=255)

                                      gCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGa                                                    <  2:3058681/93‑1 (MQ=255)

                                      gCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGa                                                    >  1:610111/1‑93 (MQ=255)

                                      gCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGa                                                    <  2:2779309/93‑1 (MQ=255)

                                        tGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                                  <  1:1081926/93‑1 (MQ=255)

                                        cgCCTAGG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCCCGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                                  <  2:477682/92‑1 (MQ=255)

                                            tAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCaa                                              >  1:3380382/1‑93 (MQ=255)

                                            tAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCaa                                              >  1:2859890/1‑93 (MQ=255)

                                               g‑cccgcccCTGCACACCCTGCGTCTATCCGCTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           <  1:1305308/93‑1 (MQ=255)

                                               g‑cccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           <  2:3859499/93‑1 (MQ=255)

                                               g‑cccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           >  2:1035469/1‑93 (MQ=255)

                                                 cccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAcc                                          <  1:2068463/93‑1 (MQ=255)

                                                  ccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCg                                         >  2:2981634/1‑93 (MQ=255)

                                                    gcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCt                                       >  2:474809/1‑93 (MQ=255)

                                                    cccccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCt                                       <  1:2721971/92‑1 (MQ=255)

                                                      cccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCGAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCttt                                     <  2:3770397/93‑1 (MQ=255)

                                                        cTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                                   <  1:2481703/93‑1 (MQ=255)

                                                         tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  <  2:3484144/93‑1 (MQ=255)

                                                         tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  >  2:62441/1‑93 (MQ=255)

                                                         tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  >  2:4187302/1‑93 (MQ=255)

                                                          gCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAAc                                 <  1:2197091/93‑1 (MQ=255)

                                                          gCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAAc                                 <  1:3831035/93‑1 (MQ=255)

                                                          gCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAAc                                 >  1:1641182/1‑93 (MQ=255)

                                                            acaACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGc                               >  1:3701716/1‑93 (MQ=255)

                                                             caACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                              >  1:3952801/1‑93 (MQ=255)

                                                             caACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                              <  1:4326417/93‑1 (MQ=255)

                                                              aaCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCtt                             <  2:4518290/93‑1 (MQ=255)

                                                              aaCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCtt                             >  1:1160943/1‑93 (MQ=255)

                                                              aaCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCtt                             <  1:235447/93‑1 (MQ=255)

                                                               aCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt                            <  2:828582/93‑1 (MQ=255)

                                                               aCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt                            >  1:2893324/1‑93 (MQ=255)

                                                               aCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt                            <  2:3956003/93‑1 (MQ=255)

                                                               aCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt                            <  2:3149764/93‑1 (MQ=255)

                                                               aCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt                            >  2:2746539/1‑93 (MQ=255)

                                                                  tGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGAt                         <  2:701355/93‑1 (MQ=255)

                                                                  tGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGAt                         <  1:4032571/93‑1 (MQ=255)

                                                                   gCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATg                        <  1:3126386/93‑1 (MQ=255)

                                                                   gCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATg                        <  1:3117113/93‑1 (MQ=255)

                                                                     gTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTa                      <  1:3593783/93‑1 (MQ=255)

                                                                      tCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAt                     >  1:1907526/1‑93 (MQ=255)

                                                                      tCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAt                     <  1:4532833/93‑1 (MQ=255)

                                                                       cTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt                    <  1:2330203/93‑1 (MQ=255)

                                                                       cTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt                    <  1:635241/93‑1 (MQ=255)

                                                                       cTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt                    <  2:15385/93‑1 (MQ=255)

                                                                       cTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt                    <  1:4457983/93‑1 (MQ=255)

                                                                        tATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTg                   >  2:3182015/1‑93 (MQ=255)

                                                                        tATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTg                   <  2:3542033/93‑1 (MQ=255)

                                                                         aTCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGc                  <  2:3238524/93‑1 (MQ=255)

                                                                           ccACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCgg                >  1:3524653/1‑93 (MQ=255)

                                                                             aCTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGAc              <  1:4174773/93‑1 (MQ=255)

                                                                             aCTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGAc              <  1:2395432/93‑1 (MQ=255)

                                                                             aCTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGAc              >  2:1514432/1‑93 (MQ=255)

                                                                              cTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGAcc             >  2:146364/1‑93 (MQ=255)

                                                                                  gcaaGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         <  1:2959108/90‑1 (MQ=255)

                                                                                  gcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         <  2:197712/93‑1 (MQ=255)

                                                                                  gcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         >  2:1111588/1‑93 (MQ=255)

                                                                                    cAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGc       <  2:330146/93‑1 (MQ=255)

                                                                                    cAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGc       >  2:2233925/1‑93 (MQ=255)

                                                                                     aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      >  1:3933216/1‑93 (MQ=255)

                                                                                     aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      <  1:4589759/93‑1 (MQ=255)

                                                                                     aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      <  2:2416416/93‑1 (MQ=255)

                                                                                       gTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGa    >  2:871492/1‑93 (MQ=255)

                                                                                        ttttCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGAt   <  2:291853/93‑1 (MQ=255)

                                                                                        ttttCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGAt   >  1:226676/1‑93 (MQ=255)

GCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAG‑CCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCT‑GCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATG  >  NC_000913/3870672‑3870851

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