Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,001,919 C→A 14.4% 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,9190CA14.4% 675.2 / 77.4 215H351N (CAT→AAT) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (75/109);  new base A (12/19);  total (87/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 217 total aligned reads displayed.

CCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCA  >  CP009273/3001771‑3002068
                                                                                                                                                    |                                                                                                                                                     
ccTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCa                                                                                                                                                      <  1:4857938/150‑1 (MQ=255)
 cTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCAt                                                                                                                                                     >  2:3988806/1‑150 (MQ=255)
      tAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCa                                                                                                                                                >  1:1047096/1‑150 (MQ=255)
      tAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCa                                                                                                                                                >  2:2650002/1‑150 (MQ=255)
      tAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCa                                                                                                                                                >  2:2109291/1‑150 (MQ=255)
           gcgcCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGc                                                                                                                                           >  2:1674249/1‑150 (MQ=255)
            cgcCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCa                                                                                                                                          <  2:3435549/150‑1 (MQ=255)
              ccACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAAc                                                                                                                                        >  2:3285114/1‑150 (MQ=255)
              ccACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAAc                                                                                                                                        <  1:3333836/150‑1 (MQ=255)
                aCCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGt                                                                                                                                      <  2:841927/150‑1 (MQ=255)
                  ccAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCa                                                                                                                                    <  1:423124/150‑1 (MQ=255)
                   cAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCaa                                                                                                                                   <  2:44664/150‑1 (MQ=255)
                     aTTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAg                                                                                                                                 >  2:3260434/1‑150 (MQ=255)
                     aTTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAg                                                                                                                                 >  2:654014/1‑150 (MQ=255)
                     aTTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAg                                                                                                                                 <  1:2841931/150‑1 (MQ=255)
                     aTTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAg                                                                                                                                 >  2:4921057/1‑150 (MQ=255)
                        cATCTCCTGCCAGGACAGCGTTGCCGCATTATGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCt                                                                                                                              <  2:1879137/150‑1 (MQ=255)
                         aTCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTg                                                                                                                             >  1:559687/1‑150 (MQ=255)
                          tctcCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGa                                                                                                                            <  2:2945126/150‑1 (MQ=255)
                           ctcCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGAt                                                                                                                           <  2:3847529/150‑1 (MQ=255)
                           ctcCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGAt                                                                                                                           >  1:2829398/1‑150 (MQ=255)
                               tGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGaa                                                                                                                       <  2:513933/150‑1 (MQ=255)
                                gCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGaaa                                                                                                                      <  1:1186313/150‑1 (MQ=255)
                                  cAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGaaaaa                                                                                                                    <  1:56481/150‑1 (MQ=255)
                                      aCAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTgg                                                                                                                >  1:360514/1‑150 (MQ=255)
                                       cAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGgc                                                                                                               <  2:3042730/150‑1 (MQ=255)
                                        aGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGgcg                                                                                                              <  1:4468039/150‑1 (MQ=255)
                                          cGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCt                                                                                                            >  1:482268/1‑150 (MQ=255)
                                          cGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCt                                                                                                            >  1:3594644/1‑150 (MQ=255)
                                           gTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTg                                                                                                           >  2:12917/1‑150 (MQ=255)
                                           gTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTg                                                                                                           <  1:3710929/150‑1 (MQ=255)
                                            ttGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCa                                                                                                                                    >  1:3477154/1‑124 (MQ=255)
                                            ttGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCa                                                                                                                                    <  2:3477154/124‑1 (MQ=255)
                                            ttGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGa                                                                                                          >  1:1782221/1‑150 (MQ=255)
                                            ttGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGa                                                                                                          <  1:923469/150‑1 (MQ=255)
                                              gCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAAc                                                                                                        <  1:3791456/150‑1 (MQ=255)
                                              gCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAAc                                                                                                        <  1:2676879/150‑1 (MQ=255)
                                                cGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACaa                                                                                                      >  2:2472351/1‑150 (MQ=255)
                                                  cATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAg                                                                                                    >  2:47493/1‑150 (MQ=255)
                                                  cATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAg                                                                                                    >  1:2415266/1‑150 (MQ=255)
                                                  cATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAg                                                                                                    <  1:2617665/150‑1 (MQ=255)
                                                    ttAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGAt                                                                                                  <  1:1989362/150‑1 (MQ=255)
                                                     tAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATg                                                                                                 <  1:203813/150‑1 (MQ=255)
                                                     tAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATg                                                                                                 <  2:1658406/150‑1 (MQ=255)
                                                        gTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCg                                                                                              >  2:1545938/1‑150 (MQ=255)
                                                          aTGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCgtg                                                                                            >  2:4970889/1‑150 (MQ=255)
                                                           tGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCgtgg                                                                                           <  1:4641285/150‑1 (MQ=255)
                                                             cgcgTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCgtggtg                                                                                         <  2:2314568/150‑1 (MQ=255)
                                                               cgTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCt                                                                                       <  1:1426834/150‑1 (MQ=255)
                                                               cgTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCt                                                                                       <  1:1426940/150‑1 (MQ=255)
                                                               cgTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCt                                                                                       <  1:5057902/150‑1 (MQ=255)
                                                                     gggTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTAt                                                                                 >  2:1945858/1‑150 (MQ=255)
                                                                       gTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCa                                                                               <  2:4872224/150‑1 (MQ=255)
                                                                       gTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCa                                                                               <  1:5018139/150‑1 (MQ=255)
                                                                         aaTCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGc                                                                             >  2:1788434/1‑150 (MQ=255)
                                                                         aaTCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGc                                                                             >  1:2538559/1‑150 (MQ=255)
                                                                                  ggCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGc                                                                    <  1:1240263/150‑1 (MQ=255)
                                                                                   gCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCg                                                                   <  1:3028084/150‑1 (MQ=255)
                                                                                         aCCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCg                                                             <  2:4703028/150‑1 (MQ=255)
                                                                                           ccGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGAc                                                           <  1:2267548/150‑1 (MQ=255)
                                                                                             gcgcATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGt                                                         <  1:2695832/150‑1 (MQ=255)
                                                                                             gcgcATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGt                                                         <  1:2378210/150‑1 (MQ=255)
                                                                                              cgcATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTg                                                        <  1:1844252/150‑1 (MQ=255)
                                                                                                catcatCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGaa                                                      <  2:667631/150‑1 (MQ=255)
                                                                                                catcatCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGaa                                                      >  1:24550/1‑150 (MQ=255)
                                                                                                 atcatcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAg                                                     <  1:3953777/150‑1 (MQ=255)
                                                                                                     tctcCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACt                                                 >  1:4909268/1‑150 (MQ=255)
                                                                                                       tcCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTAc                                               <  1:1669413/150‑1 (MQ=255)
                                                                                                       tcCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGAGCTGAACAAAGATGCCGTTGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTAc                                               <  2:4771029/150‑1 (MQ=255)
                                                                                                        ccGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTAcc                                              <  2:67596/150‑1 (MQ=255)
                                                                                                         cGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACcg                                             <  2:5037310/150‑1 (MQ=255)
                                                                                                          ggTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACcgc                                            <  1:2079373/150‑1 (MQ=255)
                                                                                                               aTCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGt                                       >  2:3328046/1‑150 (MQ=255)
                                                                                                                  cGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAACGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGt                                    <  1:1008935/150‑1 (MQ=255)
                                                                                                                  cGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGt                                    <  2:628360/150‑1 (MQ=255)
                                                                                                                  cGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGt                                    >  2:1244553/1‑150 (MQ=255)
                                                                                                                   ggCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTc                                   <  2:1620162/150‑1 (MQ=255)
                                                                                                                    gcgcTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCt                                  <  1:2998587/150‑1 (MQ=255)
                                                                                                                       cTGTCGGTTTGGGCGTTCTCGCAGCGGTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTggg                               <  2:2427227/150‑1 (MQ=255)
                                                                                                                        tGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGa                              <  2:3705223/150‑1 (MQ=255)
                                                                                                                          tCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAg                            <  2:1742307/150‑1 (MQ=255)
                                                                                                                                tGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAac                      <  1:1944958/150‑1 (MQ=255)
                                                                                                                                  ggCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAacac                    >  1:1630001/1‑150 (MQ=255)
                                                                                                                                  ggCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAacac                    <  1:3108383/150‑1 (MQ=255)
                                                                                                                                          cgtaGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGc            <  2:5098454/147‑1 (MQ=255)
                                                                                                                                          cGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGc            <  2:2207942/150‑1 (MQ=255)
                                                                                                                                            cAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCAc          <  1:2940812/150‑1 (MQ=255)
                                                                                                                                            cAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCAc          <  1:3388265/150‑1 (MQ=255)
                                                                                                                                             aGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCAcc         <  1:5114647/150‑1 (MQ=255)
                                                                                                                                             aGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCAcc         >  1:4125487/1‑150 (MQ=255)
                                                                                                                                              gCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCt        <  2:2685133/150‑1 (MQ=255)
                                                                                                                                              gCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCt        <  1:2626794/150‑1 (MQ=255)
                                                                                                                                               cGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCtt       >  1:627386/1‑150 (MQ=255)
                                                                                                                                                ggTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTCGCACCTTa      >  2:4070529/1‑150 (MQ=255)
                                                                                                                                                ggTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTa      >  1:752285/1‑150 (MQ=255)
                                                                                                                                                ggTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTa      >  2:3884333/1‑150 (MQ=255)
                                                                                                                                                ggTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTa      <  1:167157/150‑1 (MQ=255)
                                                                                                                                                ggTTAATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTa      <  2:1008709/150‑1 (MQ=255)
                                                                                                                                                 gTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTaa     <  2:4893847/150‑1 (MQ=255)
                                                                                                                                                    cATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCa  <  2:4018388/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCA  >  CP009273/3001771‑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: ‑