Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,574,820 C→A 94.6% A609E (GCG→GAG)  tktB → transketolase 2, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,574,8200CA94.6% 565.9 / 21.6 202A609E (GCG→GAG) tktBtransketolase 2, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (8/3);  new base A (105/84);  total (114/88)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 202 total aligned reads displayed.

GAAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGT  >  CP009273/2574675‑2574958
                                                                                                                                                 |                                                                                                                                          
     tACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGaa                                                                                                                                     >  1:1512775/1‑148 (MQ=255)
        ccTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCa                                                                                                                                  <  1:1563782/148‑1 (MQ=255)
         cTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAg                                                                                                                                 >  1:2732187/1‑148 (MQ=255)
          tGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAgg                                                                                                                                <  1:2256665/148‑1 (MQ=255)
          tGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAgg                                                                                                                                >  1:3424041/1‑148 (MQ=255)
            cAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTa                                                                                                                              <  1:765624/148‑1 (MQ=255)
               gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGNCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTg                                                                                                                           >  1:1809931/1‑148 (MQ=255)
               gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGNCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGTGGAAGCAGGTATTg                                                                                                                           >  1:1810622/1‑148 (MQ=255)
               gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTg                                                                                                                           >  2:621017/1‑148 (MQ=255)
                  gCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCg                                                                                                                        <  1:2322977/148‑1 (MQ=255)
                  gCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCg                                                                                                                        >  2:592485/1‑148 (MQ=255)
                    agagAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGAt                                                                                                                      <  1:28392/148‑1 (MQ=255)
                    agagAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGAt                                                                                                                      <  1:1835081/148‑1 (MQ=255)
                      agaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTa                                                                                                                    >  2:814027/1‑148 (MQ=255)
                        aaaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACt                                                                                                                  <  1:1413582/148‑1 (MQ=255)
                         aaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTg                                                                                                                 <  1:976069/148‑1 (MQ=255)
                          aTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTgg                                                                                                                >  1:392599/1‑148 (MQ=255)
                          aTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTgg                                                                                                                >  2:2025289/1‑148 (MQ=255)
                              gCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACa                                                                                                            <  1:691079/148‑1 (MQ=255)
                               cAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGACTATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACaa                                                                                                           >  2:1573322/1‑148 (MQ=255)
                               cAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACaa                                                                                                           <  2:2207415/148‑1 (MQ=255)
                                 ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGt                                                                                                         <  2:1768834/148‑1 (MQ=255)
                                 ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGt                                                                                                         <  2:5316/148‑1 (MQ=255)
                                  gagaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTa                                                                                                        <  2:3487999/148‑1 (MQ=255)
                                    gaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATg                                                                                                      >  2:1234907/1‑148 (MQ=255)
                                     aaGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGt                                                                                                     >  2:2427701/1‑148 (MQ=255)
                                     aaGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGt                                                                                                     >  2:3265812/1‑148 (MQ=255)
                                     aaGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGAGGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGt                                                                                                     >  2:637441/1‑148 (MQ=255)
                                       ggTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg                                                                                                   >  1:1173368/1‑148 (MQ=255)
                                       ggTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATAGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg                                                                                                   >  2:1169962/1‑148 (MQ=255)
                                        gTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTgg                                                                                                  <  1:2023326/148‑1 (MQ=255)
                                        gTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTgg                                                                                                  >  1:3057570/1‑148 (MQ=255)
                                        gTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTgg                                                                                                  >  1:2395887/1‑148 (MQ=255)
                                         tCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGt                                                                                                 >  1:2880940/1‑148 (MQ=255)
                                         tCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGt                                                                                                 >  1:902084/1‑148 (MQ=255)
                                            cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTg                                                                                              <  1:56890/148‑1 (MQ=255)
                                              aTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGACTATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaa                                                                                            >  2:2578542/1‑148 (MQ=255)
                                                gTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAg                                                                                          >  1:2003027/1‑148 (MQ=255)
                                                 tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgg                                                                                         <  2:1928606/148‑1 (MQ=255)
                                                  aCGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAggg                                                                                        <  2:2316512/148‑1 (MQ=255)
                                                  aCGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAggg                                                                                        <  2:3145926/148‑1 (MQ=255)
                                                   cgcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgggg                                                                                       <  1:1275530/148‑1 (MQ=255)
                                                      gTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCaa                                                                                    <  2:3327838/148‑1 (MQ=255)
                                                      gTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCaa                                                                                    <  2:1792555/148‑1 (MQ=255)
                                                      gTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCaa                                                                                    <  2:418629/148‑1 (MQ=255)
                                                      gTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCaa                                                                                    <  2:1792132/148‑1 (MQ=255)
                                                       tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGACTATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAt                                                                                   >  1:2122304/1‑148 (MQ=255)
                                                       tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAt                                                                                   >  1:2160634/1‑148 (MQ=255)
                                                       tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAt                                                                                   >  2:2957399/1‑148 (MQ=255)
                                                        aGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAtt                                                                                  <  2:690382/148‑1 (MQ=255)
                                                        aGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAtt                                                                                  >  1:591463/1‑148 (MQ=255)
                                                         gTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTg                                                                                 >  1:3032012/1‑148 (MQ=255)
                                                          tgtCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGt                                                                                <  1:1838220/146‑1 (MQ=255)
                                                           ttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTCCTGGTACAAGTATGTTGGTCTGATAGGGGCAATTGTc                                                                               >  2:1924656/1‑148 (MQ=255)
                                                             cccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCgg                                                                             >  1:825101/1‑148 (MQ=255)
                                                               cTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGa                                                                           >  1:2062115/1‑148 (MQ=255)
                                                               cTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGa                                                                           >  1:262743/1‑148 (MQ=255)
                                                                tGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGAt                                                                          <  1:1729833/148‑1 (MQ=255)
                                                                  cccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCACTTGTCGGGATGa                                                                        >  2:1880220/1‑148 (MQ=255)
                                                                   ccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGAc                                                                       >  1:1722390/1‑148 (MQ=255)
                                                                      cGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACggg                                                                    >  1:2946789/1‑148 (MQ=255)
                                                                       gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGt                                                                   <  1:3232959/148‑1 (MQ=255)
                                                                        aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTc                                                                               <  1:2657130/135‑1 (MQ=255)
                                                                        aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTc                                                                               >  2:2657130/1‑135 (MQ=255)
                                                                         ccGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTa                                                                 <  2:3424041/148‑1 (MQ=255)
                                                                           gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACg                                                               >  1:1701557/1‑148 (MQ=255)
                                                                           gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACg                                                               <  2:1411055/148‑1 (MQ=255)
                                                                            atatTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACgg                                                              >  1:2480282/1‑148 (MQ=255)
                                                                            atatTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACgg                                                              >  1:1608532/1‑148 (MQ=255)
                                                                             tatTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACggg                                                             <  2:1525900/148‑1 (MQ=255)
                                                                              atTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACgggg                                                            <  1:3535929/148‑1 (MQ=255)
                                                                               ttttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGa                                                           >  1:682112/1‑148 (MQ=255)
                                                                                tttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGaa                                                          >  1:48588/1‑148 (MQ=255)
                                                                                           aGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcggc                                               >  2:154049/1‑148 (MQ=255)
                                                                                            ggATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcggcg                                              >  1:152952/1‑148 (MQ=255)
                                                                                            ggATGAGGAATATCGGGAGTCGGTGTTGCCTTCAAACGTTGCGGCTCGCGTGGATGTGGAATCAGGTATTGACGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTTCTCcggcg                                              >  1:2514888/1‑148 (MQ=39)
                                                                                                gAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATa                                          <  2:199126/148‑1 (MQ=255)
                                                                                                  ggtatatCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCATGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATTACGGGTTACGGGGAATCTGCTCCGGCGGATAAg                                        <  1:2935107/145‑1 (MQ=255)
                                                                                                   gAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGGATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGc                                       >  1:2620058/1‑148 (MQ=255)
                                                                                                     atatCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTg                                     >  1:1158837/1‑148 (MQ=255)
                                                                                                           ggAGTCGTTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTTGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCg                               <  2:2122304/148‑1 (MQ=255)
                                                                                                             aGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGtt                             <  2:392600/148‑1 (MQ=255)
                                                                                                             aGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGtt                             <  2:392548/148‑1 (MQ=255)
                                                                                                              gTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTc                            <  1:621066/148‑1 (MQ=255)
                                                                                                               tCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                           <  2:2801877/148‑1 (MQ=255)
                                                                                                                 ggTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCttt                         >  2:1393443/1‑148 (MQ=255)
                                                                                                                     ttGCCTTCAAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATAGCAGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCt                     >  1:1170587/1‑148 (MQ=255)
                                                                                                                      tGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCtt                    >  2:2836283/1‑148 (MQ=255)
                                                                                                                        ccTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGTATAAGCTGTTCCCGTTCTTTGGCTTTa                  <  2:2999407/148‑1 (MQ=255)
                                                                                                                           tCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccg               >  1:989723/1‑148 (MQ=255)
                                                                                                                            cTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccgc              >  1:1072152/1‑148 (MQ=255)
                                                                                                                             tAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAcagcc             >  2:2081373/1‑148 (MQ=255)
                                                                                                                              aaCGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccgccg            <  1:2609943/148‑1 (MQ=255)
                                                                                                                               aCGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCt                     >  2:2296413/1‑138 (MQ=255)
                                                                                                                               aCGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCga           >  2:1307369/1‑148 (MQ=255)
                                                                                                                                cGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCgag          >  2:2189460/1‑148 (MQ=255)
                                                                                                                                  ttGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGaa        >  1:1013882/1‑148 (MQ=255)
                                                                                                                                     cGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAatat     <  1:3004372/148‑1 (MQ=255)
                                                                                                                                       gCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtg   <  1:2957399/148‑1 (MQ=255)
                                                                                                                                       gCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATtg   <  1:2315469/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
GAAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGT  >  CP009273/2574675‑2574958

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: