Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,001,919 C→A 21.5% H351N (CAT→AAT)  ygeX → 2,3‑diaminopropionate ammonia lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,001,9190CA21.5% 897.9 / 194.6 345H351N (CAT→AAT) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (118/152);  new base A (42/33);  total (160/185)
Fisher's exact test for biased strand distribution p-value = 6.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 347 total aligned reads displayed.

TCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCA  >  CP009273/3001770‑3002068
                                                                                                                                                     |                                                                                                                                                     
 ccTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCa                                                                                                                                                      <  1:1506865/150‑1 (MQ=255)
  cTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCAt                                                                                                                                                     >  1:1077742/1‑150 (MQ=255)
    aCGTTACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGGATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGATTCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTa                                                                                                                                                   <  2:5317367/150‑1 (MQ=255)
     cGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTAt                                                                                                                                                  <  1:3405255/150‑1 (MQ=255)
      gTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATc                                                                                                                                                 >  1:1139711/1‑150 (MQ=255)
         aCTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCAcc                                                                                                                                              >  1:2747660/1‑150 (MQ=255)
         aCTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCAcc                                                                                                                                              <  2:118713/150‑1 (MQ=255)
         aCTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCAcc                                                                                                                                              <  1:2761597/150‑1 (MQ=255)
               ccACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAAc                                                                                                                                        <  2:4858799/150‑1 (MQ=255)
               ccACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAAc                                                                                                                                        <  1:443567/150‑1 (MQ=255)
                cACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACg                                                                                                                                       <  2:3101966/150‑1 (MQ=255)
                  cccAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTc                                                                                                                                     <  2:1574221/150‑1 (MQ=255)
                     aaTTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCaaa                                                                                                                                  >  1:1287930/1‑150 (MQ=255)
                     aaTTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCaaa                                                                                                                                  <  2:1400728/150‑1 (MQ=255)
                     aaTTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCaaa                                                                                                                                  >  2:4301859/1‑150 (MQ=255)
                      aTTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAg                                                                                                                                 >  2:227962/1‑150 (MQ=255)
                      aTTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGAATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAg                                                                                                                                 >  2:3167947/1‑150 (MQ=255)
                         cATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCt                                                                                                                              >  2:1245273/1‑150 (MQ=255)
                         cATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCt                                                                                                                              <  2:2441658/150‑1 (MQ=255)
                         cATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCt                                                                                                                              <  1:164789/150‑1 (MQ=255)
                          aTCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTg                                                                                                                             <  1:4960062/150‑1 (MQ=255)
                          aTCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTg                                                                                                                             >  2:4455050/1‑150 (MQ=255)
                            ctcCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGAt                                                                                                                           <  1:1349668/150‑1 (MQ=255)
                             tcCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATg                                                                                                                          <  2:343180/150‑1 (MQ=255)
                             tcCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATg                                                                                                                          <  2:812024/150‑1 (MQ=255)
                              ccTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATgg                                                                                                                         >  1:4219718/1‑150 (MQ=255)
                                   cAGGACAGCGTTGCCGCATTAGGTATGCGAGTGCTGGGTAATCCGTACGTCAACGACCCGCGCATCATCTCCTGTGTATCCGGCGCTGTCGGTTTTGGCGTTCTCGCAGCGGTTCATTATCACCCGCTTCGTCTATGCCTGATGGaatat                                                                                                                    <  2:4953117/150‑4 (MQ=255)
                                     ggACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACt                                                                                                                  <  2:3572418/150‑1 (MQ=255)
                                     ggACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACt                                                                                                                  <  1:2026048/150‑1 (MQ=255)
                                      gACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTg                                                                                                                 >  1:4589930/1‑150 (MQ=255)
                                       aCAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTgg                                                                                                                >  2:2559705/1‑150 (MQ=255)
                                           cGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCt                                                                                                            <  1:970369/150‑1 (MQ=255)
                                            gTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTg                                                                                                           >  2:947899/1‑150 (MQ=255)
                                                 cGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACaa                                                                                                      >  2:3449654/1‑150 (MQ=255)
                                                  gCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACaaa                                                                                                     >  1:4126636/1‑150 (MQ=255)
                                                  gCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACaaa                                                                                                     >  1:2921446/1‑150 (MQ=255)
                                                   cATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAg                                                                                                    >  2:5178419/1‑150 (MQ=255)
                                                    aTTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGa                                                                                                   <  2:5547048/150‑1 (MQ=255)
                                                      tAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATg                                                                                                 <  2:3875261/150‑1 (MQ=255)
                                                          tATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCgt                                                                                             >  2:1489514/1‑150 (MQ=255)
                                                           aTGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCgtg                                                                                            <  1:1929442/150‑1 (MQ=255)
                                                           aTGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCgtg                                                                                            >  1:245872/1‑150 (MQ=255)
                                                           aTGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGGTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCgtg                                                                                            <  2:2742612/150‑1 (MQ=255)
                                                             gcgcGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCgtggt                                                                                          <  2:1039757/150‑1 (MQ=255)
                                                             gcgcGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCgtggt                                                                                          <  2:2879494/150‑1 (MQ=255)
                                                              cgcgTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCgtggtg                                                                                         >  1:589855/1‑150 (MQ=255)
                                                                  tGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTgg                                                                                     >  1:3263313/1‑150 (MQ=255)
                                                                     tGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTa                                                                                  >  1:3954568/1‑150 (MQ=255)
                                                                      gggTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTAt                                                                                 <  2:2778374/150‑1 (MQ=255)
                                                                       ggTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTCTCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATc                                                                                >  2:1533928/1‑150 (MQ=255)
                                                                       ggTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATc                                                                                >  2:565710/1‑150 (MQ=255)
                                                                        gTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCa                                                                               <  2:92113/150‑1 (MQ=255)
                                                                        gTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCa                                                                               >  1:5442048/1‑150 (MQ=255)
                                                                         tAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAg                                                                              >  2:3322323/1‑150 (MQ=255)
                                                                         tAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAg                                                                              >  2:536628/1‑150 (MQ=255)
                                                                            tCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCAc                                                                           <  1:1109834/150‑1 (MQ=255)
                                                                             ccGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACt                                                                          <  1:2607999/150‑1 (MQ=255)
                                                                              cGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTg                                                                         >  2:1807347/1‑150 (MQ=255)
                                                                                  cGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAgg                                                                     <  2:3900177/150‑1 (MQ=255)
                                                                                   ggCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGTTGCCGTGGTGCTGGTTATCAGCACTGAAGGc                                                                    >  1:2668304/1‑150 (MQ=255)
                                                                                      aaCGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGac                                                                 >  1:815373/1‑150 (MQ=255)
                                                                                       aCGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGaca                                                                >  1:50397/1‑150 (MQ=255)
                                                                                        cGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGacac                                                               <  1:2158344/150‑1 (MQ=255)
                                                                                         gACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACAcc                                                              <  1:3754299/150‑1 (MQ=255)
                                                                                          aCCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCg                                                             >  2:4373057/1‑150 (MQ=255)
                                                                                           cccGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGGTAATTATCACCCGCACCGGCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGa                                                            >  2:3036745/1‑150 (MQ=255)
                                                                                            ccGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGAc                                                           <  1:4026233/150‑1 (MQ=255)
                                                                                            ccGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGAc                                                           >  2:4523257/1‑150 (MQ=255)
                                                                                              gcgcATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGt                                                         >  1:4467901/1‑150 (MQ=255)
                                                                                              gcgcATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGt                                                         >  2:4659205/1‑150 (MQ=255)
                                                                                                  atcatcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAg                                                     <  2:4812232/150‑1 (MQ=255)
                                                                                                  atcatcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAg                                                     >  1:2218694/1‑150 (MQ=255)
                                                                                                    catcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCa                                                   <  1:1245277/150‑1 (MQ=255)
                                                                                                       ctcCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTCTCAGCACTGAAGGCGACACCGACGTGAAGCACTa                                                >  2:1193800/1‑150 (MQ=255)
                                                                                                        tcCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTAc                                               >  1:4823853/1‑150 (MQ=255)
                                                                                                          cGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACcg                                             <  2:4416571/150‑1 (MQ=255)
                                                                                                          cGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACcg                                             <  1:1752764/150‑1 (MQ=255)
                                                                                                           ggTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACcgc                                            <  1:5207217/150‑1 (MQ=255)
                                                                                                            gTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACcgcg                                           <  2:1564860/150‑1 (MQ=255)
                                                                                                            gTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACcgcg                                           >  1:4730860/1‑150 (MQ=255)
                                                                                                               aaTCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAg                                        <  1:4283710/150‑1 (MQ=255)
                                                                                                                aTCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGt                                       <  1:670629/150‑1 (MQ=255)
                                                                                                                aTCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGt                                       <  1:4877406/150‑1 (MQ=255)
                                                                                                                 tCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGtt                                      <  1:5295055/150‑1 (MQ=255)
                                                                                                                  ccGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTg                                     <  2:4124836/150‑1 (MQ=255)
                                                                                                                    ggCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTc                                   <  1:2976790/150‑1 (MQ=255)
                                                                                                                       gcTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTgg                                >  1:2430070/1‑150 (MQ=255)
                                                                                                                       gcTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTgg                                >  2:5286172/1‑150 (MQ=255)
                                                                                                                         tGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGa                              >  2:685549/1‑150 (MQ=255)
                                                                                                                           tCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAg                            <  1:4699780/150‑1 (MQ=255)
                                                                                                                             ggTTTGGGCGTTCTCGCAGTGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGc                          <  2:4865212/150‑1 (MQ=255)
                                                                                                                                   ggCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAacac                    >  2:1895941/1‑150 (MQ=255)
                                                                                                                                    gCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACg                   <  2:5120212/150‑1 (MQ=255)
                                                                                                                                      gTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCt                 >  2:371759/1‑150 (MQ=255)
                                                                                                                                       ttCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTg                <  2:589855/150‑1 (MQ=255)
                                                                                                                                         ctcGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTa              <  1:1546426/150‑1 (MQ=255)
                                                                                                                                          tcGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAg             <  2:1307561/150‑1 (MQ=255)
                                                                                                                                              aGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCAcc         <  1:2380551/150‑1 (MQ=255)
                                                                                                                                                   ttCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAAt    >  1:2068918/1‑150 (MQ=255)
                                                                                                                                                     cATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCa  <  2:221474/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCA  >  CP009273/3001770‑3002068

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: