Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,168,795 C→A 84.0% D336Y (GAT→TAT)  gatZ ← D‑tagatose 1,6‑bisphosphate aldolase 2, subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,168,7950CA84.0% 488.1 / 97.6 201D336Y (GAT→TAT) gatZD‑tagatose 1,6‑bisphosphate aldolase 2, subunit
Reads supporting (aligned to +/- strand):  ref base C (10/22);  new base A (80/89);  total (90/111)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 201 total aligned reads displayed.

TGGGAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTAAAG  >  CP009273/2168648‑2168941
                                                                                                                                                   |                                                                                                                                                  
tGGGAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGt                                                                                                                                                  <  1:2356738/150‑1 (MQ=255)
  ggAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTa                                                                                                                                                <  2:4478431/150‑1 (MQ=255)
  ggAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTa                                                                                                                                                >  1:480207/1‑150 (MQ=255)
    aaGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTaaa                                                                                                                                              >  2:1884974/1‑150 (MQ=255)
    aaGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCAAGTAATAACGAATACGATCCGACATGCTGTAACTAATATCCAGTAATGAATAGTTaaa                                                                                                                                              >  2:1201132/1‑150 (MQ=255)
     aGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTaaaa                                                                                                                                             <  1:4052540/150‑1 (MQ=255)
      gATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAAc                                                                                                                                            <  1:1369808/150‑1 (MQ=255)
       aTACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAAcc                                                                                                                                           <  1:3222404/150‑1 (MQ=255)
          cTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGt                                                                                                                                        >  2:2767034/1‑150 (MQ=255)
           tGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTa                                                                                                                                       >  2:672345/1‑150 (MQ=255)
            gACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTAc                                                                                                                                      <  1:2636640/150‑1 (MQ=255)
            gACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTAc                                                                                                                                      <  2:2601052/150‑1 (MQ=255)
             aCTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTCATGAATAGTTAAAACCCGTACg                                                                                                                                     >  2:435248/1‑150 (MQ=255)
             aCTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACg                                                                                                                                     <  1:395271/150‑1 (MQ=255)
             aCTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACg                                                                                                                                     <  2:199404/150‑1 (MQ=255)
              cTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATCTCCAGTAATGAATAGTTAAAACCCGTACGa                                                                                                                                    >  2:2335737/1‑150 (MQ=255)
                aaTCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGata                                                                                                                                  <  2:1522542/150‑1 (MQ=255)
                aaTCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGata                                                                                                                                  <  1:1849489/150‑1 (MQ=255)
                 aTCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGataa                                                                                                                                 <  1:4334806/150‑1 (MQ=255)
                 aTCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGataa                                                                                                                                 >  1:2073706/1‑150 (MQ=255)
                  tCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGataat                                                                                                                                <  2:1212488/150‑1 (MQ=255)
                     tGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATAtt                                                                                                                             <  1:2283144/150‑1 (MQ=255)
                      gCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATAttt                                                                                                                            <  2:835039/150‑1 (MQ=255)
                        ccAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATAttttt                                                                                                                          <  2:1540383/150‑1 (MQ=255)
                            tGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCa                                                                                                                      >  1:3586129/1‑150 (MQ=255)
                            tGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCa                                                                                                                      >  2:3325199/1‑150 (MQ=255)
                             gggATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAg                                                                                                                     <  2:3621774/150‑1 (MQ=255)
                             gggATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAg                                                                                                                     >  1:1437610/1‑150 (MQ=255)
                               gATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTa                                                                                                                   <  1:2401824/150‑1 (MQ=255)
                               gATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTa                                                                                                                   <  1:3803187/150‑1 (MQ=255)
                                   tCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGc                                                                                                               >  1:2193529/1‑150 (MQ=255)
                                     cACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCgg                                                                                                             <  2:1784432/150‑1 (MQ=255)
                                     cACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTCAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCgg                                                                                                             >  2:88107/1‑150 (MQ=255)
                                       cGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGtt                                                                                                           >  1:3842803/1‑150 (MQ=255)
                                       cGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGtt                                                                                                           >  2:2507807/1‑150 (MQ=255)
                                       cGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGtt                                                                                                           >  2:2508345/1‑150 (MQ=255)
                                           ttCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTtcgtc                                                                                                       <  2:1234257/150‑1 (MQ=255)
                                              aaGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAg                                                                                                    <  2:2962554/150‑1 (MQ=255)
                                               aGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGc                                                                                                   >  1:2420577/1‑150 (MQ=255)
                                                gATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCa                                                                                                  <  2:656889/150‑1 (MQ=255)
                                                gATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCa                                                                                                  <  2:2956830/150‑1 (MQ=255)
                                                 aTTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCAt                                                                                                 >  2:4056055/1‑150 (MQ=255)
                                                      ccATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCActt                                                                                            <  2:1311827/150‑1 (MQ=255)
                                                       catcatCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCActtc                                                                                           <  1:4466632/150‑1 (MQ=255)
                                                        atcatcGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCActtct                                                                                          >  2:2454376/1‑150 (MQ=255)
                                                        atcatcGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCActtct                                                                                          >  2:3766357/1‑150 (MQ=255)
                                                           atcGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCa                                                                                       <  2:3309400/150‑1 (MQ=255)
                                                           atcGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCa                                                                                       >  2:2131758/1‑150 (MQ=255)
                                                             cGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAAt                                                                                     <  2:3982694/150‑1 (MQ=255)
                                                                ttCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTAc                                                                                  >  2:2296216/1‑150 (MQ=255)
                                                                  cGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCg                                                                                <  2:1927062/150‑1 (MQ=255)
                                                                     cGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCa                                                                             >  2:3477535/1‑150 (MQ=255)
                                                                      gCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGATTAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAg                                                                            <  2:145889/150‑1 (MQ=255)
                                                                      gCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAg                                                                            >  1:737473/1‑150 (MQ=255)
                                                                      gCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAg                                                                            >  2:768961/1‑150 (MQ=255)
                                                                      gCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAg                                                                            >  1:2997628/1‑150 (MQ=255)
                                                                              tAATCCGACTATGCGGCCAGTTATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGc                                                                    >  2:2208402/1‑150 (MQ=255)
                                                                               aaTCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCt                                                                   <  2:3109720/150‑1 (MQ=255)
                                                                                   cGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGa                                                               >  2:2990018/1‑150 (MQ=255)
                                                                                    gACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGAt                                                              >  1:953780/1‑150 (MQ=255)
                                                                                      cTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGAttt                                                            >  2:4078087/1‑150 (MQ=255)
                                                                                      cTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGAttt                                                            <  2:1650248/150‑1 (MQ=255)
                                                                                       tATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGAtttt                                                           <  1:1207333/150‑1 (MQ=255)
                                                                                        aTGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTc                                                          >  1:1167516/1‑150 (MQ=255)
                                                                                        aTGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTc                                                          >  1:1713550/1‑150 (MQ=255)
                                                                                            ggCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAggg                                                      <  1:4184324/150‑1 (MQ=255)
                                                                                              ccAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGc                                                    <  1:3407736/150‑1 (MQ=255)
                                                                                              ccAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGc                                                    <  2:2455155/150‑1 (MQ=255)
                                                                                               cAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCg                                                   <  2:777944/150‑1 (MQ=255)
                                                                                                  taataaCGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATa                                                >  2:366690/1‑150 (MQ=255)
                                                                                                   aataaCGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATaa                                               <  2:178786/150‑1 (MQ=255)
                                                                                                       aCGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTc                                           >  1:1530387/1‑150 (MQ=255)
                                                                                                        cGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTcc                                          <  1:2453524/150‑1 (MQ=255)
                                                                                                        cGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTcc                                          <  1:3328629/150‑1 (MQ=255)
                                                                                                         gAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCt                                         >  1:3898371/1‑150 (MQ=255)
                                                                                                              cGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCa                                    <  2:3586129/150‑1 (MQ=255)
                                                                                                               gATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCaa                                   <  1:556943/150‑1 (MQ=255)
                                                                                                                aTCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAAt                                  >  2:3013668/1‑150 (MQ=255)
                                                                                                                aTCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAAt                                  <  1:2862280/150‑1 (MQ=255)
                                                                                                                 tCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAAtt                                 >  1:852358/1‑150 (MQ=255)
                                                                                                                  ccGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAAttt                                >  2:4454035/1‑150 (MQ=255)
                                                                                                                  ccGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAAttt                                <  1:2131758/150‑1 (MQ=255)
                                                                                                                   cGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTtg                               >  1:97314/1‑150 (MQ=255)
                                                                                                                   cGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTtg                               <  2:3416524/150‑1 (MQ=255)
                                                                                                                       aGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGcc                           <  2:4534765/150‑1 (MQ=255)
                                                                                                                        ggCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCa                          >  1:1614985/1‑150 (MQ=255)
                                                                                                                          cTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGt                        <  2:799612/150‑1 (MQ=255)
                                                                                                                             tAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCa                     <  1:112485/150‑1 (MQ=255)
                                                                                                                             tAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCa                     >  1:605403/1‑150 (MQ=255)
                                                                                                                              aaCGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCaa                    >  2:4398196/1‑150 (MQ=255)
                                                                                                                                  aaTATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAatat                <  2:151942/150‑1 (MQ=255)
                                                                                                                                  aaTATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAatat                <  2:3441132/150‑1 (MQ=255)
                                                                                                                                    tatCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCg              <  2:2860814/150‑1 (MQ=255)
                                                                                                                                       ccAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAGTTTGTGCCAGTGCAAATATCGCCt           <  1:2253305/150‑1 (MQ=255)
                                                                                                                                       ccAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCt           >  1:3543277/1‑150 (MQ=255)
                                                                                                                                       ccAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCt           <  2:780694/150‑1 (MQ=255)
                                                                                                                                        cAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTc          >  1:1554890/1‑150 (MQ=255)
                                                                                                                                           tAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTcgcg       >  2:2926957/1‑150 (MQ=255)
                                                                                                                                           tAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTcgcg       >  2:2054219/1‑150 (MQ=255)
                                                                                                                                                aaTAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTAAAg  >  2:3093996/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGGGAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTAAAG  >  CP009273/2168648‑2168941

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: