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,0610AG88.3% 817.4 / 86.0 317K596K (AAA→AAGb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (16/21);  new base G (141/139);  total (157/160)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

CTACGGCGACTATCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGAT  >  NC_000913/3221971‑3222151
                                                                                          |                                                                                          
  anggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAgn                                                                                        <  2:5604467/93‑2 (MQ=255)
   cggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGt                                                                                       >  2:5069650/1‑93 (MQ=255)
    ggcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGtt                                                                                      <  1:3167346/93‑1 (MQ=255)
     gcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTa                                                                                     >  2:3925699/1‑93 (MQ=255)
     gcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTa                                                                                     <  2:211062/93‑1 (MQ=255)
     gcgactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTa                                                                                     <  1:5939678/93‑1 (MQ=255)
       gactaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATc                                                                                   <  2:910476/93‑1 (MQ=255)
        actaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATcg                                                                                  <  2:3663689/93‑1 (MQ=255)
         ctaTCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATcgc                                                                                 <  2:5795554/93‑1 (MQ=255)
              ccAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCgg                                                                            <  2:1188455/93‑1 (MQ=255)
               caacaaCTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGt                                                                           >  1:1407056/1‑93 (MQ=255)
                aacaacTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTa                                                                          <  1:7398820/93‑1 (MQ=255)
                aacaacTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTa                                                                          >  1:4945585/1‑93 (MQ=255)
                aacaacTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTa                                                                          >  2:5889740/1‑93 (MQ=255)
                aacaacTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTa                                                                          >  2:1947091/1‑93 (MQ=255)
                 acaacTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTaa                                                                         >  1:5257207/1‑93 (MQ=255)
                 acaacTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTaa                                                                         >  1:5521051/1‑93 (MQ=255)
                    acTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTaaaaa                                                                      >  2:2738496/1‑93 (MQ=255)
                      tataACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATc                                                                    >  1:5819685/1‑93 (MQ=255)
                      tataACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATc                                                                    >  1:4709631/1‑93 (MQ=255)
                      tataACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATc                                                                    >  2:705771/1‑93 (MQ=255)
                      tataACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATc                                                                    <  2:6252775/93‑1 (MQ=255)
                       ataACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATcc                                                                   >  1:1724782/1‑93 (MQ=255)
                         aaCTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAc                                                                 >  1:5439606/1‑93 (MQ=255)
                          aCTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcg                                                                <  1:6394449/93‑1 (MQ=255)
                          aCTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCAcg                                                                >  2:1635099/1‑93 (MQ=255)
                           cTTCTGTCTTGATGGTTTGATCTATTCCGATCGGACGTCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCAcgc                                                               <  1:7014523/93‑1 (MQ=255)
                            ttCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcgcg                                                              >  1:2171834/1‑93 (MQ=255)
                              cTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCAcgcgcg                                                            <  2:1348512/93‑1 (MQ=255)
                               tGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCgg                                                           <  1:2992738/93‑1 (MQ=255)
                                  cTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGAt                                                        >  1:3434887/1‑93 (MQ=255)
                                   ttGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATc                                                       <  1:2561370/93‑1 (MQ=255)
                                   ttGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATc                                                       >  2:6830383/1‑93 (MQ=255)
                                   ttGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATc                                                       <  2:605269/93‑1 (MQ=255)
                                   ttGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATc                                                       >  1:3994348/1‑93 (MQ=255)
                                    tGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCt                                                      <  2:2950937/93‑1 (MQ=255)
                                     gATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTg                                                     >  2:3262971/1‑93 (MQ=255)
                                     gATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTg                                                     >  2:2613443/1‑93 (MQ=255)
                                      aTGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCCCCGGTAAAAATCCACGCGCGGGATCTGa                                                    <  2:1054821/93‑1 (MQ=255)
                                      aTGGTTTGATCTATTCCGACCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGa                                                    <  2:5002756/93‑1 (MQ=255)
                                       tGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGAc                                                   >  1:5634644/1‑93 (MQ=255)
                                        ggTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACt                                                  <  2:2832100/93‑1 (MQ=255)
                                        ggTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACt                                                  >  1:4775325/1‑93 (MQ=255)
                                         gTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTc                                                 >  1:3916201/1‑93 (MQ=255)
                                         gTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTc                                                 >  1:434715/1‑93 (MQ=255)
                                          tttGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTcg                                                >  1:3771789/1‑93 (MQ=255)
                                          tttGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTcg                                                >  1:1039219/1‑93 (MQ=255)
                                           ttGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTcgc                                               <  2:7264279/93‑1 (MQ=255)
                                           ttGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTcgc                                               <  1:1185331/93‑1 (MQ=255)
                                           ttGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTcgc                                               >  1:7179479/1‑93 (MQ=255)
                                            tGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTcgcg                                              <  1:3558337/93‑1 (MQ=255)
                                              aTCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCgcggc                                            >  2:2300636/1‑93 (MQ=255)
                                                cTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGa                                          <  1:3108404/93‑1 (MQ=255)
                                                 tATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAg                                         >  2:6002745/1‑93 (MQ=255)
                                                  aTTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGt                                        >  2:3969400/1‑93 (MQ=255)
                                                  aTTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGt                                        >  1:5900997/1‑93 (MQ=255)
                                                  aTTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGt                                        <  2:1788005/93‑1 (MQ=255)
                                                   ttCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGtt                                       <  2:4075190/93‑1 (MQ=255)
                                                   ttCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGtt                                       >  1:3700793/1‑93 (MQ=255)
                                                     ccGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACTCGCGGGATCTGACTCGCGGCGAGTTGa                                     <  2:6418528/93‑1 (MQ=255)
                                                     ccGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGa                                     >  1:5863133/1‑93 (MQ=255)
                                                      cGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGaa                                    >  1:922897/1‑93 (MQ=255)
                                                       gATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGaaa                                   <  2:1733322/93‑1 (MQ=255)
                                                       gATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGaaa                                   <  1:3338013/93‑1 (MQ=255)
                                                       gATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGaaa                                   >  1:5326764/1‑93 (MQ=255)
                                                       gATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGaaa                                   <  2:79747/93‑1 (MQ=255)
                                                        aTCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAg                                  <  2:1490850/93‑1 (MQ=255)
                                                        aTCAGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAg                                  >  2:4162056/1‑93 (MQ=255)
                                                          cAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTc                                <  2:6685668/93‑1 (MQ=255)
                                                           aGACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCg                               >  1:5992101/1‑93 (MQ=255)
                                                            gACGCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGa                              >  1:6793739/1‑93 (MQ=255)
                                                               gCCGGGACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGaaaa                           <  2:3056709/93‑1 (MQ=255)
                                                                  gggACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATaa                        >  1:6282296/1‑93 (MQ=255)
                                                                   ggACCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGCATCTGACTCGCGGCGAGTTGAAAGTCGAAAATaaa                       <  2:4525104/93‑1 (MQ=255)
                                                                   ggACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATaaa                       >  2:4782555/1‑93 (MQ=255)
                                                                    gACCGGGCCTGAAAAAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAAc                      >  2:6448210/1‑93 (MQ=255)
                                                                     aCCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACt                     <  1:2302809/93‑1 (MQ=255)
                                                                     aCCGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACt                     <  1:490613/93‑1 (MQ=255)
                                                                      ccGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtg                    <  2:4312753/93‑1 (MQ=255)
                                                                       cGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtgt                   <  2:5958672/93‑1 (MQ=255)
                                                                       cGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtgt                   <  2:6199124/93‑1 (MQ=255)
                                                                       cGGGCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtgt                   <  1:7455896/93‑1 (MQ=255)
                                                                        gggCCTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACtgtg                  <  1:203006/93‑1 (MQ=255)
                                                                           ccTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGtt               <  1:1705321/93‑1 (MQ=255)
                                                                            cTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGttt              >  2:6555116/1‑93 (MQ=255)
                                                                            cTGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGttt              >  1:7254193/1‑93 (MQ=255)
                                                                             tGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTa             >  2:5655221/1‑93 (MQ=255)
                                                                             tGAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTa             >  1:4865074/1‑93 (MQ=255)
                                                                              gAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTAc            >  1:503554/1‑93 (MQ=255)
                                                                              gAAAGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTAc            <  2:2498731/93‑1 (MQ=255)
                                                                                aaGAGTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCa          <  2:4677151/93‑1 (MQ=255)
                                                                                 agagTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCAc         >  1:4192900/1‑93 (MQ=255)
                                                                                 agagTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCAc         <  2:2179408/93‑1 (MQ=255)
                                                                                  gagTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACg        >  1:674024/1‑93 (MQ=255)
                                                                                  gagTACAAGCAGGTTATCGCGCCGCTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACg        >  2:7116055/1‑93 (MQ=255)
                                                                                   agTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGc       <  2:7431299/93‑1 (MQ=255)
                                                                                    gTACAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCt      >  1:2842357/1‑93 (MQ=255)
                                                                                       cAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTtga   <  1:1632101/93‑1 (MQ=255)
                                                                                       cAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTtga   <  2:5433462/93‑1 (MQ=255)
                                                                                       cAAGCAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTtga   <  2:5315871/93‑1 (MQ=255)
                                                                                          |                                                                                          
CTACGGCGACTATCCCAACAACTATAACTTCTGTCTTGATGGTTTGATCTATTCCGATCAGACGCCGGGACCGGGCCTGAAAGAGTACAAACAGGTTATCGCGCCGGTAAAAATCCACGCGCGGGATCTGACTCGCGGCGAGTTGAAAGTCGAAAATAAACTGTGGTTTACCACGCTTGAT  >  NC_000913/3221971‑3222151

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