Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,217,098 C→A 12.1% I369I (ATC→ATA 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,0980CA12.1% 1035.9 / 78.8 297I369I (ATC→ATAebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (126/135);  new base A (19/17);  total (145/152)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 301 total aligned reads displayed.

ATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGC  >  CP009273/3216949‑3217246
                                                                                                                                                     |                                                                                                                                                    
    aaGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACt                                                                                                                                                  <  1:3498812/150‑1 (MQ=255)
       tGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCg                                                                                                                                               >  2:3354689/1‑150 (MQ=255)
        gcgcGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGt                                                                                                                                              <  2:4543039/150‑1 (MQ=255)
         cgcgACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTg                                                                                                                                             >  1:4017836/1‑150 (MQ=255)
         cgcgACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTg                                                                                                                                             >  2:3555924/1‑150 (MQ=255)
              cGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTAc                                                                                                                                        >  1:3588766/1‑150 (MQ=255)
              cGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTAc                                                                                                                                        >  1:3113240/1‑150 (MQ=255)
               ggTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTAcc                                                                                                                                       >  2:2045266/1‑150 (MQ=255)
                    gTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCa                                                                                                                                  <  1:2730146/150‑1 (MQ=255)
                       cTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTa                                                                                                                               >  2:4455346/1‑150 (MQ=255)
                          gATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTAccc                                                                                                                            <  2:838363/150‑1 (MQ=255)
                           aTCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCg                                                                                                                           >  1:2342991/1‑150 (MQ=255)
                           aTCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCg                                                                                                                           >  2:1258156/1‑150 (MQ=255)
                           aTCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCg                                                                                                                           <  2:1721936/150‑1 (MQ=255)
                            tCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGa                                                                                                                          <  2:2367367/150‑1 (MQ=255)
                            tCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGAGGATGAAGCAGCACAATATCAACACCGTGCGTACCGCTCACTCCCCGa                                                                                                                          >  2:1432081/1‑150 (MQ=255)
                              aaTAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAAc                                                                                                                        >  2:2417275/1‑150 (MQ=255)
                               aTAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACg                                                                                                                       <  2:4328464/150‑1 (MQ=255)
                               aTAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACg                                                                                                                       <  2:2227665/150‑1 (MQ=255)
                                 aaCCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGAt                                                                                                                     <  2:4463771/150‑1 (MQ=255)
                                   ccGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATcc                                                                                                                   <  1:4537237/150‑1 (MQ=255)
                                   ccGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATcc                                                                                                                   <  2:3843327/150‑1 (MQ=255)
                                     gTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCcgc                                                                                                                 <  1:3346985/150‑1 (MQ=255)
                                       tATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGt                                                                                                               <  1:2620761/150‑1 (MQ=255)
                                       tATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGGGCGTACCGCTCACTACCCGAACGATCCGCGt                                                                                                               >  1:1205499/1‑150 (MQ=255)
                                        aTGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGtt                                                                                                              <  1:2688053/150‑1 (MQ=255)
                                        aTGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGtt                                                                                                              <  1:2503288/150‑1 (MQ=255)
                                         tgtgATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGttt                                                                                                             >  2:586911/1‑150 (MQ=255)
                                          gtgATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGtttt                                                                                                            <  2:1298495/150‑1 (MQ=255)
                                          gtgATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGtttt                                                                                                            >  1:2125693/1‑150 (MQ=255)
                                          gtgATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGtttt                                                                                                            <  1:3764576/150‑1 (MQ=255)
                                           tgATGCTGCACGGCGTCAACCATCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGttttt                                                                                                           <  2:194202/150‑1 (MQ=255)
                                           tgATGCTGCACGGCGTCAACCATCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGttttt                                                                                                           <  2:2899561/150‑1 (MQ=255)
                                            gATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTa                                                                                                          <  1:1876690/150‑1 (MQ=255)
                                            gATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTa                                                                                                          <  2:750317/150‑1 (MQ=255)
                                             aTGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTAc                                                                                                         >  2:3035010/1‑150 (MQ=255)
                                             aTGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTAc                                                                                                         <  1:575637/150‑1 (MQ=255)
                                             aTGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTAc                                                                                                         <  2:2023905/150‑1 (MQ=255)
                                              tgctgcACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGGAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACg                                                                                                        <  1:1098484/150‑1 (MQ=255)
                                                ctgcACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGaa                                                                                                      <  2:3779518/150‑1 (MQ=255)
                                                  gcACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACt                                                                                                    <  1:2536772/150‑1 (MQ=255)
                                                  gcACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACt                                                                                                    <  1:2301370/150‑1 (MQ=255)
                                                    aCGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgt                                                                                                  <  1:4025097/150‑1 (MQ=255)
                                                    aCGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgt                                                                                                  >  2:1864811/1‑150 (MQ=255)
                                                    aCGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgt                                                                                                  <  1:3056812/150‑1 (MQ=255)
                                                     cGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtg                                                                                                 >  1:2852277/1‑150 (MQ=255)
                                                     cGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtg                                                                                                 >  1:1360401/1‑150 (MQ=255)
                                                      ggCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtgt                                                                                                >  1:4496405/1‑150 (MQ=255)
                                                      ggCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtgt                                                                                                >  1:1154763/1‑150 (MQ=255)
                                                        cGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGa                                                                                              >  1:4455819/1‑150 (MQ=255)
                                                          tCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGata                                                                                            >  1:2911825/1‑150 (MQ=255)
                                                          tCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGata                                                                                            >  1:726821/1‑150 (MQ=255)
                                                          tCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGata                                                                                            >  1:3160117/1‑150 (MQ=255)
                                                           cAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGatat                                                                                           >  2:4166967/1‑150 (MQ=255)
                                                           cAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGatat                                                                                           >  2:3978895/1‑150 (MQ=255)
                                                           cAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGatat                                                                                           <  1:4144135/150‑1 (MQ=255)
                                                            aaCCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATc                                                                                          >  2:3714591/1‑150 (MQ=255)
                                                            aaCCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATc                                                                                          >  1:3208903/1‑150 (MQ=255)
                                                                 tCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACgg                                                                                     >  2:4434047/1‑150 (MQ=255)
                                                                           gATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTtgtg                                                                           <  1:3317841/150‑1 (MQ=255)
                                                                              catcGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATg                                                                        <  2:1179811/150‑1 (MQ=255)
                                                                               atcGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATgg                                                                       >  1:48799/1‑150 (MQ=255)
                                                                                   cAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATggagga                                                                   >  1:2979972/1‑150 (MQ=255)
                                                                                    aaaGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGaa                                                                  <  2:294713/150‑1 (MQ=255)
                                                                                      aGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAAc                                                                <  1:2269930/150‑1 (MQ=255)
                                                                                      aGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAAc                                                                <  1:1043203/150‑1 (MQ=255)
                                                                                       ggCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAAcc                                                               >  2:2207815/1‑150 (MQ=255)
                                                                                          cgcgCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGAc                                                            <  2:4240816/150‑1 (MQ=255)
                                                                                             gcCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTc                                                         <  2:4521824/150‑1 (MQ=255)
                                                                                             gcCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTc                                                         <  2:3741033/150‑1 (MQ=255)
                                                                                                gTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGaa                                                      <  2:468453/150‑1 (MQ=255)
                                                                                                  tGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATc                                                    <  2:4379174/150‑1 (MQ=255)
                                                                                                             gcgTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTg                                         <  2:3113240/150‑1 (MQ=255)
                                                                                                               gTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCt                                       <  1:2409392/150‑1 (MQ=255)
                                                                                                                   agaAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATg                                   <  1:2582485/150‑1 (MQ=255)
                                                                                                                     aaaGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTc                                 >  1:868684/1‑150 (MQ=255)
                                                                                                                      aataTCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCg                                <  2:663824/147‑1 (MQ=255)
                                                                                                                      aaGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCg                                <  1:3180257/150‑1 (MQ=255)
                                                                                                                      aaGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCg                                <  1:4455253/150‑1 (MQ=255)
                                                                                                                       aGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCgg                               <  1:2892220/150‑1 (MQ=255)
                                                                                                                         aTCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCg                             >  2:1365717/1‑150 (MQ=255)
                                                                                                                           ctcCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGat                           >  1:3675385/1‑150 (MQ=255)
                                                                                                                              cAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATAtt                        <  2:2826109/150‑1 (MQ=255)
                                                                                                                                  tgatgaAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGcc                    <  2:4472344/150‑1 (MQ=255)
                                                                                                                                   gatgaAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGTGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCg                   <  1:697138/150‑1 (MQ=255)
                                                                                                                                       aaGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTAtt               >  1:2772856/1‑150 (MQ=255)
                                                                                                                                        agcagcACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAAATGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTa              <  2:364468/150‑1 (MQ=255)
                                                                                                                                         gcagcaCAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGAGATTAGCCGTATTAc             <  2:1068994/150‑1 (MQ=255)
                                                                                                                                          cagcaCAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAcc            <  1:2211562/150‑1 (MQ=255)
                                                                                                                                            gcaCAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcga          >  1:1812514/1‑150 (MQ=255)
                                                                                                                                             cacaATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac         >  2:32098/1‑150 (MQ=255)
                                                                                                                                             cacaATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac         <  1:775945/150‑1 (MQ=255)
                                                                                                                                             cacaATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac         >  1:1477228/1‑150 (MQ=255)
                                                                                                                                                aaTATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGAt      >  2:1781397/1‑150 (MQ=255)
                                                                                                                                                 atatCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATc     <  2:3907929/150‑1 (MQ=255)
                                                                                                                                                 atatCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATc     >  1:1669119/1‑150 (MQ=255)
                                                                                                                                                  tataAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATcc    >  1:389225/1‑150 (MQ=255)
                                                                                                                                                  tatCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATcc    <  1:2816768/150‑1 (MQ=255)
                                                                                                                                                   atCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCg   >  1:930891/1‑150 (MQ=255)
                                                                                                                                                   atCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCg   <  1:4244257/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGC  >  CP009273/3216949‑3217246

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: