Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,016,392 G→A 6.2% intergenic (‑17/‑202) sulA ← / → sxy SOS cell division inhibitor/CRP‑S‑dependent promoter expression factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,016,3920GA 6.2% 765.9 / 25.0 212intergenic (‑17/‑202)sulA/sxySOS cell division inhibitor/CRP‑S‑dependent promoter expression factor
Reads supporting (aligned to +/- strand):  ref base G (91/107);  new base A (7/6);  total (99/113)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

TCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAA  >  CP009273/1016243‑1016540
                                                                                                                                                     |                                                                                                                                                    
 catGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCtgt                                                                                                                                                     >  2:3442234/1‑150 (MQ=255)
 catGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCtgt                                                                                                                                                     >  1:2475265/1‑150 (MQ=255)
 catGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCtgt                                                                                                                                                     >  2:1347685/1‑150 (MQ=255)
 catGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCtgt                                                                                                                                                     >  1:1190063/1‑150 (MQ=255)
   tGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGa                                                                                                                                                   >  1:1590327/1‑150 (MQ=255)
    gggCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTATGAg                                                                                                                                                  >  1:4060482/1‑150 (MQ=255)
      gCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGtt                                                                                                                                                >  1:2993440/1‑150 (MQ=255)
       cTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTa                                                                                                                                               >  1:1813624/1‑150 (MQ=255)
       cTGATATTCGCGAGAGACAACTTCACTGATTTGCCCGGGTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCAGAGCCTGAAGTGTACATAATAAATCCAGCCCCTGTGAGTTa                                                                                                                                               <  2:335144/150‑1 (MQ=255)
        tGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTAc                                                                                                                                              >  2:3781491/1‑150 (MQ=255)
            cTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTa                                                                                                                                          >  1:1635436/1‑150 (MQ=255)
               cgcgATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATgg                                                                                                                                       >  1:3369593/1‑150 (MQ=255)
               cgcgATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATgg                                                                                                                                       <  2:2239875/150‑1 (MQ=255)
                   aTAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgt                                                                                                                                   <  2:1746262/150‑1 (MQ=255)
                   aTAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgt                                                                                                                                   <  1:3418192/150‑1 (MQ=255)
                   aTAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgt                                                                                                                                   >  1:410051/1‑150 (MQ=255)
                     aGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtac                                                                                                                                 >  1:1060132/1‑150 (MQ=255)
                       aCAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacag                                                                                                                               <  1:556170/150‑1 (MQ=255)
                        cAACTTCACTGATCAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagt                                                                                                                              <  2:3692080/150‑1 (MQ=255)
                              cACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATcc                                                                                                                        >  2:65664/1‑150 (MQ=255)
                              cACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATcc                                                                                                                        >  2:65638/1‑150 (MQ=255)
                              cACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTATGAGTTACTGTATGGATGTACAGTACATcc                                                                                                                        >  2:3853554/1‑150 (MQ=255)
                               aCTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCa                                                                                                                       <  1:1520776/150‑1 (MQ=255)
                                cTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAg                                                                                                                      <  1:2380508/150‑1 (MQ=255)
                                cTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAg                                                                                                                      <  2:1571787/150‑1 (MQ=255)
                                  gATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTg                                                                                                                    >  2:1882066/1‑150 (MQ=255)
                                   aTAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGTGGATGAGAACGACGAAGACCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGa                                                                                                                   >  1:2098280/1‑150 (MQ=255)
                                         ccGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGAcaacaa                                                                                                             <  1:1579958/150‑1 (MQ=255)
                                           ggCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAg                                                                                                           >  1:412372/1‑150 (MQ=255)
                                           ggCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAg                                                                                                           <  1:1401729/150‑1 (MQ=255)
                                           ggCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAg                                                                                                           >  2:3217/1‑150 (MQ=255)
                                            gCTGTAGTGTTTTCCGTAGAGACACGCGCAATTGTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGa                                                                                                          <  1:2203740/150‑1 (MQ=255)
                                             cTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGAt                                                                                                         >  1:148007/1‑150 (MQ=255)
                                             cTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGAt                                                                                                         <  2:763406/150‑1 (MQ=255)
                                              tGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATc                                                                                                        <  1:316012/150‑1 (MQ=255)
                                              tGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATc                                                                                                        >  2:4100728/1‑150 (MQ=255)
                                                  gtgtTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAAcc                                                                                                    <  2:1193616/150‑1 (MQ=255)
                                                     ttttCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTa                                                                                                 <  2:747207/150‑1 (MQ=255)
                                                       ttCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTAtt                                                                                               <  1:2033674/150‑1 (MQ=255)
                                                          cGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTc                                                                                            >  2:2379867/1‑150 (MQ=255)
                                                           gTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCg                                                                                           >  1:43327/1‑150 (MQ=255)
                                                             agagACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGa                                                                                         >  1:3091466/1‑150 (MQ=255)
                                                              tagaCACGCGCAATTTTAATAGATGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTAAATAATCAATCAAGCCCCTCTGATTTAATTTATGGATGTACAGTACATCAAGAGACAACAAAGAACAACACTATTTTCGGaa                                                                                        >  2:1312647/2‑150 (MQ=255)
                                                              gagaCACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGaa                                                                                        >  2:2113503/1‑150 (MQ=255)
                                                                gaCACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAag                                                                                      >  2:189558/1‑150 (MQ=255)
                                                                  cacGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAagag                                                                                    >  1:1918576/1‑150 (MQ=255)
                                                                   acGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGc                                                                                   <  1:201880/150‑1 (MQ=255)
                                                                    cgcgCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGcc                                                                                  <  1:1560296/150‑1 (MQ=255)
                                                                       gcAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCg                                                                               >  1:171898/1‑150 (MQ=255)
                                                                         aaTTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCa                                                                             >  2:3371619/1‑150 (MQ=255)
                                                                                cTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                                                      >  1:4170002/1‑150 (MQ=255)
                                                                                   gctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCg                                                                   <  1:1251243/150‑1 (MQ=255)
                                                                                      gcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTg                                                                <  1:388421/150‑1 (MQ=255)
                                                                                        ggATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGt                                                              <  1:2575838/150‑1 (MQ=255)
                                                                                        ggATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGt                                                              >  1:351012/1‑150 (MQ=255)
                                                                                          aTGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGa                                                            <  1:308791/150‑1 (MQ=255)
                                                                                            gagaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACg                                                          >  2:3318301/1‑150 (MQ=255)
                                                                                             agaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACgg                                                         <  1:847392/150‑1 (MQ=255)
                                                                                             agaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACgg                                                         <  2:2676355/150‑1 (MQ=255)
                                                                                              gaACGACGAAGATCGATGTGCATAGCCTGTAGTGTACATAATCACTCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACggg                                                        <  1:178707/150‑1 (MQ=255)
                                                                                              gaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACggg                                                        <  1:3116581/150‑1 (MQ=255)
                                                                                                  gacgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaaaa                                                    >  2:16902/1‑150 (MQ=255)
                                                                                                      aagaaCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATa                                                <  1:2625184/150‑1 (MQ=255)
                                                                                                         aaCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAt                                             <  2:227636/150‑1 (MQ=255)
                                                                                                            gATGTGCATAGCCTGATGTGTACATAATCAATCCATCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGTGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTaa                                          <  2:1315691/150‑1 (MQ=255)
                                                                                                             aTGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAAt                                         <  2:3463022/150‑1 (MQ=255)
                                                                                                             aTGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAAt                                         <  2:3073315/150‑1 (MQ=255)
                                                                                                              tgtgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATc                                        <  2:867959/150‑1 (MQ=255)
                                                                                                               gtgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCt                                       <  1:3578455/150‑1 (MQ=255)
                                                                                                               gtgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCt                                       <  1:3742403/150‑1 (MQ=255)
                                                                                                                tgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCtt                                      <  2:1987033/150‑1 (MQ=255)
                                                                                                                 gCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTATGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTg                                     <  2:408694/150‑1 (MQ=255)
                                                                                                                    tAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGccc                                  <  1:2977698/150‑1 (MQ=255)
                                                                                                                     aGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGcccc                                 <  2:1285547/150‑1 (MQ=255)
                                                                                                                       ccTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCtt                               <  2:3261335/150‑1 (MQ=255)
                                                                                                                         tGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTaa                             >  2:961852/1‑150 (MQ=255)
                                                                                                                         tGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTaa                             >  2:3052611/1‑150 (MQ=255)
                                                                                                                         tGAAGTGTAAAAAATCAATCCAGCCACTGTGAGTTACTGTATGGATGTACATAACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGATACTCGCAAATTTTGTCGGTGGTGACGGGAAAACATAAATTAATCTTGCCAATTaa                             >  2:130698/1‑150 (MQ=255)
                                                                                                                          gAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAg                            >  2:2718639/1‑150 (MQ=255)
                                                                                                                           aaGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGa                           >  2:2670106/1‑150 (MQ=255)
                                                                                                                            aGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGaa                          >  1:2515565/1‑150 (MQ=255)
                                                                                                                             gtgtACATAATCAATCCAGCCCCTATGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAAt                         >  2:1140998/1‑150 (MQ=255)
                                                                                                                                   aTAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTg                   <  2:1786645/150‑1 (MQ=255)
                                                                                                                                     aatcaatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGcc                 >  2:1920380/1‑150 (MQ=255)
                                                                                                                                      atcaatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCt                >  2:4134265/1‑150 (MQ=255)
                                                                                                                                        caatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTAt              >  2:3245126/1‑150 (MQ=255)
                                                                                                                                           tcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTAtttt           <  1:642852/150‑1 (MQ=255)
                                                                                                                                           tcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACAGCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTAtttt           >  1:3938661/1‑150 (MQ=255)
                                                                                                                                            ccAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTc          <  2:4050019/150‑1 (MQ=255)
                                                                                                                                            ccAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTc          <  2:3568847/150‑1 (MQ=255)
                                                                                                                                            ccAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTc          <  2:2249856/150‑1 (MQ=255)
                                                                                                                                             cAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCg         >  2:3613666/1‑150 (MQ=255)
                                                                                                                                             cAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTCAAGAATAAGTTGCCTATTTTCg         <  2:549603/150‑1 (MQ=255)
                                                                                                                                              aGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGt        >  1:768879/1‑150 (MQ=255)
                                                                                                                                               gCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTa       <  2:3667381/150‑1 (MQ=255)
                                                                                                                                                ccccTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAg      <  1:2742875/150‑1 (MQ=255)
                                                                                                                                                 cccTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGt     >  1:2331639/1‑150 (MQ=255)
                                                                                                                                                  ccTGTGAGTTACTGTATTGATGTACAGTACATCAATTTACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATTAGTTGCCTATTTTCGTAGtt    <  2:1740619/150‑1 (MQ=255)
                                                                                                                                                    tgtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTaa  >  1:1002910/1‑150 (MQ=255)
                                                                                                                                                    tgtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTACGTAGTTaa  >  2:2204506/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAA  >  CP009273/1016243‑1016540

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: