BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,622,791	C→T	P110S (CCG→TCG)	serB →	3‑phosphoserine phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,622,791	0	C	T	100.0%	740.4 / NA	205	P110S (CCG→TCG)	serB	3‑phosphoserine phosphatase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (91/114); total (91/114)

Only 100 of 213 total aligned reads displayed.

GCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATG  >  NC_000913/4622701‑4622880

gCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCg                                                                                         >  2:669766/1‑93 (MQ=255)

 cAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCgg                                                                                        >  2:2546080/1‑93 (MQ=255)

  aGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCggg                                                                                       >  2:737247/1‑93 (MQ=255)

  aGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCggg                                                                                       >  1:3328626/1‑93 (MQ=255)

    gTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGtt                                                                                     <  1:48525/93‑1 (MQ=255)

     ttCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGttt                                                                                    >  1:1375460/1‑93 (MQ=255)

       cactcacCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGc                                                                                  >  1:1147429/1‑93 (MQ=255)

       cactcacCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGc                                                                                  <  2:3960003/93‑1 (MQ=255)

        actcacCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCt                                                                                 >  1:4504788/1‑93 (MQ=255)

         ctcacCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTg                                                                                >  1:3562198/1‑93 (MQ=255)

           cacCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTggg                                                                              >  2:1629316/1‑92 (MQ=255)

           cacCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGt                                                                              <  2:2617586/93‑1 (MQ=255)

            acCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTg                                                                             >  2:2078771/1‑93 (MQ=255)

             ccGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGa                                                                            <  2:1616566/93‑1 (MQ=255)

               gCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATg                                                                          <  1:3421271/93‑1 (MQ=255)

                 aCGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGa                                                                        >  1:617154/1‑93 (MQ=255)

                  cGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATgggt                                                                       >  2:2471712/1‑91 (MQ=255)

                  cGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGat                                                                       <  2:378304/93‑1 (MQ=255)

                   gggCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGata                                                                      <  1:3887542/93‑1 (MQ=255)

                   gggCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGata                                                                      <  1:2104906/93‑1 (MQ=255)

                    ggCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGatat                                                                     <  1:3965725/93‑1 (MQ=255)

                    ggCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGatat                                                                     <  2:1277690/93‑1 (MQ=255)

                    ggCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGatat                                                                     <  2:4124115/93‑1 (MQ=255)

                     gCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATg                                                                    <  2:3374687/93‑1 (MQ=255)

                     gCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATg                                                                    <  1:2677465/93‑1 (MQ=255)

                     gCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGATGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATg                                                                    >  1:2516495/1‑93 (MQ=255)

                       tACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGa                                                                  >  1:2858772/1‑93 (MQ=255)

                        acacGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGAc                                                                 <  2:2117896/93‑1 (MQ=255)

                           cGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTcc                                                              >  2:1946683/1‑93 (MQ=255)

                             ccTGGCCCACGAAGCGCAGGCTGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCAc                                                            <  2:598600/93‑1 (MQ=255)

                             ccTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCAc                                                            >  2:2260583/1‑93 (MQ=255)

                             ccTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCAc                                                            >  2:2267203/1‑93 (MQ=255)

                             ccTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCAc                                                            >  1:923859/1‑93 (MQ=255)

                              cTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCAcc                                                           >  2:588575/1‑93 (MQ=255)

                               tGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCg                                                          >  1:4118375/1‑93 (MQ=255)

                                ggCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGc                                                         >  1:3664940/1‑93 (MQ=255)

                                 gCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGcc                                                        >  1:578523/1‑93 (MQ=255)

                                     ccGAAGCGCAGCCGGATGTCGCCCCGCTGGGGAAAATCCCGCGCCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGCCTCCTCCGCCATcc                                                    <  2:3005366/92‑1 (MQ=255)

                                      cGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCa                                                   >  2:1432576/1‑93 (MQ=255)

                                      cGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCa                                                   <  2:4206416/93‑1 (MQ=255)

                                      cGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGACATCCa                                                   >  1:2322626/1‑93 (MQ=255)

                                       gAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAg                                                  >  1:2018596/1‑93 (MQ=255)

                                        aaGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGa                                                 <  1:4580891/93‑1 (MQ=255)

                                         aGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGAt                                                <  1:705200/93‑1 (MQ=255)

                                         aGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGAt                                                <  2:3829002/93‑1 (MQ=255)

                                         aGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGAt                                                <  1:4251643/93‑1 (MQ=255)

                                         aGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGAt                                                <  2:1898871/93‑1 (MQ=255)

                                         aGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGAt                                                <  1:3829024/93‑1 (MQ=255)

                                         aGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGAt                                                <  1:3729919/93‑1 (MQ=255)

                                          gctcAGCTGGATGTCGCCCCGCTGGGTAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGAtt                                               <  2:3369415/90‑1 (MQ=255)

                                          gcgcAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGAtt                                               <  2:3461079/93‑1 (MQ=255)

                                          gcgcAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGAtt                                               >  1:2242280/1‑93 (MQ=255)

                                          gcgcAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGAtt                                               <  1:240335/93‑1 (MQ=255)

                                          gcgcAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGAtt                                               <  1:3133599/93‑1 (MQ=255)

                                           cgcAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTg                                              <  1:4401094/93‑1 (MQ=255)

                                            gcAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGa                                             >  1:4318788/1‑93 (MQ=255)

                                             cAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGaa                                            <  1:678454/93‑1 (MQ=255)

                                             cAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGaa                                            >  1:2806453/1‑93 (MQ=255)

                                             cAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGaa                                            >  2:1460052/1‑93 (MQ=255)

                                              aGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAAt                                           <  2:2585621/93‑1 (MQ=255)

                                               gCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATg                                          >  1:824996/1‑93 (MQ=255)

                                               gCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATg                                          <  2:795091/93‑1 (MQ=255)

                                               gCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATg                                          <  2:838899/93‑1 (MQ=255)

                                                 tGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTa                                        >  2:1394645/1‑93 (MQ=255)

                                                 tGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGCTTGCTGGCGATGGATATGGACTCGACAGCCATCAAGATTGAAAGTa                                        >  2:1768606/1‑93 (MQ=255)

                                                  ggATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTAt                                       <  2:3638327/93‑1 (MQ=255)

                                                   gATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTAtt                                      >  1:3665839/1‑93 (MQ=255)

                                                   gATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTAtt                                      >  2:268149/1‑93 (MQ=255)

                                                   gATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTAtt                                      >  2:4147324/1‑93 (MQ=255)

                                                   gATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTAtt                                      >  2:4257269/1‑93 (MQ=255)

                                                    aTGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCGCCGCCATCCAGATTGAATGTATtg                                     <  1:2644790/93‑1 (MQ=255)

                                                     tGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATtga                                    <  1:4017836/93‑1 (MQ=255)

                                                     tGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATtga                                    >  2:233580/1‑93 (MQ=255)

                                                     tGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATtga                                    <  1:1986728/93‑1 (MQ=255)

                                                      gTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATtgat                                   <  2:3571784/93‑1 (MQ=255)

                                                        cGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATtgatga                                 <  1:4210650/93‑1 (MQ=255)

                                                        cGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATtgatga                                 <  1:4406119/93‑1 (MQ=255)

                                                         gCCCCGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGaa                                >  2:4494888/1‑93 (MQ=255)

                                                          ccccGCTGGGGAAAATCCCGCACCTGCGGACGTCGGGTTTGCTGGTGATGAATATGGACTCCGCCGCCATCCAGATTGAATGTATTGATGaaa                               <  2:2177765/93‑1 (MQ=255)

                                                          ccccGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGaaa                               <  2:1107856/93‑1 (MQ=255)

                                                          ccccGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGaaa                               <  2:2656793/93‑1 (MQ=255)

                                                           cccGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAAt                              <  2:2011606/93‑1 (MQ=255)

                                                           cccGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAAt                              >  1:299730/1‑93 (MQ=255)

                                                            ccGCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAAtt                             >  2:2433399/1‑93 (MQ=255)

                                                              gCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGc                           <  1:1243238/93‑1 (MQ=255)

                                                              gCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGc                           <  1:213102/93‑1 (MQ=255)

                                                              gCTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGc                           <  2:1559163/93‑1 (MQ=255)

                                                               cTGGGGAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGcc                          <  1:3179988/93‑1 (MQ=255)

                                                                   ggAAAATCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAAc                      >  2:3278133/1‑93 (MQ=255)

                                                                       aaTCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGc                  <  1:4464746/93‑1 (MQ=255)

                                                                       aaTCCCGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGc                  <  2:538246/93‑1 (MQ=255)

                                                                           ccGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGa              <  1:1532637/93‑1 (MQ=255)

                                                                            cGCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGaa             <  1:4483631/93‑1 (MQ=255)

                                                                            cGCACCTGCGCACGTCGGGTTTGCTGGAGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGaa             <  1:940252/93‑1 (MQ=255)

                                                                             gCACCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAAc            <  2:3492948/93‑1 (MQ=255)

                                                                               aCCTGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACgg          >  2:1116522/1‑93 (MQ=255)

                                                                                  tGCGCACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCg       >  1:4341262/1‑93 (MQ=255)

                                                                                   gcgcACGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCga      <  1:1427812/93‑1 (MQ=255)

                                                                                       aCGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGAtg  <  2:1289055/93‑1 (MQ=255)

                                                                                       aCGTCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGAtg  >  1:186382/1‑93 (MQ=255)

GCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATG  >  NC_000913/4622701‑4622880

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: `‑`