Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,168,795 C→A 93.3% 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,7950CA93.3% 951.0 / 53.7 309D336Y (GAT→TAT) gatZD‑tagatose 1,6‑bisphosphate aldolase 2, subunit
Reads supporting (aligned to +/- strand):  ref base C (7/14);  new base A (134/153);  total (141/168)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 310 total aligned reads displayed.

TTGGGAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTAAAG  >  CP009273/2168647‑2168941
                                                                                                                                                    |                                                                                                                                                  
   ggAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTa                                                                                                                                                >  2:3202522/1‑150 (MQ=255)
     aaGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTaaa                                                                                                                                              <  1:2428998/150‑1 (MQ=255)
      aGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTaaaa                                                                                                                                             <  1:1029298/150‑1 (MQ=255)
      aGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTaaaa                                                                                                                                             <  2:4478235/150‑1 (MQ=255)
         tACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAAccc                                                                                                                                          <  1:4345929/150‑1 (MQ=255)
         tACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAAccc                                                                                                                                          >  1:1539462/1‑150 (MQ=255)
         tACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAAccc                                                                                                                                          >  1:3846921/1‑150 (MQ=255)
          aCTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCg                                                                                                                                         >  1:2384235/1‑150 (MQ=255)
            tGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTa                                                                                                                                       >  2:769146/1‑150 (MQ=255)
             gACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTAc                                                                                                                                      >  1:518321/1‑150 (MQ=255)
               cTAATAATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGa                                                                                                                                    <  2:4088460/150‑1 (MQ=255)
                  aTCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGataa                                                                                                                                 >  1:2056319/1‑150 (MQ=255)
                   tCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGataat                                                                                                                                <  2:3903064/150‑1 (MQ=255)
                     aTGCCCAGTGGGATGTCCACTCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAatat                                                                                                                              <  1:402265/150‑1 (MQ=255)
                      tGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATAtt                                                                                                                             <  1:310385/150‑1 (MQ=255)
                       gCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATAttt                                                                                                                            <  2:4429756/150‑1 (MQ=255)
                        cccAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATAtttt                                                                                                                           >  2:1645633/1‑150 (MQ=255)
                        cccAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATAtttt                                                                                                                           <  2:1793203/150‑1 (MQ=255)
                           aGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTc                                                                                                                        >  2:4278890/1‑150 (MQ=255)
                           aGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTc                                                                                                                        >  1:5036373/1‑150 (MQ=255)
                            gTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTcc                                                                                                                       <  2:4611985/150‑1 (MQ=255)
                             tGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCa                                                                                                                      <  1:901574/150‑1 (MQ=255)
                              gggATGTCCACGCATTCAAGATTAACCATCATCGTTTCGAAGCTATTTTTAATCAGACTATGCGGCCAGTAATAACGAATACGATGCGACAGGGTTTAACGATTATCCAGTTATGAATAGTTAAAACCCTTACGATAATATTTTTTCCAg                                                                                                                     <  2:1114435/150‑1 (MQ=255)
                                  tGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTAtt                                                                                                                 <  2:1782651/150‑1 (MQ=255)
                                  tGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTAtt                                                                                                                 <  2:197316/150‑1 (MQ=255)
                                       aCGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGt                                                                                                            >  2:2685911/1‑150 (MQ=255)
                                          ccTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTtcg                                                                                                         <  1:4814025/150‑1 (MQ=255)
                                            ttCTAGATTCACCTTCATCGTTTCGTCGCTATTTTTTATCCGACTATGCGTACAGTAATAACGAATACTTTCCGACAGGCTGTAACGAATATCCATTAATTATTTGTTAATTCCCGTACGTTAATTTTTTTTCCAGTATAGCGGTtcgtc                                                                                                       <  1:1487233/150‑1 (MQ=255)
                                              cAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGa                                                                                                     <  2:898547/150‑1 (MQ=255)
                                                aGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGc                                                                                                   <  1:2873810/150‑1 (MQ=255)
                                                 gATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCa                                                                                                  >  2:2797212/1‑150 (MQ=255)
                                                 gATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCa                                                                                                  >  2:813422/1‑150 (MQ=255)
                                                     cACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCAAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCAc                                                                                              <  1:3169350/150‑1 (MQ=255)
                                                      aCCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCAct                                                                                             >  2:1291416/1‑150 (MQ=255)
                                                           catcGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACttcttc                                                                                        >  1:1917697/1‑150 (MQ=255)
                                                            atcGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCa                                                                                       >  1:16018/1‑150 (MQ=255)
                                                            atcGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCa                                                                                       >  1:15991/1‑150 (MQ=255)
                                                               gTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAAtt                                                                                    >  2:2519997/1‑150 (MQ=255)
                                                               gTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAAtt                                                                                    >  2:3673335/1‑150 (MQ=255)
                                                                  tCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTAcc                                                                                 <  1:221272/150‑1 (MQ=255)
                                                                   cGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCg                                                                                <  1:1733472/150‑1 (MQ=255)
                                                                   cGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCg                                                                                <  1:4518841/150‑1 (MQ=255)
                                                                      cGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCa                                                                             <  1:1104011/150‑1 (MQ=255)
                                                                      cGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCa                                                                             <  1:4569258/150‑1 (MQ=255)
                                                                            ttttAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAAc                                                                       <  2:1013780/150‑1 (MQ=255)
                                                                             tttAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAAcc                                                                      <  2:2024514/150‑1 (MQ=255)
                                                                               tAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGc                                                                    <  1:3241526/150‑1 (MQ=255)
                                                                                aaTCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTACTGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCt                                                                   <  2:3852208/150‑1 (MQ=255)
                                                                                aaTCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCt                                                                   >  1:2587600/1‑150 (MQ=255)
                                                                                 aTCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTg                                                                  >  2:1195705/1‑150 (MQ=255)
                                                                                  tCCGACTATGCGGACAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGc                                                                 <  1:352010/150‑1 (MQ=255)
                                                                                    cGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGa                                                               >  2:1903025/1‑150 (MQ=255)
                                                                                     gACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGAt                                                              >  1:1430317/1‑150 (MQ=255)
                                                                                     gACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGAt                                                              >  2:3223450/1‑150 (MQ=255)
                                                                                      aCTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGAtt                                                             <  1:3396836/150‑1 (MQ=255)
                                                                                      aCTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGAtt                                                             <  2:3116120/150‑1 (MQ=255)
                                                                                      aCTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGAtt                                                             <  1:2840145/150‑1 (MQ=255)
                                                                                        tATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGAtttt                                                           <  2:4544446/150‑1 (MQ=255)
                                                                                          tGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCa                                                         >  1:597435/1‑150 (MQ=255)
                                                                                               ccAGTAATAACGAAAACGATCCGACAGGCTGTAACGAATAACCAGTAATGAATAGTTAAAACCCGTACGATAATATTATTTCCAGTATTGCGGAACGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGc                                                    >  2:1483651/1‑150 (MQ=255)
                                                                                                 aGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGa                                                  >  1:509647/1‑150 (MQ=255)
                                                                                                 aGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGa                                                  >  2:624010/1‑150 (MQ=255)
                                                                                                      taaCGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGt                                             >  1:343342/1‑150 (MQ=255)
                                                                                                        aCGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTc                                           >  2:2831911/1‑150 (MQ=255)
                                                                                                         cGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTcc                                          >  1:2082471/1‑150 (MQ=255)
                                                                                                         cGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTcc                                          <  2:705259/150‑1 (MQ=255)
                                                                                                             tACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCt                                      <  1:1835687/150‑1 (MQ=255)
                                                                                                               cGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCa                                    <  2:2915713/150‑1 (MQ=255)
                                                                                                                gATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCaa                                   >  1:1655638/1‑150 (MQ=255)
                                                                                                                gATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCaa                                   <  2:3467670/150‑1 (MQ=255)
                                                                                                                 aTCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAAt                                  >  2:4799829/1‑150 (MQ=255)
                                                                                                                 aTCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAAt                                  >  2:5073862/1‑150 (MQ=255)
                                                                                                                  tCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAAtt                                 <  2:2238929/150‑1 (MQ=255)
                                                                                                                    cGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTtg                               <  1:4918842/150‑1 (MQ=255)
                                                                                                                     gACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTtgt                              >  1:2474880/1‑150 (MQ=255)
                                                                                                                       cAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGc                            >  2:3061998/1‑150 (MQ=255)
                                                                                                                         ggCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCa                          <  1:822493/150‑1 (MQ=255)
                                                                                                                           cTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGt                        <  2:2953737/150‑1 (MQ=255)
                                                                                                                             gTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATACGTTCCTGCTCAATTTGTGCCAGTGc                      >  2:4060741/1‑150 (MQ=255)
                                                                                                                             gTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGc                      <  2:892448/150‑1 (MQ=255)
                                                                                                                             gTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGc                      >  2:4899590/1‑150 (MQ=255)
                                                                                                                              tAACGAATATCCAGTAATGAATAGTTAAAACCCGTCCGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCa                     <  1:1715039/150‑1 (MQ=255)
                                                                                                                              tAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCa                     <  1:4276523/150‑1 (MQ=255)
                                                                                                                               aaCGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCaa                    >  1:1163112/1‑150 (MQ=255)
                                                                                                                               aaCGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCaa                    >  2:3620176/1‑150 (MQ=255)
                                                                                                                                aCGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAAACGCTGAGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCaaa                   >  1:2956429/1‑150 (MQ=255)
                                                                                                                                 cGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAat                  <  1:2256189/150‑1 (MQ=255)
                                                                                                                                   aaTATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAatat                >  1:3360788/1‑150 (MQ=255)
                                                                                                                                    atatCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATc               >  1:2901211/1‑150 (MQ=255)
                                                                                                                                    atatCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATc               <  2:2116297/150‑1 (MQ=255)
                                                                                                                                    atatCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATc               >  1:2901489/1‑150 (MQ=255)
                                                                                                                                     tatCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATAACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCg              >  1:787387/1‑150 (MQ=255)
                                                                                                                                      atCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGc             <  2:2413255/150‑1 (MQ=255)
                                                                                                                                        ccAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCt           <  1:3319632/150‑1 (MQ=255)
                                                                                                                                          aGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTcg         >  2:3748900/1‑150 (MQ=255)
                                                                                                                                             aaTGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGt      <  1:3701348/150‑1 (MQ=255)
                                                                                                                                             aaTGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGt      >  2:2780391/1‑150 (MQ=255)
                                                                                                                                              aTGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTa     <  1:420623/150‑1 (MQ=255)
                                                                                                                                                gAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTaaa   >  2:3394338/1‑150 (MQ=255)
                                                                                                                                                 aaTCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTAAAg  <  1:3031930/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTGGGAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTAAAG  >  CP009273/2168647‑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: