BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	16,925	A→G	R12R (CGT→CGC)	gef ←	Gef

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	16,925	0	A	G	92.4%	763.0 / 38.1	277	R12R (CGT→CGC)	gef	Gef
Reads supporting (aligned to +/- strand): ref base A (4/17); new base G (82/173); total (86/191)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 283 total aligned reads displayed.

CCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTG  >  NC_000913/16834‑17013

   gcgcCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGt                                                                                      <  1:5002946/93‑1 (MQ=255)

   gcgcCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGt                                                                                      >  2:7419425/1‑93 (MQ=255)

        gccACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTa                                                                                 <  1:273479/93‑1 (MQ=255)

        gccACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTa                                                                                 <  1:3867101/93‑1 (MQ=255)

         ccACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTaa                                                                                <  1:7199420/93‑1 (MQ=255)

         ccACTACGGCGGTGATACAGATGACGATCAGGGCAACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTaa                                                                                <  2:2860898/93‑1 (MQ=255)

          cACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAag                                                                               <  1:1038588/93‑1 (MQ=255)

           aCTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAaga                                                                              <  2:1379347/93‑1 (MQ=255)

            cTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAagag                                                                             <  2:6843172/93‑1 (MQ=255)

             tACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAgg                                                                            <  1:8793412/93‑1 (MQ=255)

             tACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAgg                                                                            >  1:178309/1‑93 (MQ=255)

             tACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAgg                                                                            <  1:6690134/93‑1 (MQ=255)

              aCGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGc                                                                           <  1:8531313/93‑1 (MQ=255)

               cggcggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCa                                                                          <  2:119621/93‑1 (MQ=255)

                ggcggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCAc                                                                         <  2:6068059/93‑1 (MQ=255)

                 gcggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCAct                                                                        >  1:519199/1‑93 (MQ=255)

                    gttaTACAGATGCGCATCCGGCCAACACTCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTTACCTTGCGGTCAGTAAGAGGCACTCTa                                                                     <  2:5540589/90‑1 (MQ=255)

                    gTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTa                                                                     <  2:7909336/93‑1 (MQ=255)

                      gATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACa                                                                   <  2:9619540/93‑1 (MQ=255)

                      gATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACa                                                                   <  1:691402/93‑1 (MQ=255)

                      gATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACa                                                                   <  1:6636546/93‑1 (MQ=255)

                      gATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACa                                                                   <  1:4209582/93‑1 (MQ=255)

                        tACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACAtg                                                                 >  1:8026315/1‑93 (MQ=255)

                        tACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACAtg                                                                 <  1:5575344/93‑1 (MQ=255)

                          cAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACAtgtg                                                               <  1:9424799/93‑1 (MQ=255)

                          cAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACAtgtg                                                               >  2:3437580/1‑93 (MQ=255)

                          cAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACAtgtg                                                               >  2:1823705/1‑93 (MQ=255)

                           aGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATgtgt                                                              >  2:6934061/1‑93 (MQ=255)

                             aTGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTc                                                            <  2:8640985/93‑1 (MQ=255)

                              tGACGATCAGGACGACAATCATCGCCTTATGCTGCTACATTGCTCTCTTCTCCGTGACCTTCCTGTCAGTAAGAGGCACTCTACATGTGTTCa                                                           <  1:57924/93‑1 (MQ=255)

                                aCGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGc                                                         <  1:4663865/93‑1 (MQ=255)

                                   aTCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCata                                                      >  2:5333921/1‑93 (MQ=255)

                                   aTCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCata                                                      <  2:6504849/93‑1 (MQ=255)

                                    tCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCatat                                                     <  2:7385441/93‑1 (MQ=255)

                                    tCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCatat                                                     <  2:5658940/93‑1 (MQ=255)

                                    tCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCatat                                                     <  2:9734097/93‑1 (MQ=255)

                                     cAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCAtata                                                    <  1:3208680/93‑1 (MQ=255)

                                     cAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCAtata                                                    <  2:3918051/93‑1 (MQ=255)

                                     cAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCAtata                                                    <  2:2945108/93‑1 (MQ=255)

                                      aGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATag                                                   <  1:9118074/93‑1 (MQ=255)

                                       gggCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATagg                                                  <  2:6231068/93‑1 (MQ=255)

                                       gggCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATagg                                                  >  2:7473741/1‑93 (MQ=255)

                                       gggCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATagg                                                  <  2:7714961/93‑1 (MQ=255)

                                       gggCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATagg                                                  <  2:77526/93‑1 (MQ=255)

                                        ggCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATagga                                                 <  2:5616667/93‑1 (MQ=255)

                                        ggCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATagga                                                 <  1:1133473/93‑1 (MQ=255)

                                        ggCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATagga                                                 >  2:7144864/1‑93 (MQ=255)

                                         gCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATaggag                                                <  1:1276200/93‑1 (MQ=255)

                                         gCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATaggag                                                <  1:1679436/93‑1 (MQ=255)

                                         gCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATaggag                                                <  2:3243713/93‑1 (MQ=255)

                                          cGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATaggagg                                               <  2:9399861/93‑1 (MQ=255)

                                            aCAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGcc                                             <  2:333362/93‑1 (MQ=255)

                                             cAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGctt                                            <  2:5625252/93‑3 (MQ=255)

                                             cAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCt                                            <  1:10006945/93‑1 (MQ=255)

                                             cAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCt                                            <  1:6068924/93‑1 (MQ=255)

                                              aaTCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTc                                           <  1:4704156/93‑1 (MQ=255)

                                              aaTCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTc                                           >  1:849732/1‑93 (MQ=255)

                                               atcatcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCg                                          <  2:2834813/93‑1 (MQ=255)

                                               atcatcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCg                                          >  2:1321809/1‑93 (MQ=255)

                                                   tcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGtt                                      >  2:8162348/1‑93 (MQ=255)

                                                     gCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGttca                                    >  1:7863955/1‑91 (MQ=255)

                                                       cTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATg                                  <  1:29569/93‑1 (MQ=255)

                                                       cTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATg                                  <  2:805834/93‑1 (MQ=255)

                                                        ttATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGGCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATgg                                 <  2:6596227/93‑1 (MQ=255)

                                                         tATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGt                                >  1:1549746/1‑93 (MQ=255)

                                                           tgctgcTTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTaa                              <  1:8311727/93‑1 (MQ=255)

                                                            gctgctTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTaaa                             <  2:4952084/93‑1 (MQ=255)

                                                                ctTCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAatat                         <  1:3511907/93‑1 (MQ=255)

                                                                ctTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAatat                         <  1:6386889/93‑1 (MQ=255)

                                                                 ttCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATc                        >  2:1559852/1‑93 (MQ=255)

                                                                  tCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCa                       <  2:6460143/93‑1 (MQ=255)

                                                                  tCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCa                       >  2:2282428/1‑93 (MQ=255)

                                                                  tCATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGGAAAATATCa                       >  2:5875945/1‑93 (MQ=255)

                                                                  tCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCa                       <  2:1522605/93‑1 (MQ=255)

                                                                   cATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCAc                      <  1:4781447/93‑1 (MQ=255)

                                                                   cATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCAc                      <  1:539093/93‑1 (MQ=255)

                                                                   cATTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCAc                      <  2:5160916/93‑1 (MQ=255)

                                                                    aTTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACt                     <  1:7693663/93‑1 (MQ=255)

                                                                    aTTGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACt                     >  1:9627/1‑93 (MQ=255)

                                                                     ttGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTc                    <  1:2164764/93‑1 (MQ=255)

                                                                     ttGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGATTGATGGTAAAATATCACTc                    <  2:3781114/93‑1 (MQ=255)

                                                                      tGCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCg                   <  2:4672656/93‑1 (MQ=255)

                                                                       gCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCgg                  <  1:5914966/93‑1 (MQ=255)

                                                                       gCTCTCTTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCgg                  <  2:8383017/93‑1 (MQ=255)

                                                                       gCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCgg                  <  1:4295306/93‑1 (MQ=255)

                                                                        ctctctTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCggg                 <  1:1901072/93‑1 (MQ=255)

                                                                        ctctctTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCggg                 <  2:6296687/93‑1 (MQ=255)

                                                                          ctctTCTCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGc               >  1:4065684/1‑93 (MQ=255)

                                                                           tcttctCCTTTACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCt              <  2:8548573/93‑1 (MQ=255)

                                                                           tcttctCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCt              <  1:3814907/93‑1 (MQ=255)

                                                                           tcttctCCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCt              >  1:259804/1‑93 (MQ=255)

                                                                             ttctCCTTGACCTTGCGGTCAGCAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCttt            >  2:902481/1‑93 (MQ=255)

                                                                              tctcCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCtttt           <  2:5488414/93‑1 (MQ=255)

                                                                               ctcCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtc          >  2:463456/1‑93 (MQ=255)

                                                                               ctcCTTGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtc          <  1:2381599/93‑1 (MQ=255)

                                                                                tcCTTGACCTTGCGGTCAGTAAGAGGCCCTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtct         <  1:3411847/93‑1 (MQ=255)

                                                                                  cTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTctct       <  1:5129850/93‑1 (MQ=255)

                                                                                   ttGACCTTGCGGGCAGTAAGAGGCACGCTACATGTGTTCAGCATATAGGAGGGCTCGGGTTGATGGTAAAATATCACTCGCGGCTTTTTTCTa      >  1:9851290/1‑93 (MQ=255)

                                                                                    tGACCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTAt     >  1:6543696/1‑93 (MQ=255)

                                                                                     taCCTTGCGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAAGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATc    <  1:1210283/92‑1 (MQ=255)

CCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTG  >  NC_000913/16834‑17013

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