Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,970,856 G→T 5.2% G407C (GGC→TGC)  wzxE → O‑antigen translocase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,970,8560GT 5.2% 1006.4 / 18.8 273G407C (GGC→TGC) wzxEO‑antigen translocase
Reads supporting (aligned to +/- strand):  ref base G (135/123);  new base T (3/12);  total (138/135)
Fisher's exact test for biased strand distribution p-value = 1.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

AAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTG  >  CP009273/3970707‑3971005
                                                                                                                                                     |                                                                                                                                                     
aaGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTtgtg                                                                                                                                                       >  2:2039727/1‑150 (MQ=255)
  gCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGc                                                                                                                                                     <  2:1250735/150‑1 (MQ=255)
   cGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCg                                                                                                                                                    >  1:2066279/1‑150 (MQ=255)
       aCTGCGGTTTTATATTCTGGCTGAAGTCAGCCAGTTCAATTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGAACTTGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGtt                                                                                                                                                <  1:420450/150‑1 (MQ=255)
       aCTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGtt                                                                                                                                                >  1:4915764/1‑150 (MQ=255)
         tGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGtttt                                                                                                                                              >  1:4152909/1‑150 (MQ=255)
          gCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGttttt                                                                                                                                             <  1:1092840/150‑1 (MQ=255)
             gTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTAct                                                                                                                                          <  1:2564607/150‑1 (MQ=255)
               tttATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTTGCGTGTTTTTActct                                                                                                                                        <  1:1240381/150‑1 (MQ=255)
                 tataTTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTgg                                                                                                                                      >  1:735009/1‑150 (MQ=255)
                     ttCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTa                                                                                                                                  >  1:2480262/1‑150 (MQ=255)
                       cTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTGTGTTGTGGCGTGTTTTTACTCTGGCGTAgg                                                                                                                                >  2:4945599/1‑150 (MQ=255)
                           cGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCCAATCCCTGCGCATGGTGCTCTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCggg                                                                                                                            <  1:4304987/150‑1 (MQ=255)
                            ggAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGc                                                                                                                           >  2:3858851/1‑150 (MQ=255)
                             gAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCa                                                                                                                          <  2:4601768/150‑1 (MQ=255)
                              aaGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTTGCGTGTTTTTACTCTGGCGTAGGCGGGCAt                                                                                                                         <  2:657288/150‑1 (MQ=255)
                              aaGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCAt                                                                                                                         >  2:245068/1‑150 (MQ=255)
                               aGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATg                                                                                                                        <  2:4464648/150‑1 (MQ=255)
                                gTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGa                                                                                                                       <  1:508319/150‑1 (MQ=255)
                                gTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGa                                                                                                                       <  1:1226487/150‑1 (MQ=255)
                                gTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGa                                                                                                                       >  1:4690525/1‑150 (MQ=255)
                                   aGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCACCTTATCTCGTCTATTTTTCTGTTGGGTGTGGCGTGGGTTTACTCTGGCGTAGGCGGGCATGACTg                                                                                                                    >  2:2156965/1‑150 (MQ=255)
                                   aGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTg                                                                                                                    <  1:5152156/150‑1 (MQ=255)
                                   aGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCAAGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTg                                                                                                                    >  2:3481167/1‑150 (MQ=255)
                                    gCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGt                                                                                                                   >  2:1582096/1‑150 (MQ=255)
                                     ccAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTa                                                                                                                  >  2:1288901/1‑150 (MQ=255)
                                     ccAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTa                                                                                                                  >  2:3221270/1‑150 (MQ=255)
                                       aGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACt                                                                                                                <  2:3188281/150‑1 (MQ=255)
                                       aGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACt                                                                                                                >  1:1487464/1‑150 (MQ=255)
                                       aGTTCACTGTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGAGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGAGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACt                                                                                                                >  2:4680669/1‑150 (MQ=255)
                                         ttCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGa                                                                                                              >  1:3621576/1‑150 (MQ=255)
                                          tCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGAt                                                                                                             <  1:2533523/150‑1 (MQ=255)
                                             cTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTNTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCa                                                                                                          <  1:2681949/150‑1 (MQ=255)
                                              tttATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCAc                                                                                                         <  2:1897336/150‑1 (MQ=255)
                                                tATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGt                                                                                                       <  1:2910160/150‑1 (MQ=255)
                                                 aTTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTa                                                                                                      <  1:2255467/150‑1 (MQ=255)
                                                  ttGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTAc                                                                                                     <  1:352087/150‑1 (MQ=255)
                                                  ttGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTAc                                                                                                     <  2:3370122/150‑1 (MQ=255)
                                                        gTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGAt                                                                                               >  1:1644533/1‑150 (MQ=255)
                                                         tATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTTGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATc                                                                                              <  2:3152019/150‑1 (MQ=255)
                                                         tATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATc                                                                                              <  1:811767/150‑1 (MQ=255)
                                                              gCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGata                                                                                         <  2:4515572/150‑1 (MQ=255)
                                                               cccACTGGCTAATCCCTGCGCAGGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGatat                                                                                        <  2:4336177/150‑1 (MQ=255)
                                                                    tGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTc                                                                                   <  2:3365512/150‑1 (MQ=255)
                                                                      gCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCAc                                                                                 <  1:3677632/150‑1 (MQ=255)
                                                                      gCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCAc                                                                                 <  2:3021957/150‑1 (MQ=255)
                                                                      gCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCAc                                                                                 <  1:4091150/150‑1 (MQ=255)
                                                                        tAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCa                                                                               >  2:2630561/1‑150 (MQ=255)
                                                                        tAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCa                                                                               >  2:694270/1‑150 (MQ=255)
                                                                        tAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCa                                                                               >  2:1646078/1‑150 (MQ=255)
                                                                          aTCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATa                                                                             <  1:495040/150‑1 (MQ=255)
                                                                          aTCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCGTTGTTGTGGCGTGTTGTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATa                                                                             >  2:382912/1‑150 (MQ=255)
                                                                           tCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATaa                                                                            <  2:1143400/150‑1 (MQ=255)
                                                                             ccTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTTGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATaacc                                                                          <  1:3309877/150‑1 (MQ=255)
                                                                             ccTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATaacc                                                                          <  1:599740/150‑1 (MQ=255)
                                                                               tGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATaaccga                                                                        >  1:159780/1‑150 (MQ=255)
                                                                                 cgcATGGTGCACTGGGCGCTGCGCATGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATaaccgaac                                                                      >  1:2830916/1‑150 (MQ=255)
                                                                                  gcATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATaaccgaacc                                                                     >  2:2843031/1‑150 (MQ=255)
                                                                                  gcATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATaaccgaacc                                                                     >  1:1875774/1‑150 (MQ=255)
                                                                                  gcATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATaaccgaacc                                                                     <  2:1762125/150‑1 (MQ=255)
                                                                                   cATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATaaccgaaccg                                                                    <  2:3003042/150‑1 (MQ=255)
                                                                                   cATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATaaccgaaccg                                                                    <  2:2920019/150‑1 (MQ=255)
                                                                                   cATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATaaccgaaccg                                                                    <  2:2919977/150‑1 (MQ=255)
                                                                                    aTGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGt                                                                   >  2:228677/1‑150 (MQ=255)
                                                                                     tGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGtt                                                                  >  1:3137240/1‑150 (MQ=255)
                                                                                     tGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGtt                                                                  <  1:3388270/150‑1 (MQ=255)
                                                                                       gTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGtttt                                                                <  1:1801049/150‑1 (MQ=255)
                                                                                           aCTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCgg                                                            >  1:4853323/1‑150 (MQ=255)
                                                                                            cTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGAACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGt                                                           >  2:2378539/1‑150 (MQ=255)
                                                                                                  gcggcgcAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCa                                                     <  2:1207735/150‑1 (MQ=255)
                                                                                                  gcggcgcAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCa                                                     >  1:785555/1‑150 (MQ=255)
                                                                                                    ggcgcAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAAt                                                   >  1:2111754/1‑150 (MQ=255)
                                                                                                    ggcgcAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAAt                                                   <  1:2135881/150‑1 (MQ=255)
                                                                                                       tcAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACAGAACCGTTTTGCGGTTTTTCAATGAc                                                >  2:4803698/2‑150 (MQ=255)
                                                                                                         aGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGAcgc                                              <  2:2078932/150‑1 (MQ=255)
                                                                                                           gCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACgcgc                                            >  1:1011798/1‑150 (MQ=255)
                                                                                                              tatatGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTgg                                         >  1:669862/1‑150 (MQ=255)
                                                                                                              tatatGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTgg                                         >  2:1383960/1‑150 (MQ=255)
                                                                                                              tatatGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTATCAATGACGCGCTgg                                         >  2:1861044/1‑150 (MQ=255)
                                                                                                                tatGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGcc                                       <  2:1218697/150‑1 (MQ=255)
                                                                                                                tatGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGcc                                       <  1:2388484/150‑1 (MQ=255)
                                                                                                                 atGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGcca                                      <  2:416931/150‑2 (MQ=255)
                                                                                                                  tGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCcgc                                     <  1:1335851/150‑1 (MQ=255)
                                                                                                                   ggCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCcgcg                                    >  2:3275233/1‑150 (MQ=255)
                                                                                                                          tataTCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGAcgagcg                             >  1:3189596/1‑150 (MQ=255)
                                                                                                                           atatCGTCTATTTTTCTCTTTGTTGTTGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGAcgagcga                            >  2:4017822/1‑150 (MQ=255)
                                                                                                                            tatCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGAcgagcgag                           <  2:887008/150‑1 (MQ=255)
                                                                                                                               cGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCTCCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCAc                        >  2:1031997/1‑150 (MQ=255)
                                                                                                                               cGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCAc                        <  1:3826999/150‑1 (MQ=255)
                                                                                                                                 tCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCAcgc                      >  1:2985922/1‑150 (MQ=255)
                                                                                                                                       tttCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACAGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGa                >  2:372151/1‑150 (MQ=255)
                                                                                                                                           tctTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGttt            <  2:2778347/150‑1 (MQ=255)
                                                                                                                                               tgttgtGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATgg        <  1:4056371/150‑1 (MQ=255)
                                                                                                                                               tgttgtGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATgg        <  1:4630353/150‑1 (MQ=255)
                                                                                                                                               tgttgtGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATgg        >  2:2550816/1‑150 (MQ=255)
                                                                                                                                                gttgtGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGgt       <  2:3674506/150‑1 (MQ=255)
                                                                                                                                                gttgtGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGgt       >  1:1983984/1‑150 (MQ=255)
                                                                                                                                                  tgtgGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGgttg     >  2:380932/1‑150 (MQ=255)
                                                                                                                                                    tgGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGgttgtt   <  2:412012/150‑1 (MQ=255)
                                                                                                                                                    tgGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGgttgtt   >  1:4427301/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATGGTGCACTGGGCGCGGCGCAGGCATATATGGCAACTTATATCGTCTATTTTTCTCTTTGTTGTGGCGTGTTTTTACTCTGGCGTAGGCGGGCATGACTGTACTGATTCACGTACTGGGATCGGATATCCCTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTG  >  CP009273/3970707‑3971005

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