Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,358,855 G→T 92.1% G123C (GGC→TGC)  puuB → gamma‑glutamylputrescine oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,358,8550GT92.1% 479.3 / 34.7 177G123C (GGC→TGC) puuBgamma‑glutamylputrescine oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (5/9);  new base T (88/75);  total (93/84)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 187 total aligned reads displayed.

TGAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACC  >  CP009273/1358709‑1358998
                                                                                                                                                  |                                                                                                                                               
  aaCTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCctgcg                                                                                                                                              >  2:188594/1‑148 (MQ=255)
    cTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCtgcggc                                                                                                                                            <  2:1137472/148‑1 (MQ=255)
    cTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCtgcggc                                                                                                                                            <  1:1645302/148‑1 (MQ=255)
    cTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCtgcggc                                                                                                                                            <  1:1645353/148‑1 (MQ=255)
       cTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTg                                                                                                                                         <  1:1198882/148‑1 (MQ=255)
       cTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTg                                                                                                                                         <  1:1744241/148‑1 (MQ=255)
        tACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGt                                                                                                                                        <  2:384212/148‑1 (MQ=255)
         aCAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGtt                                                                                                                                       <  1:2857645/148‑1 (MQ=255)
         aCAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGtt                                                                                                                                       <  1:3302634/148‑1 (MQ=255)
          cAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGttt                                                                                                                                      <  2:2091079/148‑1 (MQ=255)
          cAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGttt                                                                                                                                      <  2:2091017/148‑1 (MQ=255)
          cAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGttt                                                                                                                                      >  2:1546514/1‑148 (MQ=255)
           aGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTtg                                                                                                                                     >  2:1195185/1‑148 (MQ=255)
                 gACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGa                                                                                                                               <  2:1297296/148‑1 (MQ=255)
                     tCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGatgaa                                                                                                                           >  1:1463539/1‑148 (MQ=255)
                      cGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGatgaat                                                                                                                          >  2:594478/1‑148 (MQ=255)
                        aCGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGatgaatga                                                                                                                        <  2:2677391/148‑1 (MQ=255)
                        aCGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGatgaatga                                                                                                                        >  2:890881/1‑148 (MQ=255)
                         cGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGAt                                                                                                                       >  1:1545165/1‑148 (MQ=255)
                           tGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATaa                                                                                                                     <  2:1447741/148‑1 (MQ=255)
                           tGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATaa                                                                                                                     >  1:2363831/1‑148 (MQ=255)
                           tGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATaa                                                                                                                     >  2:948719/1‑148 (MQ=255)
                               cGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAg                                                                                                                 >  2:2307940/1‑148 (MQ=255)
                                 aaaaaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATACACAGCt                                                                                                               >  1:1072033/1‑148 (MQ=255)
                                 aaaaaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCt                                                                                                               <  2:762992/148‑1 (MQ=255)
                                  aaaaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTc                                                                                                              >  2:2806179/1‑148 (MQ=255)
                                   aaaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCg                                                                                                             >  1:892533/1‑148 (MQ=255)
                                       gCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCcac                                                                                                         >  2:1082364/1‑148 (MQ=255)
                                       gCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCcac                                                                                                         >  2:1590867/1‑148 (MQ=255)
                                       gCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCcac                                                                                                         >  1:135140/1‑148 (MQ=255)
                                       gCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCcac                                                                                                         >  2:3182261/1‑148 (MQ=255)
                                       gCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCcac                                                                                                         >  2:2126988/1‑148 (MQ=255)
                                         tACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCacac                                                                                                       >  2:2533546/1‑148 (MQ=255)
                                          aCGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACt                                                                                                      >  1:2246215/1‑148 (MQ=255)
                                          aCGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACt                                                                                                      >  1:2246186/1‑148 (MQ=255)
                                          aCGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACt                                                                                                      <  1:3007425/148‑1 (MQ=255)
                                             gCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGa                                                                                                   <  1:494503/148‑1 (MQ=255)
                                               aTGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAag                                                                                                 <  2:688688/148‑1 (MQ=255)
                                                 ggACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAagag                                                                                               <  2:1609175/148‑1 (MQ=255)
                                                  gACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGc                                                                                              >  1:3040825/1‑148 (MQ=255)
                                                      ccGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGaa                                                                                          <  2:1808109/148‑1 (MQ=255)
                                                      ccGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGaa                                                                                          <  1:686319/148‑1 (MQ=255)
                                                        gCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAag                                                                                        >  1:710277/1‑148 (MQ=255)
                                                         cccGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAaga                                                                                       <  1:2466398/148‑1 (MQ=255)
                                                           cGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAgaga                                                                                     <  1:3162243/148‑1 (MQ=255)
                                                           cGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAgaga                                                                                     <  1:1520532/148‑1 (MQ=255)
                                                             tATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAAc                                                                                   <  2:3217051/148‑1 (MQ=255)
                                                             tATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAAc                                                                                   >  2:497998/1‑148 (MQ=255)
                                                             tATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAAc                                                                                   <  1:2939490/148‑1 (MQ=255)
                                                              aTGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACt                                                                                  <  1:200940/148‑1 (MQ=255)
                                                                gCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTgg                                                                                <  1:2024052/148‑1 (MQ=255)
                                                                  tCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGa                                                                              <  2:144264/148‑1 (MQ=255)
                                                                  tCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGa                                                                              <  1:781390/148‑1 (MQ=255)
                                                                  tCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGa                                                                              <  1:248992/148‑1 (MQ=255)
                                                                   cGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGaa                                                                             <  2:2070926/148‑1 (MQ=255)
                                                                   cGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGaa                                                                             <  2:1339139/148‑1 (MQ=255)
                                                                      cagcaTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGc                                                                          <  1:166044/148‑1 (MQ=255)
                                                                           tgatgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACgg                                                                     >  2:3084287/1‑148 (MQ=255)
                                                                             atgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCa                                                                   >  1:2005208/1‑148 (MQ=255)
                                                                             atgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCa                                                                   >  1:294936/1‑148 (MQ=255)
                                                                                    aaGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACa                                                            >  1:2748715/1‑148 (MQ=255)
                                                                                     aGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAg                                                           <  2:339399/148‑1 (MQ=255)
                                                                                       gtggtgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCt                                                         <  1:2068690/148‑1 (MQ=255)
                                                                                           tgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGaa                                                     <  2:2971379/148‑1 (MQ=255)
                                                                                           tgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGaa                                                     <  1:3032977/148‑1 (MQ=255)
                                                                                            gagaTCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAAt                                                    >  2:1455673/1‑148 (MQ=255)
                                                                                            gagaTCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAAt                                                    >  2:799581/1‑148 (MQ=255)
                                                                                             agaTCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAAtt                                                   <  2:616063/148‑1 (MQ=255)
                                                                                             agaTCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAAtt                                                   >  1:3283157/1‑148 (MQ=255)
                                                                                             agaTCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAAtt                                                   <  2:226349/148‑1 (MQ=255)
                                                                                              gaTCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTg                                                  >  2:1718495/1‑148 (MQ=255)
                                                                                               atcatcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGc                                                 <  2:2147033/148‑1 (MQ=255)
                                                                                               atcatcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGc                                                 >  2:1464850/1‑148 (MQ=255)
                                                                                                 catcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTg                                               <  1:2532529/148‑1 (MQ=255)
                                                                                                  atcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTgg                                              >  2:3122087/1‑148 (MQ=255)
                                                                                                    ccGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGac                                            >  1:3058787/1‑148 (MQ=255)
                                                                                                       cgAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgcc                                         <  2:2437737/148‑1 (MQ=255)
                                                                                                       cgAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgcc                                         <  2:2437699/148‑1 (MQ=255)
                                                                                                          aCGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgccaac                                      >  1:1885732/1‑148 (MQ=255)
                                                                                                           cGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgccaacg                                     >  1:377795/1‑148 (MQ=255)
                                                                                                             tATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgccaacgcc                                   >  1:733660/1‑148 (MQ=255)
                                                                                                                   aCGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGc                             >  1:3181690/1‑148 (MQ=255)
                                                                                                                     gTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCg                           >  1:2758043/1‑148 (MQ=255)
                                                                                                                       tATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTg                         >  1:772444/1‑148 (MQ=255)
                                                                                                                       tATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTg                         >  1:772364/1‑148 (MQ=255)
                                                                                                                          cAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAg                      >  2:3079374/1‑148 (MQ=255)
                                                                                                                          cAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAg                      >  2:3012290/1‑148 (MQ=255)
                                                                                                                          cAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAg                      >  2:348046/1‑148 (MQ=255)
                                                                                                                           aGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGt                     <  1:1332146/148‑1 (MQ=255)
                                                                                                                            gATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTa                    >  1:1616999/1‑148 (MQ=255)
                                                                                                                              ttGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGc                  >  1:307387/1‑148 (MQ=255)
                                                                                                                               tGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGcc                 >  1:2903294/1‑148 (MQ=255)
                                                                                                                                       aCTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCg         <  1:2126988/148‑1 (MQ=255)
                                                                                                                                       aCTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCg         <  1:1590867/148‑1 (MQ=255)
                                                                                                                                        cTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGc        <  1:630376/148‑1 (MQ=255)
                                                                                                                                         tACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCt       >  1:2407900/1‑148 (MQ=255)
                                                                                                                                         tACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCt       <  1:1195185/148‑1 (MQ=255)
                                                                                                                                            cGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCtata    >  2:1100466/1‑148 (MQ=255)
                                                                                                                                             gCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATAc   >  2:2801479/1‑148 (MQ=255)
                                                                                                                                              ccccTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCAGTCGCCGTGAAGTAGCCAGCGATCGCTATAcc  >  2:2729221/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TGAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACC  >  CP009273/1358709‑1358998

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: