Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,984,478 G→T 44.9% 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,4780GT44.9% 175.0 / 260.6 204P208Q (CCG→CAG) yecAUPF0149 family protein
Reads supporting (aligned to +/- strand):  ref base G (64/48);  new base T (43/49);  total (107/97)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 206 total aligned reads displayed.

TGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACAT  >  CP009273/1984330‑1984626
                                                                                                                                                    |                                                                                                                                                    
tGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACgg                                                                                                                                                     <  1:2275682/150‑1 (MQ=255)
   aaTTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCa                                                                                                                                                  >  1:98180/1‑150 (MQ=255)
    aTTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCaa                                                                                                                                                 <  2:2263170/150‑1 (MQ=255)
    aTTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCaa                                                                                                                                                 >  2:633252/1‑150 (MQ=255)
       tttGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGa                                                                                                                                              <  2:3823403/150‑1 (MQ=255)
        ttGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGAt                                                                                                                                             <  1:1326708/150‑1 (MQ=255)
         tGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAAGTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATc                                                                                                                                            <  1:2944844/150‑1 (MQ=255)
           aTTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGt                                                                                                                                          <  1:4226586/150‑1 (MQ=255)
             tGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGTGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTa                                                                                                                                        <  1:3883314/150‑1 (MQ=255)
                   aGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGcccc                                                                                                                                  >  2:3057118/1‑150 (MQ=255)
                      gttgttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGt                                                                                                                               <  2:1124042/150‑1 (MQ=255)
                       ttgttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGtt                                                                                                                              <  1:759561/150‑1 (MQ=255)
                        tgttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGTATCGTTACGCCCCGGatt                                                                                                                             >  2:1539848/1‑150 (MQ=255)
                        tgttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGttt                                                                                                                             <  1:1728132/150‑1 (MQ=255)
                        tgttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGttt                                                                                                                             >  1:1356347/1‑150 (MQ=255)
                          ttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTtct                                                                                                                           >  2:1680913/1‑150 (MQ=255)
                          ttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTtct                                                                                                                           <  1:1943085/150‑1 (MQ=255)
                           tACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtctc                                                                                                                          >  2:4593855/1‑150 (MQ=255)
                           tACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtctc                                                                                                                          >  1:4323158/1‑150 (MQ=255)
                             ccAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCtctt                                                                                                                        <  2:3133495/150‑1 (MQ=255)
                              cAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCtcttc                                                                                                                       >  1:748134/1‑150 (MQ=255)
                              cAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCtcttc                                                                                                                       <  2:438498/150‑1 (MQ=255)
                               aGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCtcttct                                                                                                                      <  1:1433117/150‑1 (MQ=255)
                               aGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCtcttct                                                                                                                      >  1:3581235/1‑150 (MQ=255)
                                gAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTg                                                                                                                     >  2:1453272/1‑150 (MQ=255)
                                gAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTg                                                                                                                     >  2:2715417/1‑150 (MQ=255)
                                  aGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCt                                                                                                                   >  2:4601978/1‑150 (MQ=255)
                                   gTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCtt                                                                                                                  <  2:2363019/150‑1 (MQ=255)
                                   gTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCtt                                                                                                                  >  1:4035720/1‑150 (MQ=255)
                                      gggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGTCACCCAATTTGGCCTCTTTATGCAGGCAGCACTGCTTATTTTTCTTACCACTACCTCACTGGCAAGGATGGTTACGCCCCGGTTTGTCTTCTGCTTTGa                                                                                                               <  2:24008/150‑1 (MQ=255)
                                      gggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGa                                                                                                               <  2:2395223/150‑1 (MQ=255)
                                        gTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATc                                                                                                             <  1:1050856/150‑1 (MQ=255)
                                           aTCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGt                                                                                                          >  1:3843520/1‑150 (MQ=255)
                                                   gCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAAc                                                                                                  <  2:1639599/150‑1 (MQ=255)
                                                    ccTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACa                                                                                                 >  1:4579796/1‑150 (MQ=255)
                                                        caaTTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCtt                                                                                             >  2:2770490/2‑150 (MQ=255)
                                                         aaTTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCttt                                                                                            <  2:2621591/150‑1 (MQ=255)
                                                          aTTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCtttt                                                                                           <  2:3479039/150‑1 (MQ=255)
                                                           ttCCTTATTGGGTGCCAGAATTAACGCTGACACCCTATTTGGACTTTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTTCCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGTTCGGTTGCTTAACAGCttttt                                                                                          <  2:2053243/150‑1 (MQ=255)
                                                            tCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTc                                                                                         >  1:3037406/1‑150 (MQ=255)
                                                            tCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTc                                                                                         >  1:3036725/1‑150 (MQ=255)
                                                              cTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCt                                                                                       >  2:1922531/1‑150 (MQ=255)
                                                              cTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCt                                                                                       <  2:3869517/150‑1 (MQ=255)
                                                               ttATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTg                                                                                      >  1:3339377/1‑150 (MQ=255)
                                                                          cAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCa                                                                           >  1:1608338/1‑150 (MQ=255)
                                                                            gAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATc                                                                         >  1:39162/1‑150 (MQ=255)
                                                                             aaTTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATcc                                                                        >  2:1790915/1‑150 (MQ=255)
                                                                                 aaCGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTa                                                                    <  1:1443435/150‑1 (MQ=255)
                                                                                  aCGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTAc                                                                   <  2:1982249/150‑1 (MQ=255)
                                                                                   cGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACg                                                                  >  2:2156917/1‑150 (MQ=255)
                                                                                   cGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACg                                                                  >  2:1134702/1‑150 (MQ=255)
                                                                                      tGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCAt                                                               >  1:1747302/1‑150 (MQ=255)
                                                                                      tGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCAt                                                               >  1:123003/1‑150 (MQ=255)
                                                                                      tGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCAt                                                               >  1:2225305/1‑150 (MQ=255)
                                                                                       gACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATg                                                              <  2:1599718/150‑1 (MQ=255)
                                                                                       gACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATg                                                              >  1:3312509/1‑150 (MQ=255)
                                                                                         ctcccAATTTTGCCTCTTAATTCAGGCAGCACTTCTTAAATTTCTTACCACTAACGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGATGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTa                                                            <  2:222587/148‑1 (MQ=255)
                                                                                          acCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAg                                                           >  2:3448008/1‑150 (MQ=255)
                                                                                            ccAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGAt                                                         >  2:3528960/1‑150 (MQ=255)
                                                                                              aaTTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCa                                                       >  1:1589628/1‑150 (MQ=255)
                                                                                              aaTTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCa                                                       >  1:1729765/1‑150 (MQ=255)
                                                                                               aTTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCaa                                                      >  2:1841003/1‑150 (MQ=255)
                                                                                                    gCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGcgcc                                                 >  1:4108271/1‑150 (MQ=255)
                                                                                                     ccTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGcgccg                                                <  1:3128343/150‑1 (MQ=255)
                                                                                                      ctctTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGcgccgc                                               >  2:2243931/1‑150 (MQ=255)
                                                                                                       tctTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCgccgcc                                              <  2:2912416/150‑1 (MQ=255)
                                                                                                           aaTGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTc                                          >  1:3758560/1‑150 (MQ=255)
                                                                                                            aTGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCg                                         <  1:89918/150‑1 (MQ=255)
                                                                                                               caggcagCACTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAAt                                      >  1:477671/1‑150 (MQ=255)
                                                                                                               caggcagCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAAt                                      >  2:4052688/1‑150 (MQ=255)
                                                                                                                aggcagCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATg                                     >  1:1289449/1‑150 (MQ=255)
                                                                                                                aggcagCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATg                                     >  1:917642/1‑150 (MQ=255)
                                                                                                                   cagcaCTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCa                                  <  2:1300244/150‑1 (MQ=255)
                                                                                                                    agcaCTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCAt                                 <  2:536235/150‑1 (MQ=255)
                                                                                                                    agcaCTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCAt                                 >  2:4133127/1‑150 (MQ=255)
                                                                                                                     gcaCTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATc                                >  1:1790804/1‑150 (MQ=255)
                                                                                                                        cTGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTAc                             >  2:3190308/1‑150 (MQ=255)
                                                                                                                         tGCTTAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACg                            <  2:3967436/150‑1 (MQ=255)
                                                                                                                          gCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGc                           >  2:4410477/1‑150 (MQ=255)
                                                                                                                           cTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGct                          >  2:3004706/1‑150 (MQ=255)
                                                                                                                           cTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGct                          >  1:3707382/1‑150 (MQ=255)
                                                                                                                            ttAAATTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGctc                         <  1:1883560/150‑1 (MQ=255)
                                                                                                                             tAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGctct                        <  2:4342758/150‑1 (MQ=255)
                                                                                                                              aaaTTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGctctc                       >  2:941492/1‑150 (MQ=255)
                                                                                                                               aaTTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGctctct                      >  1:836058/1‑150 (MQ=255)
                                                                                                                                aTTTCTTACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCtt                     <  1:621413/150‑1 (MQ=255)
                                                                                                                                aTTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCtt                     <  1:3051151/150‑1 (MQ=255)
                                                                                                                                aTTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCtt                     >  2:3434377/1‑150 (MQ=255)
                                                                                                                                  ttCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCa                   <  2:582124/150‑1 (MQ=255)
                                                                                                                                     ttACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAAc                <  2:2719713/150‑1 (MQ=255)
                                                                                                                                     ttACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAAc                >  1:262956/1‑150 (MQ=255)
                                                                                                                                      tACCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACg               <  1:2420654/150‑1 (MQ=255)
                                                                                                                                      tACCACTACCGCACGGGCAAGGATCTTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTATCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCAACGCTCTCATCAAACg               >  1:2296312/1‑150 (MQ=255)
                                                                                                                                       aCCACTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGc              >  2:1320255/1‑150 (MQ=255)
                                                                                                                                           cTACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCt          <  1:2158701/150‑1 (MQ=255)
                                                                                                                                           cTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCt          <  2:3354244/150‑1 (MQ=255)
                                                                                                                                            tACCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTg         <  1:401964/150‑1 (MQ=255)
                                                                                                                                             aCCGCACTGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTgg        <  1:1539768/150‑1 (MQ=255)
                                                                                                                                                   cGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACAt  <  1:2058723/150‑1 (MQ=255)
                                                                                                                                                   cGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACAt  <  2:4473336/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCCATCCAGTACGCATGTAGATCAAGCGCCGCCAGTCGAATGGCATCTACGCTCTCTTCAAACGCTTCTGGCGACAT  >  CP009273/1984330‑1984626

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: