Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,524,056 C→A 100% A179A (GCG→GCT zipA ← FtsZ stabilizer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,524,0560CA97.7% 556.6 / 3.8 177A179A (GCG→GCTzipAFtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (93/80);  total (95/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 177 total aligned reads displayed.

TATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCCGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACGGCTGTTCAGGCTGCGG  >  CP009273/2523911‑2524199
                                                                                                                                                 |                                                                                                                                               
tatTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGc                                                                                                                                               >  2:1344867/1‑148 (MQ=255)
   tCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTAc                                                                                                                                            >  1:260490/1‑148 (MQ=255)
    cATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTAca                                                                                                                                           >  1:39369/1‑148 (MQ=255)
      tatCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTAcagg                                                                                                                                         >  2:3421943/1‑148 (MQ=255)
      tatCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTAcagg                                                                                                                                         <  2:331160/148‑1 (MQ=255)
       atCGCCAAAAATGAAGCCCGCTTGTAGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTAcaggc                                                                                                                                        >  1:2243090/1‑148 (MQ=255)
        tCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTAcaggct                                                                                                                                       >  2:3771861/1‑148 (MQ=255)
            cAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTAcaggctcagg                                                                                                                                   <  2:1402635/148‑1 (MQ=255)
               aaaTGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCtg                                                                                                                                >  2:2044799/1‑148 (MQ=255)
                   gAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTgtggt                                                                                                                            >  1:801516/1‑148 (MQ=255)
                    aaGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTgtggtg                                                                                                                           <  2:2655668/148‑1 (MQ=255)
                    aaGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTgtggtg                                                                                                                           <  2:470918/148‑1 (MQ=255)
                     aGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGc                                                                                                                          <  2:3791443/148‑1 (MQ=255)
                      gCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGcc                                                                                                                         >  1:2702136/1‑148 (MQ=255)
                       cccGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGATCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCg                                                                                                                        >  2:769840/1‑148 (MQ=255)
                       cccGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGATCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCg                                                                                                                        >  2:768898/1‑148 (MQ=255)
                           cTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGGGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTAc                                                                                                                    >  2:2206994/1‑148 (MQ=255)
                             tgttgAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACgg                                                                                                                  >  1:3265494/1‑148 (MQ=255)
                              gttgAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACggg                                                                                                                 <  2:631061/148‑1 (MQ=255)
                               ttgAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGt                                                                                                                >  2:1739126/1‑148 (MQ=255)
                                    tGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGc                                                                                                           >  2:2726769/1‑148 (MQ=255)
                                    tGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGc                                                                                                           >  2:3349587/1‑148 (MQ=255)
                                    tGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGc                                                                                                           >  2:3408551/1‑148 (MQ=255)
                                       tGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAgg                                                                                                        <  1:1842962/148‑1 (MQ=255)
                                        gTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGc                                                                                                       >  2:2606186/1‑148 (MQ=255)
                                         ttAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCt                                                                                                      <  1:1835961/148‑1 (MQ=255)
                                             gaagCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGTCGCTACGGGTTCTGCAGGCTGGaa                                                                                                  >  1:2047316/1‑148 (MQ=255)
                                             gaagCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGaa                                                                                                  <  1:1711879/148‑1 (MQ=255)
                                              aagCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGaaa                                                                                                 >  2:3156599/1‑148 (MQ=255)
                                              aagCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGaaa                                                                                                 >  2:2831334/1‑148 (MQ=255)
                                                gCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGc                                                                                               >  1:3258112/1‑148 (MQ=255)
                                                gCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGc                                                                                               <  1:1901107/148‑1 (MQ=255)
                                                 cAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGcc                                                                                              >  1:1451829/1‑148 (MQ=255)
                                                    ttCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCtgt                                                                                           >  1:3173193/1‑148 (MQ=255)
                                                    ttCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCtgt                                                                                           <  1:3272141/148‑1 (MQ=255)
                                                    ttCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCtgt                                                                                           >  2:1139953/1‑148 (MQ=255)
                                                    ttCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCtgt                                                                                           >  2:1139915/1‑148 (MQ=255)
                                                    ttCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCtgt                                                                                           >  1:916801/1‑148 (MQ=255)
                                                    ttCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCtgt                                                                                           >  1:1270415/1‑148 (MQ=255)
                                                      cACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCtgttg                                                                                         <  1:2870604/148‑1 (MQ=255)
                                                            ttAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCgg                                                                                   <  2:1505689/148‑1 (MQ=255)
                                                            ttAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCgg                                                                                   <  1:1147810/148‑1 (MQ=255)
                                                            ttAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCgg                                                                                   >  1:717318/1‑148 (MQ=255)
                                                             tAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGt                                                                                  <  1:19812/148‑1 (MQ=255)
                                                                     tACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGc                                                                          <  2:2575046/148‑1 (MQ=255)
                                                                     tACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGc                                                                          >  1:3533321/1‑148 (MQ=255)
                                                                        cGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGa                                                                       <  1:97683/148‑1 (MQ=255)
                                                                          tGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAAt                                                                     >  1:2406542/1‑148 (MQ=255)
                                                                          tGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAAt                                                                     >  1:2665923/1‑148 (MQ=255)
                                                                           gATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATg                                                                    <  1:3659031/148‑1 (MQ=255)
                                                                           gATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATg                                                                    <  1:3659030/148‑1 (MQ=255)
                                                                              gcgcCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGcac                                                                 <  1:3444072/148‑1 (MQ=255)
                                                                                    cgACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTg                                                           >  2:2661964/1‑148 (MQ=255)
                                                                                      aCGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTgcg                                                         >  1:2289076/1‑148 (MQ=255)
                                                                                       cGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGGTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTgcgg                                                        >  1:2498416/1‑148 (MQ=255)
                                                                                          tCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGgcgc                                                     >  2:2055671/1‑148 (MQ=255)
                                                                                          tCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGgcgc                                                     >  1:497763/1‑148 (MQ=255)
                                                                                           cATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCt                                                    >  1:3135833/1‑148 (MQ=255)
                                                                                           cATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCt                                                    <  1:3771828/148‑1 (MQ=255)
                                                                                            aTGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTg                                                   >  2:3668755/1‑148 (MQ=255)
                                                                                            aTGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTg                                                   <  1:3788160/148‑1 (MQ=255)
                                                                                             tGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTgg                                                  >  1:2260928/1‑148 (MQ=255)
                                                                                               aTAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGgcg                                                <  1:1433581/148‑1 (MQ=255)
                                                                                               aTAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGgcg                                                >  1:562890/1‑148 (MQ=255)
                                                                                                tAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGgcgc                                               >  1:2071452/1‑148 (MQ=255)
                                                                                                tAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGgcgc                                               >  1:174542/1‑148 (MQ=255)
                                                                                                 actCACCGCTTCTTTGCGCTTCGGTTGATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGcgcg                                              <  2:2478100/146‑1 (MQ=255)
                                                                                                 aaTCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGcgcg                                              >  2:3027729/1‑148 (MQ=255)
                                                                                                  aTCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGa                                             <  1:592861/148‑1 (MQ=255)
                                                                                                   tCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGAc                                            >  1:2152087/1‑148 (MQ=255)
                                                                                                   tCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGAc                                            >  2:642321/1‑148 (MQ=255)
                                                                                                    cACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCCGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGAcc                                           <  1:1651243/148‑1 (MQ=255)
                                                                                                    cACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGAcc                                           <  2:801516/148‑1 (MQ=255)
                                                                                                    cACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGAcc                                           <  2:1561878/148‑1 (MQ=255)
                                                                                                    cACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGGTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGAcc                                           <  2:1562555/148‑1 (MQ=255)
                                                                                                      ccGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCtg                                         <  2:239535/148‑1 (MQ=255)
                                                                                                      ccGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCtg                                         <  1:3238316/148‑1 (MQ=255)
                                                                                                       cGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCtgt                                        <  2:2418357/148‑1 (MQ=255)
                                                                                                         cttcttTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCCGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTgg                                      >  1:1734265/1‑148 (MQ=255)
                                                                                                           tcttTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCg                                    >  1:1008182/1‑148 (MQ=255)
                                                                                                             tttGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGaa                                  <  2:2191863/148‑1 (MQ=255)
                                                                                                                 cgcTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTg                              <  1:3419994/148‑1 (MQ=255)
                                                                                                                 cgcTTCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTg                              >  2:243742/1‑148 (MQ=255)
                                                                                                                    ttCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGgct                           >  1:1573611/1‑148 (MQ=255)
                                                                                                                     tCGGTTTATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGgctg                          >  2:522201/1‑148 (MQ=255)
                                                                                                                         tttATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCtgctgc                      >  2:1484519/1‑148 (MQ=255)
                                                                                                                          ttATCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCa                     <  2:3545624/148‑1 (MQ=255)
                                                                                                                            aTCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAAc                   >  1:135184/1‑148 (MQ=255)
                                                                                                                             tCCATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACg                  >  2:1165284/1‑148 (MQ=255)
                                                                                                                              ccATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACgg                 <  2:2996753/148‑1 (MQ=255)
                                                                                                                              ccATAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACgg                 >  1:1053965/1‑148 (MQ=255)
                                                                                                                                aTAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACGGCt               >  2:1273857/1‑148 (MQ=255)
                                                                                                                                 tAACTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACGGCTg              <  2:1793804/148‑1 (MQ=255)
                                                                                                                                   aCTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACGGCTGtt            <  1:3736207/148‑1 (MQ=255)
                                                                                                                                    cTGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACGGCTGTTc           <  2:308402/148‑1 (MQ=255)
                                                                                                                                     tGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACGGCTGTTCa          <  2:3122337/148‑1 (MQ=255)
                                                                                                                                     tGGAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACGGCTGTTCa          <  2:3122497/148‑1 (MQ=255)
                                                                                                                                       gAGCAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACGGCTGTTCAgg        >  2:2997556/1‑148 (MQ=255)
                                                                                                                                          cAGGTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACGGCTGTTCAGGCTg     >  1:311826/1‑148 (MQ=255)
                                                                                                                                             gTTCAGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACGGCTGTTCAGGCTGCgg  >  1:1671568/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TATTCATATCGCCAAAAATGAAGCCCGCTTGTTGAATGCTGTTAAGAAGCAGTTCACCGTTTAGCTCGCTACCGTGATGCGCCGCGACGTTCATGATAATCACCGCTTCTTTGCGCTTCGGTTTATCCATAACTGGAGCAGGTTCCGCTACAGGCTCAGGCTGTGGTGCCGCTACGGGTTCTGCAGGCTGGAAAGCCTGTTGTGCCGGTTGCGGTGCTGAATGCACAGGCTGCGGCGCTGGCGCGACCTGTGGCGAAACTGGCTGCTGCAACGGCTGTTCAGGCTGCGG  >  CP009273/2523911‑2524199

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: