BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,802	G→A	E57K (GAG→AAG)	yabB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,802	0	G	A	99.5%	728.3 / ‑5.5	202	E57K (GAG→AAG)	yabB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (120/80); total (122/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

Only 100 of 206 total aligned reads displayed.

CTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATG  >  NC_000913/89712‑89891

ctTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAg                                                                                         >  2:2944714/1‑93 (MQ=255)

 ttGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAga                                                                                        >  2:4390951/1‑92 (MQ=255)

  tGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCa                                                                                       <  1:1295247/93‑1 (MQ=255)

    agaACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCaaa                                                                                     >  1:3963033/1‑93 (MQ=255)

     gaACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCaaaa                                                                                    >  1:3449410/1‑93 (MQ=255)

      aaCGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCaaaaa                                                                                   <  2:120956/93‑1 (MQ=255)

        cGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCTCCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAAtt                                                                                 <  2:2144048/93‑1 (MQ=255)

         gCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTa                                                                                <  1:164014/93‑1 (MQ=255)

         gCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTa                                                                                >  2:3220259/1‑93 (MQ=255)

           tGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTAtt                                                                              <  1:3279699/93‑2 (MQ=255)

           tGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTAtt                                                                              >  1:3007556/1‑92 (MQ=255)

           tGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTAtt                                                                              >  2:2484016/1‑92 (MQ=255)

           tGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTAtt                                                                              <  2:780384/93‑2 (MQ=255)

            gCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTAttg                                                                             <  2:4504213/93‑3 (MQ=255)

            gCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTAttg                                                                             <  1:167170/93‑3 (MQ=255)

             ccGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATtgc                                                                            <  2:1998715/93‑1 (MQ=255)

             ccGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATtgc                                                                            <  2:3932653/93‑1 (MQ=255)

             ccGGTCAAATGGATTGCCCCATTGACATTTATCCCCCGTGCCTGCTGCTTTATCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATtgc                                                                            <  2:809664/93‑1 (MQ=255)

              cGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGTATGGGAAATTATCAAGCAAAAATTATtgcg                                                                           >  2:347018/1‑93 (MQ=255)

               ggTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGt                                                                          >  2:2534583/1‑93 (MQ=255)

               ggTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGt                                                                          >  1:3749882/1‑93 (MQ=255)

               ggTCAAATGGTTTGCACCATTGACATTTATCACCCATGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGt                                                                          >  2:4562975/1‑93 (MQ=255)

                gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                         >  2:1869095/1‑93 (MQ=255)

                   aaaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGt                                                                      >  1:3095604/1‑93 (MQ=255)

                   aaaTGATTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGt                                                                      >  2:3054703/1‑93 (MQ=255)

                    aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                     <  1:229672/93‑1 (MQ=255)

                           tGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATg                                                              <  1:1081880/93‑1 (MQ=255)

                           tGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATg                                                              <  1:4500127/93‑1 (MQ=255)

                            gCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGa                                                             <  2:3602440/93‑1 (MQ=255)

                            gCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGa                                                             >  1:4166563/1‑93 (MQ=255)

                             cACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGaa                                                            <  2:1071466/93‑1 (MQ=255)

                             cACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGaa                                                            <  1:4612575/93‑1 (MQ=255)

                               ccATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAcc                                                          <  1:1921935/93‑1 (MQ=255)

                               ccATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAcc                                                          <  2:2328163/93‑1 (MQ=255)

                               ccATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAcc                                                          >  1:4673109/1‑93 (MQ=255)

                               ccATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAcc                                                          >  1:2872911/1‑93 (MQ=255)

                               ccATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAcc                                                          >  1:4056380/1‑93 (MQ=255)

                                 aTTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCg                                                        <  2:277922/93‑1 (MQ=255)

                                  ttGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCgg                                                       >  1:2111635/1‑93 (MQ=255)

                                  ttGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCgg                                                       <  1:3481615/93‑1 (MQ=255)

                                  ttGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCgg                                                       >  1:3824404/1‑93 (MQ=255)

                                       aTTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGGATGAACCCGGTTGAg                                                  >  2:1619571/1‑93 (MQ=255)

                                        tttATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAgc                                                 >  1:924028/1‑93 (MQ=255)

                                         ttATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAgcg                                                <  1:178857/93‑1 (MQ=255)

                                              aCCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCtgtg                                           >  1:2741700/1‑93 (MQ=255)

                                               cccGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgtgt                                          >  2:3384339/1‑93 (MQ=255)

                                               cccGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgtgt                                          >  2:1935466/1‑93 (MQ=255)

                                               cccGTGCCTGCTGCTTTACCCCCTGCCTGAATGGCAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTCGAGCGCCgtgt                                          >  2:225630/1‑93 (MQ=255)

                                                 cGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGc                                        >  2:497908/1‑93 (MQ=255)

                                                     ccTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgcg                                    >  2:368191/1‑93 (MQ=255)

                                                      ctgctgCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgcgc                                   >  1:2390368/1‑93 (MQ=255)

                                                       tgctgcTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGcc                                  >  2:664964/1‑93 (MQ=255)

                                                       tgctgcTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGcc                                  >  1:4401035/1‑93 (MQ=255)

                                                        gctgctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCt                                 <  2:2025871/93‑1 (MQ=255)

                                                          tgctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTAc                               <  2:3534982/93‑1 (MQ=255)

                                                          tgctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTAc                               >  1:858680/1‑93 (MQ=255)

                                                          tgctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTAc                               >  2:1741132/1‑93 (MQ=255)

                                                           gctTTACCGCCTGCCTGAATGGGAAATTATCAAGCAAAATTTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACt                              <  2:3856809/93‑1 (MQ=255)

                                                              ttACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGtt                           >  1:2918759/1‑93 (MQ=255)

                                                               tACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTa                          <  2:2862790/93‑1 (MQ=255)

                                                                aCCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAg                         <  1:2355581/93‑1 (MQ=255)

                                                                 cccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAgg                        <  1:2687182/93‑1 (MQ=255)

                                                                  ccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGt                       <  2:3168222/93‑1 (MQ=255)

                                                                  ccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGt                       >  1:2787067/1‑93 (MQ=255)

                                                                    cctgcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCa                     <  1:2926888/93‑1 (MQ=255)

                                                                     ctgcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCAt                    >  2:1419746/1‑93 (MQ=255)

                                                                     ctgcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCAt                    <  2:4134410/93‑1 (MQ=255)

                                                                      tgcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATg                   <  1:3560806/93‑1 (MQ=255)

                                                                       gcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGc                  >  1:304055/1‑93 (MQ=255)

                                                                         ctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCa                >  1:1139202/1‑93 (MQ=255)

                                                                          tgcaTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCAAGAACCCGGATGAACGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAg               >  2:4068955/4‑93 (MQ=255)

                                                                          tgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAg               >  2:421696/1‑93 (MQ=255)

                                                                           gAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGc              <  2:4023171/93‑1 (MQ=255)

                                                                           gAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGc              <  1:1250511/93‑1 (MQ=255)

                                                                            aaTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCg             >  1:2070359/1‑93 (MQ=255)

                                                                             aTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGa            >  2:3781235/1‑93 (MQ=255)

                                                                             aTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGa            >  2:4656026/1‑93 (MQ=255)

                                                                              tGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGaa           <  1:4492972/93‑1 (MQ=255)

                                                                              tGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGaa           <  2:1297680/93‑1 (MQ=255)

                                                                              tGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGaa           <  1:647428/93‑1 (MQ=255)

                                                                               gggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAAt          >  1:3053541/1‑93 (MQ=255)

                                                                               gggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAAt          <  1:2174842/93‑1 (MQ=255)

                                                                                ggAATTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg         >  1:3981966/1‑93 (MQ=255)

                                                                                ggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg         >  2:694409/1‑93 (MQ=255)

                                                                                ggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg         <  2:318653/93‑1 (MQ=255)

                                                                                ggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg         >  2:2667891/1‑93 (MQ=255)

                                                                                ggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg         <  2:21767/93‑1 (MQ=255)

                                                                                ggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg         >  2:1493444/1‑93 (MQ=255)

                                                                                ggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg         <  2:2037424/93‑1 (MQ=255)

                                                                                ggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg         <  2:4335051/93‑1 (MQ=255)

                                                                                 gAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGt        >  2:3119755/1‑93 (MQ=255)

                                                                                   aaTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCa      >  2:2375082/1‑93 (MQ=255)

                                                                                    aTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAg     <  1:1704450/93‑1 (MQ=255)

                                                                                     ttATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAga    >  1:3163970/1‑93 (MQ=255)

                                                                                     ttATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAga    >  1:1825861/1‑93 (MQ=255)

                                                                                      tATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgat   <  2:4249960/93‑1 (MQ=255)

                                                                                      tATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgat   >  1:703387/1‑93 (MQ=255)

                                                                                      tATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgat   >  2:3419309/1‑93 (MQ=255)

                                                                                       aTCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatg  >  1:4568617/1‑93 (MQ=255)

                                                                                       aTCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatg  <  1:1482506/93‑1 (MQ=255)

CTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATG  >  NC_000913/89712‑89891

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