Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,534 C→A 29.1% T50T (ACG→ACT gatY ← D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,170,5340CA29.1% 573.8 / 209.9 296T50T (ACG→ACTgatYD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (116/94);  new base A (41/45);  total (157/139)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

Only 100 of 298 total aligned reads displayed.

AGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACAT  >  CP009273/2170386‑2170683
                                                                                                                                                    |                                                                                                                                                     
aGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCg                                                                                                                                                      >  2:4584515/1‑150 (MQ=255)
    cACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTcc                                                                                                                                                  <  2:3002246/150‑1 (MQ=255)
    cACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTcc                                                                                                                                                  >  1:2327193/1‑150 (MQ=255)
     aCGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCg                                                                                                                                                 <  2:3828002/150‑1 (MQ=255)
     aCGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCg                                                                                                                                                 >  1:3945569/1‑150 (MQ=255)
            aacgaacCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCgatga                                                                                                                                          <  2:3796414/150‑1 (MQ=255)
              cgaacCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGatgatg                                                                                                                                        >  1:5318063/1‑150 (MQ=255)
               gaacCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGAtgatga                                                                                                                                       >  1:3275757/1‑150 (MQ=255)
               gaacCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGAtgatga                                                                                                                                       >  2:3995160/1‑150 (MQ=255)
               gaacCTTCTGAGCGATATAGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGAtgatga                                                                                                                                       <  2:4742787/150‑1 (MQ=255)
                aacCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGAc                                                                                                                                      >  1:2270642/1‑150 (MQ=255)
                 acCTTCTGAGCGATATCGTCAACTGTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGAcc                                                                                                                                     <  2:1836072/150‑1 (MQ=255)
                 acCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGAcc                                                                                                                                     <  1:768386/150‑1 (MQ=255)
                       tGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGTAGTTCCGGCGATGATGACCGGCGCa                                                                                                                               >  2:893264/1‑150 (MQ=255)
                       tGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCa                                                                                                                               >  1:2158370/1‑150 (MQ=255)
                          gCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGc                                                                                                                            <  1:1259567/150‑1 (MQ=255)
                            gATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAg                                                                                                                          >  2:2191078/1‑150 (MQ=255)
                             atatCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAgg                                                                                                                         <  1:4946617/150‑1 (MQ=255)
                             atatCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAgg                                                                                                                         <  2:2685539/150‑1 (MQ=255)
                             atatCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAgg                                                                                                                         >  1:3071908/1‑150 (MQ=255)
                             atatCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAgg                                                                                                                         >  2:1241669/1‑150 (MQ=255)
                               atCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGtt                                                                                                                       >  2:1897880/1‑150 (MQ=255)
                               atCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGtt                                                                                                                       <  1:3978440/150‑1 (MQ=255)
                               atCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGtt                                                                                                                       <  1:4339078/150‑1 (MQ=255)
                               atCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGtt                                                                                                                       >  1:4518662/1‑150 (MQ=255)
                                    cAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAg                                                                                                                  >  1:2825073/1‑150 (MQ=255)
                                    cAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAg                                                                                                                  >  2:4770536/1‑150 (MQ=255)
                                     aaaTTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGc                                                                                                                 >  2:2065333/1‑150 (MQ=255)
                                     aaaTTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGc                                                                                                                 >  2:2562379/1‑150 (MQ=255)
                                       aTTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCgg                                                                                                               >  2:3642204/1‑150 (MQ=255)
                                          tCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGttt                                                                                                            <  1:5179474/150‑1 (MQ=255)
                                          tCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGttt                                                                                                            >  1:4677139/1‑150 (MQ=255)
                                              gtgATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTac                                                                                                        >  2:3888591/1‑150 (MQ=255)
                                               tgATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTacc                                                                                                       <  2:4164412/150‑1 (MQ=255)
                                                gATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTacca                                                                                                      >  2:4354656/1‑150 (MQ=255)
                                                       gagaTGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTg                                                                                               >  1:4454329/1‑150 (MQ=255)
                                                       gagaTGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTg                                                                                               >  1:272749/1‑150 (MQ=255)
                                                       gagaTGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTg                                                                                               >  1:2957052/1‑150 (MQ=255)
                                                          aTGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCAt                                                                                            <  2:1625269/150‑1 (MQ=255)
                                                          aTGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCAt                                                                                            <  2:94828/150‑1 (MQ=255)
                                                          aTGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCAt                                                                                            <  2:3479254/150‑1 (MQ=255)
                                                              aTTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGtt                                                                                        <  2:3275846/150‑1 (MQ=255)
                                                              aTTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGtt                                                                                        >  1:2127572/1‑150 (MQ=255)
                                                                tGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTc                                                                                      <  1:2049936/150‑1 (MQ=255)
                                                                tGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTc                                                                                      <  2:5001632/150‑1 (MQ=255)
                                                                tGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTc                                                                                      <  2:3160302/150‑1 (MQ=255)
                                                                 gCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCg                                                                                     >  2:3930622/1‑150 (MQ=255)
                                                                 gCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCg                                                                                     >  1:4120634/1‑150 (MQ=255)
                                                                    aGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCgaga                                                                                  >  1:620957/1‑150 (MQ=255)
                                                                    aGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCgaga                                                                                  >  1:5131267/1‑150 (MQ=255)
                                                                        gATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATtgt                                                                              >  1:368645/1‑150 (MQ=255)
                                                                         aTGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATtgtg                                                                             >  1:4873296/1‑150 (MQ=255)
                                                                          tGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGa                                                                            <  2:1484417/150‑1 (MQ=255)
                                                                           ggTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTCCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGaa                                                                           <  1:1071407/150‑1 (MQ=255)
                                                                              gATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAatat                                                                        <  1:1369185/150‑1 (MQ=255)
                                                                               atatTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATAtt                                                                       >  1:984582/1‑150 (MQ=255)
                                                                                 atTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGa                                                                     >  1:2090201/1‑150 (MQ=255)
                                                                                 atTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGa                                                                     >  2:2588161/1‑150 (MQ=255)
                                                                                 atTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATGTTCTGTACCAGCATGAGTCAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATCTTGa                                                                     >  2:4467831/1‑150 (MQ=255)
                                                                                   tGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAAt                                                                   <  1:4218225/150‑1 (MQ=255)
                                                                                   tGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAAt                                                                   <  2:1708665/150‑1 (MQ=255)
                                                                                   tGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAAt                                                                   <  1:3722969/150‑1 (MQ=255)
                                                                                    gCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATg                                                                  >  1:1053759/1‑150 (MQ=255)
                                                                                      ttCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGcc                                                                <  2:2282921/150‑1 (MQ=255)
                                                                                         gCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGa                                                             >  1:801690/1‑150 (MQ=255)
                                                                                          ccATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGaa                                                            <  2:3560222/150‑1 (MQ=255)
                                                                                           cATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAAc                                                           >  2:1185579/1‑150 (MQ=255)
                                                                                           cATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAAc                                                           >  2:490337/1‑150 (MQ=255)
                                                                                               gcgcTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTAGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCa                                                       >  1:4431494/1‑150 (MQ=255)
                                                                                                cgcTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCAt                                                      >  1:3969706/1‑150 (MQ=255)
                                                                                                cgcTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCAt                                                      <  1:3525207/150‑1 (MQ=255)
                                                                                                  cTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATaa                                                    >  1:1524811/1‑150 (MQ=255)
                                                                                                  cTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATaa                                                    >  1:317358/1‑150 (MQ=255)
                                                                                                   tGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAAc                                                   <  2:801690/150‑1 (MQ=255)
                                                                                                    gACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAAcc                                                  >  2:231807/1‑150 (MQ=255)
                                                                                                       cAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAAccgcc                                               <  2:1608517/150‑1 (MQ=255)
                                                                                                           gcCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCgcgc                                           >  1:4736784/1‑150 (MQ=255)
                                                                                                            ccAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCt                                          >  1:1105357/1‑150 (MQ=255)
                                                                                                               aCAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCtgtg                                       <  1:781645/150‑1 (MQ=255)
                                                                                                               aCAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGATTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCtgtg                                       <  2:64615/150‑1 (MQ=255)
                                                                                                                   aTTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCgtt                                   >  2:312191/1‑150 (MQ=255)
                                                                                                                   aTTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCgtt                                   <  1:2473139/150‑1 (MQ=255)
                                                                                                                   aTTTTCTGTACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCgtt                                   <  1:1139207/150‑1 (MQ=255)
                                                                                                                           tACCAGCATGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCa                           <  2:3356019/150‑1 (MQ=255)
                                                                                                                                   tGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCttt                   <  2:1978854/150‑1 (MQ=255)
                                                                                                                                   tGAGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGACTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCttt                   <  2:1546470/150‑1 (MQ=255)
                                                                                                                                     aGTAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGt                 <  1:5275255/150‑1 (MQ=255)
                                                                                                                                      gTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTc                <  2:520595/150‑1 (MQ=255)
                                                                                                                                       tAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCg               <  1:1734823/150‑1 (MQ=255)
                                                                                                                                       tAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCg               >  2:3874536/1‑150 (MQ=255)
                                                                                                                                       tAAATGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTTCGTTGTTCAGCATCTGCTTTGTCg               <  2:2877455/150‑1 (MQ=255)
                                                                                                                                        aaaTGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGa              <  1:3349652/150‑1 (MQ=255)
                                                                                                                                         aaTGTGCCAGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGAt             >  1:2659987/1‑150 (MQ=255)
                                                                                                                                           tgtgCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGACTATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATAc           >  2:1124329/1‑150 (MQ=255)
                                                                                                                                              gCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCAc        <  1:3584959/150‑1 (MQ=255)
                                                                                                                                                cAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGt      >  1:930308/1‑150 (MQ=255)
                                                                                                                                                cAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGt      >  2:1133116/1‑150 (MQ=255)
                                                                                                                                                 aGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTa     >  2:720619/1‑150 (MQ=255)
                                                                                                                                                 aGGAGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTa     <  2:5318063/150‑1 (MQ=255)
                                                                                                                                                    cGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACAt  <  2:4560152/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCGTGTGATGGTCGAGATGAATTGCCAGTGGATGGTGATATTGCTTCGCCATCGCGCTGACCAGCGCCAACAGATTTTCTGTACCAGCATGAGTAAATGTGCCAGGCGTTCCGGCGATGATGACCGGCGCATGCAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACAT  >  CP009273/2170386‑2170683

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: