Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,994,076 A→T 100% L114* (TTA→TAA)  yigE ← DUF2233 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,994,0760AT97.9% 867.2 / 5.1 246L114* (TTA→TAA) yigEDUF2233 family protein
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base T (132/108);  total (135/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 246 total aligned reads displayed.

ATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCG  >  CP009273/3993929‑3994219
                                                                                                                                                   |                                                                                                                                               
aTCAACATTGGCCCTGACTGCACCTAAAACTGAATCTCTTTACTTGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTAc                                                                                                                                               <  1:867024/150‑1 (MQ=255)
  cAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGc                                                                                                                                             >  2:1593733/1‑150 (MQ=255)
  cAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGc                                                                                                                                             <  2:3284285/150‑1 (MQ=255)
     cATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCAc                                                                                                                                          >  1:3112401/1‑150 (MQ=255)
      aTTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCAcc                                                                                                                                         <  2:2496742/150‑1 (MQ=255)
       ttGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCt                                                                                                                                        <  1:1689189/150‑1 (MQ=255)
        tGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCtt                                                                                                                                       >  1:3505624/1‑150 (MQ=255)
        tGGCCCTGACTGCACCGCAAACTGAATCTCTTAAATGGTTAAGAAGCCATCCAGACGAACGTTTACGACTTTATCTCCCGCGACATAAAACACGCAGCCAGGACCCATAAAGAAATTCCCTTCAACTGAAGAGAGATTTAACGCCACCtt                                                                                                                                       >  1:3112293/1‑150 (MQ=255)
         ggCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTc                                                                                                                                      <  1:463687/150‑1 (MQ=255)
          gCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAACACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTct                                                                                                                                     <  2:2251448/150‑1 (MQ=255)
            ccTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctgc                                                                                                                                   >  2:757186/1‑150 (MQ=255)
            ccTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctgc                                                                                                                                   >  2:757271/1‑150 (MQ=255)
             cTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctgct                                                                                                                                  <  2:1465591/150‑1 (MQ=255)
              tGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTctgctg                                                                                                                                 <  2:1918670/150‑1 (MQ=255)
               gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGa                                                                                                                                <  2:824098/150‑1 (MQ=255)
                aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGAc                                                                                                                               >  1:1729155/1‑150 (MQ=255)
                aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGAc                                                                                                                               >  1:1729175/1‑150 (MQ=255)
                 cTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGAcc                                                                                                                              <  1:1153517/150‑1 (MQ=255)
                  tGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCg                                                                                                                             <  1:2918407/150‑1 (MQ=255)
                    cACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGtt                                                                                                                           <  1:1019788/150‑1 (MQ=255)
                      ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGACCGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGtttt                                                                                                                         <  1:542975/150‑1 (MQ=255)
                       cGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTc                                                                                                                        >  1:4010206/1‑150 (MQ=255)
                                aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACaa                                                                                                               >  1:2814114/1‑150 (MQ=255)
                                    tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg                                                                                                           <  2:959705/150‑1 (MQ=255)
                                     ctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGa                                                                                                          <  1:2752799/150‑1 (MQ=255)
                                          cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc                                                                                                     <  2:1986322/150‑1 (MQ=255)
                                          cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc                                                                                                     >  1:1568186/1‑150 (MQ=255)
                                          cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc                                                                                                     <  1:2186892/150‑1 (MQ=255)
                                            ggTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGgcgc                                                                                                   <  1:543252/150‑1 (MQ=255)
                                               tttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATa                                                                                                >  2:2973546/1‑150 (MQ=255)
                                                ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAg                                                                                               >  2:2165173/1‑150 (MQ=255)
                                                ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAg                                                                                               >  2:3656966/1‑150 (MQ=255)
                                                 tGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGc                                                                                              <  1:3641869/150‑1 (MQ=255)
                                                  gAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCt                                                                                             >  1:334204/1‑150 (MQ=255)
                                                   aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAATCCGAGCGGCGCATAGCtt                                                                                            >  2:1759052/1‑150 (MQ=255)
                                                    aGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCttt                                                                                           >  1:1154165/1‑150 (MQ=255)
                                                     ggCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTtc                                                                                          >  2:3896966/1‑150 (MQ=255)
                                                        aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACAAAAAACACGCCGCCTGGACGGATAAAGAAATTCCATTCACCTGAAGCGAGATTTTACGCCACCATCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTacatc                                                                                       >  1:2038697/1‑150 (MQ=255)
                                                          ccAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTcatcat                                                                                     <  1:2428001/150‑1 (MQ=255)
                                                            aGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAg                                                                                   <  2:1436130/150‑1 (MQ=255)
                                                             gACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGa                                                                                  <  2:1746155/150‑1 (MQ=255)
                                                              acgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGAt                                                                                 <  2:3563772/150‑1 (MQ=255)
                                                              acgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGAt                                                                                 >  2:2208535/1‑150 (MQ=255)
                                                                gaacgaTGCCGACTTTATCTCCCGCGACATAAAACTCGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgc                                                                               <  2:2547766/150‑1 (MQ=255)
                                                                gaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgc                                                                               >  2:2795652/1‑150 (MQ=255)
                                                                    gaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgccgcc                                                                           >  1:4105775/1‑150 (MQ=255)
                                                                    gaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgccgcc                                                                           >  1:479652/1‑150 (MQ=255)
                                                                    gaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATgccgcc                                                                           <  2:1887044/150‑1 (MQ=255)
                                                                       gCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGtt                                                                        <  1:2228368/150‑1 (MQ=255)
                                                                       gCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGtt                                                                        >  2:3889008/1‑150 (MQ=255)
                                                                         cGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCa                                                                      >  1:2488903/1‑150 (MQ=255)
                                                                          gACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCAt                                                                     >  1:1560561/1‑150 (MQ=255)
                                                                          gACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCAt                                                                     >  2:775594/1‑150 (MQ=255)
                                                                          gACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCAt                                                                     >  1:2715095/1‑150 (MQ=255)
                                                                               tATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCa                                                                <  1:1469428/150‑1 (MQ=255)
                                                                                  ctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCt                                                             >  2:1238025/1‑150 (MQ=255)
                                                                                  ctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCt                                                             <  2:3241180/150‑1 (MQ=255)
                                                                                   tcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTg                                                            >  2:1027089/1‑150 (MQ=255)
                                                                                         gACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctg                                                      >  2:731268/1‑150 (MQ=255)
                                                                                          aCATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctga                                                     >  2:3845110/1‑150 (MQ=255)
                                                                                          aCATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctga                                                     >  2:3583789/1‑150 (MQ=255)
                                                                                           cATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgac                                                    <  1:3195822/150‑1 (MQ=255)
                                                                                           cATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgac                                                    >  2:542098/1‑150 (MQ=255)
                                                                                            aTAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgacc                                                   >  1:1834319/1‑150 (MQ=255)
                                                                                            aTAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgacc                                                   >  1:39723/1‑150 (MQ=255)
                                                                                             tAAAACACGCCGCCAGGACGGATAAAGAAATACCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccc                                                  >  2:2222725/1‑150 (MQ=255)
                                                                                              aaaaCACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccct                                                 >  1:2315967/1‑150 (MQ=255)
                                                                                              aaaaCACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccct                                                 <  2:2564371/150‑1 (MQ=255)
                                                                                               aaaCACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccctg                                                >  1:1878413/1‑150 (MQ=255)
                                                                                                 acacGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCAcctgaccctgac                                              >  1:2143274/1‑150 (MQ=255)
                                                                                                  cacGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACt                                             <  1:200601/150‑1 (MQ=255)
                                                                                                  cacGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACt                                             >  2:118944/1‑150 (MQ=255)
                                                                                                    cgccgcCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTAt                                           >  2:4099962/1‑150 (MQ=255)
                                                                                                      ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTa                                         <  1:3099312/150‑1 (MQ=255)
                                                                                                      ccgccAGGACGGATAAAGAAATGCCCTACACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGATTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTa                                         <  2:4175809/150‑1 (MQ=255)
                                                                                                          cAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAatat                                     >  2:3092186/1‑150 (MQ=255)
                                                                                                          cAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAatat                                     >  2:157198/1‑150 (MQ=255)
                                                                                                            ggACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATcc                                   <  2:342392/150‑1 (MQ=255)
                                                                                                             gACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCg                                  >  2:866812/1‑150 (MQ=255)
                                                                                                               cGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGcc                                >  1:580806/1‑150 (MQ=255)
                                                                                                                ggATAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCca                               >  1:3897519/1‑150 (MQ=255)
                                                                                                                  aTAAAGAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCcagc                             >  1:1783754/1‑150 (MQ=255)
                                                                                                                       gAAATTCCCTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGc                        >  2:2392291/1‑150 (MQ=255)
                                                                                                                             cccTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTaa                  >  1:1635844/1‑150 (MQ=255)
                                                                                                                             cccTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTaa                  <  2:2851679/150‑1 (MQ=255)
                                                                                                                              ccTTCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTGATATCCGCCAGCAGAGCATGTAAc                 >  2:3339162/1‑150 (MQ=255)
                                                                                                                               cTTCACCTGAATCGAGATTTTACGCCACCTTCTGCTGAACGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTTTTAATATCCGCCATCAGTGCATGTaaag                <  2:134395/150‑3 (MQ=255)
                                                                                                                                 tCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGtt              >  2:3640219/1‑150 (MQ=255)
                                                                                                                                 tCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGtt              <  2:4105775/150‑1 (MQ=255)
                                                                                                                                 tCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGtt              >  2:3577757/1‑150 (MQ=255)
                                                                                                                                 tCACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGtt              >  1:409594/1‑150 (MQ=255)
                                                                                                                                  cACCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTc             >  1:372909/1‑150 (MQ=255)
                                                                                                                                   aCCTGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTcc            >  1:574937/1‑150 (MQ=255)
                                                                                                                                      tGAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCa         <  1:1319493/150‑1 (MQ=255)
                                                                                                                                       gAAGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCAc        <  2:3977702/150‑1 (MQ=255)
                                                                                                                                        aaGCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCTATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCTCg       <  1:3084995/150‑1 (MQ=255)
                                                                                                                                          gCGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCt     <  2:1402254/150‑1 (MQ=255)
                                                                                                                                           cGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCtt    >  1:3592432/1‑150 (MQ=255)
                                                                                                                                           cGAGATTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCtt    >  2:109185/1‑150 (MQ=255)
                                                                                                                                             agaTTTTACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCg  >  1:304184/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ATCAACATTGGCCCTGACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCG  >  CP009273/3993929‑3994219

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: