BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,520,082	A→G	H527H (CAT→CAC)	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,520,082	0	A	G	97.7%	549.0 / 1.3	175	H527H (CAT→CAC)	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand): ref base A (2/2); new base G (92/79); total (94/81)
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.24e-01

Only 100 of 179 total aligned reads displayed.

GGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACATGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAC  >  NC_000913/4519992‑4520171

ggggCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTg                                                                                         <  2:1667034/93‑1 (MQ=255)

agggCCTCCGGCCAGCGCTTCTGGCCCAACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTg                                                                                         <  2:1667033/92‑1 (MQ=255)

  ggCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCt                                                                                       <  2:2040692/93‑1 (MQ=255)

   gCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCtt                                                                                      <  1:1218021/93‑1 (MQ=255)

     cTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTAc                                                                                    <  2:803881/93‑1 (MQ=255)

     cTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTAc                                                                                    >  1:4553195/1‑93 (MQ=255)

      tCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTAcc                                                                                   >  1:783091/1‑93 (MQ=255)

       ccGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCg                                                                                  <  2:3316404/93‑1 (MQ=255)

       ccGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCg                                                                                  >  1:2177475/1‑93 (MQ=255)

         ggCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTa                                                                                <  2:3253745/93‑1 (MQ=255)

         ggCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGGGATTAACGAAACGTGCTTACCGTa                                                                                >  1:500272/1‑93 (MQ=255)

          gCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTAt                                                                               <  2:3499347/93‑1 (MQ=255)

           ccAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATg                                                                              >  2:1021819/1‑93 (MQ=255)

            cAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGa                                                                             >  2:313095/1‑93 (MQ=255)

             aGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGAc                                                                            >  2:648345/1‑93 (MQ=255)

             aGCGCTTCTGGCCACACATGGCCGATGCACGCGCCAGTAGACACGCCCGAGCAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGAc                                                                            >  1:835839/1‑93 (MQ=255)

              gcgcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACa                                                                           <  2:2665261/93‑1 (MQ=255)

              gcgcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACa                                                                           <  1:3763867/93‑1 (MQ=255)

                gcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACAGTATGAAAGa                                                                         >  1:4553337/1‑93 (MQ=255)

                 cTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGAt                                                                        >  1:36500/1‑93 (MQ=255)

                   tCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGc                                                                      <  2:3063799/93‑1 (MQ=255)

                      ggCCCCACATGGCCGATGCACGCGCCAGTAGAACCGCCCGAGAAACGTGCATCGGTGATTAAAGAAACATGCTTACCGTATGACAGATGCTTc                                                                   <  1:2492020/93‑1 (MQ=255)

                      ggCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTc                                                                   <  1:2695410/93‑1 (MQ=255)

                       gCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCa                                                                  <  1:1721720/93‑1 (MQ=255)

                         cccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGCGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAgg                                                                >  2:3034197/1‑93 (MQ=255)

                         cccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAgg                                                                >  2:2387530/1‑93 (MQ=255)

                         cccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAgg                                                                >  2:2874303/1‑93 (MQ=255)

                          ccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAggg                                                               >  2:4456107/1‑93 (MQ=255)

                          ccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGAAAGATGATTCAggg                                                               >  1:1423008/1‑93 (MQ=255)

                           cacaTGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGc                                                              >  1:1670322/1‑93 (MQ=255)

                             caTGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCAc                                                            <  1:4135727/93‑1 (MQ=255)

                               tGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTg                                                          <  1:3451748/93‑1 (MQ=255)

                                ggCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTgg                                                         <  1:1508980/93‑1 (MQ=255)

                                 gCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGt                                                        >  2:2758164/1‑93 (MQ=255)

                                   cGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTaa                                                      <  1:199006/93‑1 (MQ=255)

                                    gATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCATGGCACTGGTAAc                                                     >  2:3558752/1‑93 (MQ=255)

                                     aTGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAAcc                                                    >  2:2112443/1‑93 (MQ=255)

                                      tGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCt                                                   <  2:2346336/93‑1 (MQ=255)

                                       gCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTg                                                  <  1:385483/93‑1 (MQ=255)

                                       gCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTg                                                  <  1:4388666/93‑1 (MQ=255)

                                       gCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTg                                                  <  1:1718789/93‑1 (MQ=255)

                                         aCGCGCCAGTAGAAACGCCCGAGAAGCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCATTGGTAACCTGGt                                                <  2:3651540/93‑1 (MQ=255)

                                         aCGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGt                                                >  1:1842157/1‑93 (MQ=255)

                                          cgcgCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTggca                                               >  1:4655845/1‑91 (MQ=255)

                                          cgcgCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTa                                               <  2:1759580/93‑1 (MQ=255)

                                           gcgcCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTCCCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAg                                              <  1:3571651/93‑1 (MQ=255)

                                             gcCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGt                                            >  2:4709008/1‑93 (MQ=255)

                                             gcCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGt                                            >  1:385890/1‑93 (MQ=255)

                                             gcCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCGTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGt                                            >  2:4212767/1‑93 (MQ=255)

                                              ccAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGtt                                           >  1:3081405/1‑93 (MQ=255)

                                              ccAGTAGAAACGCACGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGtt                                           <  1:3262102/93‑1 (MQ=255)

                                                aGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttc                                         >  1:1408219/1‑93 (MQ=255)

                                                 gTAGAAACGCCCGAGAAACGTGCATCGGTGATTAAGAAAACGGGCTTACGGAATGACAGATGCTTCAGGGAACTGGAAACCTGGAAGGttttt                                        >  2:1647441/1‑91 (MQ=255)

                                                  tAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttctt                                       <  1:1001812/93‑1 (MQ=255)

                                                    gAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTcc                                     >  1:1041220/1‑93 (MQ=255)

                                                     aaaCGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCa                                    >  1:910936/1‑93 (MQ=255)

                                                     aaaCGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCa                                    <  2:3978702/93‑1 (MQ=255)

                                                       aCGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATc                                  >  2:1797326/1‑93 (MQ=255)

                                                       aCGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATc                                  >  2:1797327/1‑93 (MQ=255)

                                                         gCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATccc                                <  1:3072028/93‑1 (MQ=255)

                                                          cccGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCt                               >  2:4722903/1‑93 (MQ=255)

                                                          cccGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCt                               <  2:4718317/93‑1 (MQ=255)

                                                           ccGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTg                              >  2:3694236/1‑93 (MQ=255)

                                                            cGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGt                             >  2:3045671/1‑93 (MQ=255)

                                                            cGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGt                             >  1:47944/1‑93 (MQ=255)

                                                            cGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGt                             <  2:4307373/93‑1 (MQ=255)

                                                            cGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGt                             >  1:262342/1‑93 (MQ=255)

                                                             gagaAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTa                            >  1:2904486/1‑93 (MQ=255)

                                                              agaAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTAc                           >  1:3049316/1‑93 (MQ=255)

                                                                  aCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGa                       >  2:1883386/1‑93 (MQ=255)

                                                                  aCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGa                       <  1:4617206/93‑1 (MQ=255)

                                                                   cGTGCATCTGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGaa                      <  2:1306814/93‑1 (MQ=255)

                                                                   cGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCGCTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGaa                      <  1:1538594/93‑1 (MQ=255)

                                                                     tGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAgg                    <  2:1390684/93‑1 (MQ=255)

                                                                        aTCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTcc                 <  1:2314192/93‑1 (MQ=255)

                                                                        aTCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTcc                 <  1:4073731/93‑1 (MQ=255)

                                                                        aTCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTcc                 <  2:1556186/93‑1 (MQ=255)

                                                                        aTCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTca                 >  1:3018767/1‑92 (MQ=255)

                                                                          cGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCaa               <  1:2078219/93‑1 (MQ=255)

                                                                          cGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCaa               >  1:3527721/1‑93 (MQ=255)

                                                                           ggTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAAc              <  1:2143674/93‑1 (MQ=255)

                                                                             tGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGc            >  2:445335/1‑93 (MQ=255)

                                                                              gATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGcc           >  1:244372/1‑93 (MQ=255)

                                                                              gATTAACGAAACATGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGcc           >  2:4106954/1‑93 (MQ=255)

                                                                               aTTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCa          >  1:1158494/1‑93 (MQ=255)

                                                                               aTTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCa          >  2:4606178/1‑93 (MQ=255)

                                                                               aTTAACGAAACATGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCa          <  2:4106171/93‑1 (MQ=255)

                                                                                ttAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaa         >  2:4350490/1‑93 (MQ=255)

                                                                                 tAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaat        >  1:3526110/1‑93 (MQ=255)

                                                                                 tAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGGTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaat        <  1:608441/93‑1 (MQ=255)

                                                                                  aaCGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaata       >  1:3940139/1‑93 (MQ=255)

                                                                                   aCGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaataa      <  2:2277292/93‑1 (MQ=255)

                                                                                   aCGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaataa      <  2:4571230/93‑1 (MQ=255)

                                                                                     gAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATg    <  2:3006896/93‑1 (MQ=255)

                                                                                     aaaaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATg    <  2:2303833/92‑1 (MQ=255)

                                                                                      aaaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGa   >  1:4465602/1‑93 (MQ=255)

                                                                                      aaaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGa   <  2:4744215/93‑1 (MQ=255)

                                                                                      aaaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGa   <  1:4325547/93‑1 (MQ=255)

                                                                                       aaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAc  <  2:1778460/93‑1 (MQ=255)

                                                                                       aaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAc  >  1:4435745/1‑93 (MQ=255)

GGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACATGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAC  >  NC_000913/4519992‑4520171

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