Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,224,895 C→T L755L (CTG→TTG)  ebgA → evolved beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,224,8950CT96.1% 1348.9 / 27.5 417L755L (CTG→TTG) ebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (11/5);  new base T (194/205);  total (206/211)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

Only 100 of 428 total aligned reads displayed.

GCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGC  >  NC_000913/3224747‑3225042
                                                                                                                                                    |                                                                                                                                                   
gCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGctgttg                                                                                                                                                   <  7:812053/151‑1 (MQ=255)
    gTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTc                                                                                                                                               >  8:653246/1‑151 (MQ=255)
     tGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTcg                                                                                                                                              >  2:1635659/1‑151 (MQ=255)
      gCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTcgc                                                                                                                                             <  7:171030/151‑1 (MQ=255)
      gCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTcgc                                                                                                                                             <  8:1524341/151‑1 (MQ=255)
          ttCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCCATCACCTTCTCAAAAATGAGTGGCAGACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGCCTCGAGAGc                                                                                                                                         >  6:1400529/1‑151 (MQ=255)
             gCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCaa                                                                                                                                      <  3:201700/151‑1 (MQ=255)
              cACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCaaa                                                                                                                                     <  1:1148140/151‑1 (MQ=255)
              cACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTATGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCAGCTTCTCAAAAATGAGTGGCAACCCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCaa                                                                                                                                      >  4:1169056/1‑150 (MQ=255)
               aCCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCCCCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCaa                                                                                                                                      >  6:29564/1‑149 (MQ=255)
               aCCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCaa                                                                                                                                      >  8:566632/1‑149 (MQ=255)
                 cAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCACACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCaa                                                                                                                                      >  1:1199563/1‑147 (MQ=255)
                   aataatGCGCGTCTGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCa                                                                                                                                <  3:752441/151‑1 (MQ=255)
                   aataatGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATc                                                                                                                                 >  1:1596250/1‑150 (MQ=255)
                      aatGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAAct                                                                                                                             >  3:580195/1‑151 (MQ=255)
                      aatGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAAct                                                                                                                             <  7:1510778/151‑1 (MQ=255)
                         gcgcGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAActtct                                                                                                                          >  6:1198047/1‑151 (MQ=255)
                             gTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCaa                                                                                                                      <  7:583550/151‑1 (MQ=255)
                               ccGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGc                                                                                                                    <  1:608475/151‑1 (MQ=255)
                               ccGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGc                                                                                                                    >  8:775331/1‑151 (MQ=255)
                                cGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGcc                                                                                                                   >  8:1047635/1‑151 (MQ=255)
                                 gCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCg                                                                                                                  >  8:272909/1‑151 (MQ=255)
                                     aCGCTGTATGACGATCGTTTGAGCTGCACCGTTCTCGGCTACTACTTGGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATGGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCgatg                                                                                                               <  6:402128/150‑1 (MQ=255)
                                      cGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGTCGAATAGCTGTTGAGTCGCGAGCCAAAGATCAACTTATTCAAGCCgatgat                                                                                                             >  6:93559/1‑151 (MQ=255)
                                      cGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCgatgat                                                                                                             >  8:1229481/1‑151 (MQ=255)
                                       gCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATc                                                                                                            >  3:464215/1‑151 (MQ=255)
                                       gCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATc                                                                                                            <  1:1305508/151‑1 (MQ=255)
                                         tGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACCTCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCgc                                                                                                          >  2:1383791/1‑150 (MQ=255)
                                          ggAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGAc                                                                                                         >  4:733137/1‑151 (MQ=255)
                                           gAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGCCTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGAc                                                                                                         >  4:6717/1‑150 (MQ=255)
                                            aaGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACa                                                                                                        >  4:1578643/1‑150 (MQ=255)
                                              gACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAAcc                                                                                                     <  2:557844/151‑1 (MQ=255)
                                               aCGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAAcc                                                                                                     >  4:397637/1‑150 (MQ=255)
                                                cGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACcac                                                                                                   <  7:467063/151‑1 (MQ=255)
                                                cGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACcac                                                                                                   <  8:1148177/151‑1 (MQ=255)
                                                 gATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACcaca                                                                                                  <  4:565970/151‑1 (MQ=255)
                                                    cGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAg                                                                                                <  3:1449817/150‑1 (MQ=255)
                                                       ttGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAgg                                                                                            >  3:125961/1‑151 (MQ=255)
                                                        tGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAgg                                                                                            >  2:391845/1‑150 (MQ=255)
                                                         gAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAg                                                                                          <  8:898880/151‑1 (MQ=255)
                                                           gCTGCACCGTTCTCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGGAAACCGACCTCCTTGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAACATCAACTTCTTCAAGGCGATGATCGACAACCACAAGCAgtagta                                                                                        <  2:331925/151‑1 (MQ=255)
                                                            cTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTAc                                                                                       >  5:1129885/1‑151 (MQ=255)
                                                              gCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTAc                                                                                       <  4:276370/149‑1 (MQ=255)
                                                               cACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGaa                                                                                    <  1:355935/151‑1 (MQ=255)
                                                               cACCGTTCGCGGCTACAACAACGCGATCACCATCTCAAAAATGAGTGGCAAACCGACATCCTGGCACGAGAATGCCGAATCGCACAAGACTCGCGAGCCAAAGATCAACTTCTTCAACCCGATGATCGACAACCACAAGCAGGACAACGa                                                                                     >  2:837618/1‑150 (MQ=255)
                                                                 ccGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAgg                                                                                  >  4:797113/1‑151 (MQ=255)
                                                                 ccGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAgg                                                                                  <  2:1241153/151‑1 (MQ=255)
                                                                     tCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCtg                                                                              >  2:1666621/1‑151 (MQ=255)
                                                                     tCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCt                                                                               <  8:44435/150‑1 (MQ=255)
                                                                       gcgGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAACCCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGCCTCGCGAGCCAAAGATCAACATCTTCAAGCCGATGATCGACAACCACACCCAGGAGTCCGAAGGGCggtg                                                                            >  5:1460395/1‑151 (MQ=255)
                                                                       gcgGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCtgtg                                                                            >  5:1536137/1‑151 (MQ=255)
                                                                       gcgGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCtgtg                                                                            >  5:1663718/1‑151 (MQ=255)
                                                                         ggCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGc                                                                          <  2:186806/151‑1 (MQ=255)
                                                                         ggCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGc                                                                          <  5:506772/151‑1 (MQ=255)
                                                                             aCAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAAcc                                                                      <  1:1520580/151‑1 (MQ=255)
                                                                               aaCTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCg                                                                     >  6:158247/1‑150 (MQ=255)
                                                                               aaCTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCg                                                                     >  7:298219/1‑150 (MQ=255)
                                                                                    cgcgATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGCTCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCAt                                                               >  1:879770/1‑151 (MQ=255)
                                                                                    cgcgATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCAt                                                               <  7:1337882/151‑1 (MQ=255)
                                                                                     gcgATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGATTCAtt                                                              >  6:1055806/1‑151 (MQ=255)
                                                                                     gcgATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCAtt                                                              >  1:1021713/1‑151 (MQ=255)
                                                                                      cgATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCAttt                                                             >  7:627573/1‑151 (MQ=255)
                                                                                      cgATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCAttt                                                             >  2:1083240/1‑151 (MQ=255)
                                                                                       gATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTg                                                            <  3:373284/151‑1 (MQ=255)
                                                                                        aTCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                           <  1:957258/151‑1 (MQ=255)
                                                                                         tCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCGAAAGATCAACTTCTTCAAGCCGATGATCGCGAACCTCAAGCAGGAGTCCGAAGGGCTGTGGCGCCCGAATCAGTTGcg                                                          >  2:1428699/1‑150 (MQ=255)
                                                                                          cACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                           >  1:240546/1‑149 (MQ=255)
                                                                                            ccTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGAt                                                       <  5:1039675/151‑1 (MQ=255)
                                                                                                ctcAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATAAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCACCCGCATCCTTTGCAGATCATg                                                   >  2:611562/1‑151 (MQ=255)
                                                                                                  cAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCAt                                                    <  8:1452747/148‑1 (MQ=255)
                                                                                                  cAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGc                                                  <  6:788236/150‑1 (MQ=255)
                                                                                                  cAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCa                                                 <  7:730400/151‑1 (MQ=255)
                                                                                                    aaaaTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGCGCCAAAGATCAACTTATTCACGCCGATGAACGCCACCCACAAGCAGGAGTACGAAGGCATGTGCCAACCGAATCATTTGACGATCATGCAgg                                               >  3:885656/1‑151 (MQ=255)
                                                                                                      aaTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGaa                                             <  6:20824/151‑1 (MQ=255)
                                                                                                       aTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAAc                                            >  8:1618312/1‑151 (MQ=255)
                                                                                                         gAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACAt                                          <  3:406880/151‑1 (MQ=255)
                                                                                                         gAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAAGCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACAt                                          <  2:971750/151‑1 (MQ=255)
                                                                                                               cAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTgcgc                                    <  1:451221/151‑1 (MQ=255)
                                                                                                                aaaCCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGcgcg                                   <  2:714299/151‑1 (MQ=255)
                                                                                                                  aCCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGAc                                 <  6:1687448/151‑1 (MQ=255)
                                                                                                                     gACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACttt                              <  1:1450480/151‑1 (MQ=255)
                                                                                                                     gACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCAAAAGATCAACTTCTTCAAGCCGATGATCGACACCCACAACCAGGAGTACGAAGGGCTGTGGCAACCCAATCATTTGCAGAACATTCAGGAACATATGCGCGACttt                              >  4:143368/1‑151 (MQ=255)
                                                                                                                       cATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGc                            >  1:765679/1‑151 (MQ=255)
                                                                                                                        aTCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGcc                           <  1:497987/151‑1 (MQ=255)
                                                                                                                         tCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCg                          >  5:304824/1‑151 (MQ=255)
                                                                                                                            tGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGt                         >  1:1060219/1‑149 (MQ=255)
                                                                                                                             ggCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAg                       <  3:1581546/150‑1 (MQ=255)
                                                                                                                             ggCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGa                      <  3:989709/151‑1 (MQ=255)
                                                                                                                             ggCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGa                      <  7:1102876/151‑1 (MQ=255)
                                                                                                                                aGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCACGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACa                   <  2:780505/151‑1 (MQ=255)
                                                                                                                                 ggTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACag                  <  1:1678425/151‑1 (MQ=255)
                                                                                                                                    gAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACagag                <  5:765514/150‑1 (MQ=255)
                                                                                                                                    gAATGGCGAATCGCTGCTGACTTGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGc               >  3:1158421/1‑151 (MQ=255)
                                                                                                                                     aaTGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGca              <  7:1029979/151‑2 (MQ=255)
                                                                                                                                      aTGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCg              >  6:49360/1‑150 (MQ=255)
                                                                                                                                       tGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGAt            >  5:1602610/1‑151 (MQ=255)
                                                                                                                                           gAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTg        >  2:154997/1‑151 (MQ=255)
                                                                                                                                            aaTCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCTCGACTTTGCCGTAGAACAGAGCGATGGTg        >  5:613557/1‑150 (MQ=255)
                                                                                                                                               cGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGt    >  6:1521475/1‑151 (MQ=255)
                                                                                                                                               cGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGt    >  5:874763/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGC  >  NC_000913/3224747‑3225042

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