Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,039,635 C→T 14.2% G11D (GGC→GAC)  etk ← tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,039,6350CT14.2% 727.1 / 92.3 254G11D (GGC→GAC) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base C (119/99);  new base T (22/14);  total (141/113)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 261 total aligned reads displayed.

GAGTATCTGCCTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACTTTTACGATAGGGATCCGGGA  >  CP009273/1039486‑1039783
                                                                                                                                                     |                                                                                                                                                    
 agTATCTGCCTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGcc                                                                                                                                                     >  1:2386621/1‑150 (MQ=255)
  gTATCTGCCTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCt                                                                                                                                                    >  2:1030373/1‑150 (MQ=255)
  gTATCTGCCTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCt                                                                                                                                                    >  2:1342013/1‑150 (MQ=255)
     tCTGCCTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGTCtggt                                                                                                                                                 <  1:184346/150‑1 (MQ=255)
      cTGCCTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCtggtg                                                                                                                                                <  1:4758821/150‑1 (MQ=255)
       tGCCTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCtggtgg                                                                                                                                               >  1:243942/1‑150 (MQ=255)
        gCCTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGc                                                                                                                                              >  2:2574512/1‑150 (MQ=255)
        gCCTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGc                                                                                                                                              >  2:3095701/1‑150 (MQ=255)
        gCCTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGc                                                                                                                                              >  1:4702796/1‑150 (MQ=255)
         ccTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCg                                                                                                                                             >  2:255269/1‑150 (MQ=255)
         ccTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCg                                                                                                                                             >  2:1920985/1‑150 (MQ=255)
            gATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTAt                                                                                                                                          >  2:1865941/1‑150 (MQ=255)
             aTAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTagt                                                                                                                                         >  2:4521338/1‑148 (MQ=255)
             aTAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTAtt                                                                                                                                         <  2:4228222/150‑1 (MQ=255)
             aTAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTAtt                                                                                                                                         >  2:2070254/1‑150 (MQ=255)
             aTAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTAtt                                                                                                                                         >  1:1260369/1‑150 (MQ=255)
             aTAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTAtt                                                                                                                                         >  1:1188221/1‑150 (MQ=255)
             aTAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGAGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTAtt                                                                                                                                         >  1:1005421/1‑150 (MQ=255)
              tAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTc                                                                                                                                        >  1:3045953/1‑150 (MQ=255)
               aaaTTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCa                                                                                                                                       <  1:3873967/150‑1 (MQ=255)
                aaTTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCat                                                                                                                                      >  2:468423/1‑150 (MQ=255)
                 aTTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCata                                                                                                                                     <  1:1006114/150‑1 (MQ=255)
                 aTTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCata                                                                                                                                     <  1:2383131/150‑1 (MQ=255)
                   tGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATAtt                                                                                                                                   <  1:292779/150‑1 (MQ=255)
                   tGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATAtt                                                                                                                                   >  1:467770/1‑150 (MQ=255)
                   tGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATAtt                                                                                                                                   <  2:238079/150‑1 (MQ=255)
                    ggTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATAttt                                                                                                                                  <  2:574050/150‑1 (MQ=255)
                     gtgtGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGTCTGGTGGCGTATTCATAtttt                                                                                                                                 >  2:2535726/1‑150 (MQ=255)
                     gtgtGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGTCTGGTGGCGTATTCATAtttt                                                                                                                                 >  1:2232907/1‑150 (MQ=255)
                      tgtgCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATAttttt                                                                                                                                >  1:32023/1‑150 (MQ=255)
                      tgtgCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATAttttt                                                                                                                                >  1:31478/1‑150 (MQ=255)
                       gtgCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTa                                                                                                                               >  2:4415810/1‑150 (MQ=255)
                        tgCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAg                                                                                                                              >  1:55073/1‑150 (MQ=255)
                         gCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGTCTGGTGGCGTATTCATATTTTTAGt                                                                                                                             >  1:3616597/1‑150 (MQ=255)
                                aGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGTCTGGTGGCGTATTCATATTTTTAGTTGTCATc                                                                                                                      <  2:3820790/150‑1 (MQ=255)
                                aGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGTCTGGTGGCGTATTCATATTTTTAGTTGTCATc                                                                                                                      <  2:3820839/150‑1 (MQ=255)
                                 gCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCg                                                                                                                     <  2:3244460/150‑1 (MQ=255)
                                    aGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTAt                                                                                                                  <  2:4581464/150‑1 (MQ=255)
                                      tAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTAttt                                                                                                                <  2:1077725/150‑1 (MQ=255)
                                       aaGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCCGACGACGCCGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTAtttt                                                                                                               >  2:927929/1‑150 (MQ=255)
                                       aaGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTAtttt                                                                                                               >  1:2894575/1‑150 (MQ=255)
                                        aGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTta                                                                                                              >  1:2252356/1‑150 (MQ=255)
                                        aGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTta                                                                                                              >  1:1546774/1‑150 (MQ=255)
                                         gCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTtat                                                                                                             >  2:2277077/1‑150 (MQ=255)
                                              aGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTCTTCATATTTTTAGTTGTCATCGTATTTTATAAACt                                                                                                        >  1:877183/1‑150 (MQ=255)
                                                     aGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTa                                                                                                 <  2:4745843/150‑1 (MQ=255)
                                                     aGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTa                                                                                                 <  2:2180651/150‑1 (MQ=255)
                                                       cGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTAcc                                                                                               <  1:4415810/150‑1 (MQ=255)
                                                       cGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTAcc                                                                                               >  2:2866917/1‑150 (MQ=255)
                                                        gTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCg                                                                                              <  2:3189007/150‑1 (MQ=255)
                                                          gAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGGTGTCATCGTATTTTCTAAACTCAAATTACCGGc                                                                                            >  1:1743757/1‑150 (MQ=255)
                                                            aTAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTg                                                                                          >  1:2481297/1‑150 (MQ=255)
                                                                         gCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGccc                                                                             >  2:3996530/1‑150 (MQ=255)
                                                                          cTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCa                                                                            >  2:667107/1‑150 (MQ=255)
                                                                             ataataAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTc                                                                         <  1:3377680/150‑1 (MQ=255)
                                                                               aataaACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCt                                                                       <  1:2432002/150‑1 (MQ=255)
                                                                                  aaaCTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGAc                                                                    <  2:539535/150‑1 (MQ=255)
                                                                                   aaCTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACt                                                                   >  1:3532650/1‑150 (MQ=255)
                                                                                     cTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTgg                                                                 <  2:3929062/150‑1 (MQ=255)
                                                                                     cTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTgg                                                                 >  2:3895934/1‑150 (MQ=255)
                                                                                      ttACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGgc                                                                <  2:3019785/150‑1 (MQ=255)
                                                                                      ttACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGgc                                                                <  2:1158321/150‑1 (MQ=255)
                                                                                       tACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGgcg                                                               <  1:728158/150‑1 (MQ=255)
                                                                                          ggTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGgcgcgc                                                            >  1:2692514/1‑150 (MQ=255)
                                                                                          ggTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGgcgcgc                                                            >  2:755689/1‑150 (MQ=255)
                                                                                           gTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCt                                                           >  1:3325273/1‑150 (MQ=255)
                                                                                               tCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATACCTGACTGGCGCGCTCCaa                                                       >  2:3264337/1‑150 (MQ=255)
                                                                                                   aTAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATa                                                   >  1:446462/1‑150 (MQ=255)
                                                                                                   aTAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTCCCTGAGTGCTGTCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATa                                                   >  2:1807607/1‑150 (MQ=255)
                                                                                                     aaCTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGTCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATAcc                                                 >  2:4384827/1‑150 (MQ=255)
                                                                                                     aaCTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGTCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATAcc                                                 >  2:2692276/1‑150 (MQ=255)
                                                                                                     aaCTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATAcc                                                 >  2:2240054/1‑150 (MQ=255)
                                                                                                        tCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTa                                              >  1:2089227/1‑150 (MQ=255)
                                                                                                         cGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAg                                             >  2:2041365/1‑150 (MQ=255)
                                                                                                         cGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAg                                             >  2:817835/1‑150 (MQ=255)
                                                                                                            cGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACa                                          >  1:2811119/1‑150 (MQ=255)
                                                                                                              aCCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGTCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATg                                        <  1:4264695/150‑1 (MQ=255)
                                                                                                                cAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGtt                                      >  2:4509468/1‑150 (MQ=255)
                                                                                                                 aGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGTCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTc                                     <  2:1275915/150‑1 (MQ=255)
                                                                                                                 aGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTc                                     >  1:3712937/1‑150 (MQ=255)
                                                                                                                      aaGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATg                                >  1:1092397/1‑150 (MQ=255)
                                                                                                                      aaGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATg                                >  2:2415795/1‑150 (MQ=255)
                                                                                                                          aGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTc                            >  1:3166862/1‑150 (MQ=255)
                                                                                                                           gatcgatcTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTcc                           <  2:306176/150‑1 (MQ=255)
                                                                                                                            atcgatcTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCt                          >  2:1604731/1‑150 (MQ=255)
                                                                                                                                atcTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACt                      <  1:69991/150‑1 (MQ=255)
                                                                                                                                atcTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACt                      <  1:70017/150‑1 (MQ=255)
                                                                                                                                 tctcATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACtt                     <  1:861761/150‑1 (MQ=255)
                                                                                                                                  ctcATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACttt                    <  2:2850029/150‑1 (MQ=255)
                                                                                                                                   tcATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACtttt                   <  2:3719534/150‑1 (MQ=255)
                                                                                                                                    cATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACTTTTa                  <  2:23202/150‑1 (MQ=255)
                                                                                                                                     aTTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACTTTTAc                 >  1:4255558/1‑150 (MQ=255)
                                                                                                                                     aTTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACTTTTAc                 >  1:1382442/1‑150 (MQ=255)
                                                                                                                                          ccTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACTTTTACGATAg            <  2:3454199/150‑1 (MQ=255)
                                                                                                                                           cTGAGTGCTGTCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACTTTTACGATAgg           >  1:4091059/1‑150 (MQ=255)
                                                                                                                                           cTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACTTTTACGATAgg           >  2:1217907/1‑150 (MQ=255)
                                                                                                                                             gAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACTTTTACGATAGGGa         <  1:3802969/150‑1 (MQ=255)
                                                                                                                                               gTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACTTTTACGATAGGGATc       <  2:243353/150‑1 (MQ=255)
                                                                                                                                               gTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACTTTTACGATAGGGATc       <  2:243973/150‑1 (MQ=255)
                                                                                                                                                 gctgcCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACTTTTACGATAGGGATCCg     >  2:2683777/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAGTATCTGCCTGATAAATTGGTGTGCTTAACAGCGAGTAAGCGACAGCGATCAGCGTGAATAACGCGGTCACGCTGATAATAAACTTACGGTGATCCCATAACTCGCCGACCAGACGAAGCAGATCGATCTCATTTTCCTGAGTGCTGCCTGGTGGCGTATTCATATTTTTAGTTGTCATCGTATTTTATAAACTCAAATTACCGGCTGAGGCGCTTCGCCCATTCCTGACTGGCGCGCTCCAACATACCGTAGACATGTTCAAATGCGTCCTGACTTTTACGATAGGGATCCGGGA  >  CP009273/1039486‑1039783

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: