Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,933,029 A→T 6.8% H375L (CAC→CTC)  pykA → pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,933,0290AT 6.8% 642.7 / 29.0 206H375L (CAC→CTC) pykApyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base A (117/75);  new base T (5/9);  total (122/84)
Fisher's exact test for biased strand distribution p-value = 8.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 217 total aligned reads displayed.

ATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTAC  >  CP009273/1932882‑1933176
                                                                                                                                                   |                                                                                                                                                   
aTCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAacca                                                                                                                                                     <  1:2653218/148‑1 (MQ=255)
aTCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAacca                                                                                                                                                     >  2:1037261/1‑148 (MQ=255)
 tCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAaccac                                                                                                                                                    <  1:3519655/148‑1 (MQ=255)
     tCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCGGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGa                                                                                                                                                >  1:1616978/1‑148 (MQ=255)
      cAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGaa                                                                                                                                               <  1:3073787/148‑1 (MQ=255)
         aaaCCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAgg                                                                                                                                            <  1:1028917/148‑1 (MQ=255)
           aCCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCg                                                                                                                                          <  1:869401/148‑1 (MQ=255)
           aCCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCg                                                                                                                                          >  2:450180/1‑148 (MQ=255)
            ccGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGt                                                                                                                                         <  1:2713205/148‑1 (MQ=255)
             cGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGtt                                                                                                                                        <  2:3500810/148‑1 (MQ=255)
             cGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGtt                                                                                                                                        >  2:1618005/1‑148 (MQ=255)
              gTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTa                                                                                                                                       <  1:1720164/148‑1 (MQ=255)
                 gCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACgg                                                                                                                                    <  1:2060461/148‑1 (MQ=255)
                    gCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGa                                                                                                                                 <  2:513141/148‑1 (MQ=255)
                      cATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGatc                                                                                                                               <  1:1646422/148‑1 (MQ=255)
                      cATGGCGCGCGATTCCCTGGGTGAGGAAAAAATCCCGAGCATAAACGTTTCTAAAAACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCAATTTCAGCAATGTACGCAGCTAACAAAATGAAAGGCGTTACGGCGagc                                                                                                                               >  2:332833/1‑146 (MQ=255)
                               cgTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATg                                                                                                                      <  2:1326381/148‑1 (MQ=255)
                                gTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGa                                                                                                                     <  2:3425709/148‑1 (MQ=255)
                                     ccTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGaa                                                                                                                <  1:246058/148‑1 (MQ=255)
                                              ggAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGt                                                                                                       <  2:1572977/148‑1 (MQ=255)
                                               gAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTa                                                                                                      <  2:304663/148‑1 (MQ=255)
                                               gAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTa                                                                                                      >  1:709543/1‑148 (MQ=255)
                                                   aaaTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACcgc                                                                                                  >  1:2412209/1‑148 (MQ=255)
                                                    aaTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACcgcg                                                                                                 <  1:2188137/148‑1 (MQ=255)
                                                     aTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCgcgc                                                                                                >  2:726223/1‑148 (MQ=255)
                                                     aTCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCgcgc                                                                                                >  2:19233/1‑148 (MQ=255)
                                                      tCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCt                                                                                               >  1:2318612/1‑148 (MQ=255)
                                                      tCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCt                                                                                               >  2:2469894/1‑148 (MQ=255)
                                                      tCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCt                                                                                               >  2:1914552/1‑148 (MQ=255)
                                                      tCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCt                                                                                               >  2:119218/1‑148 (MQ=255)
                                                       cccGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCtg                                                                                              <  2:909063/148‑1 (MQ=255)
                                                       cccGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCtg                                                                                              >  1:2877240/1‑148 (MQ=255)
                                                       cccGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCtg                                                                                              <  1:3425956/148‑1 (MQ=255)
                                                       cccGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCtg                                                                                              <  1:2763396/148‑1 (MQ=255)
                                                       cccGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCtg                                                                                              >  2:629959/1‑148 (MQ=255)
                                                               tCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTc                                                                                      >  1:2562046/1‑148 (MQ=255)
                                                                cAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTca                                                                                     >  2:2638160/1‑147 (MQ=255)
                                                                 aaCGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTccc                                                                                    >  2:608967/1‑148 (MQ=255)
                                                                   cGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGt                                                                                  >  1:2486837/1‑148 (MQ=255)
                                                                    gTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTa                                                                                 >  1:2206063/1‑148 (MQ=255)
                                                                    gTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTa                                                                                 >  2:3490507/1‑148 (MQ=255)
                                                                      ttCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATc                                                                               >  1:587583/1‑148 (MQ=255)
                                                                      ttCTAAACACCGTCTGGACGTTCAGTTCGAAAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATc                                                                               <  1:3007755/148‑1 (MQ=255)
                                                                          aaaCACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGct                                                                           <  2:1867099/148‑1 (MQ=255)
                                                                           aaCACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGctc                                                                          >  2:2833254/1‑148 (MQ=255)
                                                                           aaCACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGctc                                                                          >  2:2699138/1‑148 (MQ=255)
                                                                             cacCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTg                                                                        <  1:1756248/148‑1 (MQ=255)
                                                                             cacCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTg                                                                        <  1:534624/148‑1 (MQ=255)
                                                                             cacCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTg                                                                        <  2:2252715/148‑1 (MQ=255)
                                                                                cGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTc                                                                     <  2:1077101/148‑1 (MQ=255)
                                                                                cGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTc                                                                     >  1:818406/1‑148 (MQ=255)
                                                                                cGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTc                                                                     >  1:2554288/1‑148 (MQ=255)
                                                                                   cTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGc                                                                  >  2:3359349/1‑148 (MQ=255)
                                                                                      gACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCaa                                                               >  2:2395439/1‑148 (MQ=255)
                                                                                      gACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCaa                                                               >  1:1206042/1‑148 (MQ=255)
                                                                                       aCGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAAt                                                              >  2:1606399/1‑148 (MQ=255)
                                                                                        cGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAAtt                                                             <  2:1456387/148‑1 (MQ=255)
                                                                                        cGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAAtt                                                             >  2:2301566/1‑148 (MQ=255)
                                                                                           tCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTc                                                          <  2:2499440/148‑1 (MQ=255)
                                                                                             aGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGc                                                        >  1:622227/1‑148 (MQ=255)
                                                                                               ttCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCa                                                      <  1:324611/148‑1 (MQ=255)
                                                                                               ttCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCa                                                      >  2:300699/1‑148 (MQ=255)
                                                                                               ttCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCa                                                      >  1:2605899/1‑148 (MQ=255)
                                                                                               ttCGACAATGTGGAAGAAGCTATTGCCAAGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACTGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCa                                                      >  2:1799530/1‑148 (MQ=255)
                                                                                                tCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCAt                                                     <  2:1557638/148‑1 (MQ=255)
                                                                                                tCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCAt                                                     >  1:1963167/1‑148 (MQ=255)
                                                                                                  gACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGt                                                   <  1:2947623/148‑1 (MQ=255)
                                                                                                    cAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTcg                                                 >  2:1718753/1‑148 (MQ=255)
                                                                                                       tgtgGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCgcgc                                              <  1:1654135/148‑1 (MQ=255)
                                                                                                       tgtgGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCgcgc                                              <  1:1069740/148‑1 (MQ=255)
                                                                                                       tgtgGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCgcgc                                              >  1:3366897/1‑148 (MQ=255)
                                                                                                        gtgGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGcc                                             <  1:2087164/148‑1 (MQ=255)
                                                                                                          ggAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCAt                                           >  1:1076179/1‑148 (MQ=255)
                                                                                                            aagaagCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGa                                         <  2:549108/148‑1 (MQ=255)
                                                                                                               aagCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACg                                      <  1:2130380/148‑1 (MQ=255)
                                                                                                                agCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGt                                     <  2:1672512/148‑1 (MQ=255)
                                                                                                                agCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGt                                     <  2:2671456/148‑1 (MQ=255)
                                                                                                                 gCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTa                                    >  1:3019041/1‑148 (MQ=255)
                                                                                                                    aTTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGc                                 >  1:2154594/1‑148 (MQ=255)
                                                                                                                         cATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAAc                            >  2:944802/1‑148 (MQ=255)
                                                                                                                         cATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAAc                            <  2:2323163/148‑1 (MQ=255)
                                                                                                                          aTGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAAcc                           >  1:2982910/1‑148 (MQ=255)
                                                                                                                          aTGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAAcc                           >  2:479719/1‑148 (MQ=255)
                                                                                                                          aTGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAAcc                           >  1:2777863/1‑148 (MQ=255)
                                                                                                                                   aTGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGctc                  <  1:2469573/148‑1 (MQ=255)
                                                                                                                                   aTGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGctc                  >  1:2182519/1‑148 (MQ=255)
                                                                                                                                    tGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGctct                 >  1:2391002/1‑148 (MQ=255)
                                                                                                                                    tGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGctct                 <  1:450180/148‑1 (MQ=255)
                                                                                                                                    tGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGctct                 >  1:1305061/1‑148 (MQ=255)
                                                                                                                                    tGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGctct                 <  1:1379703/148‑1 (MQ=255)
                                                                                                                                      tACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGctctct               >  1:3098645/1‑148 (MQ=255)
                                                                                                                                       aCGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTa              >  1:1797962/1‑148 (MQ=255)
                                                                                                                                          cAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCg           <  1:1555744/148‑1 (MQ=255)
                                                                                                                                            gCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTg         >  2:800054/1‑148 (MQ=255)
                                                                                                                                            gCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTg         <  1:2817757/148‑1 (MQ=255)
                                                                                                                                             cTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTgg        >  2:596299/1‑148 (MQ=255)
                                                                                                                                               aaCCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCg      <  2:2000533/148‑1 (MQ=255)
                                                                                                                                                   accTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTAc  >  2:271919/1‑148 (MQ=255)
                                                                                                                                                   accTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTAc  >  2:1466202/1‑148 (MQ=255)
                                                                                                                                                   accTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTAc  >  2:2316053/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTAC  >  CP009273/1932882‑1933176

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: