Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,984,478 G→T 17.8% P208Q (CCG→CAG)  yecA ← UPF0149 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,984,4780GT17.8% 756.5 / 104.9 247P208Q (CCG→CAG) yecAUPF0149 family protein
Reads supporting (aligned to +/- strand):  ref base G (119/84);  new base T (22/22);  total (141/106)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 251 total aligned reads displayed.

GCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACATC  >  CP009273/1984331‑1984627
                                                                                                                                                   |                                                                                                                                                     
gCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCAAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACggg                                                                                                                                                     <  1:1413946/150‑1 (MQ=255)
  aaTTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCa                                                                                                                                                   >  1:1522245/1‑150 (MQ=255)
   aTTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCaa                                                                                                                                                  <  2:981088/150‑1 (MQ=255)
    tttttGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAg                                                                                                                                                 >  2:3927017/1‑150 (MQ=255)
      tttGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGa                                                                                                                                               >  2:3642203/1‑150 (MQ=255)
        tGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATc                                                                                                                                             >  1:1175460/1‑150 (MQ=255)
          aTTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGt                                                                                                                                           <  1:1557188/150‑1 (MQ=255)
           ttGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGtt                                                                                                                                          <  1:771243/150‑1 (MQ=255)
              cGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACg                                                                                                                                       >  2:1739776/1‑150 (MQ=255)
                  aGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGcccc                                                                                                                                   >  2:2449009/1‑150 (MQ=255)
                   gggTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCg                                                                                                                                  <  1:3199231/150‑1 (MQ=255)
                   gggTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCg                                                                                                                                  <  1:2880137/150‑1 (MQ=255)
                    ggTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCgg                                                                                                                                 >  2:913013/1‑150 (MQ=255)
                       tgttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGttt                                                                                                                              >  1:1000785/1‑150 (MQ=255)
                       tgttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGttt                                                                                                                              >  1:1983195/1‑150 (MQ=255)
                        gttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtc                                                                                                                             >  1:1163891/1‑150 (MQ=255)
                          tACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtctc                                                                                                                           >  2:1985961/1‑150 (MQ=255)
                           aCCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTctct                                                                                                                          <  1:753123/150‑1 (MQ=255)
                              aGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCtcttct                                                                                                                       <  2:1051045/150‑1 (MQ=255)
                              aGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCtcttct                                                                                                                       >  1:837189/1‑150 (MQ=255)
                              aGAAGTGGGGTGATCAGATAGCCTCAAATAACTTATTGTGTTCCAGAATTAACGCTGACACACAATTTGGACTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCtcttct                                                                                                                       >  1:2020740/1‑150 (MQ=255)
                                   tGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCttt                                                                                                                  >  1:248141/1‑150 (MQ=255)
                                    ggggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCNGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGNTNAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTNCTCTTCTGCtttn                                                                                                                 <  2:4528429/150‑2 (MQ=255)
                                    ggggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTg                                                                                                                 >  1:1778973/1‑150 (MQ=255)
                                     gggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATATCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGa                                                                                                                >  1:3825014/1‑150 (MQ=255)
                                       gTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATc                                                                                                              <  1:400600/150‑1 (MQ=255)
                                        tGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCg                                                                                                             <  2:1629822/150‑1 (MQ=255)
                                            cAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTg                                                                                                         >  2:1284708/1‑150 (MQ=255)
                                              gATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCt                                                                                                       <  2:3766216/150‑1 (MQ=255)
                                                  gCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAAc                                                                                                   >  1:3118998/1‑150 (MQ=255)
                                                  gCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAAc                                                                                                   >  1:19437/1‑150 (MQ=255)
                                                  gCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTGACCACTACCGCACGGGCACGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAAc                                                                                                   >  1:245921/1‑150 (MQ=255)
                                                   ccTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAAGTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACa                                                                                                  >  1:1655286/1‑150 (MQ=255)
                                                   acTCAAATTCCTTATTGGTTGCCATAATTAACTCTGAAACCCAATTTGGCCTCTTAATGCAGGCAGCACGGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATAGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACa                                                                                                  <  1:2464024/149‑1 (MQ=255)
                                                    cTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAg                                                                                                 >  2:771990/1‑150 (MQ=255)
                                                         aTTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCtttt                                                                                            <  2:1695990/150‑1 (MQ=255)
                                                          ttCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCttttt                                                                                           <  1:4215977/150‑1 (MQ=255)
                                                            ccTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc                                                                                         >  1:3395637/1‑150 (MQ=255)
                                                            ccTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTcc                                                                                         <  2:3493934/150‑1 (MQ=255)
                                                             cTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCt                                                                                        <  2:2251088/150‑1 (MQ=255)
                                                              ttATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTg                                                                                       >  1:4483216/1‑150 (MQ=255)
                                                                 ttGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCgg                                                                                    >  1:2737267/1‑150 (MQ=255)
                                                                   gggTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGAt                                                                                  >  1:211283/1‑150 (MQ=255)
                                                                     gTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATgc                                                                                <  2:1663127/150‑1 (MQ=255)
                                                                     gTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATgc                                                                                <  2:1175460/150‑1 (MQ=255)
                                                                      tGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATgcg                                                                               <  2:2185389/150‑1 (MQ=255)
                                                                      tGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATgcg                                                                               <  2:1925513/150‑1 (MQ=255)
                                                                        ccAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGcc                                                                             >  1:569042/1‑150 (MQ=255)
                                                                        ccAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGcc                                                                             >  2:4497596/1‑150 (MQ=255)
                                                                          aGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCAt                                                                           >  1:2364195/1‑150 (MQ=255)
                                                                                  cGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACg                                                                   >  1:1159893/1‑150 (MQ=255)
                                                                                  cGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACg                                                                   <  2:1315600/150‑1 (MQ=255)
                                                                                   gCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGc                                                                  >  2:929809/1‑150 (MQ=255)
                                                                                   gCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGc                                                                  >  1:3756697/1‑150 (MQ=255)
                                                                                   gCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGc                                                                  >  1:3757145/1‑150 (MQ=255)
                                                                                     tGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCAt                                                                >  1:177801/1‑150 (MQ=255)
                                                                                      gACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATg                                                               >  2:2428734/1‑150 (MQ=255)
                                                                                       acacCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGt                                                              >  2:1038739/1‑150 (MQ=255)
                                                                                       acacCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGt                                                              >  2:2321280/1‑150 (MQ=255)
                                                                                          cccAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGa                                                           >  2:2676752/1‑150 (MQ=255)
                                                                                            cAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATc                                                         <  1:2052754/150‑1 (MQ=255)
                                                                                            cAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATc                                                         <  2:3718879/150‑1 (MQ=255)
                                                                                               tttGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAg                                                      <  1:46424/150‑1 (MQ=255)
                                                                                                 tGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCACGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAgcg                                                    >  1:1387708/1‑150 (MQ=255)
                                                                                                       ctTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCa                                              >  2:3001509/1‑150 (MQ=255)
                                                                                                        ttAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAg                                             <  2:4153566/150‑1 (MQ=255)
                                                                                                         tAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGt                                            >  2:3613574/1‑150 (MQ=255)
                                                                                                            tGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGa                                         <  2:4346574/150‑1 (MQ=255)
                                                                                                                ggcagCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATgg                                     <  1:62665/150‑1 (MQ=255)
                                                                                                                  cagcaCTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCa                                   >  2:575542/1‑150 (MQ=255)
                                                                                                                   agcaCTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCAt                                  >  1:695072/1‑150 (MQ=255)
                                                                                                                   agcaCTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCAt                                  >  2:1569898/1‑150 (MQ=255)
                                                                                                                   agcaCTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCAt                                  >  2:3255014/1‑150 (MQ=255)
                                                                                                                     caCTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCt                                >  2:3072379/1‑150 (MQ=255)
                                                                                                                     caCTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCt                                >  2:2412799/1‑150 (MQ=255)
                                                                                                                     caCTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCt                                <  2:2737267/150‑1 (MQ=255)
                                                                                                                     caCTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCt                                >  2:4303708/1‑150 (MQ=255)
                                                                                                                       cTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTAc                              >  2:2956060/1‑150 (MQ=255)
                                                                                                                       cTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTAc                              <  1:4315671/150‑1 (MQ=255)
                                                                                                                        tGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACg                             >  1:1886999/1‑150 (MQ=255)
                                                                                                                        tGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACg                             >  1:3582150/1‑150 (MQ=255)
                                                                                                                            tAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGctct                         >  1:4144197/1‑150 (MQ=255)
                                                                                                                             aaaTTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGctctc                        <  2:740608/150‑1 (MQ=255)
                                                                                                                                 ttCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCa                    <  2:2024276/150‑1 (MQ=255)
                                                                                                                                     tACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACg                >  2:1485659/1‑150 (MQ=255)
                                                                                                                                      aCCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGc               >  1:1499898/1‑150 (MQ=255)
                                                                                                                                       ccACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCt              >  1:3027999/1‑150 (MQ=255)
                                                                                                                                        cACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCtt             <  1:3642202/150‑1 (MQ=255)
                                                                                                                                         aCTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTc            <  2:1778973/150‑1 (MQ=255)
                                                                                                                                         aCTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTc            >  1:2717907/1‑150 (MQ=255)
                                                                                                                                          cTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCt           >  2:3185013/1‑150 (MQ=255)
                                                                                                                                           tACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTg          >  2:389050/1‑150 (MQ=255)
                                                                                                                                           tACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTg          >  2:1945155/1‑150 (MQ=255)
                                                                                                                                            aCCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTgg         >  2:4613004/1‑150 (MQ=255)
                                                                                                                                               gCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGa      >  1:852055/1‑150 (MQ=255)
                                                                                                                                               gCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGa      >  1:579804/1‑150 (MQ=255)
                                                                                                                                                 aCGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACa    >  2:1385240/1‑150 (MQ=255)
                                                                                                                                                 aCGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACa    >  2:2795883/1‑150 (MQ=255)
                                                                                                                                                   gggCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACATc  <  2:2857932/150‑1 (MQ=255)
                                                                                                                                                   gggCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCAGCTACGCTCGCTTCAAACGCTTCTGGCGACATc  >  2:3337230/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACATC  >  CP009273/1984331‑1984627

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: ‑