BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,219,900	G→A	S697F (TCC→TTC)	bglX ←	beta‑D‑glucoside glucohydrolase, periplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,219,900	0	G	A	100.0%	1120.8 / NA	258	S697F (TCC→TTC)	bglX	beta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (125/133); total (125/133)

Only 100 of 261 total aligned reads displayed.

GCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACG  >  NC_000913/2219802‑2219997

gcgcCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGaaa                                                                                                 <  1:3997770/101‑1 (MQ=255)

gcgcCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGaaa                                                                                                 <  1:3998061/101‑1 (MQ=255)

  gcCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTTATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGc                                                                                               >  1:3998403/1‑101 (MQ=255)

    cTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCg                                                                                             >  1:3998407/1‑101 (MQ=255)

    cTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCg                                                                                             >  2:3998406/1‑101 (MQ=255)

      cAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACACGGCGGCTCATGAAAGCCGTc                                                                                           <  1:3998031/101‑1 (MQ=255)

        atatCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTcac                                                                                         <  2:3998022/101‑1 (MQ=255)

         tatCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTcaca                                                                                        >  1:3998411/1‑101 (MQ=255)

            cgatcgGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATcc                                                                                     >  1:3998414/1‑101 (MQ=255)

             gatcgGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCt                                                                                    >  2:3998416/1‑101 (MQ=255)

              atcgGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTg                                                                                   >  1:3998417/1‑101 (MQ=255)

              atcgGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTg                                                                                   <  1:3997888/101‑1 (MQ=255)

               tcgGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGc                                                                                  <  1:3997611/101‑1 (MQ=255)

                cgGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCa                                                                                 <  2:3997684/101‑1 (MQ=255)

                 gggAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCaa                                                                                <  1:3998119/101‑1 (MQ=255)

                  ggAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAg                                                                               >  1:3998420/1‑101 (MQ=255)

                   gAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGt                                                                              >  1:3998423/1‑101 (MQ=255)

                     aGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTAc                                                                            >  2:3998427/1‑101 (MQ=255)

                     aGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTAc                                                                            >  1:3998428/1‑101 (MQ=255)

                      gCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACa                                                                           >  2:3998431/1‑101 (MQ=255)

                       cTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAg                                                                                                <  2:3997920/79‑1 (MQ=255)

                       cTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACAt                                                                          <  1:3998010/101‑1 (MQ=255)

                       cTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACAt                                                                          >  1:3998434/1‑101 (MQ=255)

                         gACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCt                                                                        <  1:3998166/101‑1 (MQ=255)

                          aCAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTg                                                                       >  1:3998436/1‑101 (MQ=255)

                             gTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCAc                                                                    <  2:3997860/101‑1 (MQ=255)

                             gTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCAc                                                                    <  1:3997997/101‑1 (MQ=255)

                              tCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACt                                                                   >  2:3998439/1‑101 (MQ=255)

                              tCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACt                                                                   <  2:3997832/101‑1 (MQ=255)

                                 gAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGt                                                                                            <  1:3998135/73‑1 (MQ=255)

                                  aGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCg                                                               <  2:3998060/101‑1 (MQ=255)

                                   gTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGt                                                              <  1:3997893/101‑1 (MQ=255)

                                     ttCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTgg                                                            >  1:3998443/1‑101 (MQ=255)

                                         cccGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACa                                                                           >  2:3998444/1‑82 (MQ=255)

                                          ccGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCAccc                                                       <  1:3997730/101‑1 (MQ=255)

                                           cGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCt                                                      <  1:3998072/101‑1 (MQ=255)

                                            ggTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCAcc                                                        <  1:3997963/98‑1 (MQ=255)

                                            ggTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCAc                                                         <  2:3997978/97‑1 (MQ=255)

                                            ggTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTc                                                     <  1:3998240/101‑1 (MQ=255)

                                            ggTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTc                                                     <  1:3998017/101‑1 (MQ=255)

                                            ggTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTc                                                     <  2:3998019/101‑1 (MQ=255)

                                            ggTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTc                                                     >  1:3998447/1‑101 (MQ=255)

                                             gTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTcg                                                    >  1:3998449/1‑101 (MQ=255)

                                              tttCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTcgc                                                   <  1:3997840/101‑1 (MQ=255)

                                               ttCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTcgcg                                                  >  1:3998452/1‑101 (MQ=255)

                                               ttCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTcgcg                                                  >  1:3998455/1‑101 (MQ=255)

                                               ttCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTcgcg                                                  >  2:3998456/1‑101 (MQ=255)

                                                tCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCgcgc                                                 <  1:3998088/101‑1 (MQ=255)

                                                  aGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCtt                                               >  2:3998460/1‑101 (MQ=255)

                                                  aGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCtt                                               <  2:3998292/101‑1 (MQ=255)

                                                   gggTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTa                                              >  1:3998462/1‑101 (MQ=255)

                                                      tGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCg                                           <  1:3997874/101‑1 (MQ=255)

                                                      tGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCg                                           >  2:3998465/1‑101 (MQ=255)

                                                        aTTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCg                                           >  2:3998470/1‑99 (MQ=255)

                                                        aTTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGgt                                         >  1:3998468/1‑101 (MQ=255)

                                                         ttttCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATGTGCACTACCGTGGCACCCTCGCGCTTACCGgtg                                        >  1:3998474/1‑101 (MQ=255)

                                                           ttCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGtt                                      <  1:3997995/101‑1 (MQ=255)

                                                           ttCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGtt                                      <  2:3997892/101‑1 (MQ=255)

                                                            tctcAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTc                                     >  1:3998477/1‑101 (MQ=255)

                                                              tcAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGt                                   >  1:3998479/1‑101 (MQ=255)

                                                               cAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTc                                  <  2:3997939/101‑1 (MQ=255)

                                                                 aaGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCAc                                <  2:3997968/101‑1 (MQ=255)

                                                                 aaGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCAc                                >  2:3998482/1‑101 (MQ=255)

                                                                 aaGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCAc                                <  1:3998106/101‑1 (MQ=255)

                                                                 aaGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCAc                                <  1:3998197/101‑1 (MQ=255)

                                                                  aGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCAcc                               >  2:3998485/1‑101 (MQ=255)

                                                                  aGCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCAcc                               >  1:3998484/1‑101 (MQ=255)

                                                                   gCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCt                                                <  2:3998040/83‑1 (MQ=255)

                                                                   gCCTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCt                              <  1:3997912/101‑1 (MQ=255)

                                                                     cTTTCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGc                            >  1:3998490/1‑101 (MQ=255)

                                                                      tttCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCa                           >  1:3998493/1‑101 (MQ=255)

                                                                      tttCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCa                           <  2:3998079/101‑1 (MQ=255)

                                                                      tttCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCa                           >  2:3998491/1‑101 (MQ=255)

                                                                       ttCAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTTCCGGTGTTCGTCACCTGCAc                          >  2:3998495/1‑101 (MQ=255)

                                                                         cAGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGc                        >  2:3998498/1‑101 (MQ=255)

                                                                          aGCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCt                       >  1:3998501/1‑101 (MQ=255)

                                                                           gCTGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTg                      >  1:3998504/1‑101 (MQ=255)

                                                                             tGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGc                    <  1:3998056/101‑1 (MQ=255)

                                                                             tGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGc                    >  1:3998507/1‑101 (MQ=255)

                                                                             tGTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGc                    <  1:3997928/101‑1 (MQ=255)

                                                                              gTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCa                   <  2:3998184/101‑1 (MQ=255)

                                                                              gTTTCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCa                   <  1:3998300/101‑1 (MQ=255)

                                                                               tttCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAg                  <  2:3998051/101‑1 (MQ=255)

                                                                                ttCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTTGCAGt                 >  2:3998511/1‑101 (MQ=255)

                                                                                ttCACAGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGt                 >  1:3998512/1‑101 (MQ=255)

                                                                                  cacaGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCa               <  2:3998097/101‑1 (MQ=255)

                                                                                    caGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACt             >  2:3998514/1‑101 (MQ=255)

                                                                                     aGGGCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACtt            <  1:3998311/101‑1 (MQ=255)

                                                                                      gggCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACttt           <  2:3998009/101‑1 (MQ=255)

                                                                                      gggCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACttt           >  1:3998516/1‑101 (MQ=255)

                                                                                        gCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGc         <  1:3998049/101‑1 (MQ=255)

                                                                                        gCGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGc         <  2:3998115/101‑1 (MQ=255)

                                                                                         cGACTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGcc        <  1:3998161/101‑1 (MQ=255)

                                                                                           aCTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGt      <  2:3998213/101‑1 (MQ=255)

                                                                                           aCTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGt      >  1:3998521/1‑101 (MQ=255)

                                                                                           aCTCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGt      >  2:3998522/1‑101 (MQ=255)

                                                                                             tCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCa    <  1:3998193/101‑1 (MQ=255)

                                                                                             tCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCa    <  2:3998105/101‑1 (MQ=255)

                                                                                             tCATGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCa    >  1:3998524/1‑101 (MQ=255)

                                                                                               aTGAAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACg  <  2:3998044/101‑1 (MQ=255)

GCGCCTCAATATCGATCGGGAAGCTGACAGTCTGAGTTTCGCCCGGTTTCAGGGTGATTTTCTCAAAGCCTTTCAGCTGTTTCACAGGGCGACTCATGGAAGCCGTCACATCCTGCAAGTACATCTGCACTACCGTGGCACCCTCGCGCTTACCGGTGTTCGTCACCTGCACGCTGGCAGTCACTTTGCCGTCACG  >  NC_000913/2219802‑2219997

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