Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,222,944 4 bp→AGCA coding (2671‑2674/3093 nt) 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,9440CA100.0% 469.5 / NA 134Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (60/74);  total (60/74)
*NC_0009133,222,9450AG100.0% 427.2 / NA 133Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (59/74);  total (59/74)
*NC_0009133,222,9460AC100.0% 428.0 / NA 133Q891S (CAA→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (59/74);  total (59/74)
*NC_0009133,222,9470GA100.0% 458.7 / NA 133A892T (GCC→ACC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (59/74);  total (59/74)

Only 100 of 140 total aligned reads displayed.

GATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222863‑3223027
                                                                                 ||||                                                                                
gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGAt                                                                          >  1:31919/1‑93 (MQ=255)
gATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGAt                                                                          >  1:1393991/1‑93 (MQ=255)
 aTCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATg                                                                         >  1:4586229/1‑93 (MQ=255)
 aTCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATg                                                                         >  2:1723090/1‑93 (MQ=255)
    aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACATCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                      <  1:160144/93‑1 (MQ=255)
    aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                      <  2:407627/93‑1 (MQ=255)
    aGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCa                                                                      <  2:4549716/93‑1 (MQ=255)
     ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                                     >  1:498649/1‑93 (MQ=255)
     ggTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCAt                                                                     <  2:3497866/93‑1 (MQ=255)
      gTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                                    >  1:2309604/1‑93 (MQ=255)
      gTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                                    <  2:316878/93‑1 (MQ=255)
      gTGGCGTATTACGGTCGTGGACCGGGCGAAAACTAAGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATg                                                                    >  2:4700931/1‑93 (MQ=255)
       gggCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGt                                                                   <  1:2291094/92‑1 (MQ=255)
        ggCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGtt                                                                  <  2:4185237/93‑1 (MQ=255)
         gCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTc                                                                 <  2:357512/93‑1 (MQ=255)
         gCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTc                                                                 >  2:3821375/1‑93 (MQ=255)
          cGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCg                                                                <  1:4133594/93‑1 (MQ=255)
           gTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                                               >  1:2390078/1‑93 (MQ=255)
           gTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                                               >  2:2851328/1‑93 (MQ=255)
           gTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                                               <  1:1495953/93‑1 (MQ=255)
           gTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCga                                                               <  2:781707/93‑1 (MQ=255)
            tATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCgag                                                              >  1:2138683/1‑93 (MQ=255)
              ttACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGaa                                                            >  1:4427027/1‑93 (MQ=255)
               tACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAAc                                                           >  1:1810686/1‑93 (MQ=255)
               tACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAAc                                                           <  2:520824/93‑1 (MQ=255)
                 cGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTa                                                         <  2:2190314/93‑1 (MQ=255)
                   gTCGTGTACCGGGCGAAAACTACGCCGACAGCCAACAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGAGGCCATGTTCGAGTACTATc                                                       >  2:4703715/1‑93 (MQ=255)
                   gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATc                                                       >  2:611537/1‑93 (MQ=255)
                   gTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATc                                                       <  1:1025367/93‑1 (MQ=255)
                    tCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                      <  2:4560711/93‑1 (MQ=255)
                    tCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                      <  2:2132715/93‑1 (MQ=255)
                    tCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                      <  1:952277/93‑1 (MQ=255)
                    tCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCAACGTCGATGCCCTGTTCGAGAACTATcc                                                      >  1:2847129/1‑93 (MQ=255)
                    tCGGGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATcc                                                      <  2:949143/93‑1 (MQ=255)
                     cGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATccc                                                     >  1:3014488/1‑93 (MQ=255)
                      gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                                    <  2:129325/93‑1 (MQ=255)
                      gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                                    <  2:107832/93‑1 (MQ=255)
                      gTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCt                                                    <  2:3536064/93‑1 (MQ=255)
                       tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                   >  1:3713617/1‑93 (MQ=255)
                       tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                   >  1:4376785/1‑93 (MQ=255)
                       tGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCtt                                                   <  1:264231/93‑1 (MQ=255)
                        ggACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTc                                                  <  2:2404052/93‑1 (MQ=255)
                        ggACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTc                                                  >  2:4306589/1‑93 (MQ=255)
                        ggACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGAGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTc                                                  >  1:2615285/1‑93 (MQ=255)
                          aCCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTccc                                                >  2:2935707/1‑93 (MQ=255)
                           ccGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCg                                               <  1:1184680/93‑1 (MQ=255)
                              ggCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGCAAACATTCCCTTCCCgaaa                                            >  2:1746857/1‑90 (MQ=255)
                               gCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGa                                           <  2:3001780/93‑1 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                       <  1:46970/93‑1 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                       <  2:2188212/93‑1 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                       >  2:3100547/1‑93 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                       <  2:1419764/93‑1 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaa                                       <  1:3005432/93‑1 (MQ=255)
                                   aaaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCCTCCCGCAGaacaa                                       <  1:2042778/93‑1 (MQ=255)
                                    aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                      >  2:4175887/1‑93 (MQ=255)
                                    aaCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGaacaac                                      >  1:2657602/1‑93 (MQ=255)
                                     aCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                                     >  1:1374768/1‑93 (MQ=255)
                                     aCTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACg                                     >  1:1741763/1‑93 (MQ=255)
                                      cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                                    <  2:366372/93‑1 (MQ=255)
                                      cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                                    >  1:3724047/1‑93 (MQ=255)
                                      cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                                    >  2:1258874/1‑93 (MQ=255)
                                      cTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACgg                                    <  2:1092213/93‑1 (MQ=255)
                                       tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                   >  2:2826798/1‑93 (MQ=255)
                                       tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                   <  1:1306736/93‑1 (MQ=255)
                                       tACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGt                                   <  1:1306737/93‑1 (MQ=255)
                                          cccGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAAc                                <  1:791662/92‑1 (MQ=255)
                                            cGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCg                              >  2:3553203/1‑93 (MQ=255)
                                                agccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAg                          >  2:381759/1‑93 (MQ=255)
                                                  ccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCa                        <  2:4278409/93‑1 (MQ=255)
                                                  ccagcAGGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCa                        <  2:4010158/93‑1 (MQ=255)
                                                   cagcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGGAACCGTCAGCag                       >  1:1606529/1‑92 (MQ=255)
                                                     gcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGt                     >  2:454539/1‑93 (MQ=255)
                                                     gcagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGt                     >  2:3552796/1‑93 (MQ=255)
                                                      cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTc                    <  1:887083/93‑1 (MQ=255)
                                                      cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTc                    <  1:3265046/93‑1 (MQ=255)
                                                      cagGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTc                    <  2:3635900/93‑1 (MQ=255)
                                                      cagGCGAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTc                    <  2:2759251/93‑1 (MQ=255)
                                                       agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTcc                   >  1:297102/1‑93 (MQ=255)
                                                       agGCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGCAACAGTCAGCACGTcc                   >  2:4229780/1‑93 (MQ=255)
                                                        ggCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCg                  <  1:3492101/93‑1 (MQ=255)
                                                        ggCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCg                  >  2:2314505/1‑93 (MQ=255)
                                                         gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                 <  2:2319567/93‑1 (MQ=255)
                                                         gCTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGc                 <  2:3797277/93‑1 (MQ=255)
                                                          cTAACATCATCGATATCTGGCGCAGCACTGTCGATGCCATGTTCGAGAACTAACCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCt                <  1:4084049/93‑1 (MQ=255)
                                                          cTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCt                <  1:4084048/93‑1 (MQ=255)
                                                          cTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCt                >  1:2668885/1‑93 (MQ=255)
                                                          cTAACATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCt                >  1:2790167/1‑93 (MQ=255)
                                                            aaCATCATCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTCACCGTCAGCATGTCCGCTgg              >  1:1006670/1‑93 (MQ=255)
                                                              catcatCGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGAc            >  1:3146114/1‑93 (MQ=255)
                                                                tcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACgg          >  1:2529267/1‑93 (MQ=255)
                                                                tcatcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTACCCCTTCCCGCAGAACAACGGCAACCGTCAGCATGTCCGCTGGACgc          >  1:2529268/1‑92 (MQ=255)
                                                                 catcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGc         <  2:3407504/93‑1 (MQ=255)
                                                                 catcGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGc         <  2:7204/93‑1 (MQ=255)
                                                                    cGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCAcc      >  2:3344045/1‑92 (MQ=255)
                                                                    cGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACt      >  1:1749807/1‑93 (MQ=255)
                                                                    cGATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACt      <  1:601932/93‑1 (MQ=255)
                                                                     gATATCTGGCGCAGCCCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTg     <  2:2379265/93‑1 (MQ=255)
                                                                     gATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTGCGCTGGACGGCACTg     >  2:3299712/1‑93 (MQ=255)
                                                                     gATATCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTg     <  2:2379266/93‑1 (MQ=255)
                                                                       tatCTGGCGCAGCACCGTCGATGCCATGTTCGAGAACTATCCCTGCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGAc   <  1:4678722/93‑1 (MQ=255)
                                                                                 ||||                                                                                
GATCAGGTGGCGTATTACGGTCGTGGACCGGGCGAAAACTACGCCGACAGCCAGCAGGCTAACATCATCGATATCTGGCGCCAAGCCGTCGATGCCATGTTCGAGAACTATCCCTTCCCGCAGAACAACGGTAACCGTCAGCATGTCCGCTGGACGGCACTGACT  >  NC_000913/3222863‑3223027

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: ‑