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%	943.1 / NA	290	A159A (GCA→GCC)	dgoA	2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (138/152); total (138/152)

Only 100 of 317 total aligned reads displayed.

GCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCT‑G‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATG  >  NC_000913/3870672‑3870851

   gTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCT‑ggc                                                                                       <  2:2020058/93‑3 (MQ=255)

   gTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCT‑ggc                                                                                       >  2:5410323/1‑91 (MQ=255)

    tGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCT‑ggcg                                                                                      >  1:3631919/1‑90 (MQ=255)

       gcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCa                                                                                   >  1:1059739/1‑93 (MQ=255)

       gcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCa                                                                                   <  1:1862761/93‑1 (MQ=255)

        ctctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCAc                                                                                  >  1:1505718/1‑93 (MQ=255)

          ctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCAcgc                                                                                <  2:2106957/93‑1 (MQ=255)

          ctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCAcgc                                                                                <  1:5118028/93‑1 (MQ=255)

           tACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCAcgcc                                                                               >  1:5368600/1‑93 (MQ=255)

           tACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCAcgcc                                                                               >  2:3910623/1‑93 (MQ=255)

             cGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCAcgccgc                                                                             <  2:3469116/93‑1 (MQ=255)

                 ttGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAAc                                                                         <  1:4220035/93‑1 (MQ=255)

                     ccGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCa                                                                     <  1:2540303/93‑1 (MQ=255)

                       ggCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCaaa                                                                   >  2:3812654/1‑93 (MQ=255)

                        gcgcGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAg                                                                  <  1:6077300/93‑1 (MQ=255)

                         cgcgATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGa                                                                 <  2:1604593/93‑1 (MQ=255)

                         cgcgATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCGGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGa                                                                 <  1:4020092/93‑1 (MQ=255)

                          gcgATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGAc                                                                <  1:6176384/93‑1 (MQ=255)

                           cgATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACt                                                               >  1:915746/1‑93 (MQ=255)

                           cgATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACt                                                               >  2:4615723/1‑93 (MQ=255)

                             aTAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGc                                                             >  2:2470332/1‑93 (MQ=255)

                             aTAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGc                                                             >  2:1663484/1‑93 (MQ=255)

                              tAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCg                                                            >  2:3748147/1‑93 (MQ=255)

                               agagATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGa                                                           <  1:257533/93‑1 (MQ=255)

                               agagATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGa                                                           <  2:5958137/93‑1 (MQ=255)

                                 agaTCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATg                                                         <  1:5537266/93‑1 (MQ=255)

                                 agaTCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATg                                                         >  2:1984923/1‑93 (MQ=255)

                                  gaTCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGt                                                        <  1:1299589/93‑1 (MQ=255)

                                   aTCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTc                                                       >  2:5138770/1‑93 (MQ=255)

                                    tCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                      >  2:4828415/1‑93 (MQ=255)

                                    tCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                      >  2:2271483/1‑93 (MQ=255)

                                    tCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                      <  2:5178442/93‑1 (MQ=255)

                                    tCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                      >  1:3157161/1‑93 (MQ=255)

                                      gCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGa                                                    >  1:2238636/1‑93 (MQ=255)

                                      gCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGa                                                    <  2:467257/93‑1 (MQ=255)

                                      cactccTAAGGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCAAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCg                                                     <  2:4504851/89‑1 (MQ=255)

                                       cTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGAt                                                   <  1:4619975/93‑1 (MQ=255)

                                       cTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGAt                                                   <  2:6224877/93‑1 (MQ=255)

                                        tGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                                  <  2:1880436/93‑1 (MQ=255)

                                        tGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                                  <  1:3826846/93‑1 (MQ=255)

                                         gCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATgg                                                 <  1:250598/93‑1 (MQ=255)

                                         accTAA‑GCCCCACCCTGCAAAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATgg                                                 <  2:887544/92‑1 (MQ=255)

                                           cTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCa                                               >  1:1956537/1‑93 (MQ=255)

                                           cTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCa                                               >  2:6086585/1‑93 (MQ=255)

                                           cTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCa                                               >  1:20499/1‑93 (MQ=255)

                                           cTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCa                                               >  1:3549334/1‑93 (MQ=255)

                                            tAA‑ACCAGCCCCTGCACATCCTGCGTCTATCCACTGTGCCAGGTTTTCTGG‑CGTCACGCTGCCAACGGCAAAGACTTCGATGTCCGATGGCaa                                              <  1:5869935/93‑1 (MQ=255)

                                                 cccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAcc                                          <  1:3767153/93‑1 (MQ=255)

                                                  ccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCg                                         >  2:1164491/1‑93 (MQ=255)

                                                   cgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGc                                        <  2:2390576/93‑1 (MQ=255)

                                                     ccccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCtt                                      >  2:4215964/1‑93 (MQ=255)

                                                      cctggCACACCCTGCTTCTATCCGCTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCttt                                     <  2:4355433/89‑1 (MQ=255)

                                                      cccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCttt                                     <  2:1359984/93‑1 (MQ=255)

                                                       ccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCtttt                                    >  2:3603495/1‑93 (MQ=255)

                                                         tGCACAACCTGCGTCTATCCACTGCGTCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  <  2:5089414/93‑1 (MQ=255)

                                                         tGCACAACCTGCGTCTATCCACTGCGCGAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  <  1:5709092/93‑1 (MQ=255)

                                                         tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  <  2:4925213/93‑1 (MQ=255)

                                                         tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  <  2:5661544/93‑1 (MQ=255)

                                                         tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  <  2:2608380/93‑1 (MQ=255)

                                                          gCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAAc                                 >  2:1460347/1‑93 (MQ=255)

                                                           cacaACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACg                                >  2:3873182/1‑93 (MQ=255)

                                                             caACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                              <  2:1007561/93‑1 (MQ=255)

                                                             caACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                              >  2:2115350/1‑93 (MQ=255)

                                                             caACCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                              >  1:5905003/1‑93 (MQ=255)

                                                              aaCCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCtt                             >  2:3146869/1‑93 (MQ=255)

                                                              aaCCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCtt                             <  1:2759983/93‑1 (MQ=255)

                                                               aCCTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt                            <  1:1174003/93‑1 (MQ=255)

                                                                ccTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTg                           <  1:6018058/93‑1 (MQ=255)

                                                                 cTGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGa                          <  1:4453685/93‑1 (MQ=255)

                                                                  tGCTTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGAt                         <  1:5023199/93‑1 (MQ=255)

                                                                  tGCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGAt                         >  1:20730/1‑93 (MQ=255)

                                                                  cgCGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGAt                         <  1:3321308/92‑1 (MQ=255)

                                                                    cGTCTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGt                       >  2:4110652/1‑93 (MQ=255)

                                                                       cTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt                    <  1:5838064/93‑1 (MQ=255)

                                                                       cTATCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt                    <  1:1636968/93‑1 (MQ=255)

                                                                         aTCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGc                  <  2:5996477/93‑1 (MQ=255)

                                                                         aTCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGc                  >  2:721957/1‑93 (MQ=255)

                                                                         aTCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGc                  <  1:479022/93‑1 (MQ=255)

                                                                          tCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCg                 >  1:2932745/1‑93 (MQ=255)

                                                                          tCCACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCg                 >  1:1475567/1‑93 (MQ=255)

                                                                           ccACTGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCgg                <  1:799843/93‑1 (MQ=255)

                                                                               tGCGCCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCa            <  1:4507358/93‑1 (MQ=255)

                                                                                 cgcCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaa          <  2:669200/93‑1 (MQ=255)

                                                                                  gcCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         <  2:3980484/93‑1 (MQ=255)

                                                                                  gcCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         >  1:5498393/1‑93 (MQ=255)

                                                                                  gcCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         <  2:1803366/93‑1 (MQ=255)

                                                                                  gcCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         <  2:2232796/93‑1 (MQ=255)

                                                                                  gcCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         <  2:5891626/93‑1 (MQ=255)

                                                                                  gcCAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         <  1:1773375/93‑1 (MQ=255)

                                                                                    cAGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGc       <  2:5734536/93‑1 (MQ=255)

                                                                                     atgTTTTCT‑GTCGCCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      <  2:1238404/91‑1 (MQ=255)

                                                                                     aGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      <  2:5297245/93‑1 (MQ=255)

                                                                                     aGGTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      <  2:5389534/93‑1 (MQ=255)

                                                                                       gTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGa    >  1:4101669/1‑93 (MQ=255)

                                                                                       gTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGa    >  2:112928/1‑93 (MQ=255)

                                                                                       gTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGa    <  1:600130/93‑1 (MQ=255)

                                                                                       gTTTTCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGa    <  2:2793641/93‑1 (MQ=255)

                                                                                        ttttCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGAt   <  2:6185159/93‑1 (MQ=255)

                                                                                        ttttCTGG‑CGGCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGAt   <  2:1400248/93‑1 (MQ=255)

                                                                                         tttCTGG‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATg  <  1:2359796/93‑1 (MQ=255)

GCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA‑GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCT‑G‑CGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATG  >  NC_000913/3870672‑3870851

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: `‑`