BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,934,746	A→G	K281K (AAA→AAG)	rbsB →	D‑ribose periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,934,746	0	A	G	100.0%	1183.8 / NA	356	K281K (AAA→AAG)	rbsB	D‑ribose periplasmic binding protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (181/175); total (181/175)

Only 100 of 358 total aligned reads displayed.

TGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGT  >  NC_000913/3934656‑3934835

tgGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGt                                                                                         <  2:2443036/93‑1 (MQ=255)

 ggCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGtt                                                                                        >  1:3770991/1‑93 (MQ=255)

  gCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTc                                                                                       >  2:5141206/1‑93 (MQ=255)

   gaagcTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCa                                                                                      <  2:2218549/89‑1 (MQ=255)

   cAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCa                                                                                      <  2:2273958/93‑1 (MQ=255)

    aaaCTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAg                                                                                     >  2:1275153/1‑93 (MQ=255)

     aaCTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAgg                                                                                    >  2:1167316/1‑93 (MQ=255)

      aCTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGc                                                                                   <  2:2382470/93‑1 (MQ=255)

      aCTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGc                                                                                   <  2:558093/93‑1 (MQ=255)

       cTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCt                                                                                  >  1:339665/1‑93 (MQ=255)

         agcagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTaa                                                                                >  1:2978703/1‑93 (MQ=255)

          gcagcGACTATCGCGCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAg                                                                               <  2:1545900/93‑1 (MQ=255)

            agcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTa                                                                             <  2:1606573/93‑1 (MQ=255)

            agcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTa                                                                             <  2:3929450/93‑1 (MQ=255)

            agcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTa                                                                             <  2:5768498/93‑1 (MQ=255)

             gcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTAt                                                                            >  1:2806937/1‑93 (MQ=255)

             gcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTAt                                                                            >  2:4151405/1‑93 (MQ=255)

              cGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATc                                                                           <  2:911242/93‑1 (MQ=255)

               gACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATcc                                                                          <  1:1350859/93‑1 (MQ=255)

               caCTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATcc                                                                          <  1:4946541/92‑1 (MQ=255)

                 cTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCgg                                                                        <  2:2876290/93‑1 (MQ=255)

                  tATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGt                                                                       <  2:3195882/93‑1 (MQ=255)

                    tCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTg                                                                     >  1:5512801/1‑93 (MQ=255)

                       cTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATc                                                                  <  1:4259840/93‑1 (MQ=255)

                       cTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATc                                                                  >  2:4918102/1‑93 (MQ=255)

                            cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGaaa                                                             <  2:3791272/93‑1 (MQ=255)

                             tACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAAc                                                            <  1:6479373/93‑1 (MQ=255)

                                ccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTgg                                                         <  2:4500956/93‑1 (MQ=255)

                                 cGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGgt                                                        >  2:800846/1‑93 (MQ=255)

                                  gATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGgtt                                                       <  2:371466/93‑1 (MQ=255)

                                  gATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGgtt                                                       <  1:3221906/93‑1 (MQ=255)

                                   aTCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGgttg                                                      >  1:4541251/1‑93 (MQ=255)

                                    tCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGgttgt                                                     >  2:6709121/1‑93 (MQ=255)

                                       gATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTaa                                                  >  1:4369388/1‑93 (MQ=255)

                                       gATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTaa                                                  <  1:1857184/93‑1 (MQ=255)

                                        aTTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAg                                                 >  1:5370221/1‑93 (MQ=255)

                                         ttGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGc                                                >  2:4752126/1‑93 (MQ=255)

                                         ttGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTAGTTGTTAAGc                                                >  2:4942811/1‑93 (MQ=255)

                                          tGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCa                                               <  1:6260698/93‑1 (MQ=255)

                                           ggCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCag                                              >  2:1426631/1‑93 (MQ=255)

                                             cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCagta                                            <  1:5766252/93‑1 (MQ=255)

                                             cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCagta                                            <  2:1461080/93‑1 (MQ=255)

                                             cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCagta                                            <  2:6105768/93‑1 (MQ=255)

                                               ctaaaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCagtagt                                          <  1:6409236/91‑1 (MQ=255)

                                               cgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCagtagt                                          <  1:1672029/93‑1 (MQ=255)

                                                gAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGtt                                         >  1:713138/1‑93 (MQ=255)

                                                gAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGtt                                         <  2:5235683/93‑1 (MQ=255)

                                                gAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGtt                                         <  2:2627010/93‑1 (MQ=255)

                                                gAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGtt                                         >  2:3188980/1‑93 (MQ=255)

                                                  aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGtttt                                       >  1:2727794/1‑93 (MQ=255)

                                                  aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGtttt                                       >  1:5280405/1‑93 (MQ=255)

                                                   aGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTa                                      >  2:5177502/1‑93 (MQ=255)

                                                   aGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTa                                      >  1:1004247/1‑93 (MQ=255)

                                                    ggCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTaa                                     >  2:19851/1‑93 (MQ=255)

                                                       gTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCa                                  <  2:2483421/93‑1 (MQ=255)

                                                       gTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCa                                  <  1:5062205/93‑1 (MQ=255)

                                                       gTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCa                                  <  2:1858499/93‑1 (MQ=255)

                                                        tCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAg                                 >  1:20955/1‑93 (MQ=255)

                                                          gAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTCAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGt                               >  2:3564606/1‑93 (MQ=255)

                                                          gAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGt                               <  1:4445524/93‑1 (MQ=255)

                                                          gAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGt                               <  1:4004916/93‑1 (MQ=255)

                                                          gAAACCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGt                               <  2:4084010/93‑1 (MQ=255)

                                                           aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGtt                              >  2:2828446/1‑93 (MQ=255)

                                                           aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGtt                              >  1:1178453/1‑93 (MQ=255)

                                                           aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGtt                              >  2:120643/1‑93 (MQ=255)

                                                           aaaCCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGtt                              >  1:6690851/1‑93 (MQ=255)

                                                            aaCCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTg                             <  2:2601465/93‑1 (MQ=255)

                                                            aaCCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTg                             >  1:3089894/1‑93 (MQ=255)

                                                            aaCCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTg                             >  2:6455532/1‑93 (MQ=255)

                                                             aCCGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGt                            >  2:3242034/1‑93 (MQ=255)

                                                              ccGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTa                           <  1:976685/93‑1 (MQ=255)

                                                              ccGCAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTa                           <  2:1794469/93‑1 (MQ=255)

                                                                 cAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGa                        >  1:4825066/1‑93 (MQ=255)

                                                                 cAGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGa                        >  1:3570859/1‑93 (MQ=255)

                                                                  aGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGAc                       >  2:6339947/1‑93 (MQ=255)

                                                                  aGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGAc                       >  2:5073569/1‑93 (MQ=255)

                                                                  aGATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGAc                       >  1:3587464/1‑93 (MQ=255)

                                                                   gATAAAGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGAcc                      >  2:4432636/1‑93 (MQ=255)

                                                                      aaaGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGa                   <  2:2157155/93‑1 (MQ=255)

                                                                      aaaGTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGa                   >  2:5283721/1‑93 (MQ=255)

                                                                        aGTGCTGAAAGGCGAGAAGGTTCAGTCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATg                 <  1:2612814/93‑1 (MQ=255)

                                                                         gTGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATgg                <  1:585228/93‑1 (MQ=255)

                                                                         gTGCTGAAAGGCGAGAAGGTTCAGGCTAAGCATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATgg                >  1:3256754/1‑93 (MQ=255)

                                                                          tGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGt               <  2:5844617/93‑1 (MQ=255)

                                                                          tGCTGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGt               >  1:5660084/1‑93 (MQ=255)

                                                                            cTGAAAGGCGAGAAGGTTCAGGCTGAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGa             <  1:6588367/93‑1 (MQ=255)

                                                                             tGAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGAc            <  1:249585/93‑1 (MQ=255)

                                                                              gAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACa           >  1:6176282/1‑93 (MQ=255)

                                                                              gAAAGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACa           >  1:5457699/1‑93 (MQ=255)

                                                                               aaaGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACAt          >  2:5409926/1‑93 (MQ=255)

                                                                                aaGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATa         >  1:5982583/1‑93 (MQ=255)

                                                                                aaGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATa         >  2:4051795/1‑93 (MQ=255)

                                                                                aaGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATa         >  2:3143115/1‑93 (MQ=255)

                                                                                aaGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGATGTATGACCTGATGGTGACATa         >  1:3809332/1‑93 (MQ=255)

                                                                                 aGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACCTaa        >  2:2970749/1‑93 (MQ=255)

                                                                                 aGGCGAGAAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATaa        >  1:2274807/1‑93 (MQ=255)

                                                                                     gagaAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATAc    >  1:5103092/1‑93 (MQ=255)

                                                                                      agaAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACg   <  2:1684701/93‑1 (MQ=255)

                                                                                      agaAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACg   >  1:177887/1‑93 (MQ=255)

                                                                                       gaAGGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGt  >  1:846738/1‑93 (MQ=255)

TGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGT  >  NC_000913/3934656‑3934835

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