BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,594,520	A→G	S531P (TCC→CCC)	ydeK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,594,520	0	A	G	99.3%	859.8 / ‑3.2	270	S531P (TCC→CCC)	ydeK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (136/132); total (137/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

Only 100 of 280 total aligned reads displayed.

ATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAT  >  NC_000913/1594430‑1594609

aTACTCCGGCCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑Ggggg                                                                                         <  2:3056908/93‑1 (MQ=255)

 tACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGc                                                                                        <  2:3352684/93‑1 (MQ=255)

  acacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCa                                                                                       <  1:4278254/93‑1 (MQ=255)

  acacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCa                                                                                       <  2:7400602/93‑1 (MQ=255)

   cacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCAt                                                                                      <  1:6041704/93‑1 (MQ=255)

   cacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCAt                                                                                      >  2:188977/1‑93 (MQ=255)

    acCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCAtt                                                                                     >  2:6522410/1‑93 (MQ=255)

      cGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTtc                                                                                   >  1:1547751/1‑93 (MQ=255)

        aCCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTtcct                                                                                 >  1:4894416/1‑93 (MQ=255)

        aCCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTtcct                                                                                 >  2:7551378/1‑93 (MQ=255)

         ccTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTtcctc                                                                                <  1:9179911/93‑1 (MQ=255)

         ccTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTtcctc                                                                                >  2:4701246/1‑93 (MQ=255)

         ccTGTAGCAATTGTGCGTTAGTAGATGATGCCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTtcctc                                                                                <  2:2040084/93‑1 (MQ=255)

          cTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTtcctcc                                                                               <  1:6480872/93‑1 (MQ=255)

           tGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCg                                                                              >  2:7219091/1‑93 (MQ=255)

            gTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCgt                                                                             >  1:5607745/1‑92 (MQ=255)

             tAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGAt                                                                            <  2:2274609/93‑1 (MQ=255)

              aGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATa                                                                           <  1:7991132/93‑1 (MQ=255)

              aGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATa                                                                           <  2:7735909/93‑1 (MQ=255)

              aGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATa                                                                           <  2:3661/93‑1 (MQ=255)

               gCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATaa                                                                          <  2:10016652/93‑1 (MQ=255)

                cAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCACCTTACC‑GGGGGCATTTCCTCCGATaaa                                                                         >  1:1013959/1‑93 (MQ=255)

                  aTTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATc                                                                       <  2:3461020/93‑1 (MQ=255)

                    tgtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCa                                                                     <  1:543823/93‑1 (MQ=255)

                      tgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATa                                                                   <  2:2541126/93‑1 (MQ=255)

                       gCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAg                                                                  >  2:5188145/1‑93 (MQ=255)

                       gCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAg                                                                  <  1:9863342/93‑1 (MQ=255)

                         gTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGcc                                                                <  1:7491854/93‑1 (MQ=255)

                         gTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCATTACC‑GGGGGCATTTCCTCCGATAAATCCATAGcc                                                                <  1:1842991/93‑1 (MQ=255)

                          ttAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGccc                                                               >  1:6291972/1‑93 (MQ=255)

                          ttAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGccc                                                               <  2:382754/93‑1 (MQ=255)

                           tagtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCg                                                              >  2:8917483/1‑93 (MQ=255)

                            agtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGt                                                             >  2:7469992/1‑93 (MQ=255)

                               agATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGa                                                          >  2:287706/1‑93 (MQ=255)

                                gatgatGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCGCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGaa                                                         >  2:6808087/1‑93 (MQ=255)

                                gatgatGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGaa                                                         <  1:3725042/93‑1 (MQ=255)

                                 atgatgTCTTTAAATTCCAGAGACCGTCCGTTGAAACATTAACTGTCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAAt                                                        <  1:6065161/93‑1 (MQ=255)

                                  tgatgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAAtt                                                       >  2:5217329/1‑93 (MQ=255)

                                  tgatgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGGTGaagt                                                       >  2:1474084/1‑91 (MQ=255)

                                   gatgTCTTTAAATTCCAGAGACTGTCCGTTGCAATATTAACTCTCCCCTTACC‑GGGGGCATTTCCTCCGGTAAATCCAGAGGCCGTTGCATTa                                                      >  1:8153292/1‑93 (MQ=255)

                                     tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    <  2:6103464/93‑1 (MQ=255)

                                     tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    >  2:560626/1‑93 (MQ=255)

                                     tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    >  1:865629/1‑93 (MQ=255)

                                     tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    >  1:8411718/1‑93 (MQ=255)

                                     tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    >  1:7021538/1‑93 (MQ=255)

                                      gTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAgg                                                   >  1:3240207/1‑93 (MQ=255)

                                       tCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAggg                                                  <  1:2224815/93‑1 (MQ=255)

                                        cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGt                                                 <  2:9015162/93‑1 (MQ=255)

                                        cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGt                                                 <  2:870792/93‑1 (MQ=255)

                                          ttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTgg                                               >  2:567690/1‑93 (MQ=255)

                                          ttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTgg                                               <  2:603863/93‑1 (MQ=255)

                                          ttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTgg                                               >  2:1147407/1‑93 (MQ=255)

                                           tAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGc                                              <  1:9483211/93‑1 (MQ=255)

                                           tAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGc                                              >  1:9635632/1‑93 (MQ=255)

                                           tAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGc                                              >  1:3383364/1‑93 (MQ=255)

                                            aaaTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGGCCGTTGTATTAAGGGTGGcg                                             >  2:5520491/1‑92 (MQ=255)

                                            aaaTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGcg                                             <  2:9496293/93‑2 (MQ=255)

                                                 ccAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCAtt                                        <  1:1251708/93‑1 (MQ=255)

                                                   agagACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTAc                                      <  1:7884098/93‑1 (MQ=255)

                                                    gagaCTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACt                                     <  1:7672267/93‑1 (MQ=255)

                                                    gagaCTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACt                                     <  2:6298212/93‑1 (MQ=255)

                                                     agaCTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTg                                    <  2:8903047/93‑1 (MQ=255)

                                                     agaCTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTg                                    <  1:9509830/93‑1 (MQ=255)

                                                        cTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                 <  2:4312796/93‑1 (MQ=255)

                                                        cTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                 >  2:883326/1‑93 (MQ=255)

                                                        cTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                 >  2:978976/1‑93 (MQ=255)

                                                          gTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                               >  2:8247321/1‑93 (MQ=255)

                                                          gTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                               >  2:5006025/1‑93 (MQ=255)

                                                           tCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAttt                              >  1:8968838/1‑93 (MQ=255)

                                                            ccGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTa                             <  2:4658409/93‑1 (MQ=255)

                                                            ccGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTa                             <  2:1242048/93‑1 (MQ=255)

                                                            ccGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTa                             <  1:4035488/93‑1 (MQ=255)

                                                            ccGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTa                             <  1:5262657/93‑1 (MQ=255)

                                                             cGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTAc                            <  2:2195451/93‑1 (MQ=255)

                                                               ttGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                          >  1:1747394/1‑93 (MQ=255)

                                                               ttGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                          <  1:2601301/93‑1 (MQ=255)

                                                                tGAAATATTAACTATCCCCTTACCGGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                          <  1:5096257/93‑1 (MQ=255)

                                                                tGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGAt                         >  1:9325485/1‑93 (MQ=255)

                                                                 gAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATa                        >  2:2726253/1‑93 (MQ=255)

                                                                 gAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATa                        <  1:3934034/93‑1 (MQ=255)

                                                                  aaaTATTAAC‑ATCCCCTTACC‑CGGGGCATTTTCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATaac                      >  1:2100311/1‑91 (MQ=255)

                                                                     tatTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGc                    <  1:2353961/93‑1 (MQ=255)

                                                                     tatTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGc                    >  2:1580500/1‑93 (MQ=255)

                                                                        tAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTAc                 <  1:1365963/93‑1 (MQ=255)

                                                                        tAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTAc                 <  2:7892056/93‑1 (MQ=255)

                                                                        tAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTAc                 >  1:861086/1‑93 (MQ=255)

                                                                         aaCTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTAcc                <  2:6444014/93‑1 (MQ=255)

                                                                            tATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa             >  1:5476436/1‑93 (MQ=255)

                                                                             aTCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGATATTTACGATACCGCTACCaaaa            >  1:4548248/1‑93 (MQ=255)

                                                                             aTCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaaa            <  1:1187646/93‑1 (MQ=255)

                                                                              tCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAt           <  1:4641995/93‑1 (MQ=255)

                                                                              tCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAt           <  2:4205782/93‑1 (MQ=255)

                                                                                 ccTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc        <  1:8516825/93‑1 (MQ=255)

                                                                                  cTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCCAAATTAcc       >  2:8355775/1‑93 (MQ=255)

                                                                                  cTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAcc       <  2:7086861/93‑1 (MQ=255)

                                                                                   ttACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCg      >  2:36120/1‑93 (MQ=255)

                                                                                    tccc‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGCGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGa     <  1:3002175/91‑1 (MQ=255)

                                                                                    tACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGa     >  1:6643227/1‑93 (MQ=255)

                                                                                     aCC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGAc    >  2:2258755/1‑93 (MQ=255)

                                                                                      cc‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa   <  1:5897583/93‑1 (MQ=255)

ATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAT  >  NC_000913/1594430‑1594609

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