Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,767,000 G→T 100% 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,0000GT99.7% 994.4 / ‑4.5 292D454Y (GAC→TAC) mtlAfused mannitol‑specific PTS enzymes: IIA components/IIB components/IIC components
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (135/156);  total (135/157)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 297 total aligned reads displayed.

CAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTT  >  CP009273/3766854‑3767145
                                                                                                                                                  |                                                                                                                                                 
 aGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACa                                                                                                                                                >  1:1343448/1‑149 (MQ=255)
 aGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACa                                                                                                                                                >  1:2420880/1‑149 (MQ=255)
   gTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGc                                                                                                                                              >  1:1400475/1‑149 (MQ=255)
       gTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGcc                                                                                                                                          <  1:4720624/149‑1 (MQ=255)
          gCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGt                                                                                                                                       >  2:1507044/1‑149 (MQ=255)
          gCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGt                                                                                                                                       >  2:3895653/1‑149 (MQ=255)
             gATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTaca                                                                                                                                    >  1:317097/1‑149 (MQ=255)
               tttCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACAcc                                                                                                                                  >  1:46395/1‑149 (MQ=255)
               tttCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACAcc                                                                                                                                  >  2:3036740/1‑149 (MQ=255)
               tttCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACAcc                                                                                                                                  >  1:987730/1‑149 (MQ=255)
                    gTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCt                                                                                                                             <  2:4205408/149‑1 (MQ=255)
                    gTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCt                                                                                                                             <  1:4223700/149‑1 (MQ=255)
                     ttACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTg                                                                                                                            >  1:2632878/1‑149 (MQ=255)
                     ttACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTg                                                                                                                            <  2:1702633/149‑1 (MQ=255)
                      tACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGa                                                                                                                           <  1:656568/149‑1 (MQ=255)
                      tACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGa                                                                                                                           >  1:1462616/1‑149 (MQ=255)
                        cTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGAcc                                                                                                                         >  2:1980676/1‑149 (MQ=255)
                          aaCAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGa                                                                                                                       >  1:3887630/1‑149 (MQ=255)
                               cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTc                                                                                                                  <  2:216884/149‑1 (MQ=255)
                                gcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCt                                                                                                                 >  2:672076/1‑149 (MQ=255)
                                gcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCt                                                                                                                 <  2:3113159/149‑1 (MQ=255)
                                gcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCt                                                                                                                 >  2:4386294/1‑149 (MQ=255)
                                  gATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTgg                                                                                                               <  1:3016283/149‑1 (MQ=255)
                                   aTCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGt                                                                                                              <  2:3679467/149‑1 (MQ=255)
                                   aTCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGt                                                                                                              <  2:3195179/149‑1 (MQ=255)
                                   aTCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGt                                                                                                              <  1:1693888/149‑1 (MQ=255)
                                      aacaacCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtgc                                                                                                           <  2:783324/149‑1 (MQ=255)
                                       acaacCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtgct                                                                                                          <  2:4604966/149‑1 (MQ=255)
                                           ccTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGccc                                                                                                      <  2:3779591/149‑1 (MQ=255)
                                           ccTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGccc                                                                                                      >  2:317284/1‑149 (MQ=255)
                                             tGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCaa                                                                                                    <  2:57817/149‑1 (MQ=255)
                                                 gccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGcc                                                                                                >  2:2517966/1‑149 (MQ=255)
                                                  ccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCca                                                                                               >  2:1164254/1‑149 (MQ=255)
                                                  ccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCca                                                                                               >  1:40869/1‑149 (MQ=255)
                                                  ccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCca                                                                                               <  1:2314857/149‑1 (MQ=255)
                                                   cAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCcac                                                                                              >  1:4153580/1‑149 (MQ=255)
                                                     gATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCacac                                                                                            >  2:2044852/1‑149 (MQ=255)
                                                      aTGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACg                                                                                           >  1:1027541/1‑149 (MQ=255)
                                                         tgGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCa                                                                                        <  1:523695/149‑1 (MQ=255)
                                                          ggACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaa                                                                                       >  2:294833/1‑149 (MQ=255)
                                                          ggACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaa                                                                                       >  1:128832/1‑149 (MQ=255)
                                                           gACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaaa                                                                                      >  1:4474764/1‑149 (MQ=255)
                                                           gACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaaa                                                                                      <  2:2013284/149‑1 (MQ=255)
                                                           gACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaaa                                                                                      >  1:2913236/1‑149 (MQ=255)
                                                            aCCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAAc                                                                                     <  1:4749583/149‑1 (MQ=255)
                                                            aCCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAAc                                                                                     <  1:3292814/149‑1 (MQ=255)
                                                             ccTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACg                                                                                    >  2:2306545/1‑149 (MQ=255)
                                                             ccTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACg                                                                                    >  1:1621069/1‑149 (MQ=255)
                                                               tCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGaa                                                                                  <  1:2743985/149‑1 (MQ=255)
                                                               tCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGaa                                                                                  <  1:4376647/149‑1 (MQ=255)
                                                               tCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGaa                                                                                  >  2:929990/1‑149 (MQ=255)
                                                                 gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAaga                                                                                <  2:4071419/149‑1 (MQ=255)
                                                                  tcatcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAagag                                                                               >  1:1957194/1‑149 (MQ=255)
                                                                  tcatcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAagag                                                                               >  2:4041513/1‑149 (MQ=255)
                                                                         tcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaa                                                                        >  1:3123100/1‑149 (MQ=255)
                                                                           acCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaaaa                                                                      <  2:1294671/149‑1 (MQ=255)
                                                                            ccGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAg                                                                     <  2:1463012/149‑1 (MQ=255)
                                                                             cGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGa                                                                    <  2:945113/149‑1 (MQ=255)
                                                                              gTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCACACGAAGAGAAAGTAAAAGAc                                                                   >  2:2881717/1‑149 (MQ=255)
                                                                              gTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGAc                                                                   >  2:1261852/1‑149 (MQ=255)
                                                                               tgacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACa                                                                  <  2:1326755/149‑1 (MQ=255)
                                                                               tgacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACa                                                                  <  1:1838097/149‑1 (MQ=255)
                                                                                gacctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAg                                                                 <  1:1558727/149‑1 (MQ=255)
                                                                                     gaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGa                                                            <  1:2164177/149‑1 (MQ=255)
                                                                                          aaCGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGAc                                                       <  2:3305884/149‑1 (MQ=255)
                                                                                            cgcgCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAg                                                     <  1:1172650/149‑1 (MQ=255)
                                                                                             gcgcTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGc                                                    <  2:4703566/149‑1 (MQ=255)
                                                                                               gcTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCtt                                                  >  1:4305959/1‑149 (MQ=255)
                                                                                                   tGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGAc                                              >  1:2687801/1‑149 (MQ=255)
                                                                                                       ccAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGAtt                                          >  1:2471761/1‑149 (MQ=255)
                                                                                                          ggTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCa                                       <  2:4336271/149‑1 (MQ=255)
                                                                                                           gTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAg                                      >  2:4490653/1‑149 (MQ=255)
                                                                                                              ccGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGc                                   >  2:2873015/1‑149 (MQ=255)
                                                                                                               ctcAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCt                                  >  1:3747288/3‑149 (MQ=255)
                                                                                                               cGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCt                                  <  2:1373207/149‑1 (MQ=255)
                                                                                                                gCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTa                                 >  1:3937231/1‑149 (MQ=255)
                                                                                                                gCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTa                                 >  1:3091063/1‑149 (MQ=255)
                                                                                                                  aGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAAc                               <  2:4669360/149‑1 (MQ=255)
                                                                                                                   ggCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAAcc                              >  2:4107604/1‑149 (MQ=255)
                                                                                                                    gCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCt                             >  1:4328341/1‑149 (MQ=255)
                                                                                                                     cacaGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTg                            >  1:916659/1‑149 (MQ=255)
                                                                                                                      acaGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACTCCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt                           <  2:2103803/149‑1 (MQ=255)
                                                                                                                        aGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTc                         <  2:712559/149‑1 (MQ=255)
                                                                                                                         gCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCa                        <  2:4561771/149‑1 (MQ=255)
                                                                                                                          cATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCaa                       <  2:275422/149‑1 (MQ=255)
                                                                                                                             atTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCt                    >  1:4661064/1‑149 (MQ=255)
                                                                                                                                tCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAgg                 <  2:3490154/149‑1 (MQ=255)
                                                                                                                                 cGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGgc                <  2:1066956/149‑1 (MQ=255)
                                                                                                                                 cGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGgc                <  1:2460055/149‑1 (MQ=255)
                                                                                                                                  gCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGgcg               >  1:4526227/1‑149 (MQ=255)
                                                                                                                                     gACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCGTTGACGATTCCAGTGCTAACCTGTTCAAGCTCGGCGCgg            >  2:3569829/1‑149 (MQ=255)
                                                                                                                                      aCCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGga           <  2:1834800/149‑1 (MQ=255)
                                                                                                                                        cAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGgaga         <  1:4163220/149‑1 (MQ=255)
                                                                                                                                         aaCTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGaa        >  1:3849860/1‑149 (MQ=255)
                                                                                                                                          aCTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAAc       <  1:2676537/149‑1 (MQ=255)
                                                                                                                                           cTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACa      >  2:2657960/1‑149 (MQ=255)
                                                                                                                                           cTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACa      >  1:1881388/1‑149 (MQ=255)
                                                                                                                                            ttCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACAt     >  2:2732340/1‑149 (MQ=255)
                                                                                                                                              ccTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCt   <  1:3532892/149‑1 (MQ=255)
                                                                                                                                               ctctACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCtt  >  2:2713441/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTT  >  CP009273/3766854‑3767145

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