Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 546,782 G→T 6.7% S296S (TCG→TCT ybcF → putative carbonate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273546,7820GT 6.7% 950.7 / 27.7 246S296S (TCG→TCTybcFputative carbonate kinase
Reads supporting (aligned to +/- strand):  ref base G (123/106);  new base T (8/9);  total (131/115)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

TTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATC  >  CP009273/546633‑546931
                                                                                                                                                     |                                                                                                                                                     
  aGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCt                                                                                                                                                     >  1:1509065/1‑150 (MQ=255)
  aGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCt                                                                                                                                                     >  1:1509289/1‑150 (MQ=255)
    cgcCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGt                                                                                                                                                   >  1:1061901/1‑150 (MQ=255)
       cATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGt                                                                                                                                                <  1:1126274/150‑1 (MQ=255)
          ttGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGt                                                                                                                                             <  2:3793670/150‑1 (MQ=255)
               aaaGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCa                                                                                                                                        >  1:533036/1‑150 (MQ=255)
               aaaGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCa                                                                                                                                        >  1:2119004/1‑150 (MQ=255)
                 aGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCAtt                                                                                                                                      <  2:3049868/150‑1 (MQ=255)
                     gATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGAc                                                                                                                                  >  2:1451137/1‑150 (MQ=255)
                       tGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACAt                                                                                                                                <  1:2075203/150‑1 (MQ=255)
                            cGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTtgtg                                                                                                                           >  1:3336216/1‑150 (MQ=255)
                             gATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGc                                                                                                                          >  1:1004824/1‑150 (MQ=255)
                              aTGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGcc                                                                                                                         >  2:570254/1‑150 (MQ=255)
                               tGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCt                                                                                                                        <  2:603740/150‑1 (MQ=255)
                                ggggCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCTCTGTAGTCGTAGGCATTAGACATTTGTGCCTg                                                                                                                       <  2:3751665/150‑1 (MQ=255)
                                ggggCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTg                                                                                                                       >  1:859912/1‑150 (MQ=255)
                                 gggCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGa                                                                                                                      <  2:1018072/150‑1 (MQ=255)
                                  ggCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGAt                                                                                                                     <  2:543758/150‑1 (MQ=255)
                                    ccGAATGTATCGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGc                                                                                                                   <  2:1546322/150‑1 (MQ=255)
                                     cGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCg                                                                                                                  >  2:4064857/1‑150 (MQ=255)
                                       aaTGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGAc                                                                                                                <  1:3141724/150‑1 (MQ=255)
                                         tGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCTCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGc                                                                                                              >  1:1596740/1‑150 (MQ=255)
                                              cggcggTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGAc                                                                                                         >  1:1131471/1‑150 (MQ=255)
                                                  ggTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGt                                                                                                     <  1:1040616/150‑1 (MQ=255)
                                                  ggTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGt                                                                                                     >  2:1089766/1‑150 (MQ=255)
                                                   gTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTc                                                                                                    >  1:2249618/1‑150 (MQ=255)
                                                    tGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCt                                                                                                   >  1:2766209/1‑150 (MQ=255)
                                                       gTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTAt                                                                                                >  2:1760244/1‑150 (MQ=255)
                                                       gTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTAt                                                                                                <  2:796636/150‑1 (MQ=255)
                                                          gTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAg                                                                                             >  1:1106217/1‑150 (MQ=255)
                                                           ttATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAgg                                                                                            <  1:1648520/150‑1 (MQ=255)
                                                            tATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGc                                                                                           >  1:3942306/1‑150 (MQ=255)
                                                             aTGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                                                          <  1:3580061/150‑1 (MQ=255)
                                                               gTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGTAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTa                                                                                        <  2:729986/150‑1 (MQ=255)
                                                               gTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTa                                                                                        <  2:3850049/150‑1 (MQ=255)
                                                                tCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTAc                                                                                       >  2:3004128/1‑150 (MQ=255)
                                                                tCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTAc                                                                                       <  2:1692790/150‑1 (MQ=255)
                                                                 cAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACa                                                                                      >  2:822706/1‑150 (MQ=255)
                                                                 cAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACa                                                                                      <  2:2598228/150‑1 (MQ=255)
                                                                 cAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACa                                                                                      <  1:4154240/150‑1 (MQ=255)
                                                                  aGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCTCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACaa                                                                                     <  1:3816424/150‑1 (MQ=255)
                                                                  aGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACaa                                                                                     >  1:3481408/1‑150 (MQ=255)
                                                                    aaGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAcc                                                                                   >  1:806995/1‑150 (MQ=255)
                                                                    aaGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAcc                                                                                   >  2:2055354/1‑150 (MQ=255)
                                                                       ccGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGt                                                                                <  2:1168756/150‑1 (MQ=255)
                                                                       ccGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGt                                                                                >  2:398304/1‑150 (MQ=255)
                                                                       ccGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGt                                                                                <  2:1299801/150‑1 (MQ=255)
                                                                         gTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGc                                                                              >  1:3420962/1‑150 (MQ=255)
                                                                          tGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGcc                                                                             <  1:601317/150‑1 (MQ=255)
                                                                               aaCCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATc                                                                        >  1:1687310/1‑150 (MQ=255)
                                                                                aCCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATcc                                                                       >  1:3091791/1‑150 (MQ=255)
                                                                                  ccGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCgg                                                                     <  1:2176412/150‑1 (MQ=255)
                                                                                      gTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAAt                                                                 >  1:1851819/1‑150 (MQ=255)
                                                                                        ggATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCTCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTg                                                               <  2:641395/150‑1 (MQ=255)
                                                                                        ggATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTg                                                               >  1:3333146/1‑150 (MQ=255)
                                                                                          aTTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGt                                                             >  1:901489/1‑150 (MQ=255)
                                                                                          aTTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGt                                                             >  2:4169230/1‑150 (MQ=255)
                                                                                           ttGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTg                                                            <  2:2985758/150‑1 (MQ=255)
                                                                                            tGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGc                                                           <  2:3693242/150‑1 (MQ=255)
                                                                                            tGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGc                                                           >  2:362167/1‑150 (MQ=255)
                                                                                               ggCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGcaca                                                        >  1:2500069/1‑150 (MQ=255)
                                                                                                 cGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCTCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAAt                                                      >  2:3170707/1‑150 (MQ=255)
                                                                                                 cGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAAt                                                      >  1:2842379/1‑150 (MQ=255)
                                                                                                 cGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAAt                                                      >  1:3250963/1‑150 (MQ=255)
                                                                                                 cGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAAt                                                      >  2:3786556/1‑150 (MQ=255)
                                                                                                        gcgAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgat                                               <  2:2914325/150‑1 (MQ=255)
                                                                                                          gAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgc                                             <  1:3391731/150‑1 (MQ=255)
                                                                                                            agtGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCTCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgcga                                           <  2:2275654/148‑1 (MQ=255)
                                                                                                             ttGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgcgat                                          <  1:3057260/150‑1 (MQ=255)
                                                                                                               gAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgcgatgc                                        <  2:3483503/150‑1 (MQ=255)
                                                                                                                aaGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCt                                       >  1:1717070/1‑150 (MQ=255)
                                                                                                                aaGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCt                                       >  1:3794841/1‑150 (MQ=255)
                                                                                                                         tGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATTCCTGATGCGATGCTTGACGCATc                              <  2:1620488/150‑1 (MQ=255)
                                                                                                                           gcgggcgAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCtt                            <  1:2698125/150‑1 (MQ=255)
                                                                                                                           gcgggcgAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCtt                            <  1:4015386/150‑1 (MQ=255)
                                                                                                                           gcgggcgAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCtt                            <  2:1626090/150‑1 (MQ=255)
                                                                                                                            cgggcgAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTa                           >  1:3156642/1‑150 (MQ=255)
                                                                                                                             gggcgAAGCGGGGACCTGTATGTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt                          <  2:890193/150‑1 (MQ=255)
                                                                                                                               gcgAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCGTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCa                        <  2:134809/150‑1 (MQ=255)
                                                                                                                                 gAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAgg                      >  2:2773998/1‑150 (MQ=255)
                                                                                                                                 gAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAgg                      >  1:3313550/1‑150 (MQ=255)
                                                                                                                                  aaGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGc                     >  2:3986015/1‑150 (MQ=255)
                                                                                                                                   aGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGcc                    >  1:4051048/1‑150 (MQ=255)
                                                                                                                                     cGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTa                  <  1:3371281/150‑1 (MQ=255)
                                                                                                                                     cGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTa                  <  2:3369541/150‑1 (MQ=255)
                                                                                                                                     cGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTa                  >  1:2855782/1‑150 (MQ=255)
                                                                                                                                      ggggACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTAc                 >  1:545815/1‑150 (MQ=255)
                                                                                                                                       gggACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACa                >  1:190385/1‑150 (MQ=255)
                                                                                                                                        ggACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACaa               >  2:354803/1‑150 (MQ=255)
                                                                                                                                         gACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt              <  2:1962778/150‑1 (MQ=255)
                                                                                                                                         gACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAAt              >  1:3243917/1‑150 (MQ=255)
                                                                                                                                           ccTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATgg            >  2:2212265/1‑150 (MQ=255)
                                                                                                                                                aTTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTAcc       >  2:261892/1‑150 (MQ=255)
                                                                                                                                                aTTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTAcc       >  2:3337327/1‑150 (MQ=255)
                                                                                                                                                 tttCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCg      <  1:605040/150‑1 (MQ=255)
                                                                                                                                                 tttCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCg      <  2:2482163/150‑1 (MQ=255)
                                                                                                                                                  ttCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCgg     <  1:2640409/150‑1 (MQ=255)
                                                                                                                                                   tCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGa    >  2:2139008/1‑150 (MQ=255)
                                                                                                                                                    cGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGAt   <  2:3664710/150‑1 (MQ=255)
                                                                                                                                                     gCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATc  <  2:915724/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATC  >  CP009273/546633‑546931

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: