BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	916,045	A→G	E117G (GAG→GGG)	ybjD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	916,045	0	A	G	96.3%	397.7 / 3.2	137	E117G (GAG→GGG)	ybjD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/2); new base G (43/88); total (46/91)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 145 total aligned reads displayed.

GACATCGGGTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAG  >  NC_000913/915955‑916133

 aCATCGGGTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCCTCTGGAAGGGGggag                                                                                       <  2:2708645/93‑1 (MQ=255)

  cATCGGGTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGt                                                                                      <  1:1854622/93‑1 (MQ=255)

   aTCGGGTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTg                                                                                     >  2:2298936/1‑93 (MQ=255)

   aTCGGGTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTg                                                                                     >  1:4294660/1‑93 (MQ=255)

   aTCGGGTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTg                                                                                     >  2:4308564/1‑93 (MQ=255)

   aTCGGGTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTg                                                                                     >  2:467523/1‑93 (MQ=255)

      gggTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCgg                                                                                  <  2:532050/93‑1 (MQ=255)

       ggTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGa                                                                                 >  2:3453111/1‑93 (MQ=255)

       ggTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGa                                                                                 >  1:466627/1‑93 (MQ=255)

        gTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGaa                                                                                <  1:1797508/93‑1 (MQ=255)

        gTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGaa                                                                                <  2:3814133/93‑1 (MQ=255)

        gTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGaa                                                                                <  1:506254/93‑1 (MQ=255)

         ttCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAg                                                                               >  1:2173127/1‑93 (MQ=255)

          tCGCGGTTACCGACCGCTGGAAGCGTGCTGGGCGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCagaaga                                                                              <  2:4238373/93‑1 (MQ=255)

          tCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGa                                                                              <  1:2076554/93‑1 (MQ=255)

           cGCTGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGAc                                                                             <  2:3409475/93‑1 (MQ=255)

            cccGTTATCGGCCGGTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACg                                                                            <  2:4032775/92‑1 (MQ=255)

             ccGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACgg                                                                           <  2:171479/93‑1 (MQ=255)

              cGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGc                                                                          >  2:1064524/1‑93 (MQ=255)

              cGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGc                                                                          <  1:1080280/93‑1 (MQ=255)

               gTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCa                                                                         <  2:172886/93‑1 (MQ=255)

                ttATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAg                                                                        <  1:1497358/93‑1 (MQ=255)

                 tATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGGGCGGAAGACGGCAGc                                                                       >  1:2204584/1‑93 (MQ=255)

                     ggCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGtga                                                                   <  1:1584605/93‑1 (MQ=255)

                     ggCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGtga                                                                   <  2:3380651/93‑1 (MQ=255)

                     ggCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGtga                                                                   <  2:1440912/93‑1 (MQ=255)

                      gCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGtgat                                                                  <  2:1884525/93‑1 (MQ=255)

                        cGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGtgatga                                                                <  1:2568442/93‑1 (MQ=255)

                        cGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGtgatga                                                                <  1:1206472/93‑1 (MQ=255)

                         gCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGac                                                               <  1:599370/93‑1 (MQ=255)

                         gCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAGCGTGATGac                                                               >  1:2601128/1‑93 (MQ=255)

                            ggAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACt                                                            <  1:4141889/93‑1 (MQ=255)

                            ggAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACt                                                            >  2:875583/1‑93 (MQ=255)

                                gCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGGGCGGAAGACGGCAGCGTGATGACACTgcgc                                                        >  2:3402739/1‑93 (MQ=255)

                                gCGTGCTGGACGACATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGACGGGGGGATGGCGGAAGACGGCAGCGTGAGGACACTgcgc                                                        >  2:4643216/1‑93 (MQ=255)

                                  gTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAg                                                      >  1:252572/1‑93 (MQ=255)

                                    gCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGtt                                                    <  2:3083949/93‑1 (MQ=255)

                                      tGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAGCGCGATGACACTGCGCAGtttt                                                  >  2:1425556/1‑93 (MQ=255)

                                      cgcaCGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGCGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGtttt                                                  <  1:1629725/90‑1 (MQ=255)

                                        gACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTtct                                                <  1:4644398/93‑1 (MQ=255)

                                         gcGCAATGCACCGATGCCTATCACCGTATTCTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTGTtctc                                               <  2:2150/92‑1 (MQ=255)

                                          cGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCg                                              >  2:2787382/1‑93 (MQ=255)

                                           gCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGa                                             <  1:2865596/93‑1 (MQ=255)

                                           gCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGa                                             <  1:1810003/93‑1 (MQ=255)

                                             cATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATa                                           <  1:2963390/93‑1 (MQ=255)

                                              aTGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATaa                                          >  2:321468/1‑93 (MQ=255)

                                              aTGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGAAGCGTGATGACACTGCGCAGTTTTCTCGATaa                                          >  2:491754/1‑93 (MQ=255)

                                               tGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATaaa                                         <  2:2901373/93‑1 (MQ=255)

                                               tGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATaaa                                         <  2:3897098/93‑1 (MQ=255)

                                               tGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATaaa                                         <  1:3766287/93‑1 (MQ=255)

                                                gCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAg                                        <  2:1867187/93‑1 (MQ=255)

                                                  aCCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGAc                                      >  1:2175000/1‑93 (MQ=255)

                                                   ccGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAAACg                                     >  1:676322/1‑93 (MQ=255)

                                                     gATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGa                                   <  1:15660/93‑1 (MQ=255)

                                                     gATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGa                                   >  2:3846262/1‑93 (MQ=255)

                                                     gATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGa                                   >  1:1626154/1‑93 (MQ=255)

                                                     gATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGa                                   >  1:983148/1‑93 (MQ=255)

                                                       tGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACa                                 <  2:1233716/93‑1 (MQ=255)

                                                         gCTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATc                               <  1:2555851/93‑1 (MQ=255)

                                                          cTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATcc                              <  1:4404917/93‑1 (MQ=255)

                                                          cTATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATcc                              <  2:4052833/93‑1 (MQ=255)

                                                          cTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATcc                              <  2:2771284/93‑1 (MQ=255)

                                                           tATCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCg                             >  2:2649131/1‑93 (MQ=255)

                                                            aTCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGa                            <  2:3933537/93‑1 (MQ=255)

                                                            aTCACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGa                            <  1:486476/93‑1 (MQ=255)

                                                              cACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGAtt                          <  2:1474683/93‑1 (MQ=255)

                                                              cACCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGAtt                          <  2:2035607/93‑1 (MQ=255)

                                                               aCCGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTg                         >  2:1377168/1‑93 (MQ=255)

                                                                 cGTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGAt                       <  1:3284733/93‑1 (MQ=255)

                                                                  gTATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATg                      <  2:2921725/93‑1 (MQ=255)

                                                                  gTATTTTTTATCGTCTGGAAGGGGAGAGGGCGGAAGACGGAAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGAAGGACATCCGATTGATg                      >  1:1209593/1‑93 (MQ=255)

                                                                   tATTTTTTATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGt                     >  1:1905966/1‑93 (MQ=255)

                                                                   tATTTTTTATCGTCTGGAAGGGGGGAGTGCAGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGt                     <  1:2256251/93‑1 (MQ=255)

                                                                     ttttttATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGCTGTCg                   >  1:1139968/1‑93 (MQ=255)

                                                                     ttttttATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCg                   <  1:647060/93‑1 (MQ=255)

                                                                     ttttttATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCg                   >  1:1295220/1‑93 (MQ=255)

                                                                       ttttATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAg                 <  2:2302243/93‑1 (MQ=255)

                                                                       tgttATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAg                 <  1:437970/91‑1 (MQ=255)

                                                                          tATCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGat              <  1:3382080/93‑1 (MQ=255)

                                                                          tATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGat              <  2:4541028/93‑1 (MQ=255)

                                                                           aTCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGata             <  2:3071393/93‑1 (MQ=255)

                                                                           aTCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGata             <  1:1586547/93‑1 (MQ=255)

                                                                            tCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGatat            >  2:2855408/1‑93 (MQ=255)

                                                                            tCGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGatat            <  1:2199762/93‑1 (MQ=255)

                                                                            tCGTCTGGAAGGGGGGAGTGAGGAAGACGGCATCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGatat            <  1:1525775/93‑1 (MQ=255)

                                                                             cGTCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATAtt           <  1:244501/93‑1 (MQ=255)

                                                                               tCTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTaa         <  1:199117/93‑1 (MQ=255)

                                                                               gcTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTaa         <  1:322352/92‑1 (MQ=255)

                                                                                cTGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAAc        <  1:4210866/93‑1 (MQ=255)

                                                                                 tGGAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACg       >  1:2338019/1‑93 (MQ=255)

                                                                                  ggAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGa      <  1:1125869/93‑1 (MQ=255)

                                                                                  ggAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGa      >  2:3398908/1‑93 (MQ=255)

                                                                                  ggAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGa      <  2:1075621/93‑1 (MQ=255)

                                                                                  ggAAGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGa      <  2:1601063/93‑1 (MQ=255)

                                                                                    aaGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATc    >  1:3459648/1‑93 (MQ=255)

                                                                                    aaGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATc    <  1:1266296/93‑1 (MQ=255)

                                                                                     aGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCa   <  1:2738794/93‑1 (MQ=255)

                                                                                     aGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCa   <  1:1956805/93‑1 (MQ=255)

                                                                                     aGGGGGGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCa   <  2:3282898/93‑1 (MQ=255)

                                                                                     aGGGGGGAGTGCGGAAGACGGAAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCa   >  2:436807/1‑93 (MQ=255)

GACATCGGGTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAGCGTGATGACACTGCGCAGTTTTCTCGATAAAGACGGACATCCGATTGATGTCGAGGATATTAACGATCAG  >  NC_000913/915955‑916133

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