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.6% 848.2 / ‑4.7 244K596K (AAA→AAGb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (138/105);  total (138/106)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

Only 100 of 247 total aligned reads displayed.

ACTACGGCGACTATCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTT  >  NC_000913/3221970‑3222148
                                                                                           |                                                                                       
actacggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAAc                                                                                        <  1:1374590/93‑1 (MQ=255)
  tacggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAg                                                                                      >  1:124172/1‑93 (MQ=255)
   acggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAgg                                                                                     >  1:576573/1‑93 (MQ=255)
   acggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAgg                                                                                     >  1:5934891/1‑93 (MQ=255)
     ggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGtt                                                                                   <  2:3339075/93‑1 (MQ=255)
     ggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGtt                                                                                   >  2:2469323/1‑93 (MQ=255)
     ggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGtt                                                                                   >  1:299205/1‑93 (MQ=255)
      gcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTa                                                                                  >  2:879733/1‑93 (MQ=255)
       cgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTAt                                                                                 <  1:5816045/93‑1 (MQ=255)
        gactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATc                                                                                <  1:6263037/93‑1 (MQ=255)
        gactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATc                                                                                <  2:2812606/93‑1 (MQ=255)
         actaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATcg                                                                               <  1:2604370/93‑1 (MQ=255)
            aTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCgcgc                                                                            >  2:5640093/1‑93 (MQ=255)
            aTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGCTATCgcgc                                                                            >  1:4110/1‑93 (MQ=255)
              cccAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCg                                                                          <  2:1259056/93‑1 (MQ=255)
               ccAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCgg                                                                         >  1:1910200/1‑93 (MQ=255)
               ccAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCgg                                                                         >  2:6041493/1‑93 (MQ=255)
                 aacaacTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTa                                                                       >  2:4655370/1‑93 (MQ=255)
                  acaacTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTaa                                                                      <  2:1873633/93‑1 (MQ=255)
                   caacTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTaaa                                                                     >  2:3393806/1‑93 (MQ=255)
                     acTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTataaa                                                                   >  2:6125493/1‑93 (MQ=255)
                       tataACTTCTGTCTTGATGGGTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATc                                                                 >  1:6104959/1‑93 (MQ=255)
                        ataACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATcc                                                                >  1:956720/1‑93 (MQ=255)
                         taACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCa                                                               >  2:4151150/1‑93 (MQ=255)
                          aaCTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAc                                                              >  2:5451256/1‑93 (MQ=255)
                           aCTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcg                                                             <  1:710168/93‑1 (MQ=255)
                           aCTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcg                                                             >  2:6070999/1‑93 (MQ=255)
                            cTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcgc                                                            >  1:3852825/1‑93 (MQ=255)
                            cTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcgc                                                            <  1:2152676/93‑1 (MQ=255)
                               cTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcgcgcg                                                         <  1:1723407/93‑1 (MQ=255)
                               cTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcgcgcg                                                         >  2:4469994/1‑93 (MQ=255)
                                tGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCgg                                                        >  2:2759431/1‑93 (MQ=255)
                                 gTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCggg                                                       >  1:4645464/1‑93 (MQ=255)
                                 gTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCggg                                                       <  1:1718095/93‑1 (MQ=255)
                                 gTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCggg                                                       >  1:2222552/1‑93 (MQ=255)
                                 gTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATACACGCGCggg                                                       >  1:4383990/1‑93 (MQ=255)
                                 gTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGGTATCGCGCCGGTAAAAATCCACGCGCggg                                                       >  2:5824733/1‑93 (MQ=255)
                                    ttGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATc                                                    <  2:4792669/93‑1 (MQ=255)
                                    ttGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATc                                                    <  1:845100/93‑1 (MQ=255)
                                     tGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCt                                                   <  2:1644458/93‑1 (MQ=255)
                                      gATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTg                                                  >  2:3846299/1‑93 (MQ=255)
                                       aTGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGa                                                 >  1:3452695/1‑93 (MQ=255)
                                       aTGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGa                                                 <  2:616291/93‑1 (MQ=255)
                                        tGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGAc                                                >  2:5362171/1‑93 (MQ=255)
                                         ggTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACt                                               <  2:3690444/93‑1 (MQ=255)
                                          gTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTc                                              <  2:2131196/93‑1 (MQ=255)
                                           tttGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTcg                                             >  1:2008062/1‑93 (MQ=255)
                                            ttGATCTAGTCCGATCAGACGCCGGCACCGGCCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTcgc                                            <  1:5510750/93‑1 (MQ=255)
                                              gATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCgcgg                                          >  2:5510068/1‑93 (MQ=255)
                                              gATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCgcgg                                          >  2:1119382/1‑93 (MQ=255)
                                                 cTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGa                                       <  2:285070/93‑1 (MQ=255)
                                                 cTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGa                                       <  1:6235534/93‑1 (MQ=255)
                                                 cTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGa                                       <  2:1819854/93‑1 (MQ=255)
                                                 atATTCCGAGCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGa                                       <  1:5326824/92‑1 (MQ=255)
                                                  tATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCgcggcgcg                                      >  1:4957222/1‑91 (MQ=255)
                                                   aTTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGt                                     <  1:1149125/93‑1 (MQ=255)
                                                    ttCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGtt                                    <  1:1577017/93‑1 (MQ=255)
                                                    ttCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGtt                                    >  2:5766690/1‑93 (MQ=255)
                                                     tCCGATCAGACGCCGGGACCGGCCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTg                                   <  2:2014653/93‑1 (MQ=255)
                                                      ccGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGa                                  >  1:3380056/1‑93 (MQ=255)
                                                        gATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGaaa                                <  1:1853890/93‑1 (MQ=255)
                                                        gATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGaaa                                <  1:2931260/93‑1 (MQ=255)
                                                         aTCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAg                               >  2:6255672/1‑93 (MQ=255)
                                                         aTCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAg                               >  1:1808413/1‑93 (MQ=255)
                                                         aTCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAg                               >  1:3244727/1‑93 (MQ=255)
                                                         aTCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAg                               >  2:2912575/1‑93 (MQ=255)
                                                         aTCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGCTGAAAg                               >  1:4998321/1‑93 (MQ=255)
                                                           cAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAACATCCACGCGCGGGATCTGACTCGCCGCGAGTTGAAAGTc                             >  1:6011807/1‑93 (MQ=255)
                                                             gACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGa                           >  1:3197301/1‑93 (MQ=255)
                                                             gACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGa                           >  1:3572804/1‑93 (MQ=255)
                                                              aCGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGaa                          >  1:5140268/1‑93 (MQ=255)
                                                               cGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGaaa                         >  1:906807/1‑93 (MQ=255)
                                                                gCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGaaaa                        >  2:839246/1‑93 (MQ=255)
                                                                 ccGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAAt                       >  2:1686929/1‑93 (MQ=255)
                                                                 ccGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAAt                       <  1:3787485/93‑1 (MQ=255)
                                                                 ccGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAAt                       >  1:6357227/1‑93 (MQ=255)
                                                                  cGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGATGAAAGTCGAAAATa                      >  1:2376947/1‑93 (MQ=255)
                                                                   gggACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATaa                     <  2:1046657/93‑1 (MQ=255)
                                                                   gggACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATaa                     >  1:3164843/1‑93 (MQ=255)
                                                                    ggACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATaaa                    <  2:3548692/93‑1 (MQ=255)
                                                                    ggACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATaaa                    >  1:4882764/1‑93 (MQ=255)
                                                                    ggACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATaaa                    >  2:2335246/1‑93 (MQ=255)
                                                                    ggACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGACGAAAATaaa                    >  1:4477317/1‑93 (MQ=255)
                                                                      aCCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACt                  >  1:2054199/1‑93 (MQ=255)
                                                                      aCCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACt                  >  2:5374226/1‑93 (MQ=255)
                                                                       ccGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtg                 >  2:3729957/1‑93 (MQ=255)
                                                                       ccGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtg                 <  2:5115049/93‑1 (MQ=255)
                                                                        cGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtgt                >  2:3136851/1‑93 (MQ=255)
                                                                        cGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtgt                <  2:1946536/93‑1 (MQ=255)
                                                                         gggCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtgtg               >  2:6207395/1‑93 (MQ=255)
                                                                         gggCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtgtg               >  2:1496048/1‑93 (MQ=255)
                                                                            ccTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGtt            <  1:1044213/93‑1 (MQ=255)
                                                                             cTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGttt           <  1:4147720/93‑1 (MQ=255)
                                                                               gAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTAc         <  2:3159670/93‑1 (MQ=255)
                                                                               gAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTAc         <  2:501208/93‑1 (MQ=255)
                                                                                aaaGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTAcc        >  1:4287405/1‑93 (MQ=255)
                                                                                   gagTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACg     >  2:5382373/1‑93 (MQ=255)
                                                                                    agTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGc    <  1:5843815/93‑1 (MQ=255)
                                                                                    agTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGc    <  2:5617327/93‑1 (MQ=255)
                                                                                      tACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGGGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCtt  <  1:5609833/93‑1 (MQ=255)
                                                                                           |                                                                                       
ACTACGGCGACTATCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTT  >  NC_000913/3221970‑3222148

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