BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,968,566	A→G	L287L (TTA→TTG)	wecB →	UDP‑N‑acetyl glucosamine ‑2‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,968,566	0	A	G	100.0%	868.4 / NA	273	L287L (TTA→TTG)	wecB	UDP‑N‑acetyl glucosamine ‑2‑epimerase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (159/114); total (159/114)

Only 100 of 282 total aligned reads displayed.

GTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGT  >  NC_000913/3968477‑3968655

gTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCg                                                                                        <  2:1782796/93‑1 (MQ=255)

gTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATACCCAGGAGTATTTGCCg                                                                                        >  1:2749828/1‑93 (MQ=255)

 tGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGt                                                                                       <  2:294429/93‑1 (MQ=255)

 tGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGt                                                                                       <  2:2790935/93‑1 (MQ=255)

   cATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGttt                                                                                     >  1:357424/1‑93 (MQ=255)

   cATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGttt                                                                                     >  2:3375674/1‑93 (MQ=255)

    aTCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTg                                                                                    >  1:4474716/1‑93 (MQ=255)

    aTCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTg                                                                                    <  1:577575/93‑1 (MQ=255)

     tctcAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGt                                                                                   >  1:1855687/1‑93 (MQ=255)

        cAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTg                                                                                >  1:801647/1‑93 (MQ=255)

         aaCCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTgg                                                                               >  1:2041516/1‑93 (MQ=255)

         aaCCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTgg                                                                               >  1:1215704/1‑93 (MQ=255)

          aCCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGc                                                                              <  1:957238/93‑1 (MQ=255)

          aCCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGc                                                                              >  2:4408433/1‑93 (MQ=255)

          aCCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGc                                                                              >  2:234921/1‑93 (MQ=255)

          aCCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGc                                                                              >  1:1470683/1‑93 (MQ=255)

           cccGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGCATTTGCCGTTTGTCTGGCt                                                                             >  2:2730975/1‑93 (MQ=255)

            ccGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtg                                                                            >  2:2269507/1‑93 (MQ=255)

             cGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtga                                                                           <  2:2686045/93‑1 (MQ=255)

                aCGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtgatga                                                                        >  2:3987673/1‑93 (MQ=255)

                aCGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtgatga                                                                        <  2:1633120/93‑1 (MQ=255)

                aCGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtgatga                                                                        >  2:1249528/1‑93 (MQ=255)

                aCGTCAGAGAACCGATCAATCGCATTCTGGGGAATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtgatga                                                                        <  2:2894788/93‑1 (MQ=255)

                  gTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAAc                                                                      <  1:3851736/93‑1 (MQ=255)

                  gTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAAc                                                                      >  2:1853364/1‑93 (MQ=255)

                   tCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAAcc                                                                     >  1:3736544/1‑93 (MQ=255)

                    cAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCa                                                                    <  1:3416959/93‑1 (MQ=255)

                     agagAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCAc                                                                   >  1:4270013/1‑93 (MQ=255)

                     agagAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCAc                                                                   >  1:4449197/1‑93 (MQ=255)

                     agagAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCAc                                                                   <  1:2118940/93‑1 (MQ=255)

                      gagaACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACg                                                                  >  1:3498340/1‑93 (MQ=255)

                        gaACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGcc                                                                <  2:4004442/93‑1 (MQ=255)

                         aaCCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCt                                                               >  2:1479426/1‑93 (MQ=255)

                         aaCCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCt                                                               <  1:1094672/93‑1 (MQ=255)

                          aCCGGTCAATCGCATTCTGGGGTATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTg                                                              <  2:3131976/93‑1 (MQ=255)

                          aCCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTg                                                              >  1:1079855/1‑93 (MQ=255)

                          aCCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTg                                                              >  1:2838050/1‑93 (MQ=255)

                           ccGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTgg                                                             <  1:878765/93‑1 (MQ=255)

                              gTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTg                                                          >  1:4187965/1‑93 (MQ=255)

                              gTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTg                                                          <  1:2329677/93‑1 (MQ=255)

                                cAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAt                                                        <  1:708262/93‑1 (MQ=255)

                                cAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAt                                                        <  1:333999/93‑1 (MQ=255)

                                 aaTCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAtt                                                       >  1:2983999/1‑93 (MQ=255)

                                 aaTCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAtt                                                       >  1:1477149/1‑93 (MQ=255)

                                  aTCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAttt                                                      <  1:2883312/93‑1 (MQ=255)

                                  aTCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAttt                                                      <  1:133364/93‑1 (MQ=255)

                                   tCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAtttt                                                     >  1:2198371/1‑93 (MQ=255)

                                   gcGCAATCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTGGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAtttt                                                     <  2:3867384/92‑1 (MQ=255)

                                       aTTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGAcc                                                 >  2:2060212/1‑93 (MQ=255)

                                       aTTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGAcc                                                 <  1:1949423/93‑1 (MQ=255)

                                       aTTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGAcc                                                 <  2:824531/93‑1 (MQ=255)

                                        ttCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCg                                                <  2:4031708/93‑1 (MQ=255)

                                             gggCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCa                                           >  2:849691/1‑93 (MQ=255)

                                             gggCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCa                                           <  2:2012614/93‑1 (MQ=255)

                                             gggCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCa                                           <  1:1724270/93‑1 (MQ=255)

                                             gggCATGTAAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCa                                           <  2:4483571/93‑1 (MQ=255)

                                               gCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAgg                                         >  2:4198942/1‑93 (MQ=255)

                                                cATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAggc                                        <  1:2593207/93‑1 (MQ=255)

                                                  tgtgAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAggcgg                                      >  2:3151480/1‑93 (MQ=255)

                                                    tgAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGTCGGCa                                    >  2:4522919/1‑93 (MQ=255)

                                                       aaaaTGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTc                                 <  2:3581087/93‑1 (MQ=255)

                                                       aaaaTGTCATTCTGATCGATCCCCAGGAGTATTTGCCGATTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTc                                 >  2:4632308/1‑93 (MQ=255)

                                                         aaTGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAg                               <  1:3978322/93‑1 (MQ=255)

                                                           tGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGga                             <  2:2165675/93‑1 (MQ=255)

                                                            gTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaa                            >  1:4609631/1‑93 (MQ=255)

                                                            gTCATTCTGATCGATCCCCAGGAGGATTTGCCGTTTGGCTGGCTGATGAACCACGCCTGGCTGATTTTGACTGACTCAGGCGGCATTCAgtca                            >  1:1627946/1‑90 (MQ=255)

                                                            gTCATTCTGATCGATACCCAGGAGTATTTGCCGTTTGTCGGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaa                            <  1:64973/93‑1 (MQ=255)

                                                             tCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaag                           >  1:4636904/1‑93 (MQ=255)

                                                             tCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaag                           <  2:1168841/93‑1 (MQ=255)

                                                              cATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaaga                          >  1:1796092/1‑93 (MQ=255)

                                                               aTTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaagaa                         <  1:544721/93‑1 (MQ=255)

                                                               aTTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaagaa                         <  2:1116570/93‑1 (MQ=255)

                                                                ttCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGaagaag                        >  1:1580223/1‑93 (MQ=255)

                                                                ttCTGATCGATCCCCAGAAGTATTTGCCGTTTGTCTGGCTGATGAACCGCGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGaagaag                        <  2:3192397/93‑1 (MQ=255)

                                                                 tCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAgc                       <  2:2423089/93‑1 (MQ=255)

                                                                    gatcgatcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGcc                    >  2:2551970/1‑93 (MQ=255)

                                                                    gatcgatcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGcc                    <  2:466958/93‑1 (MQ=255)

                                                                    gaacgaccCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGcc                    <  2:321781/90‑1 (MQ=255)

                                                                     atcgatcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCt                   <  1:3395983/93‑1 (MQ=255)

                                                                        gatcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCg                >  2:3783054/1‑93 (MQ=255)

                                                                        gatcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCg                >  1:1832303/1‑93 (MQ=255)

                                                                        gatcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCg                >  2:3505202/1‑93 (MQ=255)

                                                                         atcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGc               >  1:3162777/1‑93 (MQ=255)

                                                                         atcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCCCCTTCGc               <  2:3339848/93‑1 (MQ=255)

                                                                         atcCCCAGGAGTATTTGCCGTGTGTCTTGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGc               <  1:2653712/93‑1 (MQ=255)

                                                                          tcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCt              >  2:4093246/1‑93 (MQ=255)

                                                                           ccccAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTg             <  2:4149453/93‑1 (MQ=255)

                                                                           ccccAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTg             <  2:647362/93‑1 (MQ=255)

                                                                            cccAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTgg            <  2:545127/93‑1 (MQ=255)

                                                                              cAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTgggg          >  2:2093949/1‑93 (MQ=255)

                                                                               aGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGa         >  1:1998232/1‑93 (MQ=255)

                                                                                ggAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGaa        >  2:2961834/1‑93 (MQ=255)

                                                                                ggAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGaa        >  2:1085080/1‑93 (MQ=255)

                                                                                ggAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGaa        <  1:1293376/93‑1 (MQ=255)

                                                                                 gattATTTGGCGTTTGTCTGGCTGATGAACCACGACTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGaaa       <  1:386842/90‑1 (MQ=255)

                                                                                 gAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGaaa       >  1:2820970/1‑93 (MQ=255)

                                                                                  aGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAAc      <  2:2588857/93‑1 (MQ=255)

                                                                                     aTTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCtg   >  2:4312026/1‑93 (MQ=255)

                                                                                      tttGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCtgt  >  2:3712922/1‑93 (MQ=255)

                                                                                      tttGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCtgt  <  2:3311441/93‑1 (MQ=255)

GTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGT  >  NC_000913/3968477‑3968655

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