Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,218,296 G→T 17.9% D769Y (GAC→TAC)  ebgA → evolved beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,218,2960GT17.9% 551.2 / 69.3 189D769Y (GAC→TAC) ebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (54/101);  new base T (11/23);  total (65/124)
Fisher's exact test for biased strand distribution p-value = 8.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 190 total aligned reads displayed.

GAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGAC  >  CP009273/3218149‑3218442
                                                                                                                                                   |                                                                                                                                                  
gAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCg                                                                                                                                                    >  2:1399005/1‑148 (MQ=255)
 aaGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGa                                                                                                                                                   <  1:4231/148‑1 (MQ=255)
   gACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGAcg                                                                                                                                                 <  2:2311992/148‑2 (MQ=255)
   gACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACa                                                                                                                                                 >  2:1080387/1‑148 (MQ=255)
    aCGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACaa                                                                                                                                                <  2:2391351/148‑1 (MQ=255)
     cGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAAc                                                                                                                                               <  1:2986542/148‑1 (MQ=255)
       aTCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACca                                                                                                                                             <  2:1671816/148‑1 (MQ=255)
        tCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACcac                                                                                                                                            >  1:2973283/1‑148 (MQ=255)
            ttGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGc                                                                                                                                        >  1:2140695/1‑148 (MQ=255)
               aGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAgg                                                                                                                                     <  2:3145763/148‑1 (MQ=255)
               aGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAgg                                                                                                                                     >  2:1226326/1‑148 (MQ=255)
               aGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAgg                                                                                                                                     <  1:2589079/148‑1 (MQ=255)
                gCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGa                                                                                                                                    <  2:450337/148‑1 (MQ=255)
                  tGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGt                                                                                                                                  >  2:3260257/1‑148 (MQ=255)
                      ccGTTCGCGTCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGa                                                                                                                              <  2:3724323/148‑1 (MQ=255)
                      ccGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGa                                                                                                                              >  2:3597474/1‑148 (MQ=255)
                       cGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGaa                                                                                                                             <  1:120745/148‑1 (MQ=255)
                        gTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAg                                                                                                                            <  2:241948/148‑1 (MQ=255)
                        gTTCGCGGCTACAACTTCGAGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTTGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAg                                                                                                                            <  2:1614213/148‑1 (MQ=255)
                          tCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAggg                                                                                                                          >  2:2359449/1‑148 (MQ=255)
                          tCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAggg                                                                                                                          <  1:3359139/148‑1 (MQ=255)
                            gcgGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCt                                                                                                                        <  1:1672345/148‑1 (MQ=255)
                                  aCAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCa                                                                                                                  <  1:959909/148‑1 (MQ=255)
                                   cAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCaa                                                                                                                 <  2:512724/148‑1 (MQ=255)
                                      cTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCg                                                                                                              >  2:1882204/1‑148 (MQ=255)
                                        tCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGaa                                                                                                            >  1:3059852/1‑148 (MQ=255)
                                           cgATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCa                                                                                                         <  1:1844334/148‑1 (MQ=255)
                                            gATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCAt                                                                                                        <  2:2174282/148‑1 (MQ=255)
                                            gATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCAt                                                                                                        <  2:3416515/148‑1 (MQ=255)
                                            gATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCAt                                                                                                        <  2:1328451/148‑1 (MQ=255)
                                            gATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCAt                                                                                                        <  1:2515754/148‑1 (MQ=255)
                                             aTCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCagt                                                                                                       >  2:1990593/1‑146 (MQ=255)
                                              tCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCAttt                                                                                                      >  1:1402806/1‑148 (MQ=255)
                                               cACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTg                                                                                                     >  1:3260956/1‑148 (MQ=255)
                                                aCCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                                                                    >  2:3814202/1‑148 (MQ=255)
                                                  cTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAg                                                                                                  >  1:3428855/1‑148 (MQ=255)
                                                   ttCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGa                                                                                                 <  2:1785923/148‑1 (MQ=255)
                                                     ctcAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATc                                                                                               <  2:2225751/148‑1 (MQ=255)
                                                      tcAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCa                                                                                              <  2:3313416/148‑1 (MQ=255)
                                                      tcAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCa                                                                                              <  1:2491386/148‑1 (MQ=255)
                                                         aaaaTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGc                                                                                           <  1:2754693/148‑1 (MQ=255)
                                                          aaaTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCa                                                                                          <  2:574277/148‑1 (MQ=255)
                                                            aTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAgg                                                                                        >  2:688408/1‑148 (MQ=255)
                                                                gTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACa                                                                                    <  1:1082140/148‑1 (MQ=255)
                                                                 tGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACAt                                                                                   <  2:3510641/148‑1 (MQ=255)
                                                                  ggCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATc                                                                                  <  2:1738035/148‑1 (MQ=255)
                                                                  ggCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATc                                                                                  <  2:641933/148‑1 (MQ=255)
                                                                          gACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGAc                                                                          <  1:527549/148‑1 (MQ=255)
                                                                           aCATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACt                                                                         <  1:2128180/148‑1 (MQ=255)
                                                                            cATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACtt                                                                        >  2:2889737/1‑148 (MQ=255)
                                                                            cATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACtt                                                                        >  2:3374941/1‑148 (MQ=255)
                                                                              tCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTg                                                                      >  1:776903/1‑148 (MQ=255)
                                                                              tCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTg                                                                      >  1:776874/1‑148 (MQ=255)
                                                                                  ggCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGt                                                                  <  2:3068311/148‑1 (MQ=255)
                                                                                      ggTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGaa                                                              <  1:2712505/148‑1 (MQ=255)
                                                                                      ggTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGaa                                                              <  2:2378658/148‑1 (MQ=255)
                                                                                      ggTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGaa                                                              <  1:1954537/148‑1 (MQ=255)
                                                                                       gTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAAc                                                             >  2:116553/1‑148 (MQ=255)
                                                                                        tGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACa                                                            <  2:1573875/148‑1 (MQ=255)
                                                                                         gAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACag                                                           <  1:1341597/148‑1 (MQ=255)
                                                                                         gAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACag                                                           <  1:758592/148‑1 (MQ=255)
                                                                                          aaTGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACaga                                                          <  2:3457340/148‑1 (MQ=255)
                                                                                          aaTGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACaga                                                          <  2:3219328/148‑1 (MQ=255)
                                                                                            tGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGc                                                        <  2:2140695/148‑1 (MQ=255)
                                                                                              gCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGa                                                      <  2:1696057/148‑1 (MQ=255)
                                                                                               cGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGAt                                                     <  2:2609620/148‑1 (MQ=255)
                                                                                               cGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGAt                                                     <  1:240687/148‑1 (MQ=255)
                                                                                                 aaTCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACCAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATgg                                                   <  2:2899817/148‑1 (MQ=255)
                                                                                                    cGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGt                                                  <  2:2353417/146‑1 (MQ=255)
                                                                                                    cGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGa                                                <  1:3566875/148‑1 (MQ=255)
                                                                                                     gctgctGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGaa                                               <  2:1461983/148‑1 (MQ=255)
                                                                                                     gctgctGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGaa                                               <  2:1186429/148‑1 (MQ=255)
                                                                                                        gctgACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTg                                            <  1:3712430/148‑1 (MQ=255)
                                                                                                             cTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGat                                       >  2:872792/1‑148 (MQ=255)
                                                                                                               cgcgAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGatca                                     <  1:689085/148‑1 (MQ=255)
                                                                                                                     ccAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGcc                               >  1:1765805/1‑148 (MQ=255)
                                                                                                                      cAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCg                              >  1:262556/1‑148 (MQ=255)
                                                                                                                        aaGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGca                            >  1:892797/1‑148 (MQ=255)
                                                                                                                         aGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGcac                           >  1:2847893/1‑148 (MQ=255)
                                                                                                                          gATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGcaca                          >  2:3677703/1‑148 (MQ=255)
                                                                                                                           aTCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAg                         <  1:1034594/148‑1 (MQ=255)
                                                                                                                           aTCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAg                         <  2:2741231/148‑1 (MQ=255)
                                                                                                                           aTCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAg                         <  1:1033901/148‑1 (MQ=255)
                                                                                                                             cAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGtt                       <  1:1098080/148‑1 (MQ=255)
                                                                                                                             cAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGtt                       >  1:2337562/1‑148 (MQ=255)
                                                                                                                              aaCTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTa                      <  1:3506067/148‑1 (MQ=255)
                                                                                                                                 ttcttcAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTg                   <  1:189157/148‑1 (MQ=255)
                                                                                                                                   cttcAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGcc                 <  1:3330643/148‑1 (MQ=255)
                                                                                                                                   cttcAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGcc                 <  1:3329726/148‑1 (MQ=255)
                                                                                                                                    ttcAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGccc                <  1:598193/148‑1 (MQ=255)
                                                                                                                                       aaGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCccgc             >  2:9238/1‑148 (MQ=255)
                                                                                                                                       aaGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCccgc             >  2:1940189/1‑148 (MQ=255)
                                                                                                                                         gCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCccgccg           >  1:1589013/1‑148 (MQ=255)
                                                                                                                                            gatgatCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGgtg        <  1:3277447/148‑1 (MQ=255)
                                                                                                                                              tgatCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGtt      >  1:1815404/1‑148 (MQ=255)
                                                                                                                                                atCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTg    <  1:436261/148‑1 (MQ=255)
                                                                                                                                                atCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTg    <  1:1882157/148‑1 (MQ=255)
                                                                                                                                                atCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTg    >  2:2819556/1‑148 (MQ=255)
                                                                                                                                                 tCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGa   <  1:2359449/148‑1 (MQ=255)
                                                                                                                                                  cGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGAc  <  2:2229217/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGAC  >  CP009273/3218149‑3218442

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