Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 997,965 G→T 100% 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,9650GT96.6% 638.5 / 11.9 209D235Y (GAT→TAT) elfGputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base T (107/95);  total (109/100)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 209 total aligned reads displayed.

GATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCA  >  CP009273/997819‑998109
                                                                                                                                                  |                                                                                                                                                
gATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGa                                                                                                                                                 <  2:2002313/148‑1 (MQ=255)
gATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGa                                                                                                                                                 <  2:176837/148‑1 (MQ=255)
  ttATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATg                                                                                                                                               <  1:2358047/148‑1 (MQ=255)
   tataTGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGt                                                                                                                                              <  2:3172134/148‑1 (MQ=255)
    atatGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTc                                                                                                                                             <  1:1031566/148‑1 (MQ=255)
    atatGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTc                                                                                                                                             >  1:135327/1‑148 (MQ=255)
     tatgtatgCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGCTACACCGGTTTTAGTGTATTACCTGTCTGTTACGATCAATTCACTTGGCTCCTGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCt                                                                                                                                            <  2:3263639/148‑1 (MQ=255)
     tatGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCt                                                                                                                                            <  1:2635884/148‑1 (MQ=255)
      atGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCtt                                                                                                                                           >  1:2602317/1‑148 (MQ=255)
      atGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCtt                                                                                                                                           >  2:773237/1‑148 (MQ=255)
        gAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCg                                                                                                                                         <  2:1892278/148‑1 (MQ=255)
           tGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTg                                                                                                                                      <  2:157534/148‑1 (MQ=255)
           tGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTg                                                                                                                                      <  1:1922400/148‑1 (MQ=255)
            gCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGc                                                                                                                                     <  2:1347859/148‑1 (MQ=255)
             cTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGcc                                                                                                                                    >  2:652618/1‑148 (MQ=255)
                caacaTACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAAt                                                                                                                                 <  2:373774/148‑1 (MQ=255)
                caacaTACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAAt                                                                                                                                 <  1:2206661/148‑1 (MQ=255)
                  acaTACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAAttt                                                                                                                               >  1:193691/1‑148 (MQ=255)
                   caTACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTc                                                                                                                              >  2:1679807/1‑148 (MQ=255)
                   caTACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTc                                                                                                                              <  1:909401/148‑1 (MQ=255)
                       cccGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCgtg                                                                                                                          <  2:1091910/148‑1 (MQ=255)
                         cGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGtgtt                                                                                                                        >  1:2452080/1‑148 (MQ=255)
                         cGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGtgtt                                                                                                                        <  1:660396/148‑1 (MQ=255)
                          gCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCACTTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGtgttg                                                                                                                       <  2:1875911/148‑1 (MQ=255)
                           cAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGtgttgt                                                                                                                      <  1:1075735/148‑1 (MQ=255)
                           cAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGtgttgt                                                                                                                      >  1:1651818/1‑148 (MQ=255)
                              gTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAgg                                                                                                                   >  2:266854/1‑148 (MQ=255)
                              gTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAgg                                                                                                                   >  2:2786401/1‑148 (MQ=255)
                               tGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAggg                                                                                                                  >  1:3552078/1‑148 (MQ=255)
                                 aTTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCa                                                                                                                >  1:853924/1‑148 (MQ=255)
                                   tCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATa                                                                                                              <  1:3184230/148‑1 (MQ=255)
                                    ccTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATaa                                                                                                             <  1:50986/148‑1 (MQ=255)
                                        cACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACct                                                                                                         <  2:2464903/148‑1 (MQ=255)
                                              tACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTggg                                                                                                   >  2:3518775/1‑148 (MQ=255)
                                                cAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGc                                                                                                 <  1:1784579/148‑1 (MQ=255)
                                                 aGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGcc                                                                                                <  1:2619829/148‑1 (MQ=255)
                                                    tACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGa                                                                                             <  1:2803364/148‑1 (MQ=255)
                                                       acCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACg                                                                                          <  1:3416451/148‑1 (MQ=255)
                                                        ccGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACgg                                                                                         >  1:3325483/1‑148 (MQ=255)
                                                        ccGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACgg                                                                                         >  1:1153192/1‑148 (MQ=255)
                                                         cGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGc                                                                                        >  2:569594/1‑148 (MQ=255)
                                                                  tgtATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGc                                                                               >  2:905558/1‑148 (MQ=255)
                                                                  tgtATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGc                                                                               <  2:1535387/148‑1 (MQ=255)
                                                                  tgtATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGc                                                                               >  1:1963217/1‑148 (MQ=255)
                                                                     aTTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAAt                                                                            <  1:980542/148‑1 (MQ=255)
                                                                     aTTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAAt                                                                            <  1:2756348/148‑1 (MQ=255)
                                                                      ttACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAAtt                                                                           >  2:2042221/1‑148 (MQ=255)
                                                                      ttACCTGTCTGGTACGATCAAGTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAAtt                                                                           <  1:3124433/148‑1 (MQ=255)
                                                                       tACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTc                                                                          >  1:1730799/1‑148 (MQ=255)
                                                                       tACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTc                                                                          <  2:3439631/148‑1 (MQ=255)
                                                                       tACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTc                                                                          >  2:2330265/1‑148 (MQ=255)
                                                                         ccTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCa                                                                        >  1:162627/1‑148 (MQ=255)
                                                                          cTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAg                                                                       <  1:2276702/148‑1 (MQ=255)
                                                                           tGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGt                                                                      >  1:3208769/1‑148 (MQ=255)
                                                                           tGTCTGGTACGATCAATTCACTTAGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGt                                                                      <  1:2158717/148‑1 (MQ=255)
                                                                            ttCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATAGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTc                                                                     >  2:2218391/2‑148 (MQ=255)
                                                                              cTGGTACGATCAATTCACTTGGCGCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAg                                                                   >  1:2393730/1‑148 (MQ=255)
                                                                                ggTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGgt                                                                 <  1:2042221/148‑1 (MQ=255)
                                                                                   aCGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTa                                                              >  2:1232062/1‑148 (MQ=255)
                                                                                   aCGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTa                                                              >  2:2489445/1‑148 (MQ=255)
                                                                                    cGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTaa                                                             >  1:569337/1‑148 (MQ=255)
                                                                                     gATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAac                                                            >  1:2532658/1‑148 (MQ=255)
                                                                                      aTCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAaca                                                           >  2:525282/1‑148 (MQ=255)
                                                                                      aTCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAaca                                                           >  2:3502862/1‑148 (MQ=255)
                                                                                        cAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACg                                                         <  1:533841/148‑1 (MQ=255)
                                                                                           ttCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGa                                                      >  2:2479699/1‑148 (MQ=255)
                                                                                            tCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGaa                                                     <  1:2355428/148‑1 (MQ=255)
                                                                                              aCTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTATTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTGGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAAc                                                   >  2:829303/1‑148 (MQ=255)
                                                                                                 tGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAAcgcg                                                >  1:771276/1‑148 (MQ=255)
                                                                                                 tGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAAcgcg                                                >  1:844623/1‑148 (MQ=255)
                                                                                                  ggCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCgg                                               <  1:2575377/148‑1 (MQ=255)
                                                                                                      cATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGAtt                                           >  1:233425/1‑148 (MQ=255)
                                                                                                       aTGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTa                                          >  1:371851/1‑148 (MQ=255)
                                                                                                        tGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAg                                         >  1:1597160/1‑148 (MQ=255)
                                                                                                         gTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGt                                        >  2:1932873/1‑148 (MQ=255)
                                                                                                         gTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGt                                        >  1:3120142/1‑148 (MQ=255)
                                                                                                            ccGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTaa                                     <  1:780641/148‑1 (MQ=255)
                                                                                                            ccGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTaa                                     <  1:3114149/148‑1 (MQ=255)
                                                                                                            ccGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATACCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTaa                                     <  1:1221725/148‑1 (MQ=255)
                                                                                                               tCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGACACGCGGGATTAGTTAATGt                                  >  2:1816372/1‑148 (MQ=255)
                                                                                                               tCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGt                                  <  1:3238702/148‑1 (MQ=255)
                                                                                                                  aTGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAACTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCaa                               <  1:1630662/148‑1 (MQ=255)
                                                                                                                   tGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAAc                              >  2:1489639/1‑148 (MQ=255)
                                                                                                                       ggAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGa                          >  1:642115/1‑148 (MQ=255)
                                                                                                                        gagagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAg                         <  2:2107838/148‑1 (MQ=255)
                                                                                                                         agagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGt                        >  2:1254893/1‑148 (MQ=255)
                                                                                                                         agagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGt                        <  2:3195394/148‑1 (MQ=255)
                                                                                                                           agaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGAGTTGTATGGCATAAACCTCTTGGGGCCAGAACGGCAGAACATGCAATTCCAGTCAGGTGTAACACCGGAAACGCGGGATTAGTTAAAGTCAACCTTAGTCt                      >  1:817758/1‑148 (MQ=255)
                                                                                                                             aCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGa                    <  2:431810/148‑1 (MQ=255)
                                                                                                                              cAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGAc                   >  1:2807962/1‑148 (MQ=255)
                                                                                                                               aaTTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACg                  >  1:2936002/1‑148 (MQ=255)
                                                                                                                                aTTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACgg                 >  2:2236856/1‑148 (MQ=255)
                                                                                                                                 ttGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGc                >  1:328040/1‑148 (MQ=255)
                                                                                                                                  tGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCa               >  1:2667679/1‑148 (MQ=255)
                                                                                                                                    aaGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAAc             >  1:139798/1‑148 (MQ=255)
                                                                                                                                     aGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAAcc            >  1:799047/1‑148 (MQ=255)
                                                                                                                                          aTTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGa       >  1:3526640/1‑148 (MQ=255)
                                                                                                                                           tttAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGAc      >  1:2063881/1‑148 (MQ=255)
                                                                                                                                            ttAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGAcc     <  2:2844583/148‑1 (MQ=255)
                                                                                                                                               ggTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCa  >  2:2920610/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCA  >  CP009273/997819‑998109

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: