Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,948,854 G→T 10.5% R55L (CGC→CTC)  ilvA → L‑threonine dehydratase, biosynthetic; also known as threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,948,8540GT10.5% 1077.2 / 64.7 294R55L (CGC→CTC) ilvAL‑threonine dehydratase, biosynthetic; also known as threonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (122/141);  new base T (19/12);  total (141/153)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 311 total aligned reads displayed.

ACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGG  >  CP009273/3948705‑3949003
                                                                                                                                                     |                                                                                                                                                     
aCCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCg                                                                                                                                                       >  1:3318329/1‑150 (MQ=255)
    cTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCgccag                                                                                                                                                   <  2:197426/150‑1 (MQ=255)
       tCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCgccagcca                                                                                                                                                <  1:4438975/150‑1 (MQ=255)
        ccGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGccagccag                                                                                                                                               >  2:572290/1‑150 (MQ=255)
        ccGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGccagccag                                                                                                                                               <  2:1673159/150‑1 (MQ=255)
        ccGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGccagccag                                                                                                                                               >  1:541938/1‑150 (MQ=255)
         cGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGt                                                                                                                                              <  2:258773/150‑1 (MQ=255)
              cTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGcaca                                                                                                                                         >  1:37473/1‑150 (MQ=255)
              cTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGcaca                                                                                                                                         <  2:4076075/150‑1 (MQ=255)
                 cGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCt                                                                                                                                      >  2:2393803/1‑150 (MQ=255)
                  ggAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCtt                                                                                                                                     >  1:2773696/1‑150 (MQ=255)
                  ggAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGACATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCtt                                                                                                                                     >  1:2711895/1‑150 (MQ=255)
                   gAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCttt                                                                                                                                    >  1:1038281/1‑150 (MQ=255)
                    aaGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTa                                                                                                                                   <  1:469543/150‑1 (MQ=255)
                         gCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTg                                                                                                                              <  1:1672823/150‑1 (MQ=255)
                          ccGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTgc                                                                                                                             >  1:1047026/1‑150 (MQ=255)
                           cGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTgcg                                                                                                                            >  2:1747650/1‑150 (MQ=255)
                           cGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTgcg                                                                                                                            >  1:1935689/1‑150 (MQ=255)
                              atatTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCgcgg                                                                                                                         >  1:581468/1‑150 (MQ=255)
                              atatTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCgcgg                                                                                                                         >  1:4113386/1‑150 (MQ=255)
                               tatTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCgcggc                                                                                                                        <  1:1811773/150‑1 (MQ=255)
                                atTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCgcggcg                                                                                                                       >  1:3997214/1‑150 (MQ=255)
                                atTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCgcggcg                                                                                                                       >  2:77090/1‑150 (MQ=255)
                                    aaGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATa                                                                                                                   >  2:977313/1‑150 (MQ=255)
                                    aaGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATa                                                                                                                   >  1:1552680/1‑150 (MQ=255)
                                        gCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGcc                                                                                                               >  2:1916574/1‑150 (MQ=255)
                                         cAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCa                                                                                                              <  2:3881679/150‑1 (MQ=255)
                                          aGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCat                                                                                                             <  1:2880461/150‑1 (MQ=255)
                                            tgctgcGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCatga                                                                                                           <  2:2467727/150‑1 (MQ=255)
                                              ctgcGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCatgatg                                                                                                         >  2:3605744/1‑150 (MQ=255)
                                               tgcGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATgg                                                                                                        <  2:2224343/150‑1 (MQ=255)
                                                gcgcgcgcCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGc                                                                                                       <  2:1576272/150‑1 (MQ=255)
                                                     cgcCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGcc                                                                                                  <  1:1988810/150‑1 (MQ=255)
                                                         ggTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGAc                                                                                              >  1:1165475/1‑150 (MQ=255)
                                                            ttACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGga                                                                                           <  1:894820/150‑1 (MQ=255)
                                                             tACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGgaa                                                                                          >  1:4061583/1‑150 (MQ=255)
                                                             tACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGgaa                                                                                          <  2:857609/150‑1 (MQ=255)
                                                                gAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGgaagaa                                                                                       >  2:1935766/1‑150 (MQ=255)
                                                                gAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGgaagaa                                                                                       >  1:3126300/1‑150 (MQ=255)
                                                                 aGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAAc                                                                                      >  1:2942248/1‑150 (MQ=255)
                                                                 aGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAAc                                                                                      <  2:1327371/150‑1 (MQ=255)
                                                                      gcgcAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGaaa                                                                                 <  1:2413214/150‑1 (MQ=255)
                                                                        gcAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAgc                                                                               >  1:3981090/1‑150 (MQ=255)
                                                                            gTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCAc                                                                           <  2:4538036/150‑1 (MQ=255)
                                                                            gTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCAc                                                                           >  2:1469356/1‑150 (MQ=255)
                                                                             ttACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACg                                                                          >  1:278396/1‑150 (MQ=255)
                                                                             ttACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACg                                                                          >  1:349799/1‑150 (MQ=255)
                                                                                 gccgcTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGt                                                                      >  2:1824786/1‑150 (MQ=255)
                                                                                  ccgcTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTg                                                                     <  2:3499444/150‑1 (MQ=255)
                                                                                   cgcTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGa                                                                    >  2:1138194/1‑150 (MQ=255)
                                                                                      tACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCa                                                                 <  2:3318274/150‑1 (MQ=255)
                                                                                      tACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCa                                                                 >  2:122849/1‑150 (MQ=255)
                                                                                         aaaaaaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTg                                                              <  2:3790651/150‑1 (MQ=255)
                                                                                          aaaaaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGc                                                             >  2:3375469/1‑150 (MQ=255)
                                                                                          aaaaaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGc                                                             <  2:3116951/150‑1 (MQ=255)
                                                                                          aaaaaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGc                                                             <  2:2608327/150‑1 (MQ=255)
                                                                                            aaaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCtt                                                           <  2:1988772/150‑1 (MQ=255)
                                                                                             aaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTc                                                          >  2:1526585/1‑150 (MQ=255)
                                                                                             aaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTc                                                          >  1:3456874/1‑150 (MQ=255)
                                                                                             aaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTc                                                          >  2:908007/1‑150 (MQ=255)
                                                                                             aaTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTc                                                          >  1:3596285/1‑150 (MQ=255)
                                                                                              aTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCt                                                         >  2:2319541/1‑150 (MQ=255)
                                                                                              aTGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCt                                                         <  2:4054137/150‑1 (MQ=255)
                                                                                                ggAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGc                                                       <  1:2336508/150‑1 (MQ=255)
                                                                                                ggAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGc                                                       >  1:1949243/1‑150 (MQ=255)
                                                                                                  aaaaaCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCgg                                                     >  1:4327172/1‑150 (MQ=255)
                                                                                                      aCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTaa                                                 >  2:482572/1‑150 (MQ=255)
                                                                                                      aCTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTaa                                                 <  1:4188970/150‑1 (MQ=255)
                                                                                                            gtcgCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCAcgc                                           >  2:3662488/1‑150 (MQ=255)
                                                                                                             tcgCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCAcgcg                                          <  2:664356/150‑1 (MQ=255)
                                                                                                             tcgCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCAcgcg                                          <  1:697130/150‑1 (MQ=255)
                                                                                                                cgTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAg                                       <  2:2949792/150‑1 (MQ=255)
                                                                                                                cgTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAg                                       <  2:4157829/150‑1 (MQ=255)
                                                                                                                cgTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAg                                       >  1:2422454/1‑150 (MQ=255)
                                                                                                                 gTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAgg                                      >  2:11604/1‑150 (MQ=255)
                                                                                                                    ttGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCg                                   >  2:1012548/1‑150 (MQ=255)
                                                                                                                      gATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTc                                 <  1:2155351/150‑1 (MQ=255)
                                                                                                                      gATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTc                                 >  1:3907903/1‑150 (MQ=255)
                                                                                                                      gATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTc                                 >  2:1075712/1‑150 (MQ=255)
                                                                                                                       aTAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTcg                                <  2:4113386/150‑1 (MQ=255)
                                                                                                                          aCGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGt                             <  1:2539668/150‑1 (MQ=255)
                                                                                                                           cGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGtt                            <  1:3072162/150‑1 (MQ=255)
                                                                                                                            gTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGttt                           <  2:2734535/150‑1 (MQ=255)
                                                                                                                            gTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGttt                           >  1:590803/1‑150 (MQ=255)
                                                                                                                            gTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGttt                           >  2:632139/1‑150 (MQ=255)
                                                                                                                             tCATTCTGGTGAAGCGCGAAGATCTCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGtttt                          >  2:2538243/1‑150 (MQ=255)
                                                                                                                               aTTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTttct                        >  1:2101709/1‑150 (MQ=255)
                                                                                                                                ttCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTttctt                       <  2:2773696/150‑1 (MQ=255)
                                                                                                                                   tGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTg                    <  1:3369090/150‑1 (MQ=255)
                                                                                                                                     gTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTgcg                  <  2:3709381/150‑1 (MQ=255)
                                                                                                                                         aGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGt              >  2:426125/1‑150 (MQ=255)
                                                                                                                                         aGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGt              >  2:1286877/1‑150 (MQ=255)
                                                                                                                                             cgAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAgg          >  2:2504649/1‑150 (MQ=255)
                                                                                                                                              gAAGATCTCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGc         >  2:43056/1‑150 (MQ=255)
                                                                                                                                              gAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGc         >  1:1872043/1‑150 (MQ=255)
                                                                                                                                               aaGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCg        <  2:3930004/150‑1 (MQ=255)
                                                                                                                                               aaGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCg        <  1:784645/150‑1 (MQ=255)
                                                                                                                                                aGATCTCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGt       >  2:2065234/1‑150 (MQ=255)
                                                                                                                                                 gATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTg      >  1:84253/1‑150 (MQ=255)
                                                                                                                                                     gccagccaGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAgg  <  1:96782/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCCCTGTCCGGTGCTCCGGAAGGTGCCGAATATTTAAGAGCAGTGCTGCGCGCGCCGGTTTACGAGGCGGCGCAGGTTACGCCGCTACAAAAAATGGAAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTGCTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGG  >  CP009273/3948705‑3949003

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