Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,222,061 A→G K596K (AAA→AAG)  b3076 → evolved beta‑D‑galactosidase alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,222,0610AG99.2% 1268.6 / ‑3.9 371K596K (AAA→AAG) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (204/164);  total (205/166)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 384 total aligned reads displayed.

CTACGGCGACTATCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGA  >  NC_000913/3221971‑3222153
                                                                                          |                                                                                            
ctacggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCa                                                                                            >  1:5368541/1‑93 (MQ=255)
ctacggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCa                                                                                            >  1:8253029/1‑93 (MQ=255)
ctacggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCa                                                                                            <  2:5691988/93‑1 (MQ=255)
  acggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAgg                                                                                          <  2:7200272/93‑1 (MQ=255)
   cggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGt                                                                                         <  1:984616/93‑1 (MQ=255)
      cgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTAt                                                                                      >  2:4850755/1‑93 (MQ=255)
      cgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTAt                                                                                      >  1:7737751/1‑93 (MQ=255)
       cactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGCCGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATc                                                                                     <  1:3175642/92‑1 (MQ=255)
        actaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATcg                                                                                    >  2:7713658/1‑93 (MQ=255)
        actaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATcg                                                                                    >  1:2466509/1‑93 (MQ=255)
        actaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCAGGACCGGGCCTGAAAGAGAACAAGCAGGTTATcg                                                                                    >  2:8692560/1‑93 (MQ=255)
          taTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATcgcg                                                                                  >  1:4982764/1‑93 (MQ=255)
           aTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCgcgc                                                                                 >  1:7961896/1‑93 (MQ=255)
            tCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGcc                                                                                >  1:7584019/1‑93 (MQ=255)
            tCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGcc                                                                                <  1:1873294/93‑1 (MQ=255)
             cccAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCg                                                                               <  2:2906337/93‑1 (MQ=255)
             cccAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCg                                                                               >  2:1467202/1‑93 (MQ=255)
              ccAACAACTTTAACTTCTGTTTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCTGGCCTGAAATTGTACAAGCCGGTTATCGCGCCgt                                                                              >  2:8072491/1‑92 (MQ=255)
               caacagCTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGt                                                                             <  2:3577592/93‑1 (MQ=255)
               caacaaCTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGt                                                                             >  1:3266644/1‑93 (MQ=255)
                aacaacTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTa                                                                            >  1:2645999/1‑93 (MQ=255)
                    acTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTaaaaa                                                                        >  2:1053043/1‑93 (MQ=255)
                       ataACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATcc                                                                     <  1:716553/93‑1 (MQ=255)
                       ataACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATcc                                                                     >  1:7928681/1‑93 (MQ=255)
                        taACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCa                                                                    >  2:3432591/1‑93 (MQ=255)
                        taACTTCTGTCTTGATGGGTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCa                                                                    >  1:3568650/1‑93 (MQ=255)
                         aaCTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAc                                                                   >  2:8776935/1‑93 (MQ=255)
                          aCTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcg                                                                  >  1:5408238/1‑93 (MQ=255)
                          aCTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcg                                                                  >  1:6265296/1‑93 (MQ=255)
                           cTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcgc                                                                 <  2:1579483/93‑1 (MQ=255)
                             tCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcgcgc                                                               >  1:5265094/1‑93 (MQ=255)
                             tCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcgcgc                                                               <  1:487572/93‑1 (MQ=255)
                                gTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCggg                                                            >  2:3664553/1‑93 (MQ=255)
                                gTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCggg                                                            <  1:1113981/93‑1 (MQ=255)
                                gTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCggg                                                            >  1:587359/1‑93 (MQ=255)
                                gTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCggg                                                            <  1:4219713/93‑1 (MQ=255)
                                 tCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGa                                                           >  1:5504526/1‑93 (MQ=255)
                                  cTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGAt                                                          <  1:3107423/93‑1 (MQ=255)
                                  cTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGAt                                                          <  1:1414904/93‑1 (MQ=255)
                                  cTTGATGGTGTGATCTATTCCGATCAGACGCCTGGACCGGGCCTGAAAGAGTACAAGCAGGTTAGCGCGCCGGTAAAAATCCACGCGCGGGAt                                                          <  2:5813705/93‑1 (MQ=255)
                                   ttGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATc                                                         >  2:4544007/1‑93 (MQ=255)
                                   ttGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATc                                                         <  2:2460368/93‑1 (MQ=255)
                                   ttGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATc                                                         >  2:189820/1‑93 (MQ=255)
                                   ttGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATc                                                         <  1:4714928/93‑1 (MQ=255)
                                   ttGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATc                                                         <  2:6462759/93‑1 (MQ=255)
                                    tGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCt                                                        >  1:614774/1‑93 (MQ=255)
                                    tGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCt                                                        <  1:4332339/93‑1 (MQ=255)
                                     gATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTg                                                       <  1:5482044/93‑1 (MQ=255)
                                      aTGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGa                                                      <  1:4591394/93‑1 (MQ=255)
                                      aTGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGa                                                      >  1:2130476/1‑93 (MQ=255)
                                       tGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGAc                                                     <  1:5961250/93‑1 (MQ=255)
                                        ggTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACt                                                    >  1:811322/1‑93 (MQ=255)
                                        ggTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACt                                                    <  2:2640558/93‑1 (MQ=255)
                                         gTTTGATCTATACCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCTACGCGCGGGCTCTGACTc                                                   <  2:5052002/93‑1 (MQ=255)
                                            tGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTcgcg                                                >  1:6783566/1‑93 (MQ=255)
                                              aTCTATTCCGATCAGACGCCGGCACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCgcggc                                              <  1:6427258/93‑1 (MQ=255)
                                               tCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCgcggcg                                             >  1:3343903/1‑93 (MQ=255)
                                                cTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGa                                            <  1:8637350/93‑1 (MQ=255)
                                                 tATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAg                                           >  2:7006158/1‑93 (MQ=255)
                                                 tATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAg                                           >  2:772958/1‑93 (MQ=255)
                                                 tATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAg                                           >  1:2812617/1‑93 (MQ=255)
                                                  aTTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGt                                          <  2:7649329/93‑1 (MQ=255)
                                                   ttCCGATCAGACGCCGGGACCGGGCCTGATAGAGTACAAGCAGGTTATCGCTCCGGTAAAAATCCACGCGCGGGATCCGACTCGCGGCGAGtt                                         >  1:7839837/1‑93 (MQ=255)
                                                   ttCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGtt                                         >  2:4797599/1‑93 (MQ=255)
                                                    tCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTg                                        >  2:88781/1‑93 (MQ=255)
                                                        aTCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAg                                    <  2:6367390/93‑1 (MQ=255)
                                                         tCAGCCGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGt                                   <  2:6661236/93‑1 (MQ=255)
                                                         tCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGt                                   >  1:8498196/1‑93 (MQ=255)
                                                           aGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCg                                 <  2:1164375/93‑1 (MQ=255)
                                                            gACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGa                                >  1:7253059/1‑93 (MQ=255)
                                                            gACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGa                                >  2:8114439/1‑93 (MQ=255)
                                                             aCGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGaa                               >  2:7454651/1‑93 (MQ=255)
                                                                ccGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGANAAt                            >  2:2001106/1‑93 (MQ=255)
                                                                  gggACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATaa                          >  2:215953/1‑93 (MQ=255)
                                                                   ggACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATaaa                         <  1:5780678/93‑1 (MQ=255)
                                                                    gACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAAc                        >  2:5500297/1‑93 (MQ=255)
                                                                    gACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAAc                        >  1:1042534/1‑93 (MQ=255)
                                                                     aCCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACt                       <  1:1222986/93‑1 (MQ=255)
                                                                      ccGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtg                      <  2:166834/93‑1 (MQ=255)
                                                                        gggCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtgtg                    >  2:2213197/1‑93 (MQ=255)
                                                                        gggCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtgtg                    >  1:174167/1‑93 (MQ=255)
                                                                        gggCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtgtg                    >  1:2095730/1‑93 (MQ=255)
                                                                         ggCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTgg                   >  2:35962/1‑93 (MQ=255)
                                                                         ggCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTgg                   <  2:4371272/93‑1 (MQ=255)
                                                                          gCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGt                  >  2:5942932/1‑93 (MQ=255)
                                                                              gAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTAc              <  1:3691166/93‑1 (MQ=255)
                                                                              gAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTAc              >  1:7376858/1‑93 (MQ=255)
                                                                               aaaGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTAcc             <  2:4201969/93‑1 (MQ=255)
                                                                                aaGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCa            >  2:6312804/1‑93 (MQ=255)
                                                                                 agagTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCAc           <  2:3765368/93‑1 (MQ=255)
                                                                                 agagTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCAc           <  1:3069878/93‑1 (MQ=255)
                                                                                  gagTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACg          <  2:535802/93‑1 (MQ=255)
                                                                                   agTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGc         <  2:72690/93‑1 (MQ=255)
                                                                                   agTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGc         <  1:2234351/93‑1 (MQ=255)
                                                                                   agTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGc         <  2:2766328/93‑1 (MQ=255)
                                                                                    gTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCt        <  2:6808801/93‑1 (MQ=255)
                                                                                    gTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCt        <  2:3047725/93‑1 (MQ=255)
                                                                                     tACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCtt       <  2:8272187/93‑1 (MQ=255)
                                                                                       cAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTtga     <  2:8469387/93‑1 (MQ=255)
                                                                                       cAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTtga     >  2:3172570/1‑93 (MQ=255)
                                                                                          |                                                                                            
CTACGGCGACTATCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGATGA  >  NC_000913/3221971‑3222153

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