BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,607,710	A→G	L332P (CTT→CCT)	uxaB ←	altronate oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,607,710	0	A	G	100.0%	1128.7 / NA	322	L332P (CTT→CCT)	uxaB	altronate oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (147/175); total (147/175)

Only 100 of 327 total aligned reads displayed.

TTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAAGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACGGTA  >  NC_000913/1607620‑1607799

ttGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGc                                                                                         <  1:8321305/93‑1 (MQ=255)

ttGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGc                                                                                         <  2:8479490/93‑1 (MQ=255)

ttGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGc                                                                                         >  2:371114/1‑93 (MQ=255)

   aGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaaa                                                                                      <  1:6303667/93‑1 (MQ=255)

   aGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaaa                                                                                      >  1:5173609/1‑93 (MQ=255)

     cgcgATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATc                                                                                    <  2:4094975/93‑1 (MQ=255)

     cgcgATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATc                                                                                    <  2:7164782/93‑1 (MQ=255)

      gcgATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATcc                                                                                   >  1:5606451/1‑93 (MQ=255)

       cgATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCa                                                                                  <  2:1655090/93‑1 (MQ=255)

       cgATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCa                                                                                  >  2:5794843/1‑93 (MQ=255)

         aTAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGt                                                                                >  2:133324/1‑93 (MQ=255)

           aGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAc                                                                              >  2:6617236/1‑93 (MQ=255)

            gACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAcc                                                                             >  2:4267891/1‑93 (MQ=255)

            gACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAcc                                                                             >  2:6303016/1‑93 (MQ=255)

            gACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAcc                                                                             <  2:2877882/93‑1 (MQ=255)

              cagcagCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCgg                                                                           <  1:3590438/93‑1 (MQ=255)

                gcagcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaa                                                                         <  1:3881887/93‑1 (MQ=255)

                gcagcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaa                                                                         >  1:3498972/1‑93 (MQ=255)

                gcagcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaa                                                                         >  2:7097718/1‑93 (MQ=255)

                   gcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataa                                                                      >  2:8354938/1‑93 (MQ=255)

                   gcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataa                                                                      >  1:5008563/1‑93 (MQ=255)

                   gcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataa                                                                      >  2:5561721/1‑93 (MQ=255)

                    cTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataat                                                                     <  2:587197/93‑1 (MQ=255)

                     tGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAGTCCAGTACCGGAATAAtt                                                                    <  2:357274/93‑1 (MQ=255)

                     tGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAAtt                                                                    >  1:8455341/1‑93 (MQ=255)

                      gATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAAttt                                                                   >  1:1962582/1‑93 (MQ=255)

                        tGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttct                                                                 <  2:2018074/93‑1 (MQ=255)

                         tcTTAATGTACGGGTTACGTAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttctt                                                                <  2:3099208/92‑1 (MQ=255)

                         gCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttctt                                                                <  2:2260974/93‑1 (MQ=255)

                            tAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGt                                                             <  2:7559555/93‑1 (MQ=255)

                             aaTGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTa                                                            <  1:7912155/93‑1 (MQ=255)

                              aTGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAg                                                           <  2:1490595/93‑1 (MQ=255)

                               tGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGa                                                          <  1:2584792/93‑1 (MQ=255)

                               tGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGa                                                          >  2:1373618/1‑93 (MQ=255)

                               tGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGa                                                          <  2:2440766/93‑1 (MQ=255)

                               tGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGa                                                          <  1:1581103/93‑1 (MQ=255)

                               tGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGa                                                          >  1:7662648/1‑93 (MQ=255)

                                gTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGAt                                                         <  1:6028676/93‑1 (MQ=255)

                                 tACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATa                                                        <  2:3953962/93‑1 (MQ=255)

                                 tACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATa                                                        <  1:3998107/93‑1 (MQ=255)

                                 tACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATa                                                        >  2:1856685/1‑93 (MQ=255)

                                 tACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATa                                                        <  1:1364118/93‑1 (MQ=255)

                                 tACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATa                                                        >  1:2434926/1‑93 (MQ=255)

                                   cGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGc                                                      >  1:7475625/1‑93 (MQ=255)

                                    gggTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCt                                                     <  1:2953636/93‑1 (MQ=255)

                                       ttACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCtttt                                                  >  1:7303143/1‑93 (MQ=255)

                                       ttACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCtttt                                                  <  2:4800508/93‑1 (MQ=255)

                                         aCGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTc                                                >  1:6966800/1‑93 (MQ=255)

                                          gggAAACGCCACGTGACTGCACTGGCGGAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAAGTTCTTCGTAGATAGCTTTTTCa                                               <  1:7236811/92‑1 (MQ=255)

                                           ggAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCaa                                              <  1:6808295/93‑1 (MQ=255)

                                           ggAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCaa                                              <  1:4277673/93‑1 (MQ=255)

                                            gAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAAc                                             <  1:6622266/93‑1 (MQ=255)

                                              agcGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCGCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGa                                           <  1:5842958/91‑1 (MQ=255)

                                               aCGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGaa                                          <  1:1749197/93‑1 (MQ=255)

                                               aCGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGaa                                          >  1:7997991/1‑93 (MQ=255)

                                               aCGCCCGGTGACTGCACTCGCGAAACATGCCAGTTCATCACGAGGCAAATCCAGTGCCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGaa                                          <  1:7776853/93‑1 (MQ=255)

                                                cGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAAt                                         <  2:8682497/93‑1 (MQ=255)

                                                   ccGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATgca                                      >  1:2787565/1‑92 (MQ=255)

                                                    cGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATgcgc                                     >  1:5915792/1‑93 (MQ=255)

                                                     ggTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCa                                    >  2:7910185/1‑93 (MQ=255)

                                                      gTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCaa                                   <  1:6084075/93‑1 (MQ=255)

                                                       tGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCaat                                  >  1:4348359/1‑92 (MQ=255)

                                                          cTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAAttt                               >  1:8148251/1‑93 (MQ=255)

                                                          cTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAAttt                               <  2:4929683/93‑1 (MQ=255)

                                                           tGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTc                              >  2:5949263/1‑93 (MQ=255)

                                                            gCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCa                             <  2:6212688/93‑1 (MQ=255)

                                                              aCTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGc                           <  1:2279126/93‑1 (MQ=255)

                                                              aCTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGc                           <  1:709077/93‑1 (MQ=255)

                                                              aCTGACGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGc                           <  1:4467163/93‑1 (MQ=255)

                                                               cTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCa                          <  1:98125/93‑1 (MQ=255)

                                                                tGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCAt                         <  1:3749636/93‑1 (MQ=255)

                                                                tGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCAt                         <  2:2674583/93‑1 (MQ=255)

                                                                 ggCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATc                        <  2:6816001/93‑1 (MQ=255)

                                                                 ggCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATc                        <  2:4596148/93‑1 (MQ=255)

                                                                  gCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCg                       >  2:5012899/1‑93 (MQ=255)

                                                                  gCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCg                       >  2:8740013/1‑93 (MQ=255)

                                                                   ggAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGt                      <  2:381695/92‑1 (MQ=255)

                                                                   cGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGt                      >  2:330973/1‑93 (MQ=255)

                                                                   cGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAAGGCGCAAATTTCAGCATCGt                      <  1:4621147/93‑1 (MQ=255)

                                                                     aaaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGtta                    >  2:7971750/1‑92 (MQ=255)

                                                                     aaaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTc                    <  2:4196437/93‑1 (MQ=255)

                                                                     aaaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTc                    >  2:1220049/1‑93 (MQ=255)

                                                                     aaaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTc                    <  2:1033581/93‑1 (MQ=255)

                                                                      aaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCa                   <  1:2130256/93‑1 (MQ=255)

                                                                      aaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCa                   <  2:2088846/93‑1 (MQ=255)

                                                                      aaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCa                   >  1:2658997/1‑93 (MQ=255)

                                                                      aaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCa                   >  1:3647360/1‑93 (MQ=255)

                                                                       aGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCAt                  >  2:5072133/1‑93 (MQ=255)

                                                                        gATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATc                 >  1:6400008/1‑93 (MQ=255)

                                                                         aTTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCg                <  1:1836780/93‑1 (MQ=255)

                                                                          ttCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGc               >  2:4705532/1‑93 (MQ=255)

                                                                             cAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTc            >  2:1951784/1‑93 (MQ=255)

                                                                                ttCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCAcc         >  2:7449245/1‑93 (MQ=255)

                                                                                 tcatcaCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCt        >  1:8695996/1‑93 (MQ=255)

                                                                                  catcaCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTa       >  1:310587/1‑93 (MQ=255)

                                                                                  catcaCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTa       >  2:8234802/1‑93 (MQ=255)

                                                                                  catcaCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTa       <  1:1115094/93‑1 (MQ=255)

                                                                                  catcaCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCAGCGTTCATCGCTTCACCTa       <  1:1188515/93‑1 (MQ=255)

                                                                                   atcaCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTAc      >  1:8075245/1‑93 (MQ=255)

                                                                                     caCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACgg    <  2:3506632/93‑1 (MQ=255)

TTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAAGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACGGTA  >  NC_000913/1607620‑1607799

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