BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,787,263	C→G	G149A (GGC→GCC)	yibD ←	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,787,263	0	C	G	98.4%	597.4 / 0.7	190	G149A (GGC→GCC)	yibD	putative regulator
Reads supporting (aligned to +/- strand): ref base C (2/1); new base G (104/83); total (106/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01

Only 100 of 201 total aligned reads displayed.

TTTTTAACAATAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCAC  >  NC_000913/3787173‑3787352

 ttttAACAATAACATCACTGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCgg                                                                                        >  2:2977667/1‑93 (MQ=255)

 ttttAACAATAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGCAGAAAGCCCCATCCGCAGCCAGTCCGGGGCgg                                                                                        >  1:933014/1‑93 (MQ=255)

  tttAACAATAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGt                                                                                       >  1:2733918/1‑93 (MQ=255)

  tttAACAATAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGt                                                                                       <  1:4088500/93‑1 (MQ=255)

  tttAACAATAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGt                                                                                       >  1:4550425/1‑93 (MQ=255)

   ttAACAATAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGtt                                                                                      >  1:237823/1‑93 (MQ=255)

    tAACAATAACATCACGGCGATAAA‑CCCCCATCCAGCCAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCAGCAGTCAGTCCGGGGCGGATa                                                                                     >  2:1707963/1‑93 (MQ=255)

     aaCAATAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTaa                                                                                    <  2:2734260/93‑1 (MQ=255)

      aCAATAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAAt                                                                                   >  2:674339/1‑93 (MQ=255)

      aCAATAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAAt                                                                                   >  2:3816811/1‑93 (MQ=255)

        aaTAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGCAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAAtcc                                                                                 >  1:3191715/1‑91 (MQ=255)

        aaTAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATAc                                                                                 >  1:2201500/1‑93 (MQ=255)

        aaTAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATAc                                                                                 >  2:4648696/1‑93 (MQ=255)

         aTAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACg                                                                                >  2:1502692/1‑93 (MQ=255)

         aTAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACg                                                                                >  1:3459994/1‑93 (MQ=255)

          tAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGc                                                                               >  1:970188/1‑93 (MQ=255)

          tAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGc                                                                               >  2:2333885/1‑93 (MQ=255)

          tAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGc                                                                               >  2:2509387/1‑93 (MQ=255)

          tAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGc                                                                               >  1:3901079/1‑93 (MQ=255)

          tAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGc                                                                               >  1:131587/1‑93 (MQ=255)

           aaCATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGcc                                                                              <  1:2520565/93‑1 (MQ=255)

           aaCATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGcc                                                                              >  2:3577872/1‑93 (MQ=255)

            aCATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCg                                                                             >  2:4322684/1‑93 (MQ=255)

            aCATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCg                                                                             <  2:4485891/93‑1 (MQ=255)

             cATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCgg                                                                            >  2:2894120/1‑93 (MQ=255)

              aTCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGt                                                                           >  1:1575307/1‑93 (MQ=255)

              aTCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGATAATACGCCGGt                                                                           >  1:193956/1‑93 (MQ=255)

               tCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGtt                                                                          <  1:1994294/93‑1 (MQ=255)

               ccgcGGCAATAAACCCCCCGTCGAGACAACGGGAGTCCAACGGCGCGAAAAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGt                                                                           <  1:3281550/90‑1 (MQ=255)

                cACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTg                                                                         <  1:827971/93‑1 (MQ=255)

                 aCGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGa                                                                        <  2:4622190/93‑1 (MQ=255)

                 aCGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGa                                                                        <  2:2619991/93‑1 (MQ=255)

                 aCGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGa                                                                        >  2:3935633/1‑93 (MQ=255)

                 aCGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGa                                                                        <  2:2433101/93‑1 (MQ=255)

                  cGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATGCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAg                                                                       >  1:1794521/1‑93 (MQ=255)

                     cGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGa                                                                    >  2:1067263/1‑93 (MQ=255)

                     cGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGa                                                                    <  2:775701/93‑1 (MQ=255)

                      gATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGaa                                                                   <  1:1611944/93‑1 (MQ=255)

                        tAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAgg                                                                 >  2:2217616/1‑93 (MQ=255)

                        tAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAgg                                                                 >  2:4267992/1‑93 (MQ=255)

                         aaa‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGc                                                                >  1:3602480/1‑93 (MQ=255)

                         aaa‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGc                                                                >  1:1860225/1‑93 (MQ=255)

                          aa‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCg                                                               >  2:2171440/1‑93 (MQ=255)

                          aa‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCg                                                               <  2:3113297/93‑1 (MQ=255)

                           a‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGa                                                              >  2:274112/1‑93 (MQ=255)

                             cccccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGAt                                                             >  1:2074651/1‑93 (MQ=255)

                             cccccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGAt                                                             <  2:3520840/93‑1 (MQ=255)

                             cccccACCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGAt                                                             <  1:2784798/93‑1 (MQ=255)

                              ccccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATc                                                            <  2:2288174/93‑1 (MQ=255)

                              ccccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATc                                                            <  1:3163009/93‑1 (MQ=255)

                              ccccATCCAAACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCAAGCCCGGCGCGGTTAATACGCCGGTTGAGCGAAGGCGATc                                                            >  2:4059322/1‑93 (MQ=255)

                               cccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCg                                                           >  2:354586/1‑93 (MQ=255)

                               cccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCg                                                           >  1:1414786/1‑93 (MQ=255)

                                ccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCgg                                                          <  2:2698137/93‑1 (MQ=255)

                                ccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCgg                                                          <  1:673530/93‑1 (MQ=255)

                                  aTCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTg                                                        >  2:4183513/1‑93 (MQ=255)

                                  aTCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTg                                                        <  1:2010113/93‑1 (MQ=255)

                                   tCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTgg                                                       >  2:2395210/1‑93 (MQ=255)

                                    ccAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTggg                                                      >  1:1172464/1‑93 (MQ=255)

                                     cAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTgggg                                                     <  1:4039571/93‑1 (MQ=255)

                                      aGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGgga                                                    <  1:3355133/93‑1 (MQ=255)

                                       gACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCTGTTGAGCGAAGGCGATCGGTGGggat                                                   >  1:235451/1‑93 (MQ=255)

                                       gACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGggat                                                   <  2:494550/93‑1 (MQ=255)

                                          aaCGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGggatgga                                                >  2:3204609/1‑93 (MQ=255)

                                          aaCGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGACCGAAGGCCATCGGTGGggatggt                                                >  2:1245796/1‑92 (MQ=255)

                                           aCGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGggatggat                                               >  2:4105643/1‑93 (MQ=255)

                                           aCGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGggatggat                                               <  2:3654620/93‑1 (MQ=255)

                                           aCGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGggatggat                                               >  1:604650/1‑93 (MQ=255)

                                           aCGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGggatggat                                               <  1:4644788/93‑1 (MQ=255)

                                            cGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGAtt                                              <  2:136779/93‑1 (MQ=255)

                                               gAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGcc                                           <  1:2125757/93‑1 (MQ=255)

                                                aGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCa                                          >  1:182054/1‑93 (MQ=255)

                                                 gTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAg                                         <  1:3747744/93‑1 (MQ=255)

                                                 gTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGACCGGTGGGGATGGATTGCCAg                                         <  1:2772913/93‑1 (MQ=255)

                                                  tCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAgg                                        >  1:3531148/1‑93 (MQ=255)

                                                   ccAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGt                                       <  1:4291600/93‑1 (MQ=255)

                                                    cAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGtt                                      <  1:1752708/93‑1 (MQ=255)

                                                     aaCGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGttt                                     >  1:2905727/1‑93 (MQ=255)

                                                       cGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTtct                                   <  2:2021968/93‑1 (MQ=255)

                                                        ggCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTtctc                                  >  1:3806604/1‑93 (MQ=255)

                                                        ggCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTtctc                                  >  1:2230133/1‑93 (MQ=255)

                                                         gcgcGAAGAAAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTcc                                 <  1:4459240/93‑1 (MQ=255)

                                                              aagaaAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGttt                            >  1:4671830/1‑93 (MQ=255)

                                                                gaaAGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCa                          >  1:1376495/1‑93 (MQ=255)

                                                                 aaaGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCAc                         >  2:3070983/1‑93 (MQ=255)

                                                                  aaGCCCCATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACg                        <  2:814345/93‑1 (MQ=255)

                                                                      cccATCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGaaaa                    >  2:4236888/1‑93 (MQ=255)

                                                                         aTCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAacac                 >  1:1701838/1‑93 (MQ=255)

                                                                         aTCCGCAGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAacac                 <  1:2803804/93‑1 (MQ=255)

                                                                            cGCAGCCAGTCCGGGGCGGGTAATACGCCGGTGGAGCGAAGGCGATCGGTGGGGATGGAGTGCCAGGTGTCTCCCGTCTCACGAAAACACCAg              >  1:4200075/1‑93 (MQ=255)

                                                                               aGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCa           <  2:2063176/93‑1 (MQ=255)

                                                                               aGCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCa           >  2:192996/1‑93 (MQ=255)

                                                                                gCCAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAg          >  2:3442143/1‑93 (MQ=255)

                                                                                  cAGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCg        >  2:718247/1‑93 (MQ=255)

                                                                                   aGTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGt       <  2:1303024/93‑1 (MQ=255)

                                                                                    gTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGtt      >  1:1453903/1‑93 (MQ=255)

                                                                                    gTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGtt      >  1:649067/1‑93 (MQ=255)

                                                                                    gTCCGGGGCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCAAGAAAAAACCAGTGAGCGtt      >  2:4433400/1‑93 (MQ=255)

                                                                                      ccGGGGTGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGc    <  2:4541210/93‑1 (MQ=255)

                                                                                        ggggCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCAc  >  1:3969096/1‑93 (MQ=255)

TTTTTAACAATAACATCACGGCGATAAA‑CCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCAC  >  NC_000913/3787173‑3787352

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