BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,725,518	(G)_6→5	coding (836/1098 nt)	nemA →	N‑ethylmaleimide reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,725,513	0	G	.	100.0%	570.3 / NA	125	coding (831/1098 nt)	nemA	N‑ethylmaleimide reductase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base . (57/68); total (57/68)

Only 100 of 141 total aligned reads displayed.

GGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGT  >  NC_000913/1725421‑1725599

ggCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCg                                                                                        >  1:411956/1‑93 (MQ=255)

 gCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgg                                                                                       <  2:3097826/93‑1 (MQ=255)

    cGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCggggg                                                                                    >  2:2130086/1‑93 (MQ=255)

      aaTGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCCTATGTCAGAACCAGATTGGGCgggggtg                                                                                  >  1:1303966/1‑91 (MQ=255)

      aaTGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggtg                                                                                  >  1:2982936/1‑91 (MQ=255)

      aaTGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggtg                                                                                  >  1:1727337/1‑91 (MQ=255)

      aaTGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGAATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggtg                                                                                  >  1:1303967/1‑91 (MQ=255)

      aaTGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGAGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCggggggg                                                                                  >  2:1562514/1‑92 (MQ=255)

       aTGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggtga                                                                                 <  1:2429590/93‑4 (MQ=255)

        tGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggtgaa                                                                                <  1:2321600/93‑5 (MQ=255)

        tGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACAAGATTGGGCgggggtgaa                                                                                >  2:470078/1‑89 (MQ=255)

         gaagaaGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACTCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAAc                                                                              >  1:2677342/1‑93 (MQ=255)

           agaagCCGATGCACTGTATCTGATTGATCAACTGAGTAAACGCGGCATTGCTTATCTGCATATGTCACAGCCAGATTGGGC‑GGGGGTGAACCg                                                                            <  2:2344046/93‑1 (MQ=255)

            aaagCCGATGCACTGTATCTGATTGAACAACTGGGTAATCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGt                                                                           <  1:1001574/92‑1 (MQ=255)

            aaagCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGt                                                                           <  2:1077230/92‑1 (MQ=255)

              agCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGtat                                                                         >  1:477690/1‑93 (MQ=255)

              agCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGtat                                                                         >  2:1373930/1‑93 (MQ=255)

               gCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGtata                                                                        >  2:4320139/1‑93 (MQ=255)

               gCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGtata                                                                        <  1:2196538/93‑1 (MQ=255)

                 cGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACt                                                                      <  2:3489889/93‑1 (MQ=255)

                 cGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACt                                                                      >  2:553750/1‑93 (MQ=255)

                  gATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTg                                                                     <  1:1035869/93‑1 (MQ=255)

                   aTGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGa                                                                    >  1:2020156/1‑93 (MQ=255)

                   aTGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGa                                                                    <  2:4575810/93‑1 (MQ=255)

                   aTGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGa                                                                    <  1:3636588/93‑1 (MQ=255)

                    tGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGAt                                                                   <  1:3545316/93‑1 (MQ=255)

                     gCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATg                                                                  <  2:3991434/93‑1 (MQ=255)

                     gCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATg                                                                  >  2:2559213/1‑93 (MQ=255)

                     gCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATg                                                                  <  1:2271407/93‑1 (MQ=255)

                     gCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGCGAACCGTATACTGATg                                                                  >  2:2559214/1‑93 (MQ=255)

                      cACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGc                                                                 <  1:2094136/93‑1 (MQ=255)

                       aCTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCg                                                                >  2:4002711/1‑93 (MQ=255)

                       aCTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCg                                                                <  2:2761981/93‑1 (MQ=255)

                       aCTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCg                                                                <  2:2616393/93‑1 (MQ=255)

                        cTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGt                                                               >  1:1635164/1‑93 (MQ=255)

                        cTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGt                                                               >  1:1635163/1‑93 (MQ=255)

                          gTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTc                                                             >  1:3646399/1‑93 (MQ=255)

                           tATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAATCCTATACTGATGCGTTcc                                                            >  2:4427363/1‑93 (MQ=255)

                            aTCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCcg                                                           <  1:1840044/93‑1 (MQ=255)

                            aTCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCcg                                                           <  1:476282/93‑1 (MQ=255)

                             tCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCcgc                                                          >  1:3423064/1‑93 (MQ=255)

                              cTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCcgcg                                                         >  1:3592117/1‑93 (MQ=255)

                                  ttGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGaaaa                                                     >  2:3502990/1‑93 (MQ=255)

                                   tGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGaaaaa                                                    <  1:1337607/93‑1 (MQ=255)

                                   tGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGaaaaa                                                    <  2:1440175/93‑1 (MQ=255)

                                      acaacTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTa                                                 <  1:1306836/93‑1 (MQ=255)

                                      acaacTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTa                                                 <  1:1989358/93‑1 (MQ=255)

                                      acaacTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTa                                                 <  1:3009227/93‑1 (MQ=255)

                                        aacTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAcg                                               >  2:3290519/1‑93 (MQ=255)

                                        aacTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAcg                                               <  2:2901/93‑1 (MQ=255)

                                        aacTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAcg                                               <  2:2693035/93‑1 (MQ=255)

                                         acTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAcgc                                              <  2:2301534/93‑1 (MQ=255)

                                         acTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAcgc                                              <  2:1857534/93‑1 (MQ=255)

                                         acTGGGCAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAcgc                                              >  2:3446960/1‑93 (MQ=255)

                                           tGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACgcgc                                            >  2:1191091/1‑93 (MQ=255)

                                           tGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTACGCGAAAAAGTACgcgc                                            >  2:2972770/1‑93 (MQ=255)

                                             ggTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGccc                                          >  1:325104/1‑93 (MQ=255)

                                             ggTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGccc                                          <  1:2553863/93‑1 (MQ=255)

                                               tAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGt                                        >  1:2614315/1‑93 (MQ=255)

                                                aaaCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGtt                                       <  1:2009410/93‑1 (MQ=255)

                                                aaaCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGtt                                       <  2:1709575/93‑1 (MQ=255)

                                                aaaCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGtt                                       <  2:4186294/93‑1 (MQ=255)

                                                  aCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTc                                     <  2:1270310/93‑1 (MQ=255)

                                                  aCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTc                                     <  2:2532158/93‑1 (MQ=255)

                                                    gcgGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTccc                                   >  2:3632735/1‑92 (MQ=255)

                                                    gcgGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCa                                   <  1:448811/93‑1 (MQ=255)

                                                    gcgGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGCTTCCa                                   >  1:300218/1‑93 (MQ=255)

                                                     cgGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCAc                                  <  2:1410741/93‑1 (MQ=255)

                                                      ggCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACg                                 >  2:1917630/1‑93 (MQ=255)

                                                      ggCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACg                                 <  1:2681759/93‑1 (MQ=255)

                                                       gaaTTGCTTATCTGCATATGTCAGAACCAGATGGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACgg                                <  1:3433385/91‑1 (MQ=255)

                                                       gCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACgg                                <  1:4072594/93‑1 (MQ=255)

                                                       gCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACgg                                <  2:2897606/93‑1 (MQ=255)

                                                        cATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGt                               <  1:1200261/93‑1 (MQ=255)

                                                         aTTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTc                              <  1:1506490/93‑1 (MQ=255)

                                                         aTTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTc                              >  2:3305948/1‑93 (MQ=255)

                                                            gCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGa                           >  1:4514876/1‑93 (MQ=255)

                                                               tATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTa                        <  2:1175548/93‑1 (MQ=255)

                                                                  cTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCg                     <  1:1903179/93‑1 (MQ=255)

                                                                   tGCATATGTCGGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGGACGCGCCCGTTTCCACGCTCCGATTATCgg                    >  2:162504/1‑93 (MQ=255)

                                                                   tGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAGCCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCgg                    >  2:162372/1‑93 (MQ=255)

                                                                    gCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGgc                   >  1:3153368/1‑93 (MQ=255)

                                                                     cATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTATGAGCTTGTTTCCACGGTCCGATTATCGgcg                  >  2:3958529/1‑93 (MQ=255)

                                                                      atatGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGgcgc                 >  1:4739838/1‑93 (MQ=255)

                                                                        atGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAg               <  1:1079922/93‑1 (MQ=255)

                                                                         tGTCAGAACCAGATTGGGC‑GTGGGTGAACCGTATACTGATTCGTTCCGTGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAgg              <  1:2725559/93‑1 (MQ=255)

                                                                         tGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAgg              >  1:667332/1‑93 (MQ=255)

                                                                          gTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAggg             >  1:1945728/1‑92 (MQ=255)

                                                                          gTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGt             <  2:4389285/93‑1 (MQ=255)

                                                                           tCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTg            <  2:4566187/93‑1 (MQ=255)

                                                                             aGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCa          <  2:1808718/93‑1 (MQ=255)

                                                                             aGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCa          <  1:4185811/93‑1 (MQ=255)

                                                                              gAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCAt         >  1:1644333/1‑93 (MQ=255)

                                                                                aCCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATac       >  1:1187098/1‑93 (MQ=255)

                                                                                 ccAGATTGGGCGGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATac       <  1:2731681/93‑1 (MQ=255)

                                                                                 ccAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATaca      <  2:4231525/93‑1 (MQ=255)

                                                                                    gATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAg   >  2:3986305/1‑93 (MQ=255)

                                                                                    gATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAg   <  1:2639224/93‑1 (MQ=255)

                                                                                    gATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAg   >  1:3230090/1‑93 (MQ=255)

                                                                                     aTTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGt  >  1:967306/1‑93 (MQ=255)

GGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGT  >  NC_000913/1725421‑1725599

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