Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,217,289 A→T 89.0% H433L (CAC→CTC)  ebgA → evolved beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,217,2890AT89.0% 869.8 / 76.8 278H433L (CAC→CTC) ebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (13/18);  new base T (122/125);  total (135/143)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

ACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGAAGTGG  >  CP009273/3217142‑3217438
                                                                                                                                                   |                                                                                                                                                     
aCGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTcg                                                                                                                                                     <  1:2045266/150‑1 (MQ=255)
    aCTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCa                                                                                                                                                 <  2:4114726/150‑1 (MQ=255)
       gtgtgATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGaa                                                                                                                                              >  2:469612/1‑150 (MQ=255)
                 tACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCg                                                                                                                                    <  1:770954/150‑1 (MQ=255)
                  aCGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGt                                                                                                                                   >  2:3619952/1‑150 (MQ=255)
                    ggCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCg                                                                                                                                 >  2:4283184/1‑150 (MQ=255)
                     gCCTGTTTGTGATGGCGGAAACCGACGTCTAATCGAACGGCTTTGCTAATGTCGGCTATATTAGCAGTATTACCGACGATCCTCAGTGGTAAAAAGTATACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGa                                                                                                                                <  2:1609374/150‑1 (MQ=255)
                     gCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGa                                                                                                                                <  1:932037/150‑1 (MQ=255)
                           ttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatca                                                                                                                          <  1:2618962/150‑1 (MQ=255)
                           ttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatca                                                                                                                          <  2:3588766/150‑1 (MQ=255)
                             gtgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatcatc                                                                                                                        >  2:1236935/1‑150 (MQ=255)
                               gATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTg                                                                                                                      <  2:2837542/150‑1 (MQ=255)
                               gATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTg                                                                                                                      <  2:3160268/150‑1 (MQ=255)
                                  ggCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTc                                                                                                                   <  2:2342991/150‑1 (MQ=255)
                                     ggAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCt                                                                                                                >  1:415604/1‑150 (MQ=255)
                                     ggAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCt                                                                                                                >  2:522592/1‑150 (MQ=255)
                                       aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTgg                                                                                                              <  1:1567845/150‑1 (MQ=255)
                                          ccGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCa                                                                                                           <  1:4051604/150‑1 (MQ=255)
                                          ccGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCa                                                                                                           <  2:1087680/150‑1 (MQ=255)
                                          ccGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCa                                                                                                           <  2:2648552/150‑1 (MQ=255)
                                           cGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCaa                                                                                                          <  2:2125693/150‑1 (MQ=255)
                                            gACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAAt                                                                                                         >  2:832367/1‑150 (MQ=255)
                                             aCGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATg                                                                                                        <  1:47008/150‑1 (MQ=255)
                                             aCGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATg                                                                                                        <  1:1077174/150‑1 (MQ=255)
                                                  gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATcc                                                                                                   <  1:484500/150‑1 (MQ=255)
                                                  gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGTTCATCATCTGGTCGCTGGGCAATGAATcc                                                                                                   <  2:4000878/150‑1 (MQ=255)
                                                       gCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTa                                                                                              >  2:3291544/1‑150 (MQ=255)
                                                       gCACGGCTTTGCTAATGTCGGCGATATTAGCCATATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTa                                                                                              >  2:4478420/1‑150 (MQ=255)
                                                          cGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATgg                                                                                           >  1:3374175/1‑150 (MQ=255)
                                                          cGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATgg                                                                                           >  2:2204838/1‑150 (MQ=255)
                                                           ggCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGc                                                                                          >  1:747754/1‑150 (MQ=255)
                                                           ggCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGc                                                                                          >  2:2980214/1‑150 (MQ=255)
                                                           ggCGTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGc                                                                                          >  1:1028297/1‑150 (MQ=255)
                                                             cTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTg                                                                                        >  2:3922038/1‑150 (MQ=255)
                                                             cTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTg                                                                                        >  2:3713054/1‑150 (MQ=255)
                                                               ttGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTa                                                                                      <  2:2773601/150‑1 (MQ=255)
                                                                tGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTaa                                                                                     <  1:1938230/150‑1 (MQ=255)
                                                                 gCTAATGTCGGCGATTTTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGATAAACCATCCTTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAAc                                                                                    <  1:1864811/150‑1 (MQ=255)
                                                                 gCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAAc                                                                                    <  1:4192294/150‑1 (MQ=255)
                                                                  cTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATAATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACa                                                                                   <  2:703130/150‑1 (MQ=255)
                                                                   tAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTTGGAAAAAGTCTACGTCGAGAGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATTAATCCGGCTATGGCTGTAACAt                                                                                  >  2:3167780/1‑150 (MQ=255)
                                                                    aaTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATc                                                                                 >  1:3739565/1‑150 (MQ=255)
                                                                     aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATcc                                                                                >  2:73874/1‑150 (MQ=255)
                                                                     aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATcc                                                                                >  2:2381341/1‑150 (MQ=255)
                                                                      tGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCcg                                                                               <  1:1744954/150‑1 (MQ=255)
                                                                       gTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCcgc                                                                              <  2:2439439/150‑1 (MQ=255)
                                                                        tCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCcgcg                                                                             <  1:4462250/150‑1 (MQ=255)
                                                                            cGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGAt                                                                         >  2:2921624/1‑150 (MQ=255)
                                                                             gATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATg                                                                        >  2:2009449/1‑150 (MQ=255)
                                                                              atatTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGt                                                                       >  2:490655/1‑150 (MQ=255)
                                                                               tatTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTa                                                                      >  1:2443520/1‑150 (MQ=255)
                                                                                atTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTAc                                                                     >  1:411961/1‑150 (MQ=255)
                                                                                 ttAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTAcc                                                                    >  1:638854/1‑150 (MQ=255)
                                                                                  tAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCa                                                                   >  2:2936676/1‑150 (MQ=255)
                                                                                   aGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCAt                                                                  >  1:2408760/1‑150 (MQ=255)
                                                                                   aGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCAt                                                                  >  1:401084/1‑150 (MQ=255)
                                                                                     ccGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATgc                                                                <  1:1655083/150‑1 (MQ=255)
                                                                                         aTTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATgcggcg                                                            >  2:3508172/1‑150 (MQ=255)
                                                                                         aTTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATgcggcg                                                            >  1:4495521/1‑150 (MQ=255)
                                                                                         aTTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATgcggcg                                                            <  1:3035010/150‑1 (MQ=255)
                                                                                           tACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGaa                                                          <  1:3714591/150‑1 (MQ=255)
                                                                                             ccGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                        >  1:3224081/1‑150 (MQ=255)
                                                                                             ccGACGATCCGCAGGGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCAACCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGACAg                                                        >  1:4347938/1‑150 (MQ=255)
                                                                                                  gaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTg                                                   <  1:1432081/150‑1 (MQ=255)
                                                                                                   aTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTgg                                                  >  2:3299933/1‑150 (MQ=255)
                                                                                                     ccGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGAt                                                >  2:384478/1‑150 (MQ=255)
                                                                                                      cGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGACAGCGCTGGATg                                               >  2:3096888/1‑150 (MQ=255)
                                                                                                      cGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATg                                               <  1:1258173/150‑1 (MQ=255)
                                                                                                      cGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATg                                               <  2:1875292/150‑1 (MQ=255)
                                                                                                      cGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATg                                               <  1:2207815/150‑1 (MQ=255)
                                                                                                         aGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGaca                                            >  1:1979200/1‑150 (MQ=255)
                                                                                                         aGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGaca                                            <  2:1360458/150‑1 (MQ=255)
                                                                                                          gTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGacac                                           <  2:4311005/150‑1 (MQ=255)
                                                                                                           tGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACAcg                                          <  2:1812514/150‑1 (MQ=255)
                                                                                                              gAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGa                                       >  1:4546284/1‑150 (MQ=255)
                                                                                                                 aaaGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTg                                    >  1:1022936/1‑150 (MQ=255)
                                                                                                                 aaaGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTg                                    <  1:1791541/150‑1 (MQ=255)
                                                                                                                    gTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTg                                 <  2:2855831/150‑1 (MQ=255)
                                                                                                                    gTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTg                                 >  1:676702/1‑150 (MQ=255)
                                                                                                                         cGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTa                            >  1:123252/1‑150 (MQ=255)
                                                                                                                          gTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTAc                           >  2:921571/1‑150 (MQ=255)
                                                                                                                           tCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACg                          <  1:2908277/150‑1 (MQ=255)
                                                                                                                            cGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACga                         >  1:2288604/1‑150 (MQ=255)
                                                                                                                            cGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACga                         >  1:2618170/1‑150 (MQ=255)
                                                                                                                             gAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACgaa                        >  1:852366/1‑150 (MQ=255)
                                                                                                                               gcgcATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACgaaga                      <  1:2859802/150‑1 (MQ=255)
                                                                                                                                 gcATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGaagaag                    <  2:2619006/150‑1 (MQ=255)
                                                                                                                                  cATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAagaaga                   >  2:3945295/1‑150 (MQ=255)
                                                                                                                                      gTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcgc               >  1:4626520/1‑150 (MQ=255)
                                                                                                                                       ttCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcgcg              <  2:1396669/150‑1 (MQ=255)
                                                                                                                                       ttCGCCATATCCTCGCGCAGAAAAAACATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcgcg              <  2:4455819/150‑1 (MQ=255)
                                                                                                                                         cGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGAt            <  1:3431518/150‑1 (MQ=255)
                                                                                                                                           ccATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGc          <  1:3613632/150‑1 (MQ=255)
                                                                                                                                             atatCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTg        <  2:3118537/150‑1 (MQ=255)
                                                                                                                                             atatCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTg        <  2:2193930/150‑1 (MQ=255)
                                                                                                                                             atatCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTg        >  1:870153/1‑150 (MQ=255)
                                                                                                                                             atatCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTg        >  1:4471145/1‑150 (MQ=255)
                                                                                                                                             atatCCTCGCGCAGAAAAACAATCCGTCGATCATCATCTGGTAGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATTACACGAGACTGGTGCATTACGAAGAAGATCGCGATGCtc        >  1:3095365/1‑149 (MQ=255)
                                                                                                                                                   aCGCGCAGAACAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGAAGTgg  <  2:4496405/150‑1 (MQ=255)
                                                                                                                                                   aCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGAAGTgg  <  2:1154763/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGAAGTGG  >  CP009273/3217142‑3217438

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: