Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,303,617 G→T 5.8% intergenic (+1/‑68) rpiB → / → yjdP ribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase/putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,303,6170GT 5.8% 1083.8 / 24.7 282intergenic (+1/‑68)rpiB/yjdPribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase/putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (99/166);  new base T (7/10);  total (106/176)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 303 total aligned reads displayed.

CACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACA  >  CP009273/4303469‑4303766
                                                                                                                                                    |                                                                                                                                                     
 acCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGAt                                                                                                                                                     >  1:2435706/1‑150 (MQ=255)
  ccAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGAtt                                                                                                                                                    <  2:4401241/150‑1 (MQ=255)
   cAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTc                                                                                                                                                   <  1:1131817/150‑1 (MQ=255)
   cAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTc                                                                                                                                                   >  2:4572294/1‑150 (MQ=255)
    aaCGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCa                                                                                                                                                  <  2:3029303/150‑1 (MQ=255)
     aCGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCAt                                                                                                                                                 <  1:1189157/150‑1 (MQ=255)
      cGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATc                                                                                                                                                >  2:528666/1‑150 (MQ=255)
       gTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGNTTCATcc                                                                                                                                               <  1:2548974/150‑1 (MQ=255)
       gTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATcc                                                                                                                                               <  1:4391990/150‑1 (MQ=255)
       gTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATcc                                                                                                                                               >  1:70578/1‑150 (MQ=255)
         gCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCac                                                                                                                                             >  1:450492/1‑150 (MQ=255)
         gCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCac                                                                                                                                             >  2:1609917/1‑150 (MQ=255)
         gCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCac                                                                                                                                             <  2:2413612/150‑1 (MQ=255)
           tGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCacta                                                                                                                                           >  1:541466/1‑150 (MQ=255)
             gCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCactact                                                                                                                                         <  1:1264006/150‑1 (MQ=255)
                 ttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCa                                                                                                                                     <  2:3306987/150‑1 (MQ=255)
                   ggTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATg                                                                                                                                   <  2:2589111/150‑1 (MQ=255)
                      tCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGAt                                                                                                                                <  2:2374952/150‑1 (MQ=255)
                      tCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGAt                                                                                                                                >  2:2119331/1‑150 (MQ=255)
                       cACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATg                                                                                                                               <  2:3360023/150‑1 (MQ=255)
                        aCGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGa                                                                                                                              <  2:3729030/150‑1 (MQ=255)
                        aCGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGa                                                                                                                              <  2:193005/150‑1 (MQ=255)
                           aGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTa                                                                                                                           <  1:2352780/150‑1 (MQ=255)
                             tGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAAt                                                                                                                         <  1:4555856/150‑1 (MQ=255)
                               gTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGa                                                                                                                       <  2:3878088/150‑1 (MQ=255)
                                ttGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGAt                                                                                                                      >  1:174553/1‑150 (MQ=255)
                                 tGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGAtt                                                                                                                     >  2:3917256/1‑150 (MQ=255)
                                    ccTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAt                                                                                                                  <  1:2532973/150‑1 (MQ=255)
                                          aCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACATTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCTATTACGGCAATAGAGAAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATTATTAATGTGGAt                                                                                                            <  1:2778540/150‑1 (MQ=255)
                                            tGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATag                                                                                                          <  2:1810370/150‑1 (MQ=255)
                                            tGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATag                                                                                                          <  2:2263413/150‑1 (MQ=255)
                                             ggCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATaga                                                                                                         >  1:1330523/1‑150 (MQ=255)
                                               cAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGt                                                                                                       <  1:4581342/150‑1 (MQ=255)
                                                  aaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTc                                                                                                    >  1:111645/1‑150 (MQ=255)
                                                  aaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTc                                                                                                    <  2:1836568/150‑1 (MQ=255)
                                                   aaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCt                                                                                                   <  2:1998888/150‑1 (MQ=255)
                                                    aTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCtt                                                                                                  <  2:825447/150‑1 (MQ=255)
                                                    aTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCtt                                                                                                  <  2:903691/150‑1 (MQ=255)
                                                      gATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCtttt                                                                                                <  1:2331994/150‑1 (MQ=255)
                                                       aTTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCttttt                                                                                               <  1:1620529/150‑1 (MQ=255)
                                                       aTTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCttttt                                                                                               <  2:3631853/150‑1 (MQ=255)
                                                        ttGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTg                                                                                              <  1:4301077/150‑1 (MQ=255)
                                                             gATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGtt                                                                                         >  2:2199428/1‑150 (MQ=255)
                                                               tGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTgg                                                                                       <  1:4170801/150‑1 (MQ=255)
                                                                gCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGc                                                                                      <  1:3219508/150‑1 (MQ=255)
                                                                    ggCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAAt                                                                                  <  1:3449033/150‑1 (MQ=255)
                                                                      cTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGa                                                                                <  2:3123516/150‑1 (MQ=255)
                                                                      cTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGa                                                                                <  2:2582231/150‑1 (MQ=255)
                                                                       tGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGaa                                                                               >  1:4155797/1‑150 (MQ=255)
                                                                       tGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGaa                                                                               >  1:651390/1‑150 (MQ=255)
                                                                        gggCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGaaa                                                                              <  1:3110272/150‑1 (MQ=255)
                                                                             cacaGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCtt                                                                         <  2:3828302/150‑1 (MQ=255)
                                                                              acaGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCttt                                                                        <  1:1038568/150‑1 (MQ=255)
                                                                                aGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTcc                                                                      >  1:2133383/1‑150 (MQ=255)
                                                                                     gAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGATATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAt                                                                 <  2:1662088/150‑1 (MQ=255)
                                                                                     gAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAt                                                                 <  2:712312/150‑1 (MQ=255)
                                                                                      aaGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGATATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAtt                                                                <  2:574716/150‑1 (MQ=255)
                                                                                      aaGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCAACCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTACTTTTTGAGGTTGGCTAATGAAACGCTTTCCGcgagt                                                                >  1:4268708/1‑146 (MQ=255)
                                                                                       aGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGCGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAttt                                                               <  2:4111622/150‑1 (MQ=255)
                                                                                         gCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCt                                                             <  1:2234687/150‑1 (MQ=255)
                                                                                         gCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCt                                                             <  1:1899569/150‑1 (MQ=255)
                                                                                          cGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCtt                                                            >  2:1465257/1‑150 (MQ=255)
                                                                                          cGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCtt                                                            <  1:2694083/150‑1 (MQ=255)
                                                                                          cGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCtt                                                            <  2:4207968/150‑1 (MQ=255)
                                                                                            gtcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGATATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCtttt                                                          >  2:136662/1‑150 (MQ=255)
                                                                                            gtcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTGGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCtttt                                                          <  2:1768586/150‑1 (MQ=255)
                                                                                             tcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTa                                                         <  1:523574/150‑1 (MQ=255)
                                                                                             tcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTa                                                         <  1:3653837/150‑1 (MQ=255)
                                                                                               gtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTAtt                                                       >  1:852326/1‑150 (MQ=255)
                                                                                                tcatcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTc                                                      >  1:979613/1‑150 (MQ=255)
                                                                                                tcatcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTc                                                      >  1:3606159/1‑150 (MQ=255)
                                                                                                        aacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCt                                              <  1:2005814/150‑1 (MQ=255)
                                                                                                        aacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCt                                              >  1:3632631/1‑150 (MQ=255)
                                                                                                         acGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTc                                             >  1:953532/1‑150 (MQ=255)
                                                                                                            cgTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACg                                          >  2:1048316/1‑150 (MQ=255)
                                                                                                               ggAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTg                                       >  2:3987989/1‑150 (MQ=255)
                                                                                                                     gATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCAcc                                 <  1:513530/150‑1 (MQ=255)
                                                                                                                      aTTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCg                                >  2:3906671/1‑150 (MQ=255)
                                                                                                                      aTTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCg                                >  2:3462127/1‑150 (MQ=255)
                                                                                                                          cGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTcc                            >  1:2647312/1‑150 (MQ=255)
                                                                                                                           ggCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCg                           <  1:2119331/150‑1 (MQ=255)
                                                                                                                                tAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCa                      <  2:2547941/150‑1 (MQ=255)
                                                                                                                                      agcGGAGAAATAGATATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGa                >  2:832257/1‑150 (MQ=255)
                                                                                                                                      agcGGAGAAATAGATATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGa                >  1:3900297/1‑150 (MQ=255)
                                                                                                                                      agcGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGa                >  1:3049167/1‑150 (MQ=255)
                                                                                                                                        cGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGata              <  1:270422/150‑1 (MQ=255)
                                                                                                                                          gagaAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtatt            <  1:3825032/150‑1 (MQ=255)
                                                                                                                                           agaAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtatta           <  1:3189896/150‑1 (MQ=255)
                                                                                                                                            gaAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtattat          <  1:4025128/150‑1 (MQ=255)
                                                                                                                                            gaAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtattat          <  2:1738193/150‑1 (MQ=255)
                                                                                                                                            gaAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtattat          <  2:1094823/150‑1 (MQ=255)
                                                                                                                                             aaaTAGATATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATc         <  2:2969314/150‑1 (MQ=255)
                                                                                                                                             aaaTAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATc         >  2:282339/1‑150 (MQ=255)
                                                                                                                                                tagagaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAGCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgac      <  1:700332/148‑1 (MQ=255)
                                                                                                                                                tagagaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgac      >  1:3207878/3‑150 (MQ=255)
                                                                                                                                                tagagaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgac      <  1:1454517/148‑1 (MQ=255)
                                                                                                                                                 agagaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgacg     <  2:4256832/149‑1 (MQ=255)
                                                                                                                                                 agagaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgacg     >  1:258258/2‑150 (MQ=255)
                                                                                                                                                   agaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCgacgac   >  2:3955794/1‑150 (MQ=255)
                                                                                                                                                   agaTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCgacgac   <  2:2141667/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACA  >  CP009273/4303469‑4303766

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: