BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,088	+C	coding (2/624 nt)	yhgKJ ←	RNA phosphage cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,554,088	1	.	C	100.0%	705.7 / NA	220	S133W (TCG→TGG)	yhgKJ	RNA phosphage cyclase
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (119/101); total (119/101)

Only 100 of 239 total aligned reads displayed.

GGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT  >  NC_000913/3553997‑3554174

 gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑Cgca                                                                                        <  1:476150/93‑3 (MQ=255)

     ggTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAggcg                                                                                    <  2:2891196/93‑1 (MQ=255)

      gTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAggcgg                                                                                   <  2:1643803/93‑1 (MQ=255)

       tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAggcggc                                                                                  >  2:5259769/1‑93 (MQ=255)

       tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAggcggc                                                                                  <  2:2610462/93‑1 (MQ=255)

       tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAggcggc                                                                                  <  2:2758271/93‑1 (MQ=255)

        aaaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGgcggcg                                                                                 >  2:3553539/1‑93 (MQ=255)

         aaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGgcgc                                                                                <  2:3243161/93‑1 (MQ=255)

         aaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGgcgc                                                                                >  1:923809/1‑93 (MQ=255)

         aaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGgcgc                                                                                <  1:3706881/93‑1 (MQ=255)

         aaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTACTATTTTCGCCAGCAGCGGCTCCAGCACCCAGCGGATAAAATC‑CGCAGGCGgcgc                                                                                >  1:2038302/1‑93 (MQ=255)

            aCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCgc                                                                             >  2:1983036/1‑92 (MQ=255)

            aCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGa                                                                             >  2:2681335/1‑93 (MQ=255)

            aCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGa                                                                             >  2:2489279/1‑93 (MQ=255)

             ccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGAc                                                                            <  2:5579086/93‑1 (MQ=255)

               gTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACgg                                                                          >  2:3259467/1‑93 (MQ=255)

               gTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACgg                                                                          >  1:5817563/1‑93 (MQ=255)

                tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACggg                                                                         >  2:3477576/1‑93 (MQ=255)

                 ggCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGt                                                                        >  2:1441727/1‑93 (MQ=255)

                  gCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGtt                                                                       <  2:5360943/93‑1 (MQ=255)

                  gCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGtt                                                                       <  1:2751197/93‑1 (MQ=255)

                   cGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGttt                                                                      >  1:5147281/1‑92 (MQ=255)

                   cGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCAACGGGTTa                                                                      <  1:6518325/93‑1 (MQ=255)

                    gTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTAt                                                                     <  2:2972273/93‑1 (MQ=255)

                    gTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGGTAt                                                                     >  1:349620/1‑93 (MQ=255)

                    gTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAAGCGGCGCCGACGGGTTAt                                                                     >  1:6674027/1‑93 (MQ=255)

                    ctAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTAt                                                                     <  1:3267901/92‑1 (MQ=255)

                     tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATc                                                                    <  2:3446892/93‑1 (MQ=255)

                     tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATc                                                                    <  1:5597910/93‑1 (MQ=255)

                      agcAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGGTATCg                                                                   >  1:5312610/3‑93 (MQ=255)

                      aaCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCg                                                                   >  2:2432798/1‑93 (MQ=255)

                       aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCTGCGGATAAAATC‑CGCAGGCGGCGCCGACCGGTTATCgg                                                                  >  2:397013/1‑93 (MQ=255)

                         aGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGtg                                                                >  2:5842539/1‑92 (MQ=255)

                         aGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGtg                                                                <  1:5984011/93‑2 (MQ=255)

                          gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGtgc                                                               <  2:5498561/93‑3 (MQ=255)

                          gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATA‑CGCAGGCGGCGCCGACGGGTTATCGGtgc                                                               >  1:4865348/1‑91 (MQ=255)

                            gTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGtgccg                                                             >  1:3233908/1‑89 (MQ=255)

                             ttgTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                          <  1:1886217/91‑1 (MQ=255)

                             tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                          >  1:1184793/1‑93 (MQ=255)

                             tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                          <  2:4663685/93‑1 (MQ=255)

                              ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                         <  1:634892/93‑2 (MQ=255)

                                tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccgc                                                       <  2:2087357/93‑4 (MQ=255)

                                   gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCa                                                     <  2:3648075/93‑1 (MQ=255)

                                    cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                                    >  1:2875726/1‑93 (MQ=255)

                                      gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACtt                                                  >  1:1456706/1‑93 (MQ=255)

                                        tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCa                                                >  1:1438275/1‑93 (MQ=255)

                                         gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCaa                                               >  2:6104735/1‑93 (MQ=255)

                                           aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                             >  2:1523296/1‑93 (MQ=255)

                                            ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAacac                                            >  2:5004053/1‑93 (MQ=255)

                                              ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgc                                          <  2:6510905/93‑1 (MQ=255)

                                               cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                         >  1:2171326/1‑93 (MQ=255)

                                               cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                         >  1:3142799/1‑93 (MQ=255)

                                               cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                         >  1:5036827/1‑93 (MQ=255)

                                               cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                         <  1:5884196/93‑1 (MQ=255)

                                                 aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGACGCTCACTTCATCACGCGaa                                       <  1:6567894/93‑1 (MQ=255)

                                                   tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAgg                                     <  1:6119632/93‑1 (MQ=255)

                                                    ttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGt                                    <  2:1313045/93‑1 (MQ=255)

                                                     tCGCCAGCAGCGCCGCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTTAACACGCGAAGGTc                                   <  1:4064809/93‑1 (MQ=255)

                                                       gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCa                                 <  1:4254765/93‑1 (MQ=255)

                                                          agcagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                              >  2:6269557/1‑93 (MQ=255)

                                                          agcagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGCTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCCTCg                              >  1:2230111/1‑93 (MQ=255)

                                                           gcagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                             >  2:4557130/1‑93 (MQ=255)

                                                           gcagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                             >  1:5257897/1‑93 (MQ=255)

                                                           gcagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGG‑CGCCGCCGCTCACCTCAACACGCGAAGGTCCATCgg                             >  2:2572462/1‑93 (MQ=255)

                                                            cagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                            <  2:5096698/93‑1 (MQ=255)

                                                            cagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                            >  2:4190803/1‑93 (MQ=255)

                                                            cagcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                            >  2:4428300/1‑93 (MQ=255)

                                                              gcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                          >  1:3386405/1‑93 (MQ=255)

                                                              gcGGCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                          >  1:4402191/1‑93 (MQ=255)

                                                                ggCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                        >  1:6224836/1‑93 (MQ=255)

                                                                 gCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                       >  2:3133858/1‑93 (MQ=255)

                                                                 gCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                       <  2:6706/93‑1 (MQ=255)

                                                                 gCTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGGCCATCGGCAAAcc                       >  1:228725/1‑93 (MQ=255)

                                                                  cTCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                      >  1:1619799/1‑93 (MQ=255)

                                                                   tCCAGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                     >  1:5183524/1‑93 (MQ=255)

                                                                    ctaGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGTTCCATCGGCAAACcaca                    <  1:2374253/91‑1 (MQ=255)

                                                                      aGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                  >  2:929899/1‑93 (MQ=255)

                                                                      aGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                  >  1:6454007/1‑93 (MQ=255)

                                                                      aGCACCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                  <  2:2349757/93‑1 (MQ=255)

                                                                         aCCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg               >  2:3888818/1‑93 (MQ=255)

                                                                         aCCCGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg               >  2:4937950/1‑93 (MQ=255)

                                                                          cccGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg              >  2:1215372/1‑93 (MQ=255)

                                                                          cccGGCGGATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg              >  1:696881/1‑93 (MQ=255)

                                                                            cggcggATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc            >  1:3981924/1‑93 (MQ=255)

                                                                            cggcggATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc            <  1:3649741/93‑1 (MQ=255)

                                                                             ggcggAGAAAAGC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca           <  1:2423795/93‑1 (MQ=255)

                                                                             cgcggATAAAATCGAC‑AGGCGCGGCGGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca           <  1:4351551/92‑1 (MQ=255)

                                                                                ggATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca        >  2:5610986/1‑93 (MQ=255)

                                                                                 gATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       >  2:4785525/1‑93 (MQ=255)

                                                                                 gATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       >  2:3850952/1‑93 (MQ=255)

                                                                                 gATAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       <  1:5533053/93‑1 (MQ=255)

                                                                                  aTAAAATC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc      >  1:1244931/1‑93 (MQ=255)

                                                                                    aaaaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt    >  2:3739162/1‑93 (MQ=255)

                                                                                     aaaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCACCACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt   <  2:771253/93‑1 (MQ=255)

                                                                                     aaaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGTGCGGGCAGCACCGtt   >  2:3816963/1‑93 (MQ=255)

                                                                                     aaaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt   >  2:3927438/1‑93 (MQ=255)

                                                                                     aaaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt   <  1:2136774/93‑1 (MQ=255)

                                                                                     aaaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt   <  1:3870242/93‑1 (MQ=255)

                                                                                     aaaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt   <  2:1572357/93‑1 (MQ=255)

                                                                                      aaTC‑CGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt  >  1:1102627/1‑93 (MQ=255)

GGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATC‑CG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT  >  NC_000913/3553997‑3554174

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