Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,767,000 G→T 90.2% D454Y (GAC→TAC)  mtlA → fused mannitol‑specific PTS enzymes: IIA components/IIB components/IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,767,0000GT90.2% 781.6 / 78.9 282D454Y (GAC→TAC) mtlAfused mannitol‑specific PTS enzymes: IIA components/IIB components/IIC components
Reads supporting (aligned to +/- strand):  ref base G (17/11);  new base T (128/126);  total (145/137)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

GCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTC  >  CP009273/3766853‑3767149
                                                                                                                                                   |                                                                                                                                                     
gCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTAc                                                                                                                                                     <  1:294986/150‑1 (MQ=255)
 cAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACa                                                                                                                                                    >  1:1991742/1‑150 (MQ=255)
 cAGGTCTGTCGCAGATTTCTGATACTAACCGCGCGATCAACAAACTGCCGCCAGATGTGGACCTCGTCATCACACACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACa                                                                                                                                                    >  2:3446876/1‑150 (MQ=255)
  aGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAg                                                                                                                                                   >  2:1148419/1‑150 (MQ=255)
    gTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCg                                                                                                                                                 <  2:2970400/150‑1 (MQ=255)
         tCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTg                                                                                                                                            >  1:813374/1‑150 (MQ=255)
         tCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTg                                                                                                                                            >  1:4268662/1‑150 (MQ=255)
           gCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTa                                                                                                                                          >  2:823997/1‑150 (MQ=255)
            cAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTac                                                                                                                                         >  2:3155275/1‑150 (MQ=255)
               aTTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACAcc                                                                                                                                      >  2:1689695/1‑150 (MQ=255)
                tttCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCa                                                                                                                                     <  2:542867/150‑1 (MQ=255)
                       tACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGAc                                                                                                                              <  1:1034642/150‑1 (MQ=255)
                       tACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGAc                                                                                                                              >  1:3355162/1‑150 (MQ=255)
                        aCTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGAcc                                                                                                                             >  2:412613/1‑150 (MQ=255)
                             cAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACg                                                                                                                        <  2:3349378/150‑1 (MQ=255)
                              aGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGt                                                                                                                       >  1:748164/1‑150 (MQ=255)
                                cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCt                                                                                                                     <  2:612045/150‑1 (MQ=255)
                                cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCt                                                                                                                     >  1:2712876/1‑150 (MQ=255)
                                  cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTgg                                                                                                                   >  2:4370018/1‑150 (MQ=255)
                                  cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTgg                                                                                                                   >  1:2823573/1‑150 (MQ=255)
                                  cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTgg                                                                                                                   <  1:3017835/150‑1 (MQ=255)
                                   gATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGt                                                                                                                  <  1:3721469/150‑1 (MQ=255)
                                    aTCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGtt                                                                                                                 <  2:3431464/150‑1 (MQ=255)
                                     tCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtg                                                                                                                <  2:617044/150‑1 (MQ=255)
                                       aacaacCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtgct                                                                                                              <  1:1550698/150‑1 (MQ=255)
                                       aacaacCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtgct                                                                                                              <  1:3182569/150‑1 (MQ=255)
                                           acCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGccc                                                                                                          <  2:2524333/150‑1 (MQ=255)
                                           acCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGccc                                                                                                          <  1:819798/150‑1 (MQ=255)
                                            ccTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCGTATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCa                                                                                                         <  1:852031/150‑1 (MQ=255)
                                              tGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGACCGTCTGGTTGCTGCCCAAc                                                                                                       <  2:4195314/150‑1 (MQ=255)
                                               gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACg                                                                                                      >  2:723170/1‑150 (MQ=255)
                                                 cgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGcc                                                                                                    <  1:4215883/150‑1 (MQ=255)
                                                  gccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCca                                                                                                   >  1:1974058/1‑150 (MQ=255)
                                                     aGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCacac                                                                                                <  2:3878811/150‑1 (MQ=255)
                                                     aGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCacac                                                                                                >  2:1013925/1‑150 (MQ=255)
                                                     aGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCacac                                                                                                <  2:1562041/150‑1 (MQ=255)
                                                     aGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCacac                                                                                                >  1:312754/1‑150 (MQ=255)
                                                      gATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACg                                                                                               <  1:2518029/150‑1 (MQ=255)
                                                       aTGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACgg                                                                                              <  1:3215149/150‑1 (MQ=255)
                                                         gtgGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCa                                                                                            <  1:379467/150‑1 (MQ=255)
                                                          tgGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaa                                                                                           >  1:3557222/1‑150 (MQ=255)
                                                          tgGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaa                                                                                           >  2:152023/1‑150 (MQ=255)
                                                           ggACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaaa                                                                                          <  2:2732940/150‑1 (MQ=255)
                                                             aCCTCGTCATCACTCACCGTGACCTGACCGAACGCGCGATGCGCCAGGTTCCGCAGGCCCAGCATGTTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCGGCCTGACAGAAAGACTGGTAGCTGCCCAACGCCACACGGAAAACg                                                                                        >  2:899516/1‑150 (MQ=255)
                                                                  gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAagag                                                                                   <  2:4494688/150‑1 (MQ=255)
                                                                   tcatcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAgaga                                                                                  >  1:2926042/1‑150 (MQ=255)
                                                                    catcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGaa                                                                                 >  1:1682768/1‑150 (MQ=255)
                                                                     atcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCTCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGaaa                                                                                >  2:1770338/1‑150 (MQ=255)
                                                                       cactcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGt                                                                              <  2:502630/150‑1 (MQ=255)
                                                                        actcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTa                                                                             <  2:2076081/150‑1 (MQ=255)
                                                                          tcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaaa                                                                           >  1:436223/1‑150 (MQ=255)
                                                                            acCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAg                                                                         <  2:2377105/150‑1 (MQ=255)
                                                                              cGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAATGTAAAAGAc                                                                       <  1:358621/150‑1 (MQ=255)
                                                                              cGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGAc                                                                       <  2:2260306/150‑1 (MQ=255)
                                                                              cGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGAc                                                                       >  2:3275187/1‑150 (MQ=255)
                                                                              cGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGAc                                                                       <  2:3599132/150‑1 (MQ=255)
                                                                               gTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACa                                                                      >  1:1068566/1‑150 (MQ=255)
                                                                               gTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACa                                                                      >  1:1839133/1‑150 (MQ=255)
                                                                                tgacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAg                                                                     >  1:2015950/1‑150 (MQ=255)
                                                                                tgacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAg                                                                     <  1:3406193/150‑1 (MQ=255)
                                                                                 gacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGc                                                                    <  1:385840/150‑1 (MQ=255)
                                                                                 gacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGc                                                                    <  1:121126/150‑1 (MQ=255)
                                                                                 gacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGc                                                                    >  1:2334769/1‑150 (MQ=255)
                                                                                   cctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCt                                                                  <  2:4062183/150‑1 (MQ=255)
                                                                                    ctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTg                                                                 <  1:3248331/150‑1 (MQ=255)
                                                                                      gaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGaa                                                               >  1:1904919/1‑150 (MQ=255)
                                                                                      gaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGaa                                                               >  1:1539314/1‑150 (MQ=255)
                                                                                      gaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGaa                                                               >  1:4582272/1‑150 (MQ=255)
                                                                                      gaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACATCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGaa                                                               >  1:4404206/1‑150 (MQ=255)
                                                                                          gAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGAc                                                           >  1:280565/1‑150 (MQ=255)
                                                                                          gAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGAc                                                           >  1:2413519/1‑150 (MQ=255)
                                                                                           aaCGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACa                                                          <  1:2087461/150‑1 (MQ=255)
                                                                                           aaCGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACa                                                          <  2:1099088/150‑1 (MQ=255)
                                                                                            aCGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAg                                                         <  1:2332846/150‑1 (MQ=255)
                                                                                              gcgcTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCt                                                       >  1:644959/1‑150 (MQ=255)
                                                                                               cgcTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCtt                                                      <  1:921848/150‑1 (MQ=255)
                                                                                                 cTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTg                                                    <  1:2126410/150‑1 (MQ=255)
                                                                                                    tGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACg                                                 <  1:1803779/150‑1 (MQ=255)
                                                                                                      cgcCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGAt                                               >  2:1669837/1‑150 (MQ=255)
                                                                                                        ccAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTc                                             >  2:3234591/1‑150 (MQ=255)
                                                                                                        ccAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTc                                             >  1:2568713/1‑150 (MQ=255)
                                                                                                          aGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCa                                           <  2:3217517/150‑1 (MQ=255)
                                                                                                          aGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCa                                           >  1:4439600/1‑150 (MQ=255)
                                                                                                                cGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTa                                     <  2:3049325/150‑1 (MQ=255)
                                                                                                                   aGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAAcc                                  >  2:2476627/1‑150 (MQ=255)
                                                                                                                    ggCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCt                                 <  1:2397137/150‑1 (MQ=255)
                                                                                                                       acaGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGtt                              <  2:3731731/150‑1 (MQ=255)
                                                                                                                       acaGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGtt                              <  1:3610087/150‑1 (MQ=255)
                                                                                                                       acaGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGtt                              <  2:2162833/150‑1 (MQ=255)
                                                                                                                        caGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTc                             <  2:1209153/150‑1 (MQ=255)
                                                                                                                              atTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTa                       >  1:680431/1‑150 (MQ=255)
                                                                                                                                ttCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAgg                     <  1:2929249/150‑1 (MQ=255)
                                                                                                                                  cGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGgcg                   >  1:4630332/1‑150 (MQ=255)
                                                                                                                                         ccaCTTCCTCTACAGCGGCCTGTACCCCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGaa            <  1:3948552/148‑1 (MQ=255)
                                                                                                                                         cAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGaa            >  2:4566705/1‑150 (MQ=255)
                                                                                                                                         cAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCGGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGaa            >  2:1825798/1‑150 (MQ=255)
                                                                                                                                          aaCTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAAc           >  2:65350/1‑150 (MQ=255)
                                                                                                                                           aCTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACa          <  2:2352282/150‑1 (MQ=255)
                                                                                                                                            cTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACAt         >  2:2556649/1‑150 (MQ=255)
                                                                                                                                             ttCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATc        <  1:880496/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCCTC  >  CP009273/3766853‑3767149

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: