Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,713 C→A 20.0% intergenic (‑30/+278) gatY ← / ← fbaB D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,170,7130CA20.0% 357.8 / 28.2 175intergenic (‑30/+278)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base C (45/84);  new base A (23/17);  total (74/101)
Fisher's exact test for biased strand distribution p-value = 1.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

GTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGG  >  CP009273/2170564‑2170862
                                                                                                                                                     |                                                                                                                                                       
gTTGGCAGCGGTTTCTACCACCACGTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTc                                                                                                                                                         <  2:1716522/150‑1 (MQ=255)
 tttgCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCa                                                                                                                                                        <  2:3626671/147‑1 (MQ=255)
   ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTTAATGCCGGAACCGCATAACCTCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaa                                                                                                                                                      <  2:263344/150‑1 (MQ=255)
   ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaa                                                                                                                                                      <  2:2529188/150‑1 (MQ=255)
      aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAata                                                                                                                                                  >  2:250161/1‑150 (MQ=255)
       gCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATTTTGAATGCCGGAACCGCATAACCTCAGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACAACGTACTTTTTCATATCCTGTCGTTTTTTTTCTATTTCAAAatat                                                                                                                                                 <  1:2371116/150‑1 (MQ=255)
        cGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAtata                                                                                                                                                <  1:3319999/150‑1 (MQ=255)
         ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATaa                                                                                                                                               <  2:1437622/150‑1 (MQ=255)
         ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATaa                                                                                                                                               >  2:2273501/1‑150 (MQ=255)
         ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTCAAAATATaa                                                                                                                                               >  2:3717154/1‑150 (MQ=255)
         ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAaaaaaaaaa                                                                                                                                                >  1:354402/1‑145 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCCTATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                              <  2:3952104/150‑1 (MQ=255)
           tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATg                                                                                                                                             <  1:3315080/150‑1 (MQ=255)
             tCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaa                                                                                                                                           <  2:1733193/150‑1 (MQ=255)
             tCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGACGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaa                                                                                                                                           <  2:1301142/150‑1 (MQ=255)
              cTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                          <  2:2992921/150‑1 (MQ=255)
              cTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                          <  1:659054/150‑1 (MQ=255)
               tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                         <  2:522704/150‑1 (MQ=255)
                   accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                                     <  2:658668/150‑1 (MQ=255)
                        ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttt                                                                                                                                <  2:822304/150‑1 (MQ=255)
                          gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTTGATTTAAAAATAAAAATAAAATAAtttttttt                                                                                                                              >  1:3812385/1‑150 (MQ=255)
                          gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTATATTTAAAAATATAATTAAAtttttttttttt                                                                                                                              >  1:2929647/1‑150 (MQ=255)
                          gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttttt                                                                                                                              <  2:1897199/150‑1 (MQ=255)
                          gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGACGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttttt                                                                                                                              <  2:1898706/150‑1 (MQ=255)
                           cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTa                                                                                                                             <  1:2630114/150‑1 (MQ=255)
                           cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTa                                                                                                                             <  1:1924206/150‑1 (MQ=255)
                              cGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAat                                                                                                                          <  1:1331306/150‑1 (MQ=255)
                              cGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTTTATTTAAAAAAAAAAATAAAATATTTTTTTTTAAat                                                                                                                          >  1:580944/1‑150 (MQ=255)
                              cGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGGTTTTTATTTACAAAAATAAAGAAAATATtttttttttttt                                                                                                                          >  1:4321439/1‑146 (MQ=255)
                               gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAata                                                                                                                         <  2:70245/150‑1 (MQ=255)
                               gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAata                                                                                                                         <  1:3621866/150‑1 (MQ=255)
                                 ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATc                                                                                                                       <  2:55250/150‑1 (MQ=255)
                                        ttGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATTAAAATAATTTTTTTTAAAAATAGATATaa                                                                                                                >  2:145915/1‑150 (MQ=255)
                                          gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACg                                                                                                              <  1:1254511/150‑1 (MQ=255)
                                          gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACg                                                                                                              <  2:3340727/150‑1 (MQ=255)
                                          gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACg                                                                                                              <  2:2347619/150‑1 (MQ=255)
                                           tgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGAATTCAAAATATAATTAAAATATTTGTTTTTAAAAATAGAGATaaaca                                                                                                             >  2:867191/1‑147 (MQ=255)
                                            gAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGAt                                                                                                            <  1:248067/150‑1 (MQ=255)
                                             aaTATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTGAAAATATTATTATATTATTTTTTTTTTTATATAGAGTTAActtta                                                                                                           >  1:2256525/1‑146 (MQ=255)
                                             aaTATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGGTTTCGCTGTCAAAAAATAAAGAAAATATTTTTTTTTAAAAATAGAGATAACCATc                                                                                                           >  2:3273379/1‑150 (MQ=255)
                                                atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATcaca                                                                                                        <  2:3284570/150‑1 (MQ=255)
                                                    aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAAAATAATGAAAATATTTTTTTTTTAATATAGAGTTAAAGATAACaaaaa                                                                                                    >  2:593552/1‑150 (MQ=255)
                                                        ccGGAACCGCATAACCGCCGCGCTGTGCTTTGTTCAGCATCTGCTTTGACGATACCACGTACATTTTAATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGAc                                                                                                <  2:3861171/150‑1 (MQ=255)
                                                        ccGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGAc                                                                                                <  2:2073883/150‑1 (MQ=255)
                                                         cGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACa                                                                                               <  2:633768/150‑1 (MQ=255)
                                                         cGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTGGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACa                                                                                               <  2:1430404/150‑1 (MQ=255)
                                                         cGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTGAATATATAAATAAATTATTTTTTTTTTAATATAGATATAAAGaaaaaaaaaaaaaaa                                                                                               >  1:494834/1‑145 (MQ=255)
                                                          ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACaa                                                                                              <  2:1535003/150‑1 (MQ=255)
                                                             aCCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGCTTTAAAAATAATAATAAATTTTTTTTTTTTAAAAATCGAGATAAAGATCACAAAAACGACAatat                                                                                           >  2:3107621/1‑150 (MQ=255)
                                                                 cATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACAGTTTCATATCCTGTCGTGTGTTTTTGATTTAAAAATATAATTAAATTTTTTGTTTTTTAATATCGGGTTAAAGATGAAAAAAAAGACAATATGtaa                                                                                       >  1:2565289/1‑150 (MQ=255)
                                                                  aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAATAAAAATAAAATATTTTTTTTTTAATATATAGATAACGATAACAAAAAACACAAAATGtaaa                                                                                      >  2:3422621/1‑150 (MQ=255)
                                                                  aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACAGTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAACTATTTGTTTTTTAATATCGAGATAACGATCACCAAAACCACAAGATGaaaa                                                                                      >  2:2823021/1‑150 (MQ=255)
                                                                    aaCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTTGATTTAAAAAAAAAAATTAAATATTTTTTTTTTAAAATAGAGATAAACATAAAAAAAACCACAaaatataaaat                                                                                    >  2:775494/1‑143 (MQ=255)
                                                                              cTGTGCGTTTTTCAGCATCTGCTTTGCCGATACCCCGTACATTTTAATATCCTGTCTTTTTTTTTCGATTTCAAAATATAATCAAATTTTTTTTTTTTAAATATCGATATACCGATCAAAAAAACGACAATATGAAAATTATTCgaggag                                                                          <  1:1363169/150‑1 (MQ=255)
                                                                              cTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCgaggag                                                                          <  2:1861334/150‑1 (MQ=255)
                                                                               tgtgCGTTGTTCAGAATCTGCTTTGTCGATCCCACGTAATTTTTCATATCGTGTCGTTTGTTTTCTATTACAAAATATATTGTAATTATTTTTTTTTAAATATCGAGATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGt                                                                         <  2:4178842/150‑1 (MQ=255)
                                                                                gtgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATGTAATTATTTTTTTTTTAATATCGAGTTTAAGATCAAAAAAACAAAAAAATAAAAATTATTAGAGGAGtt                                                                        >  1:1806585/1‑149 (MQ=255)
                                                                                 tgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGAATTAAAAAAAAAATTAAAATAATTTTTTTTTAAAAAAGAGATAAACATAACAAAAAAAACAAAAAGAAAAAAATTAGAGGCGtta                                                                       >  2:3808012/1‑148 (MQ=255)
                                                                                  gCGTTGTACAGCATCTGCTTTGTCGATACCCCGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGaa                                                                     <  2:1791827/150‑1 (MQ=255)
                                                                                     ttgttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATTTCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAgg                                                                  <  2:4311337/150‑1 (MQ=255)
                                                                                      tgtaCAGCATCTGCTTGGTCGATACCACGTAAATTTTCATATCCTGTCTTTTTTTTGCGATTTCAAAATATATTTAATTTTTTTTTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGc                                                                 <  2:1916923/150‑1 (MQ=255)
                                                                                       gttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAAATAAATTATTTTTTTTTAAATATAGAGATAAAGATAACAAAAAACACAAAATTAAAATTATTAGAGGAGTGAAATGCa                                                                >  2:206028/1‑150 (MQ=255)
                                                                                           aGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTGAAAAAATAAATAAAATATTTTTTTTTTAATATAGAGATAAACAACACAAAAAAGAAAAAATTTAAAAAATTAGTGGAGTGTGAAGGCaaaa                                                             >  2:2679043/1‑150 (MQ=255)
                                                                                              aTCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTTAATATAGAGATAAAGATAACAAAAAAGAAAAAATGAAAATTATTAGAGTTGTGTATTGCaaaaaaaa                                                         >  2:644798/1‑149 (MQ=255)
                                                                                                cTGTTTTGTCGATCCCAATTAATTTTTAATATCTTTTCTTTTTTTTTCTATTAAAAAATATAATTAATTTTTTTTTTTTTAAATATCTATAAAACGAAAAAAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACgg                                                       <  1:867863/150‑1 (MQ=255)
                                                                                                cTGCTTTGTCGATACCACTTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACgg                                                       <  1:1840524/150‑1 (MQ=255)
                                                                                                cTGATTTTTCGATACCACTTAAATTTTCATATCCTGTCGTTTGTTTTCTATTACAAAATATAATGAATTTATTTTTTTTTAATTATCGATATAACGACCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACgg                                                       <  2:1559435/150‑1 (MQ=255)
                                                                                                 tGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTATATATCGAGATAAAGATAACAAAAACGACAATATGAAAATTATTCGCGGAGTGAAAGGCCAAAATACGGc                                                      >  2:4148979/1‑150 (MQ=255)
                                                                                                 tGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATGAAAATAATTTTTTTTTAAAATAGAGATAACGATAACAAAAACAAAAAAATGAAAATTTATAGAGGAGTGTAATGCAAAAAAAAGGc                                                      >  2:4010831/1‑150 (MQ=255)
                                                                                                      tGTCGATCCCACTTACTTTTAAATAACTTGTCGTTTTTTTTCTTTTTAAAAATTAATTTAATTTTTTTTTTTTTAATTATCTATATAACTAAAACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTccc                                                 <  2:3785456/150‑1 (MQ=255)
                                                                                                         ctaTACCACTTACTTTTTCATATCCTTTCTTTTTTTTTCTATTTAAAAATATAATTAAATTTTTTTTTTTTAAATTTCTAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGAt                                              <  2:1353530/148‑1 (MQ=255)
                                                                                                         cGATACCACGTACATTTTCATATCCTTTCGTTTGTTTTCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGAt                                              <  2:1352224/150‑1 (MQ=255)
                                                                                                          gATACCACGTACAGTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTAAAAATAGAGATAACGATCACAAAAACGCCAATATTAAAATTATTGGGGGAGTGAAAGGCAAAAAAAAGGCCTCCCGaga                                             >  2:4022462/1‑148 (MQ=255)
                                                                                                           aTACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAg                                            <  2:2163261/150‑1 (MQ=255)
                                                                                                           aTACCACGTACATTTTCATATCCTGTCGTGTGTTTTTGATTTCAAAATATAAATAAAATATTTTTTTTTTTAAATAGAGATAAAGAGAACAAAAAACACCAAATGTAAATTAATAGAGGGGTGAAAGGCAAAAAAAAAGCCAcccccgag                                            >  2:455974/1‑145 (MQ=255)
                                                                                                            aaCCACTTAAATTTTAATATCCTGTCTTTTTTTTTCTATTTAAAAATATATTGAAATTATTTTTTTTTAAATATCTATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAgg                                           <  1:3023716/149‑1 (MQ=255)
                                                                                                             acaaCGTAAATTTTCATATCCTGTCGTTTTTTTTGTATTTCAAAATATAATTAATTTATTTTTTTTTAAATATCTTTATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAggg                                          <  1:2961047/147‑1 (MQ=255)
                                                                                                                     ctttttCAGAGCCTGTCGTTTGTTTTCGATTACAAAATATAATGAATTTTTTTGTTTTTAAATATCGAGATAACTCTCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                  <  2:3082406/148‑1 (MQ=255)
                                                                                                                     cATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                  <  2:4488734/150‑1 (MQ=255)
                                                                                                                     atttttCATATCTTGTGTTTTTTTTTCTATTAAAAAATATAATTAAATTATTTTTTTTTAAATATCTATATAACGATCACAACCACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                  <  1:973928/148‑1 (MQ=255)
                                                                                                                     atttttCATACCATGTGTTTTTTTTTCTATTAAAAAATATAATTAATTTTTTTTTTTTTAAATTTCGATAAAAAGATAACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                  <  1:3238521/148‑1 (MQ=255)
                                                                                                                          tCATATCCTGTCGTGTGTTTTGGATTTAAAAAAAAAAATAAAATATTTTTTTTTTAATAAAGACAAAAAGAACACAAAAAAGACAAAAATAAAATTATTAGAGGAGTGAAAGGCAAAAAAAAAGCCACCCCAAAAGGAAGCCCGAGCaaa                             >  2:152040/1‑150 (MQ=255)
                                                                                                                             tatCCTTTTTTTTTTTTTCTTTTAAAAAATATAAATAAATTTTTTTTTTTTAAATATCTATAAAAGTAAAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTg                          <  1:858411/150‑1 (MQ=255)
                                                                                                                              tcccTGTCGTTTTTTTTCGATTAAAAAATAAAATGAAATTATTTTTTTTTAAATATCTATATAACGAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTga                         <  1:4412155/148‑2 (MQ=255)
                                                                                                                                cgtGTCTTTTGTTTTCTATTTCAAAATATAATGAATTTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcgc                       <  2:3882431/148‑1 (MQ=255)
                                                                                                                                   gTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGCGGAGTGAAAGGCAAAAAAACGGCCGCCCGATAGGGAAGCCGTCGCAAAGTGCGCgtg                    >  1:3904991/1‑150 (MQ=255)
                                                                                                                                   gTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                    >  2:3221570/1‑150 (MQ=255)
                                                                                                                                         tgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa              >  1:195279/1‑150 (MQ=255)
                                                                                                                                          gtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt             >  1:2409653/1‑150 (MQ=255)
                                                                                                                                          gtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt             >  1:676077/1‑150 (MQ=255)
                                                                                                                                          gtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt             >  1:676087/1‑150 (MQ=255)
                                                                                                                                           ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg            >  1:2217388/1‑150 (MQ=255)
                                                                                                                                            tttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc           >  1:766888/1‑150 (MQ=255)
                                                                                                                                                taTTTAAAAATATAATGAAATTTTTTTTTTTTAAATATCTAGATAACGAACAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa       <  1:1994752/149‑1 (MQ=255)
                                                                                                                                                gATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa       >  1:764121/1‑150 (MQ=255)
                                                                                                                                                 aTTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt      >  1:3449752/1‑150 (MQ=255)
                                                                                                                                                 aTTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt      >  1:4445950/1‑150 (MQ=255)
                                                                                                                                                     cAAAATATAATTAAATTATTTTTTTTTAATTATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg  <  2:3915226/150‑1 (MQ=255)
                                                                                                                                                     cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGACGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg  >  1:3662286/1‑150 (MQ=255)
                                                                                                                                                     cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg  >  1:3769135/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
GTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGG  >  CP009273/2170564‑2170862

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