BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,808,078	A→G	intergenic (+6/‑157)	yniC → / → b1728	putative phosphatase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,808,078	0	A	G	100.0%	596.1 / NA	182	intergenic (+6/‑157)	yniC/b1728	putative phosphatase/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (98/84); total (98/84)

Only 100 of 195 total aligned reads displayed.

AGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACAC  >  NC_000913/1807988‑1808166

aGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCggg                                                                                        >  2:3091604/1‑93 (MQ=255)

  cgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGgggg                                                                                      <  1:1079018/93‑1 (MQ=255)

  cgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGgggg                                                                                      <  2:2994962/93‑1 (MQ=255)

   gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGc                                                                                     >  1:999564/1‑93 (MQ=255)

    cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCa                                                                                    >  2:2932736/1‑93 (MQ=255)

     aaaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAg                                                                                   <  1:1992576/93‑1 (MQ=255)

     aaaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGCGGGCAg                                                                                   >  1:1016555/1‑93 (MQ=255)

      aaaTGCTCCACGTTTTGTATTAGCAGACGTCAAACTATCATCACTGCCAGAACTCACCGCAAAAGACCTTCTCGGATAGTGACCGGGGGGAGt                                                                                  >  2:2434887/1‑93 (MQ=255)

      aaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGt                                                                                  >  2:1018892/1‑93 (MQ=255)

      aaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGt                                                                                  <  2:4419367/93‑1 (MQ=255)

       aaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTg                                                                                 >  1:2610025/1‑93 (MQ=255)

        aTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGa                                                                                >  2:4498869/1‑93 (MQ=255)

           aTCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtct                                                                             >  2:591986/1‑93 (MQ=255)

            tCCACGTTTTGTATTAGCAGACTTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtctc                                                                            <  1:1019403/93‑1 (MQ=255)

            tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtctc                                                                            >  2:2849399/1‑93 (MQ=255)

            tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtctc                                                                            <  2:47640/93‑1 (MQ=255)

            tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtctc                                                                            >  1:1400772/1‑93 (MQ=255)

              cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGc                                                                          <  2:4537536/93‑1 (MQ=255)

              cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGc                                                                          >  2:1685467/1‑93 (MQ=255)

              cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGc                                                                          <  1:4645203/93‑1 (MQ=255)

               aCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCt                                                                         >  1:2165523/1‑93 (MQ=255)

               aCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCt                                                                         <  2:3082857/93‑1 (MQ=255)

                cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTg                                                                        >  1:4120959/1‑93 (MQ=255)

                     tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCt                                                                   <  2:912736/93‑1 (MQ=255)

                     tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCt                                                                   <  2:2206762/93‑1 (MQ=255)

                       tATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGAAGTGATCTCGCTGCCCCTgg                                                                 >  1:2441252/1‑93 (MQ=255)

                        aTTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGt                                                                >  1:4735786/1‑93 (MQ=255)

                         ttAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGtt                                                               >  2:226686/1‑93 (MQ=255)

                          tAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGCTCTCGCTGCCCCTGGttt                                                              >  2:2846872/1‑92 (MQ=255)

                            gCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                            >  2:3451591/1‑93 (MQ=255)

                             cAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCttt                                                           <  1:618415/93‑1 (MQ=255)

                             cAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCttt                                                           >  1:3838693/1‑93 (MQ=255)

                              aGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                          <  1:556007/93‑1 (MQ=255)

                              aGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                          <  1:1294224/93‑1 (MQ=255)

                               gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                         >  1:1663273/1‑93 (MQ=255)

                                  gTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTg                                                      >  1:1173168/1‑93 (MQ=255)

                                   tCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGa                                                     >  1:4418546/1‑93 (MQ=255)

                                    cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTgca                                                    >  2:3160695/1‑91 (MQ=255)

                                     aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaaa                                                   <  2:2381422/93‑1 (MQ=255)

                                     aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaaa                                                   >  2:1550972/1‑93 (MQ=255)

                                     aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaaa                                                   >  2:3880616/1‑93 (MQ=255)

                                      aaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGCAAt                                                  >  1:122365/1‑93 (MQ=255)

                                      aaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                  <  1:3223168/93‑1 (MQ=255)

                                        cTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTg                                                >  1:842495/1‑93 (MQ=255)

                                          ttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCa                                              >  2:935348/1‑93 (MQ=255)

                                          ttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCa                                              >  2:442519/1‑93 (MQ=255)

                                           tCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAt                                             >  1:1382317/1‑93 (MQ=255)

                                            cATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAtt                                            <  2:3673138/93‑1 (MQ=255)

                                             aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTc                                           >  1:1689863/1‑93 (MQ=255)

                                             aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTc                                           >  2:3769240/1‑93 (MQ=255)

                                              tCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCa                                          >  1:3451434/1‑93 (MQ=255)

                                               cGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCaa                                         <  2:2671248/93‑1 (MQ=255)

                                               cGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCaa                                         <  2:1354705/93‑1 (MQ=255)

                                                gCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                        >  1:499422/1‑93 (MQ=255)

                                                ccTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                        <  1:3555473/92‑1 (MQ=255)

                                                  tGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTg                                      <  1:2765509/93‑1 (MQ=255)

                                                   gACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGa                                     >  1:1852533/1‑93 (MQ=255)

                                                    gcAGAACTCACCGCAAAATACCTTCGCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                    <  2:1624804/92‑1 (MQ=255)

                                                    aCAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                    <  2:2367392/93‑1 (MQ=255)

                                                     cAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACg                                   >  2:2046013/1‑93 (MQ=255)

                                                     cAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACg                                   <  1:2861568/93‑1 (MQ=255)

                                                      aGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACgg                                  >  1:403873/1‑93 (MQ=255)

                                                      aGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACgg                                  >  2:2449500/1‑93 (MQ=255)

                                                      aGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACgg                                  <  2:2667909/93‑1 (MQ=255)

                                                       gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGa                                 >  1:2262957/1‑93 (MQ=255)

                                                       gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGa                                 <  1:3686747/93‑1 (MQ=255)

                                                          cTCCCCGCAAAAGACCTTCTCGGTTAATGACCGGATGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTa                              <  2:2039895/93‑1 (MQ=255)

                                                           tcgccGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTaa                             >  2:1071834/4‑93 (MQ=255)

                                                           tCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTaa                             <  1:3475405/93‑1 (MQ=255)

                                                            cACCGCAAAAGACCTTCTCGGTTAATTACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAt                            <  1:4631865/93‑1 (MQ=255)

                                                            cACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAt                            >  1:2101901/1‑93 (MQ=255)

                                                              ccGCAAAAGACCTCCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCg                          <  1:4716126/93‑1 (MQ=255)

                                                               cGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGt                         <  1:2538621/93‑1 (MQ=255)

                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                        <  1:3036312/93‑1 (MQ=255)

                                                                 cAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                       >  1:2841167/1‑93 (MQ=255)

                                                                  aaaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCaa                      <  1:4429764/93‑1 (MQ=255)

                                                                  aaaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCaa                      <  1:3518002/93‑1 (MQ=255)

                                                                   aaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                     >  2:2168862/1‑93 (MQ=255)

                                                                   aaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                     <  2:2030614/93‑1 (MQ=255)

                                                                    aaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                    >  2:2082991/1‑93 (MQ=255)

                                                                    aaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                    >  2:770916/1‑93 (MQ=255)

                                                                        ccTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                <  2:368390/93‑1 (MQ=255)

                                                                        ccTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                <  2:4116888/93‑1 (MQ=255)

                                                                        ccTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                <  1:42762/93‑1 (MQ=255)

                                                                         cTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga               <  1:1785776/93‑1 (MQ=255)

                                                                          ttCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAagag              >  1:2776954/1‑93 (MQ=255)

                                                                          ttCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAagag              <  2:1113641/93‑1 (MQ=255)

                                                                             tcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaaa           <  1:2949956/93‑1 (MQ=255)

                                                                             tcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaaa           >  1:3676915/1‑93 (MQ=255)

                                                                             tcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaaa           >  2:1985132/1‑93 (MQ=255)

                                                                              cGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAag          <  2:3197092/93‑1 (MQ=255)

                                                                              cGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAag          <  2:4148094/93‑1 (MQ=255)

                                                                               ggTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga         <  2:1160591/93‑1 (MQ=255)

                                                                                ggtAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag        >  2:4640718/3‑93 (MQ=255)

                                                                                 ttAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGt       <  2:4718238/93‑1 (MQ=255)

                                                                                 ttAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGt       >  2:4074023/1‑93 (MQ=255)

                                                                                 ttAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGt       >  1:2804711/1‑93 (MQ=255)

                                                                                  tAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGtt      <  1:1965368/93‑1 (MQ=255)

                                                                                     tGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTaca   <  2:4476747/93‑1 (MQ=255)

                                                                                      gACCGGGGGCAGTGATCTCGCTGCCCGTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac  <  2:4421849/93‑1 (MQ=255)

AGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACAC  >  NC_000913/1807988‑1808166

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