Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 997,965 G→T 92.1% D235Y (GAT→TAT)  elfG → putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273997,9650GT92.1% 750.6 / 51.7 257D235Y (GAT→TAT) elfGputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (13/8);  new base T (124/112);  total (137/120)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

CGCGATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCAGC  >  CP009273/997816‑998111
                                                                                                                                                     |                                                                                                                                                  
  cgATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTAt                                                                                                                                                  <  2:2642382/150‑1 (MQ=255)
  cgATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGAt                                                                                                                                                  >  2:4347202/1‑150 (MQ=255)
  cgATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGAACACCGGATTTAGTGTATTACCTGTCTGGAAAGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAAATGAAGTTGATTTAGGTTAt                                                                                                                                                  >  1:2396606/1‑150 (MQ=255)
   gATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATg                                                                                                                                                 >  1:4232346/1‑150 (MQ=255)
    aTTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGt                                                                                                                                                <  1:4331335/150‑1 (MQ=255)
     ttATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTc                                                                                                                                               <  1:2079050/150‑1 (MQ=255)
       atatGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCtt                                                                                                                                             <  2:2977082/150‑1 (MQ=255)
       atatGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCtt                                                                                                                                             >  1:508016/1‑150 (MQ=255)
       atatGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCtt                                                                                                                                             <  2:3692820/150‑1 (MQ=255)
       atatGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCtt                                                                                                                                             >  1:3666683/1‑150 (MQ=255)
        tatGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTc                                                                                                                                            <  2:1314036/150‑1 (MQ=255)
         atGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCg                                                                                                                                           <  2:1624893/150‑1 (MQ=255)
           gAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCt                                                                                                                                         >  1:407633/1‑150 (MQ=255)
             aTGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGc                                                                                                                                       <  2:1990143/150‑1 (MQ=255)
                 tACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAAt                                                                                                                                   <  1:2262312/150‑1 (MQ=255)
                 tACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAAt                                                                                                                                   >  2:1477928/1‑150 (MQ=255)
                 tACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAAt                                                                                                                                   <  1:3111534/150‑1 (MQ=255)
                    aacaTACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTc                                                                                                                                <  1:3937508/150‑1 (MQ=255)
                    aacaTACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTc                                                                                                                                <  1:3935091/150‑1 (MQ=255)
                    aacaTACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTc                                                                                                                                >  1:1318712/1‑150 (MQ=255)
                      caTACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCg                                                                                                                              >  2:4193203/1‑150 (MQ=255)
                       aTACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCgt                                                                                                                             >  1:193655/1‑150 (MQ=255)
                           ccGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGtgttg                                                                                                                         >  2:4424387/1‑150 (MQ=255)
                           ccGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGtgttg                                                                                                                         <  2:1754309/150‑1 (MQ=255)
                             gCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTa                                                                                                                       <  1:1394686/150‑1 (MQ=255)
                              cAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAg                                                                                                                      >  1:2685900/1‑150 (MQ=255)
                               aGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAgg                                                                                                                     >  2:1367341/1‑150 (MQ=255)
                                ggTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAggg                                                                                                                    <  1:2721973/150‑1 (MQ=255)
                                 gTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGc                                                                                                                   >  2:1771220/1‑150 (MQ=255)
                                 gTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGc                                                                                                                   >  1:2873951/1‑150 (MQ=255)
                                     ttCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATaa                                                                                                               >  1:321488/1‑150 (MQ=255)
                                      tCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATaaa                                                                                                              <  2:4218416/150‑1 (MQ=255)
                                       ccTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAAc                                                                                                             <  2:792284/150‑1 (MQ=255)
                                         tGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACct                                                                                                           <  2:2800851/150‑1 (MQ=255)
                                          gCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACctc                                                                                                          <  1:3813683/150‑1 (MQ=255)
                                            aCGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCtt                                                                                                        >  1:1815185/1‑150 (MQ=255)
                                              gACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTgg                                                                                                      <  1:312708/150‑1 (MQ=255)
                                               aCTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTggg                                                                                                     <  1:938356/150‑1 (MQ=255)
                                               aCTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTggg                                                                                                     <  2:1950644/150‑1 (MQ=255)
                                                  aCAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGcc                                                                                                  <  1:1960731/150‑1 (MQ=255)
                                                   cAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCa                                                                                                 <  1:817852/150‑1 (MQ=255)
                                                    aGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAg                                                                                                >  1:4593768/1‑150 (MQ=255)
                                                       tACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAAc                                                                                             <  2:845398/150‑1 (MQ=255)
                                                         ctccGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACgg                                                                                           <  2:2371052/148‑1 (MQ=255)
                                                         cacCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACgg                                                                                           >  2:2481685/1‑150 (MQ=255)
                                                         cacCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACgg                                                                                           >  1:4625970/1‑150 (MQ=255)
                                                          acCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGc                                                                                          >  2:938051/1‑150 (MQ=255)
                                                          acCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGc                                                                                          >  2:3567405/1‑150 (MQ=255)
                                                             ggTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGa                                                                                       <  1:3303860/150‑1 (MQ=255)
                                                             ggTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGa                                                                                       >  1:1983546/1‑150 (MQ=255)
                                                              gTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGaa                                                                                      >  1:4114427/1‑150 (MQ=255)
                                                              gTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGaa                                                                                      >  2:4034472/1‑150 (MQ=255)
                                                               ttttAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAAc                                                                                     <  2:617183/150‑1 (MQ=255)
                                                                    gtgtATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCa                                                                                <  2:3673213/150‑1 (MQ=255)
                                                                    gtgtATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCa                                                                                <  2:3989538/150‑1 (MQ=255)
                                                                      gtATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAAt                                                                              <  2:3338429/150‑1 (MQ=255)
                                                                        aTTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTc                                                                            >  1:1735867/1‑150 (MQ=255)
                                                                         ttACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTcc                                                                           <  2:3116243/150‑1 (MQ=255)
                                                                         ttACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTcc                                                                           <  1:376533/150‑1 (MQ=255)
                                                                             cTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTc                                                                       >  1:2107783/1‑150 (MQ=255)
                                                                              tGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCa                                                                      >  2:2848298/1‑150 (MQ=255)
                                                                               gTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAg                                                                     >  2:238809/1‑150 (MQ=255)
                                                                                tCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAgg                                                                    >  2:2252490/1‑150 (MQ=255)
                                                                                 cTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGgt                                                                   <  1:1678499/150‑1 (MQ=255)
                                                                                  tGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGgtg                                                                  >  1:360241/1‑150 (MQ=255)
                                                                                  tGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGgtg                                                                  <  2:2559344/150‑1 (MQ=255)
                                                                                        gATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAacac                                                            >  2:1102439/1‑150 (MQ=255)
                                                                                         aTCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACg                                                           <  2:1217346/150‑1 (MQ=255)
                                                                                         aTCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACg                                                           <  2:2156841/150‑1 (MQ=255)
                                                                                           cAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACggg                                                         >  1:1227223/1‑150 (MQ=255)
                                                                                            aaTTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAAGTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGa                                                        >  1:2937734/1‑150 (MQ=255)
                                                                                             aTTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGaa                                                       <  1:21923/150‑1 (MQ=255)
                                                                                               tCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAAc                                                     >  2:875177/1‑150 (MQ=255)
                                                                                                  cTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAAcgcg                                                  <  2:2071754/150‑1 (MQ=255)
                                                                                                  cTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAAcgcg                                                  <  2:685642/150‑1 (MQ=255)
                                                                                                   ttGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCgg                                                 <  2:328830/150‑1 (MQ=255)
                                                                                                     ggCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGa                                               <  2:182219/150‑1 (MQ=255)
                                                                                                        tCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTa                                            <  2:2391359/150‑1 (MQ=255)
                                                                                                        tCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTa                                            <  2:4322079/150‑1 (MQ=255)
                                                                                                         cATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAg                                           <  2:1240432/150‑1 (MQ=255)
                                                                                                          aTGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGt                                          <  1:863308/150‑1 (MQ=255)
                                                                                                          aTGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGt                                          <  2:272522/150‑1 (MQ=255)
                                                                                                           tggtCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGtt                                         <  1:3544030/147‑1 (MQ=255)
                                                                                                                            agagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTTCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCt                        >  1:2795014/1‑150 (MQ=255)
                                                                                                                            agagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCt                        >  1:1022974/1‑150 (MQ=255)
                                                                                                                              agaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGa                      <  2:3814008/150‑1 (MQ=255)
                                                                                                                               gaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGAc                     >  2:4282638/1‑150 (MQ=255)
                                                                                                                               gaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGAc                     <  1:4204525/150‑1 (MQ=255)
                                                                                                                               gaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGAc                     >  2:1142998/1‑150 (MQ=255)
                                                                                                                                 cAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACgg                   >  1:3757667/1‑150 (MQ=255)
                                                                                                                                   aTTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCa                 >  1:728254/1‑150 (MQ=255)
                                                                                                                                   aTTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCa                 <  2:1863029/150‑1 (MQ=255)
                                                                                                                                     tGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAAc               >  2:3807970/1‑150 (MQ=255)
                                                                                                                                     tGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAAc               >  2:4552916/1‑150 (MQ=255)
                                                                                                                                      gAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAAcc              >  1:2556055/1‑150 (MQ=255)
                                                                                                                                       aaGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACca             <  2:3140058/150‑1 (MQ=255)
                                                                                                                                            gATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGAc        >  1:283548/1‑150 (MQ=255)
                                                                                                                                             aTTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGAcc       >  2:1485970/1‑150 (MQ=255)
                                                                                                                                                tAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCa    >  1:4573914/1‑150 (MQ=255)
                                                                                                                                                 aGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCAg   <  1:163966/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGCGATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCAGC  >  CP009273/997816‑998111

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: