Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,767,000 G→T 51.6% 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,0000GT51.6% 80.4 / inf 279D454Y (GAC→TAC) mtlAfused mannitol‑specific PTS enzymes: IIA components/IIB components/IIC components
Reads supporting (aligned to +/- strand):  ref base G (62/73);  new base T (76/68);  total (138/141)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 281 total aligned reads displayed.

ATGCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCC  >  CP009273/3766851‑3767147
                                                                                                                                                     |                                                                                                                                                   
aTGCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCg                                                                                                                                                     <  1:1473295/150‑1 (MQ=255)
   cAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACa                                                                                                                                                  >  2:3505202/1‑150 (MQ=255)
    aGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAg                                                                                                                                                 >  1:742779/1‑150 (MQ=255)
      gTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                                                                               <  2:1783872/150‑1 (MQ=255)
       tCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCgg                                                                                                                                              <  2:3099345/150‑1 (MQ=255)
        cTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGc                                                                                                                                             <  2:614385/150‑1 (MQ=255)
        cTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGc                                                                                                                                             >  1:1684035/1‑150 (MQ=255)
         tGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGcc                                                                                                                                            <  1:791429/150‑1 (MQ=255)
          gTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCt                                                                                                                                           <  1:2585249/150‑1 (MQ=255)
          gTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCt                                                                                                                                           <  1:1826754/150‑1 (MQ=255)
          gTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCt                                                                                                                                           >  1:3304025/1‑150 (MQ=255)
           tCGCAGATTTCTGTTACTAACAGCGCGATCAACTACCTGCCGAGAGTTGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGAGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGTCCATCTTCCTCGACAGCGGCCTg                                                                                                                                          <  2:3829796/150‑1 (MQ=255)
             gCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTa                                                                                                                                        <  1:195945/150‑1 (MQ=255)
              cAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTac                                                                                                                                       >  2:2416153/1‑150 (MQ=255)
                  tttCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCa                                                                                                                                   >  2:4046446/1‑150 (MQ=255)
                      tGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCt                                                                                                                               <  2:657991/150‑1 (MQ=255)
                      tGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCt                                                                                                                               <  1:1386228/150‑1 (MQ=255)
                              aCAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAAc                                                                                                                       <  1:1382082/150‑1 (MQ=255)
                                   gcgATCAACAACCTACCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTg                                                                                                                  <  1:597744/150‑1 (MQ=255)
                                     gATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGt                                                                                                                >  1:7222/1‑150 (MQ=255)
                                     gATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGt                                                                                                                <  1:2974558/150‑1 (MQ=255)
                                      aTCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGtt                                                                                                               >  1:2948047/1‑150 (MQ=255)
                                      aTCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGtt                                                                                                               >  1:2004505/1‑150 (MQ=255)
                                      aTCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGtt                                                                                                               <  2:1394112/150‑1 (MQ=255)
                                       tCAACAACCTGCCGCCAGATGTGGACCTCGACATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtg                                                                                                              >  1:3861737/1‑150 (MQ=255)
                                        caacaaCCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtgc                                                                                                             >  1:3996919/1‑150 (MQ=255)
                                         aacaacCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtgct                                                                                                            <  2:2733572/150‑1 (MQ=255)
                                         aacaacCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtgct                                                                                                            >  2:84013/1‑150 (MQ=255)
                                           caacCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtgctgc                                                                                                          >  1:1818782/1‑150 (MQ=255)
                                            aacCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGcc                                                                                                         >  2:1917957/1‑150 (MQ=255)
                                               cTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCaa                                                                                                      >  1:2849824/1‑150 (MQ=255)
                                               cTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCaa                                                                                                      >  2:156358/1‑150 (MQ=255)
                                               cTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCaa                                                                                                      <  2:1248447/150‑1 (MQ=255)
                                               cTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCaa                                                                                                      >  2:4003675/1‑150 (MQ=255)
                                                tGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAAc                                                                                                     <  1:1764540/150‑1 (MQ=255)
                                                tGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAAc                                                                                                     <  1:3255554/150‑1 (MQ=255)
                                                tGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAAc                                                                                                     >  1:3032413/1‑150 (MQ=255)
                                                 gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACg                                                                                                    <  2:2610822/150‑1 (MQ=255)
                                                     ccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCcac                                                                                                <  1:2388840/150‑1 (MQ=255)
                                                        gATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACg                                                                                             >  1:3634958/1‑150 (MQ=255)
                                                          tgtgGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGc                                                                                           <  2:2654918/150‑1 (MQ=255)
                                                          tgtgGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGc                                                                                           <  2:3493887/150‑1 (MQ=255)
                                                            tgGACCTCGTCATCAAACACCGTGACCTGACCGAACGCGCTATGCTCCAGGTTCCGAAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGaaa                                                                                         >  2:697257/1‑150 (MQ=255)
                                                             ggACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCaaa                                                                                        >  2:3295083/1‑150 (MQ=255)
                                                              gACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAAc                                                                                       <  2:549260/150‑1 (MQ=255)
                                                                 cTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGaa                                                                                    <  1:3430028/150‑1 (MQ=255)
                                                                   cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAaga                                                                                  <  1:1775294/150‑1 (MQ=255)
                                                                     tcatcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAgaga                                                                                <  1:3543140/150‑1 (MQ=255)
                                                                      catcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGaa                                                                               <  1:689533/150‑1 (MQ=255)
                                                                       atcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGaaa                                                                              >  1:3661185/1‑150 (MQ=255)
                                                                       atcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGaaa                                                                              >  1:3662547/1‑150 (MQ=255)
                                                                        tcactcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAg                                                                             >  1:3609878/1‑150 (MQ=255)
                                                                           ctcacCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTaa                                                                          <  2:1506924/150‑1 (MQ=255)
                                                                              acCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAg                                                                       >  1:711093/1‑150 (MQ=255)
                                                                                cGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGAc                                                                     <  2:3115522/150‑1 (MQ=255)
                                                                                      ctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTg                                                               >  2:3420345/1‑150 (MQ=255)
                                                                                      ctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTg                                                               <  2:2015474/150‑1 (MQ=255)
                                                                                      ctgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTg                                                               <  2:3537208/150‑1 (MQ=255)
                                                                                       tgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGa                                                              <  1:812092/150‑1 (MQ=255)
                                                                                       tgaccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGa                                                              >  1:1514652/1‑150 (MQ=255)
                                                                                         accGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGACAGTAAAAGACAGCCTGaaa                                                            >  2:2078336/1‑150 (MQ=255)
                                                                                         accGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGaaa                                                            >  1:4172007/1‑150 (MQ=255)
                                                                                          ccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAg                                                           <  2:3682381/150‑1 (MQ=255)
                                                                                          ccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAg                                                           >  1:4120512/1‑150 (MQ=255)
                                                                                          ccGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAg                                                           <  1:3109623/150‑1 (MQ=255)
                                                                                            gAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGAc                                                         >  2:3951580/1‑150 (MQ=255)
                                                                                               cgcgCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGc                                                      >  1:2473068/1‑150 (MQ=255)
                                                                                               cgcgCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGc                                                      <  1:1095470/150‑1 (MQ=255)
                                                                                                 cgcTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCtt                                                    >  1:3470791/1‑150 (MQ=255)
                                                                                                 cgcTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCtt                                                    <  2:3454695/150‑1 (MQ=255)
                                                                                                  gcTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCttt                                                   <  1:839626/150‑1 (MQ=255)
                                                                                                   cTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGGCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTg                                                  >  2:933169/1‑150 (MQ=255)
                                                                                                     aTGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGAc                                                >  2:852993/1‑150 (MQ=255)
                                                                                                      tGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACg                                               <  2:2836724/150‑1 (MQ=255)
                                                                                                       gcgcCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGa                                              >  2:3236150/1‑150 (MQ=255)
                                                                                                        cgcCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGAt                                             <  1:2912883/150‑1 (MQ=255)
                                                                                                          ccAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTc                                           >  1:578244/1‑150 (MQ=255)
                                                                                                          ccAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTc                                           >  2:2152276/1‑150 (MQ=255)
                                                                                                           cAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTcc                                          <  2:1007683/150‑1 (MQ=255)
                                                                                                             ggTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAg                                        <  1:936592/150‑1 (MQ=255)
                                                                                                                 ccGCAGGCACAGCATATTTCGCTGAACAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGATGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCt                                    <  1:1306200/150‑1 (MQ=255)
                                                                                                                     aGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAAcc                                >  1:1710099/1‑150 (MQ=255)
                                                                                                                     aGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAAcc                                >  1:2018188/1‑150 (MQ=255)
                                                                                                                        cacaGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt                             >  1:2121710/1‑150 (MQ=255)
                                                                                                                        cacaGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGt                             >  2:1639527/1‑150 (MQ=255)
                                                                                                                         acaGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGtt                            <  2:1686294/150‑1 (MQ=255)
                                                                                                                         acaGCATATTTCGCTGACCAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGtt                            >  2:3931121/1‑150 (MQ=255)
                                                                                                                          caGCATATTTCGCTGACCAACTTCCTCGACGGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTc                           <  1:126758/150‑1 (MQ=255)
                                                                                                                                atTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTa                     >  2:1542236/1‑150 (MQ=255)
                                                                                                                                  ttCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAgg                   <  1:2597072/150‑1 (MQ=255)
                                                                                                                                    cGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGgcg                 <  2:1090804/150‑1 (MQ=255)
                                                                                                                                    cGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGgcg                 >  2:3052314/1‑150 (MQ=255)
                                                                                                                                           cAACTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGaa          <  2:3367757/150‑1 (MQ=255)
                                                                                                                                           cAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGaa          <  2:4130863/150‑1 (MQ=255)
                                                                                                                                           cAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGaa          >  1:1631131/1‑150 (MQ=255)
                                                                                                                                            aaCTTCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAAc         >  2:3012504/1‑150 (MQ=255)
                                                                                                                                               ttCCTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATc      >  2:3177807/1‑150 (MQ=255)
                                                                                                                                                 ccTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCtt    <  2:1834758/150‑1 (MQ=255)
                                                                                                                                                 ccTCTACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCtt    >  2:2262954/1‑150 (MQ=255)
                                                                                                                                                  ctctACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTc   <  1:1178902/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATGCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCTTTGACGATTCCAGTGCTAACCTGTTCAAGCTAGGCGCGGAGAACATCTTCC  >  CP009273/3766851‑3767147

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: