Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,754,840 T→A 16.9% E687V (GAA→GTA)  yfjK ← DEAD/H helicase‑like protein, CP4‑57 putative defective prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,754,8400TA16.9% 681.8 / 122.3 248E687V (GAA→GTA) yfjKDEAD/H helicase‑like protein, CP4‑57 putative defective prophage
Reads supporting (aligned to +/- strand):  ref base T (114/92);  new base A (16/26);  total (130/118)
Fisher's exact test for biased strand distribution p-value = 4.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

GTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGG  >  CP009273/2754691‑2754989
                                                                                                                                                     |                                                                                                                                                     
  aGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTcc                                                                                                                                                     <  2:3881430/150‑1 (MQ=255)
      aaaCGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCa                                                                                                                                                 <  2:3329376/150‑1 (MQ=255)
        aCGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAg                                                                                                                                               >  1:4015377/1‑150 (MQ=255)
          gTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTc                                                                                                                                             >  2:2926745/1‑150 (MQ=255)
           tACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCt                                                                                                                                            >  1:1718168/1‑150 (MQ=255)
            aCCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTa                                                                                                                                           >  1:1053175/1‑150 (MQ=255)
             ccTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTAc                                                                                                                                          <  1:95009/150‑1 (MQ=255)
               tAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCa                                                                                                                                        <  2:4381586/150‑1 (MQ=255)
                   tCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCg                                                                                                                                    <  2:4344682/150‑1 (MQ=255)
                   tCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCg                                                                                                                                    >  1:761313/1‑150 (MQ=255)
                      aTGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGgtg                                                                                                                                 <  2:4426226/150‑1 (MQ=255)
                      aTGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGgtg                                                                                                                                 <  2:2095355/150‑1 (MQ=255)
                        gtcgtcTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTa                                                                                                                               <  1:1978875/150‑1 (MQ=255)
                             cTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTc                                                                                                                          <  2:2230756/150‑1 (MQ=255)
                             cTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTc                                                                                                                          <  1:164678/150‑1 (MQ=255)
                              tAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCt                                                                                                                         <  1:3253690/150‑1 (MQ=255)
                              tAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCt                                                                                                                         >  1:318566/1‑150 (MQ=255)
                               aaTACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTg                                                                                                                        >  2:2285337/1‑150 (MQ=255)
                                  aCCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAAt                                                                                                                     <  1:280648/150‑1 (MQ=255)
                                    cGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATgg                                                                                                                   >  1:3810688/1‑150 (MQ=255)
                                    cGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATgg                                                                                                                   <  2:1786684/150‑1 (MQ=255)
                                        cgcgcgcCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAAt                                                                                                               >  1:3063529/1‑150 (MQ=255)
                                        cgcgcgcCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAAt                                                                                                               <  2:3210487/150‑1 (MQ=255)
                                          cgcgcCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGc                                                                                                             >  2:3707690/1‑150 (MQ=255)
                                             gcCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCaa                                                                                                          <  2:1462652/150‑1 (MQ=255)
                                              ccTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAAt                                                                                                         >  2:4239464/1‑150 (MQ=255)
                                                  aTAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATctct                                                                                                     >  1:3854320/1‑150 (MQ=255)
                                                     aaCATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCtcttct                                                                                                  <  1:2493955/150‑1 (MQ=255)
                                                      aCATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTa                                                                                                 <  1:3410633/150‑1 (MQ=255)
                                                      aCATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTa                                                                                                 <  1:1301923/150‑1 (MQ=255)
                                                        aTCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAg                                                                                               >  2:764101/1‑150 (MQ=255)
                                                          cTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGc                                                                                             >  1:2083950/1‑150 (MQ=255)
                                                            gATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTATGGCAt                                                                                           <  1:1371552/150‑1 (MQ=255)
                                                             aTCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCAtt                                                                                          <  1:100299/150‑1 (MQ=255)
                                                             aTCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCAtt                                                                                          <  1:3689138/150‑1 (MQ=255)
                                                              tCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTg                                                                                         <  2:1717576/150‑1 (MQ=255)
                                                              tCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTg                                                                                         <  2:1717900/150‑1 (MQ=255)
                                                              tCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTg                                                                                         <  1:2282156/150‑1 (MQ=255)
                                                               cAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGa                                                                                        >  1:4228007/1‑150 (MQ=255)
                                                                aaCGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGac                                                                                       >  1:414714/1‑150 (MQ=255)
                                                                  cGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGacac                                                                                     >  2:854496/1‑150 (MQ=255)
                                                                   gTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACAcc                                                                                    >  2:2188971/1‑150 (MQ=255)
                                                                   gTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACAcc                                                                                    <  2:517163/150‑1 (MQ=255)
                                                                     aCCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCtg                                                                                  <  1:567644/150‑1 (MQ=255)
                                                                     aCCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCtg                                                                                  <  1:3747269/150‑1 (MQ=255)
                                                                        cAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCtggtg                                                                               <  2:679119/150‑1 (MQ=255)
                                                                        cAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCtggtg                                                                               >  2:304461/1‑150 (MQ=255)
                                                                          aTCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTggtggt                                                                             >  2:2164710/1‑150 (MQ=255)
                                                                            cGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAg                                                                           >  2:3589355/1‑150 (MQ=255)
                                                                              aTCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGAt                                                                         >  1:1578518/1‑150 (MQ=255)
                                                                               tCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATg                                                                        >  1:591067/1‑150 (MQ=255)
                                                                                 tACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGa                                                                      <  1:4338212/150‑1 (MQ=255)
                                                                                  aCTCCTAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACCTCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGaa                                                                     <  2:2666795/150‑1 (MQ=255)
                                                                                   cTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGaaa                                                                    <  2:4564161/150‑1 (MQ=255)
                                                                                   cTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGaaa                                                                    <  1:3957079/150‑1 (MQ=255)
                                                                                    tCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAAc                                                                   <  2:371077/150‑1 (MQ=255)
                                                                                      cAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACtt                                                                 <  1:4064204/150‑1 (MQ=255)
                                                                                       aaaTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACttt                                                                >  1:3353852/1‑150 (MQ=255)
                                                                                       aaaTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACttt                                                                >  1:50708/1‑150 (MQ=255)
                                                                                       aaaTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACttt                                                                <  1:935404/150‑1 (MQ=255)
                                                                                        aaTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTc                                                               <  2:4094846/150‑1 (MQ=255)
                                                                                         aTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCg                                                              >  1:465205/1‑150 (MQ=255)
                                                                                         aTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCg                                                              >  2:2160009/1‑150 (MQ=255)
                                                                                         aTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCg                                                              <  2:3936147/150‑1 (MQ=255)
                                                                                         aTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCg                                                              >  1:353722/1‑150 (MQ=255)
                                                                                          tGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGa                                                             >  2:1091776/1‑150 (MQ=255)
                                                                                          tGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGa                                                             >  2:1529924/1‑150 (MQ=255)
                                                                                           gTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGaa                                                            >  1:3291348/1‑150 (MQ=255)
                                                                                           gTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGaa                                                            >  1:2846249/1‑150 (MQ=255)
                                                                                             ccTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACg                                                          >  1:2096459/1‑150 (MQ=255)
                                                                                             ccTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACg                                                          <  1:3637304/150‑1 (MQ=255)
                                                                                                gAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCa                                                       >  2:2971983/1‑150 (MQ=255)
                                                                                                 aGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCaa                                                      <  1:1217207/150‑1 (MQ=255)
                                                                                                     tCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTg                                                  <  2:1624297/150‑1 (MQ=255)
                                                                                                         gAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACAt                                              >  2:1162355/1‑150 (MQ=255)
                                                                                                         gAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACAt                                              >  1:4402407/1‑150 (MQ=255)
                                                                                                          aaGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACAtg                                             >  1:2988250/1‑150 (MQ=255)
                                                                                                           aGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACAtgt                                            >  2:1570958/1‑150 (MQ=255)
                                                                                                            gATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACAtgtg                                           <  2:3079637/150‑1 (MQ=255)
                                                                                                             aTACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACAtgtgt                                          >  1:229563/1‑150 (MQ=255)
                                                                                                                cTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTa                                       >  2:2702893/1‑150 (MQ=255)
                                                                                                                cTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTa                                       <  2:620881/150‑1 (MQ=255)
                                                                                                                cTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTa                                       >  1:2949833/1‑150 (MQ=255)
                                                                                                                  ggCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAAt                                     >  2:2007510/1‑150 (MQ=255)
                                                                                                                   gCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATc                                    <  2:3032104/150‑1 (MQ=255)
                                                                                                                   gCTCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATc                                    <  1:2363890/150‑1 (MQ=255)
                                                                                                                     tCAACTCATCCACATCAGCCTCGTCGGAAAATACCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGa                                  >  1:3828098/1‑150 (MQ=255)
                                                                                                                        aCTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACtt                               <  1:2674471/150‑1 (MQ=255)
                                                                                                                           cATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTtgc                            >  1:770724/1‑150 (MQ=255)
                                                                                                                                acaTCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGt                       >  2:4158775/1‑150 (MQ=255)
                                                                                                                                acaTCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGt                       <  2:3502537/150‑1 (MQ=255)
                                                                                                                                   tCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGtata                    >  1:1311355/1‑150 (MQ=255)
                                                                                                                                   tCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGtata                    >  1:2909219/1‑150 (MQ=255)
                                                                                                                                    cAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATAc                   >  1:1427728/1‑150 (MQ=255)
                                                                                                                                      gCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCa                 >  2:3386950/1‑150 (MQ=255)
                                                                                                                                        cTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACg               <  1:2223782/150‑1 (MQ=255)
                                                                                                                                        cTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACg               <  2:2383814/150‑1 (MQ=255)
                                                                                                                                          cgtcgGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCt             >  2:505592/1‑150 (MQ=255)
                                                                                                                                              ggAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTAc         >  2:1899266/1‑150 (MQ=255)
                                                                                                                                                     tCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtgg  <  2:3964151/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAGACAAACGTACCTAAATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGG  >  CP009273/2754691‑2754989

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: