BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,815	C→T	S61L (TCG→TTG)	yabB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,815	0	C	T	99.4%	1135.5 / ‑4.6	310	S61L (TCG→TTG)	yabB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (180/128); total (181/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 326 total aligned reads displayed.

CCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTC  >  NC_000913/89725‑89904

   gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                                      <  1:2120649/93‑1 (MQ=255)

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

     cAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTg                                                                                    <  1:1711698/93‑1 (MQ=255)

     cAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTg                                                                                    <  1:5352784/93‑1 (MQ=255)

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

       aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTCCCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                                  <  2:6320176/93‑1 (MQ=255)

       aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                                  >  2:2125998/1‑93 (MQ=255)

        aTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCg                                                                                 >  2:5146921/1‑93 (MQ=255)

         tGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGa                                                                                <  1:6015471/93‑1 (MQ=255)

             ttGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCAt                                                                            >  1:7451074/1‑93 (MQ=255)

              tGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATg                                                                           >  2:3254678/1‑93 (MQ=255)

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

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

                 aCCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAc                                                                        >  2:7696553/1‑93 (MQ=255)

                  ccATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCTGAAATTATCGCGTCTGTCGATCATGAAcc                                                                       >  2:4265625/1‑93 (MQ=255)

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

                   cATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAccc                                                                      >  1:4291604/1‑93 (MQ=255)

                   cATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAccc                                                                      >  2:49475/1‑93 (MQ=255)

                    aTTGACATTTATCACCCGTGCCTGCTGCTTTACACCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCg                                                                     >  1:925693/1‑93 (MQ=255)

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

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

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

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

                       gACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGtt                                                                  >  1:3944762/1‑93 (MQ=255)

                        aCATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTg                                                                 >  1:44310/1‑93 (MQ=255)

                         cATTTATCACCCGTGCCTGCTGCTTTCCCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGa                                                                <  2:1313865/93‑1 (MQ=255)

                         cATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGa                                                                <  2:3623011/93‑1 (MQ=255)

                         cATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGa                                                                >  1:5048699/1‑93 (MQ=255)

                          aTTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAg                                                               >  2:2460812/1‑93 (MQ=255)

                           tttATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGc                                                                              <  1:2862082/77‑1 (MQ=255)

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

                                cACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgt                                                         <  2:3030423/93‑1 (MQ=255)

                                cACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgt                                                         <  1:7599466/93‑1 (MQ=255)

                                cACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgt                                                         >  2:3397109/1‑93 (MQ=255)

                                cACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgt                                                         <  2:5070796/93‑1 (MQ=255)

                                cACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgt                                                         >  2:3755260/1‑93 (MQ=255)

                                 tcccGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgtg                                                        <  1:365009/92‑1 (MQ=255)

                                 ccccGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgtg                                                        <  2:5719966/92‑1 (MQ=255)

                                 aCCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTTAGCGCCgtg                                                        >  2:2901069/1‑93 (MQ=255)

                                 aCCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgtg                                                        >  1:193801/1‑93 (MQ=255)

                                    cGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGCGCATGAACCCGGTTGCGCGCCGTGTGc                                                     >  1:3589272/1‑93 (MQ=255)

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

                                     gTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCa                                                    >  2:603412/1‑93 (MQ=255)

                                      tGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAg                                                   >  1:6338404/1‑93 (MQ=255)

                                      tGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAg                                                   >  2:1639819/1‑93 (MQ=255)

                                       gCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgc                                                  <  2:593234/93‑1 (MQ=255)

                                       gCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgc                                                  <  2:7055067/93‑1 (MQ=255)

                                       gCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgc                                                  >  1:2290807/1‑93 (MQ=255)

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

                                          tgctgcTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGcc                                               <  2:8786/93‑1 (MQ=255)

                                          tgctgcTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGcc                                               <  1:7818679/93‑1 (MQ=255)

                                           gctgctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCt                                              >  1:4583035/1‑93 (MQ=255)

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

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

                                               ctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTg                                          >  2:6809602/1‑93 (MQ=255)

                                                tttACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGt                                         >  2:2554457/1‑93 (MQ=255)

                                                tttACCCCCTGACTGAATGAAAAATTAACAAAAAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGt                                         <  2:4611654/93‑1 (MQ=255)

                                                  tACGCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTa                                       <  2:7210116/93‑1 (MQ=255)

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

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

                                                    cccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAgg                                     >  2:7868277/1‑93 (MQ=255)

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

                                                      cccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTCCTGTTAGGTc                                   >  1:7066825/1‑93 (MQ=255)

                                                      cccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTc                                   >  2:7383145/1‑93 (MQ=255)

                                                      cccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTc                                   >  2:4437971/1‑93 (MQ=255)

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

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

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

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

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

                                                                 tGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGaa                        >  2:2527144/1‑93 (MQ=255)

                                                                 tGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGaa                        >  2:408449/1‑93 (MQ=255)

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

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

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

                                                                     aaaTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTc                    <  2:1448326/93‑1 (MQ=255)

                                                                     aaaTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTc                    <  1:5882388/93‑1 (MQ=255)

                                                                       aTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCCTGCCAGCGAATGTCAg                  >  1:1106646/1‑93 (MQ=255)

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

                                                                       aTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAg                  >  2:1088179/1‑93 (MQ=255)

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

                                                                        atATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAga                 <  1:2218989/92‑1 (MQ=255)

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

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

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

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

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

                                                                            caaGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatgga             >  1:7706965/3‑93 (MQ=255)

                                                                             aaGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatggat            <  1:6628754/92‑1 (MQ=255)

                                                                             aaGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatggat            >  2:4661785/2‑93 (MQ=255)

                                                                               gCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGatggatgg          <  1:6526952/93‑1 (MQ=255)

                                                                               gCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGatggatgg          >  1:839969/1‑93 (MQ=255)

                                                                               gCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGatggatgg          >  2:7557766/1‑93 (MQ=255)

                                                                                cAAAAATTATTGCGTCTGTCGAGCATGAACCCGGATGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGgc         <  1:6081835/93‑1 (MQ=255)

                                                                                  aaaaTTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATATCAGATGGATGgcgc       >  2:5436425/1‑93 (MQ=255)

                                                                                    aaTTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAg     <  2:1788408/93‑1 (MQ=255)

                                                                                    aaTTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAg     <  1:5573960/93‑1 (MQ=255)

                                                                                     attattGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAgg    >  2:7807526/1‑93 (MQ=255)

                                                                                      atATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTTTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGAGGGATGGCGCAGGt   <  1:3242386/92‑1 (MQ=255)

                                                                                       tATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTc  <  1:473387/93‑1 (MQ=255)

CCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTC  >  NC_000913/89725‑89904

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