Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,218,296 G→T 8.1% 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,2960GT 8.1% 712.9 / 24.8 196D769Y (GAC→TAC) ebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (75/105);  new base T (3/13);  total (78/118)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 196 total aligned reads displayed.

GAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACT  >  CP009273/3218149‑3218443
                                                                                                                                                   |                                                                                                                                                   
 aaGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGAc                                                                                                                                                   <  2:3309702/149‑1 (MQ=255)
    aCGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAAc                                                                                                                                                <  2:52051/149‑1 (MQ=255)
      gATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACca                                                                                                                                              <  1:1132082/149‑1 (MQ=255)
       aTCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACcac                                                                                                                                             >  2:3952042/1‑149 (MQ=255)
         cGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACaa                                                                                                                                           >  2:1670933/1‑149 (MQ=255)
         cGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACaa                                                                                                                                           <  2:1120312/149‑1 (MQ=255)
           tttGAGCTGCACCGTTCGCGGCTACAGCTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCACGCCGATGATCGACAACCACAAGc                                                                                                                                         >  1:364604/1‑149 (MQ=255)
           tttGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGc                                                                                                                                         >  1:364759/1‑149 (MQ=255)
            ttGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCa                                                                                                                                        <  1:1672478/149‑1 (MQ=255)
              gAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAgg                                                                                                                                      <  1:1038691/149‑1 (MQ=255)
                gCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAg                                                                                                                                    >  1:3147831/1‑149 (MQ=255)
                      ccGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGaa                                                                                                                              <  1:2985002/149‑1 (MQ=255)
                      ccGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGaa                                                                                                                              <  1:537380/149‑1 (MQ=255)
                          tCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGc                                                                                                                          <  1:3656920/149‑1 (MQ=255)
                                  aCAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCaa                                                                                                                  >  1:3907582/1‑149 (MQ=255)
                                   cAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAAc                                                                                                                 >  1:748098/1‑149 (MQ=255)
                                    aaCTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAGCCACAAGCAGGAGTACGAAGGGCTGTGGCAAcc                                                                                                                >  1:3045039/1‑149 (MQ=255)
                                     aCTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCg                                                                                                               >  1:825108/1‑149 (MQ=255)
                                        tCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAAt                                                                                                            <  1:957762/149‑1 (MQ=255)
                                         cgcgATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATc                                                                                                           >  1:2594171/1‑149 (MQ=255)
                                               cACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                                                                     <  2:2579527/149‑1 (MQ=255)
                                               cACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                                                                     <  1:1946305/149‑1 (MQ=255)
                                                aCCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCa                                                                                                    <  2:564371/149‑1 (MQ=255)
                                                 ccTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAg                                                                                                   >  2:3824154/1‑149 (MQ=255)
                                                   ttCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGAt                                                                                                 <  1:2470269/149‑1 (MQ=255)
                                                   ttCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGAt                                                                                                 <  1:1608797/149‑1 (MQ=255)
                                                     ctcAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAAGCATTTGCAGATCa                                                                                               >  2:3443490/1‑149 (MQ=255)
                                                      tcAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCAt                                                                                              <  1:1610709/149‑1 (MQ=255)
                                                        aaaaaTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGc                                                                                            <  1:3580845/149‑1 (MQ=255)
                                                         aaaaTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCa                                                                                           <  2:1570320/149‑1 (MQ=255)
                                                         aaaaTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCa                                                                                           <  1:1002946/149‑1 (MQ=255)
                                                         aaaaTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCa                                                                                           <  1:1826689/149‑1 (MQ=255)
                                                          aaaTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAg                                                                                          >  2:3921161/1‑149 (MQ=255)
                                                           aaTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAgg                                                                                         >  1:2654822/1‑149 (MQ=255)
                                                             tGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGaa                                                                                       >  2:801504/1‑149 (MQ=255)
                                                              gAGTGACAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAAc                                                                                      >  2:3117660/1‑149 (MQ=255)
                                                               aGTGGCAAACCGACATCCTGGCAGGTGAATGGCGTATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACa                                                                                     <  2:3094118/149‑1 (MQ=255)
                                                               aGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACa                                                                                     <  2:632229/149‑1 (MQ=255)
                                                               aGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACa                                                                                     <  2:177568/149‑1 (MQ=255)
                                                                 tGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATc                                                                                   <  1:1333718/149‑1 (MQ=255)
                                                                  ggCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCt                                                                                  >  2:1477118/1‑149 (MQ=255)
                                                                    cAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTgc                                                                                <  1:2673948/149‑1 (MQ=255)
                                                                    cAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTgc                                                                                <  1:232949/149‑1 (MQ=255)
                                                                        ccGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGa                                                                            <  1:1287955/149‑1 (MQ=255)
                                                                         cGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGAc                                                                           <  1:1765857/149‑1 (MQ=255)
                                                                          gACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACt                                                                          <  1:2654747/149‑1 (MQ=255)
                                                                          gACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACt                                                                          <  1:2654698/149‑1 (MQ=255)
                                                                           aCATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACtt                                                                         <  1:3610728/149‑1 (MQ=255)
                                                                             ntCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTg                                                                       >  1:711661/2‑149 (MQ=255)
                                                                              tCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGc                                                                      >  2:1934750/1‑149 (MQ=255)
                                                                                   gCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAg                                                                 >  1:1137971/1‑149 (MQ=255)
                                                                                     aGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGaa                                                               >  2:449258/1‑149 (MQ=255)
                                                                                       gTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACa                                                             >  1:3540874/1‑149 (MQ=255)
                                                                                       gTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACa                                                             <  1:324849/149‑1 (MQ=255)
                                                                                         gAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACaga                                                           <  1:425891/149‑1 (MQ=255)
                                                                                         gAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACaga                                                           >  2:3814019/1‑149 (MQ=255)
                                                                                         gAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACaga                                                           <  1:3129778/149‑1 (MQ=255)
                                                                                          aaTGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACagag                                                          >  2:1325858/1‑149 (MQ=255)
                                                                                             ggCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGa                                                       <  1:2064149/149‑1 (MQ=255)
                                                                                             ggCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGa                                                       <  1:2523597/149‑1 (MQ=255)
                                                                                              gCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGAt                                                      <  1:1532756/149‑1 (MQ=255)
                                                                                               cGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATg                                                     <  1:506403/149‑1 (MQ=255)
                                                                                                gAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATgg                                                    <  2:432297/149‑1 (MQ=255)
                                                                                                 aaTCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGt                                                   >  1:2237149/1‑149 (MQ=255)
                                                                                                 aaTCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGt                                                   <  1:2541415/149‑1 (MQ=255)
                                                                                                 aaTCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGt                                                   >  2:2516400/1‑149 (MQ=255)
                                                                                                 aaTCGCTGCTGACTCGCGAGCAAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGt                                                   >  2:2670975/1‑149 (MQ=255)
                                                                                                 aaTCGCTGCTGACTCGCGAGCAAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGt                                                   >  2:2670920/1‑149 (MQ=255)
                                                                                                     gctgctGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGTGCGACTTTGCCGTAGAACAGAGCGATGGTGAAg                                               <  1:3033894/149‑1 (MQ=255)
                                                                                                       tgctgACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTg                                             >  1:163450/1‑149 (MQ=255)
                                                                                                        gctgACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGc                                            <  1:785025/149‑1 (MQ=255)
                                                                                                        gctgACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGc                                            <  1:3932665/149‑1 (MQ=255)
                                                                                                        gctgACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGc                                            <  1:294612/149‑1 (MQ=255)
                                                                                                        gctgACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGc                                            <  2:2909167/149‑1 (MQ=255)
                                                                                                           gACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGa                                         >  2:2340018/1‑149 (MQ=255)
                                                                                                            aCTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGat                                        <  2:253583/149‑1 (MQ=255)
                                                                                                             cTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGatc                                       <  2:102334/149‑1 (MQ=255)
                                                                                                              tCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGatca                                      <  1:449258/149‑1 (MQ=255)
                                                                                                               cgcgAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGatcat                                     <  2:992060/149‑1 (MQ=255)
                                                                                                                gcgAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGatcatc                                    <  1:1182376/149‑1 (MQ=255)
                                                                                                                gcgAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGatcatc                                    <  2:3152760/149‑1 (MQ=255)
                                                                                                                  gAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAg                                  >  2:1621589/1‑149 (MQ=255)
                                                                                                                      cAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGc                              >  2:904245/1‑149 (MQ=255)
                                                                                                                        aaGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGTGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGcac                            >  2:3664447/1‑149 (MQ=255)
                                                                                                                        aaGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGTGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGcac                            >  2:3670637/1‑149 (MQ=255)
                                                                                                                         aGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGcaca                           >  1:2201253/1‑149 (MQ=255)
                                                                                                                           aTCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGAACAGt                         >  2:3825829/1‑149 (MQ=255)
                                                                                                                              aaCTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTAt                      <  1:3409841/149‑1 (MQ=255)
                                                                                                                               aCTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTAtt                     >  1:2758510/1‑149 (MQ=255)
                                                                                                                               aCTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTAtt                     <  1:1537320/149‑1 (MQ=255)
                                                                                                                               aCTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTAtt                     >  2:671879/1‑149 (MQ=255)
                                                                                                                               aCTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTAtt                     <  1:3893685/149‑1 (MQ=255)
                                                                                                                               aCTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTAtt                     <  1:3548815/149‑1 (MQ=255)
                                                                                                                                 ttcttcAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGc                   <  2:1179763/149‑1 (MQ=255)
                                                                                                                                  tcttcAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGcc                  >  1:3776394/1‑149 (MQ=255)
                                                                                                                                    ttcAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGcccc                <  1:1233873/149‑1 (MQ=255)
                                                                                                                                      cAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCccgc              >  2:4089397/1‑149 (MQ=255)
                                                                                                                                       aaGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCccgcc             >  1:4075582/1‑149 (MQ=255)
                                                                                                                                         gCCGATGATCTACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCgg           <  1:643876/149‑1 (MQ=255)
                                                                                                                                                  cGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACt  <  2:1115209/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACT  >  CP009273/3218149‑3218443

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: