Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,713 C→A 21.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,7130CA21.0% 459.9 / 19.3 214intergenic (‑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 (47/106);  new base A (37/14);  total (94/120)
Fisher's exact test for biased strand distribution p-value = 2.03e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

GTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGG  >  CP009273/2170564‑2170862
                                                                                                                                                      |                                                                                                                                                     
gTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTc                                                                                                                                                       <  1:5231273/150‑1 (MQ=255)
 ttGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCa                                                                                                                                                      <  1:1837872/150‑1 (MQ=255)
   ggCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCTATTTCaaa                                                                                                                                                    >  2:3035529/1‑150 (MQ=255)
    gCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaaa                                                                                                                                                   <  1:2508406/150‑1 (MQ=255)
    gCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCaaaa                                                                                                                                                   <  2:1620602/150‑1 (MQ=255)
     cAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAat                                                                                                                                                  <  1:4136557/150‑1 (MQ=255)
     cAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAat                                                                                                                                                  <  1:3929804/150‑1 (MQ=255)
      aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAata                                                                                                                                                 <  2:2898402/150‑1 (MQ=255)
      aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAata                                                                                                                                                 <  1:4609404/150‑1 (MQ=255)
      aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAata                                                                                                                                                 <  2:1883823/150‑1 (MQ=255)
      aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAata                                                                                                                                                 >  2:5231073/1‑150 (MQ=255)
      aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTTGTTTTTaaaaaa                                                                                                                                                 >  1:956429/1‑148 (MQ=255)
      aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTTaaaaaa                                                                                                                                                 >  1:1966044/1‑148 (MQ=255)
         ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATaa                                                                                                                                              <  2:4609871/150‑1 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                             <  2:171184/150‑1 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                             <  1:1896052/150‑1 (MQ=255)
          gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                             <  2:1974340/150‑1 (MQ=255)
           tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATg                                                                                                                                            <  2:4396549/150‑1 (MQ=255)
             tCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaa                                                                                                                                          <  2:195097/150‑1 (MQ=255)
              cTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGaaa                                                                                                                                         <  1:2563418/150‑1 (MQ=255)
               tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                        <  2:311195/150‑1 (MQ=255)
               tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                        <  1:5114587/150‑1 (MQ=255)
                accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATAAAATGaaaat                                                                                                                                       >  2:4694890/1‑148 (MQ=255)
                   accacTTGCATCGTTTCGAGGTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAAATAAattatt                                                                                                                                    >  2:549064/1‑150 (MQ=255)
                   accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTAAAAAAAAAAATAAaatatt                                                                                                                                    >  1:4485764/1‑150 (MQ=255)
                   accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTAAattatt                                                                                                                                    >  2:1252666/1‑150 (MQ=255)
                     cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTg                                                                                                                                  <  1:1181620/150‑1 (MQ=255)
                     cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTg                                                                                                                                  <  1:2042808/150‑1 (MQ=255)
                      acTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGt                                                                                                                                 <  1:2623744/150‑1 (MQ=255)
                       cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGtt                                                                                                                                <  2:2641319/150‑1 (MQ=255)
                       cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTAAAAAAATAATTTAAATAtttttt                                                                                                                                >  1:1369848/1‑150 (MQ=255)
                        ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttt                                                                                                                               <  2:4309435/150‑1 (MQ=255)
                        ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttt                                                                                                                               <  1:2916438/150‑1 (MQ=255)
                         tGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGtttt                                                                                                                              <  1:2083433/150‑1 (MQ=255)
                          gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttttt                                                                                                                             <  2:282762/150‑1 (MQ=255)
                           cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTa                                                                                                                            <  1:4375446/150‑1 (MQ=255)
                            aTCGTTTCGAGATTGTGAATATTGAATGCAGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAAAATAATGAAAATATTTTTTTTTaa                                                                                                                           >  2:2075386/1‑150 (MQ=255)
                             tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaaa                                                                                                                          <  1:1904588/150‑1 (MQ=255)
                             tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAATGAAAATAATTTTTTTTaaa                                                                                                                          >  1:1741365/1‑150 (MQ=255)
                               gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAata                                                                                                                        <  2:3869312/150‑1 (MQ=255)
                               gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTATAGTTAAAAAAAAAATTAAAATATTTTTTTTTaaaaa                                                                                                                        >  1:5375816/1‑148 (MQ=255)
                               gTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGGTTTGGATTTCAAAAAAAAAATAAAAAATTTTTTTTTTAata                                                                                                                        >  1:1472005/1‑150 (MQ=255)
                                tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                       >  2:3685752/1‑150 (MQ=255)
                                tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                       <  2:2698901/150‑1 (MQ=255)
                                tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                       <  2:3540870/150‑1 (MQ=255)
                                 ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATc                                                                                                                      <  1:3510800/150‑1 (MQ=255)
                                    gagaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgag                                                                                                                   >  2:3903056/1‑150 (MQ=255)
                                    gagaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgag                                                                                                                   <  1:3694327/150‑1 (MQ=255)
                                     agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGTAAATATTTTTTTTTTAATATAgata                                                                                                                  >  2:727723/1‑148 (MQ=255)
                                     agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGTATTTAAAAAAAAAAATAAAAAAATTTTTTTTAAATAAAgaga                                                                                                                  >  1:2626910/1‑150 (MQ=255)
                                      gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGAt                                                                                                                 >  2:351946/1‑150 (MQ=255)
                                      gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGAt                                                                                                                 <  2:5207346/150‑1 (MQ=255)
                                          gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTAAAAAAATAATGAAATTATTTTTTTTTTAAAATATAGATaaac                                                                                                             >  1:2269375/1‑148 (MQ=255)
                                             aaTATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAATATAATTAAATTATTTGTTTTTAAATATAGAGATAAAGATc                                                                                                          >  2:4059908/1‑150 (MQ=255)
                                             aaTATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATc                                                                                                          <  2:1676653/150‑1 (MQ=255)
                                                atTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATcaca                                                                                                       <  1:2230691/150‑1 (MQ=255)
                                                    aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGGGTTTTAGAATTCAAAAAAAAAAGAAAAAATATTTTTTTAAATAAAGAGATAACCATAACaaaaa                                                                                                   >  1:2888698/1‑150 (MQ=255)
                                                      tGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACg                                                                                                 <  2:1818957/150‑1 (MQ=255)
                                                        ccGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTGTTTTGATTTAAAAAAATAAATAAATTATTTTTTTTTTAATTTAGATATAAAGATaaaaaaaaaaaa                                                                                               >  1:3410832/1‑146 (MQ=255)
                                                        ccGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGATTTAAAAATATAAAGAAAATATTTTTTTTTTAAAATAGAGATAAAGAACACaaaaaacac                                                                                               >  2:1553060/1‑146 (MQ=255)
                                                               cGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGATTTTAAAATATAATTAAAATATTTTTTTTTTAATATAGAGATAAAAATAAAAAAAAAAACAaaattt                                                                                        >  1:1618820/1‑148 (MQ=255)
                                                                  aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACTTAAATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGaaaa                                                                                     <  2:2239353/150‑1 (MQ=255)
                                                                   tAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGTGATAACGATCACAAAAACGACAATATGAAAat                                                                                    <  2:3213467/150‑1 (MQ=255)
                                                                    acccgccgCGGTGTGCTTTTTTCAGCATCTGCTTTGTCTATACCACTTAAATTTTAATATCCTTTCTTTTTTTTTCTTTTTCAAAATATAATTAATTTTTTTTTTTTTAAATATCGAGATAACTATCAAAAAAACGACAATATGAAAatt                                                                                   <  1:3002752/148‑1 (MQ=255)
                                                                      ccgccgCGCTGGGCGTTGTTCACAATCTGCTTTGTCGATACCACTTAAATTTTCATACCCTGTCGTTTTTTTTCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAACACGACAATATGAAAattat                                                                                 <  1:1845748/150‑1 (MQ=255)
                                                                       cgccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTAAAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAattatt                                                                                <  2:3823816/150‑1 (MQ=255)
                                                                        gccgcGCTGTGAGTTGTTCAGCATGTGCTTTGTCTATACCACTTACATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATTAAATTTTTTTTTTTTAAATATCGAGATAACGATCAAAAAAACGACAATATGAAAATTATTc                                                                               <  2:4407338/150‑1 (MQ=255)
                                                                         ccgcGCTGTGCTTTGTTCCGAATCGGCTTTGTCTATACCACGTACATTTTCATAACCTGTCGTTTTTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGATATAACGAACACAAACACGACAATATGAAAATTATTCg                                                                              <  2:5180038/150‑1 (MQ=255)
                                                                         ccgcGCTGTGCGTTTTTCAGCATCTGCTTTGGCGATACCACGTACATTTTCATATCCTGTCTTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCTAGATAACGATCACAAAAACGACAATATGAAAATTATTCg                                                                              <  2:5180739/150‑1 (MQ=255)
                                                                          cgcgCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAAAATAATTAAATAATTTTTTTTTAAATAAAGAGAGAAACAAAACAAAAAAAACAAAAAGAAAaataatagc                                                                             >  2:948574/1‑145 (MQ=255)
                                                                              cTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGAATTAAAAAAAAAATGTAAATATTTTTTTTTTAAAAAAGATAAAACCATAACAAAAACAAAAAAAATAAAAATATTCgatgag                                                                         >  2:17348/1‑150 (MQ=255)
                                                                                  ggtttgttCAGCATCGGATTTGTCGATACCACCTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGaa                                                                     <  2:2450769/147‑1 (MQ=255)
                                                                                      tgttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAAAATAATGAAATTATTTTTTTTTAAATATAGAGATAACGATCACAAAAAAGACAAAATGAAAATTATTCGAGGAGTGAAAGGc                                                                 >  2:4891340/1‑150 (MQ=255)
                                                                                       gttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATAGAGATAACGATCACAAAAAAGACAATATTAAAATTATTAGAGGAGTGAAAGGCa                                                                >  2:866031/1‑150 (MQ=255)
                                                                                            gCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaaaaaa                                                           <  2:4900168/150‑1 (MQ=255)
                                                                                                             aCCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAggg                                          <  2:3610167/150‑1 (MQ=255)
                                                                                                                ccctaaaTTTTAATATCCTGTGTTTTTTTTCCTATTACAAAATATAATTAATTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAg                                       <  2:956429/144‑1 (MQ=255)
                                                                                                                  ttaatttttCATATCGTGTCTTTTTTTTTCTTTTTAAAAATATATTTAATTTTTTTTTTTTTAAATATCTATAAAACTATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGcc                                     <  1:226387/145‑1 (MQ=255)
                                                                                                                  gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGcc                                     <  2:4534478/150‑1 (MQ=255)
                                                                                                                     cATTTTCATATCCTGTCGTGGGTTTTGGATTTCAAAATATAATGTAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACCATATGAAAAGTATTGGAGGAGTGCAAGGCCAAAAAACGGCCTCCCGCGAGGGAAGCCGTa                                  >  2:4519027/1‑150 (MQ=255)
                                                                                                                     aaTTTTCAAACCCTGTCGTTTTTTTACGAATAAAAAATATAATTAAATTATTTTTTTTTAAATATCGATATAACGAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                  <  1:4522585/149‑1 (MQ=255)
                                                                                                                        tttCATATCCGGTCGTTTGTTTTCTATTTCAAAATATAATGAATTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCa                               <  2:4944734/150‑1 (MQ=255)
                                                                                                                              atCCTGTCGTTTGTTTTCGATTTAAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgc                         <  2:796783/150‑1 (MQ=255)
                                                                                                                                ccTGTCGTTTGTTTTCGATTTCAAAATATAATGAAAATATTTGTTTTTAAATATCGAGATAACGATAACAAAAACGACAATATGAAAATTATTCGAGGAGTGTAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTATCAAAGTgcgc                       >  2:1370904/1‑150 (MQ=255)
                                                                                                                                    tCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtgt                   <  2:2623540/150‑1 (MQ=255)
                                                                                                                                     cGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                  >  2:2293306/1‑150 (MQ=255)
                                                                                                                                       tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                >  1:207585/1‑150 (MQ=255)
                                                                                                                                         ttgtttTCGTTTTCAAAATATAATGAAATTTTTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt               <  1:5534083/150‑1 (MQ=255)
                                                                                                                                         ttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt               >  1:1507961/1‑150 (MQ=255)
                                                                                                                                          tgtttTCTATTTCAAAATATAATGAATTTTTTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa              <  2:2228871/150‑1 (MQ=255)
                                                                                                                                          tgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAATCGGCCTCCCGATCGGGAAGCCGTAGCAACGTGCGCGTGTTTTTa              >  2:1474376/1‑150 (MQ=255)
                                                                                                                                          tgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa              >  2:5481987/1‑150 (MQ=255)
                                                                                                                                           gtttTCGATTTCAAAATATAATGAATTTATTTGTTTTTAATTATCGAGATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt             <  2:1392851/150‑1 (MQ=255)
                                                                                                                                            ttttCTTTTTAAAAATATAATTAAATTATTTTTTTTTAAATATCTATATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg            <  1:866031/150‑1 (MQ=255)
                                                                                                                                            ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg            <  2:4455534/150‑1 (MQ=255)
                                                                                                                                              ttCTATTTAAAAATATAATTAAATTATTTTTTTTTAAATATCGATATAACGAACACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc          <  1:5132615/150‑1 (MQ=255)
                                                                                                                                              ttCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc          <  2:3587016/150‑1 (MQ=255)
                                                                                                                                                cGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCgg        >  1:657973/1‑150 (MQ=255)
                                                                                                                                                  aTTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt      >  2:2798032/1‑150 (MQ=255)
                                                                                                                                                      cAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg  >  2:1445046/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 < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: