BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,506,541	G→A	D302D (GAC→GAT)	glk ←	glucokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,506,541	0	G	A	82.9%	774.4 / 110.1	293	D302D (GAC→GAT)	glk	glucokinase
Reads supporting (aligned to +/- strand): ref base G (20/30); new base A (100/143); total (120/173)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

Only 100 of 294 total aligned reads displayed.

AGAGTTACCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAA  >  NC_000913/2506451‑2506633

agagTTACCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATc                                                                                            <  2:3987300/93‑1 (MQ=255)

 gagTTACCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCa                                                                                           >  1:6406533/1‑93 (MQ=255)

 gagTTACCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCa                                                                                           <  1:6007145/93‑1 (MQ=255)

 gagTTACCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCa                                                                                           <  2:717755/93‑1 (MQ=255)

 gagTTACCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCa                                                                                           <  1:283966/93‑1 (MQ=255)

   gTTACCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATg                                                                                         >  2:5297121/1‑93 (MQ=255)

   gTTACCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATg                                                                                         >  2:4683639/1‑93 (MQ=255)

     tACCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGa                                                                                       >  2:5440223/1‑93 (MQ=255)

     tACCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGa                                                                                       <  2:450776/93‑1 (MQ=255)

     tACCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGa                                                                                       <  1:3865641/93‑1 (MQ=255)

      aCCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGAc                                                                                      <  2:362948/93‑1 (MQ=255)

       ccTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACg                                                                                     >  2:3776144/1‑93 (MQ=255)

        cTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGa                                                                                    <  2:2729640/93‑1 (MQ=255)

        cTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGa                                                                                    <  1:1550933/93‑1 (MQ=255)

        cTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGa                                                                                    <  2:5837927/93‑1 (MQ=255)

         tCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGAt                                                                                   <  2:2799451/93‑1 (MQ=255)

         tCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGAt                                                                                   >  1:4293013/1‑93 (MQ=255)

          cccGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATg                                                                                  >  1:6702638/1‑93 (MQ=255)

           ccGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATga                                                                                 >  1:4564437/1‑93 (MQ=255)

           ccGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATga                                                                                 >  1:6623226/1‑93 (MQ=255)

             gATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATgaga                                                                               <  1:6521866/93‑1 (MQ=255)

             gATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATgaga                                                                               <  2:3516773/93‑1 (MQ=255)

             gATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATgaga                                                                               >  2:2219163/1‑93 (MQ=255)

              atatAAAAGTAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGAt                                                                              >  2:4025551/1‑93 (MQ=255)

              atatAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGTACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGAt                                                                              <  2:381377/93‑1 (MQ=255)

              atatAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGAt                                                                              <  1:1974994/93‑1 (MQ=255)

               tataAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATa                                                                             <  2:3368197/93‑1 (MQ=255)

                ataAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATac                                                                            >  1:7279754/1‑93 (MQ=255)

                ataAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATac                                                                            <  2:7105563/93‑1 (MQ=255)

                 taAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATaca                                                                           >  2:217068/1‑93 (MQ=255)

                  aaaaGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATacac                                                                          <  2:1660594/93‑1 (MQ=255)

                  aaaaGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATacac                                                                          <  2:2450615/93‑1 (MQ=255)

                  aaaaGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATacac                                                                          >  1:2545558/1‑93 (MQ=255)

                   aaaGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACAcc                                                                         <  2:66498/93‑1 (MQ=255)

                   aaaGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACAcc                                                                         <  1:1013070/93‑1 (MQ=255)

                   aaaGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACAcc                                                                         <  2:6215236/93‑1 (MQ=255)

                    aaggaaggATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGTACGATGAGATACACCg                                                                        >  1:222994/1‑93 (MQ=255)

                    aaggaaggATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCg                                                                        >  1:5605655/1‑93 (MQ=255)

                     aggaaggaTTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCgg                                                                       >  1:4484161/1‑93 (MQ=255)

                     aggaaggaTTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCgg                                                                       <  1:2770101/93‑1 (MQ=255)

                      ggaaggaTTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGa                                                                      >  1:1215759/1‑93 (MQ=255)

                      ggaaggaTTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGa                                                                      >  2:2964565/1‑93 (MQ=255)

                       gaaggaTTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGaa                                                                     <  2:4260521/93‑1 (MQ=255)

                         aggaTTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAata                                                                   >  2:5953411/1‑93 (MQ=255)

                           gaTTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATc                                                                 >  1:7477907/1‑93 (MQ=255)

                           gaTTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATc                                                                 <  1:2395003/93‑1 (MQ=255)

                             tttACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCAt                                                               <  1:5684538/93‑1 (MQ=255)

                             tttACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCAt                                                               <  2:6795227/93‑1 (MQ=255)

                             tttACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCAt                                                               >  1:3972738/1‑93 (MQ=255)

                               tACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATgg                                                             <  2:4236791/93‑1 (MQ=255)

                               tACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATgg                                                             <  2:6976773/93‑1 (MQ=255)

                                aCAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGa                                                            >  2:962982/1‑93 (MQ=255)

                                 cAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGAc                                                           >  1:1466747/1‑93 (MQ=255)

                                   gAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACat                                                         <  1:5812666/93‑1 (MQ=255)

                                   gAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACat                                                         >  1:6858384/1‑93 (MQ=255)

                                    aaTGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACata                                                        <  2:3904700/93‑1 (MQ=255)

                                        tgACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCt                                                    <  1:1614537/93‑1 (MQ=255)

                                            cTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTaa                                                <  2:7514768/93‑1 (MQ=255)

                                             tAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTaaa                                               <  2:5646266/93‑1 (MQ=255)

                                             tAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTaaa                                               <  1:6242582/93‑1 (MQ=255)

                                              aaGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAg                                              <  2:1537194/93‑1 (MQ=255)

                                              aaGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGACTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAg                                              <  2:859373/93‑1 (MQ=255)

                                               aGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAgc                                             <  1:4015864/93‑1 (MQ=255)

                                               aGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAgc                                             <  1:4835079/93‑1 (MQ=255)

                                                 gTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAgcgc                                           <  2:2631042/93‑1 (MQ=255)

                                                 gTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAgcgc                                           <  1:5342941/93‑1 (MQ=255)

                                                  tCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGcc                                          >  2:3101362/1‑93 (MQ=255)

                                                  tCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGcc                                          >  1:1803587/1‑93 (MQ=255)

                                                  tCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGcc                                          <  2:4894688/93‑1 (MQ=255)

                                                   cTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGccc                                         <  1:472275/93‑1 (MQ=255)

                                                    tGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCt                                        >  2:5015486/1‑93 (MQ=255)

                                                    tGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCt                                        >  2:6851196/1‑93 (MQ=255)

                                                     ggCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCtt                                       <  1:342257/93‑1 (MQ=255)

                                                       cGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTa                                     <  1:3400560/93‑1 (MQ=255)

                                                        gTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTAt                                    >  2:4091065/1‑93 (MQ=255)

                                                         tAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATc                                   >  1:3010847/1‑93 (MQ=255)

                                                         tAAATGTGCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATc                                   <  2:1952139/93‑1 (MQ=255)

                                                              gtgCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCaa                              <  1:1389234/93‑1 (MQ=255)

                                                               tgCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCaaa                             <  2:6045557/93‑1 (MQ=255)

                                                                gCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAAt                            <  2:1214508/93‑1 (MQ=255)

                                                                gCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAAt                            >  1:6068467/1‑93 (MQ=255)

                                                                gCACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAAt                            <  2:6464909/93‑1 (MQ=255)

                                                                 cACCGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATg                           <  1:5544789/93‑1 (MQ=255)

                                                                   ccGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCg                         >  2:4806294/1‑93 (MQ=255)

                                                                    cGGAACCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCgg                        <  1:972900/93‑1 (MQ=255)

                                                                      gAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCa                      >  2:6384398/1‑93 (MQ=255)

                                                                        aCCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACg                    <  2:6014558/93‑1 (MQ=255)

                                                                        aCCGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACg                    <  1:3770976/93‑1 (MQ=255)

                                                                          cGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGa                  <  2:1598621/93‑1 (MQ=255)

                                                                          cGAGAAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGa                  <  1:3323691/93‑1 (MQ=255)

                                                                           gagaAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGaa                 >  2:6268828/1‑93 (MQ=255)

                                                                            agaAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGaaa                >  1:3473862/1‑93 (MQ=255)

                                                                            agaAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCCAATGCGGCACGGaaa                >  1:6188102/1‑93 (MQ=255)

                                                                             gaAGGCCCGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAAc               <  1:2619786/93‑1 (MQ=255)

                                                                                  cccGGATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAg          >  1:4933890/1‑93 (MQ=255)

                                                                                    cGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGc        <  1:3641930/93‑1 (MQ=255)

                                                                                      gATTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCtt      >  2:7462879/1‑93 (MQ=255)

                                                                                       aTTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCttt     >  2:144074/1‑93 (MQ=255)

                                                                                       aTTATCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCttt     >  2:5623129/1‑93 (MQ=255)

                                                                                         tGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTga   <  2:7240434/93‑1 (MQ=255)

AGAGTTACCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAATGCGGCACGGAAACCGGAGGCTTTGAA  >  NC_000913/2506451‑2506633

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