Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,038,316 G→T 100% R451S (CGT→AGT)  etk ← tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,038,3160GT97.2% 552.7 / 6.1 179R451S (CGT→AGT) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base T (93/80);  total (99/80)
Fisher's exact test for biased strand distribution p-value = 6.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 179 total aligned reads displayed.

TTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGC  >  CP009273/1038175‑1038460
                                                                                                                                             |                                                                                                                                                
ttttAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCaa                                                                                                                                            >  1:1989340/1‑148 (MQ=255)
ttttAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCaa                                                                                                                                            >  1:1989322/1‑148 (MQ=255)
   tAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACAt                                                                                                                                         <  2:2369782/148‑1 (MQ=255)
   tAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACAt                                                                                                                                         <  1:714728/148‑1 (MQ=255)
   tAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTCGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACAACAACGCAACAt                                                                                                                                         >  2:3431035/1‑148 (MQ=255)
     gTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCg                                                                                                                                       <  2:2152369/148‑1 (MQ=255)
       aCGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCa                                                                                                                                     <  2:1465364/148‑1 (MQ=255)
       aCGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCa                                                                                                                                     <  2:1003470/148‑1 (MQ=255)
       aCGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCa                                                                                                                                     <  2:3235742/148‑1 (MQ=255)
          aTGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCAcgc                                                                                                                                  >  2:2359342/1‑148 (MQ=255)
            ggCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACgcgc                                                                                                                                >  1:1599892/1‑148 (MQ=255)
            ggCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACgcgc                                                                                                                                >  1:1365457/1‑148 (MQ=255)
            ggCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACgcgc                                                                                                                                >  1:3479462/1‑148 (MQ=255)
                 tgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCa                                                                                                                           >  2:3538182/1‑148 (MQ=255)
                 tgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCa                                                                                                                           >  2:1776470/1‑148 (MQ=255)
                    tgATTAGAAAATAAATTTGTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACgg                                                                                                                        <  2:2284668/148‑1 (MQ=255)
                      aTTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCa                                                                                                                      >  1:3349904/1‑148 (MQ=255)
                      aTTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCa                                                                                                                      <  1:589193/148‑1 (MQ=255)
                         aGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCAccc                                                                                                                   >  1:842758/1‑148 (MQ=255)
                            aaataaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCcaca                                                                                                                >  1:3289445/1‑148 (MQ=255)
                            aaataaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCcaca                                                                                                                >  1:1995866/1‑148 (MQ=255)
                              ataaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGa                                                                                                              >  1:1622161/1‑148 (MQ=255)
                               taaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGaa                                                                                                             >  1:3290077/1‑148 (MQ=255)
                                aaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACTCGCCAGCACGGCACCCACAGaaa                                                                                                            >  1:1870032/1‑148 (MQ=255)
                                aaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGaac                                                                                                            >  2:3294430/1‑147 (MQ=255)
                                aaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGaaa                                                                                                            >  1:1372105/1‑148 (MQ=255)
                                aaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGaaa                                                                                                            >  2:2516425/1‑148 (MQ=255)
                                aaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGaaa                                                                                                            >  2:1592518/1‑148 (MQ=255)
                                aaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGaaa                                                                                                            >  1:2753265/1‑148 (MQ=255)
                                  atTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATa                                                                                                          <  2:2793317/148‑1 (MQ=255)
                                    ttttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATaaa                                                                                                        >  2:1077764/1‑148 (MQ=255)
                                    ttttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATaaa                                                                                                        >  2:1077485/1‑148 (MQ=255)
                                          tACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGcc                                                                                                  >  2:1951435/1‑148 (MQ=255)
                                             gCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAg                                                                                               >  1:81121/1‑148 (MQ=255)
                                             gCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAg                                                                                               >  1:517279/1‑148 (MQ=255)
                                              cAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGa                                                                                              <  1:2879892/148‑1 (MQ=255)
                                                gACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAAt                                                                                            >  1:3139787/1‑148 (MQ=255)
                                                gACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAAt                                                                                            >  1:1420764/1‑148 (MQ=255)
                                                gACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAAt                                                                                            >  2:531642/1‑148 (MQ=255)
                                                 aCGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATa                                                                                           >  2:2612816/1‑148 (MQ=255)
                                                   gggTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATaaa                                                                                         >  2:1396301/1‑148 (MQ=255)
                                                    ggTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATaaaa                                                                                        <  2:672431/148‑1 (MQ=255)
                                                    ggTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATaaaa                                                                                        <  1:2531074/148‑1 (MQ=255)
                                                       gCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCa                                                                                     >  1:955887/1‑148 (MQ=255)
                                                         gTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAg                                                                                   <  1:1076951/148‑1 (MQ=255)
                                                         gTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAg                                                                                   <  2:894678/148‑1 (MQ=255)
                                                          tttATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGc                                                                                  >  2:2834682/1‑148 (MQ=255)
                                                          tttATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGc                                                                                  >  2:1642947/1‑148 (MQ=255)
                                                                cAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACg                                                                            <  1:1579717/148‑1 (MQ=255)
                                                                  gCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGtt                                                                          >  2:2436163/1‑148 (MQ=255)
                                                                  gCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGtt                                                                          <  2:3549102/148‑1 (MQ=255)
                                                                  gCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGtt                                                                          >  2:2436294/1‑148 (MQ=255)
                                                                     aCTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCaa                                                                       <  2:734098/148‑1 (MQ=255)
                                                                          gACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtt                                                                  <  1:1011245/148‑1 (MQ=255)
                                                                          gACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtt                                                                  <  1:2197932/148‑1 (MQ=255)
                                                                           aCATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCttt                                                                 >  1:3252299/1‑148 (MQ=255)
                                                                            cATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtttt                                                                <  2:1232048/148‑1 (MQ=255)
                                                                            cATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtttt                                                                <  2:2133168/148‑1 (MQ=255)
                                                                              ttGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtttttt                                                              <  2:384966/148‑1 (MQ=255)
                                                                               tGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTc                                                             <  2:757811/148‑1 (MQ=255)
                                                                                gggATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCg                                                            <  1:2923027/148‑1 (MQ=255)
                                                                                  gATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGt                                                          <  1:953624/148‑1 (MQ=255)
                                                                                       tGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCAc                                                     >  1:311138/1‑148 (MQ=255)
                                                                                       tGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCAc                                                     <  2:1201830/148‑1 (MQ=255)
                                                                                       tGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCAc                                                     <  2:1202022/148‑1 (MQ=255)
                                                                                       tGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCAc                                                     <  1:2735944/148‑1 (MQ=255)
                                                                                           aTAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGc                                                 <  2:3139787/148‑1 (MQ=255)
                                                                                                cGCTGATGCCGTGCTCTGCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCgg                                            <  2:1413504/148‑1 (MQ=255)
                                                                                                   tGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTg                                         <  2:3290077/148‑1 (MQ=255)
                                                                                                    gATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGa                                        >  1:2232322/1‑148 (MQ=255)
                                                                                                    gATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGa                                        >  1:2026869/1‑148 (MQ=255)
                                                                                                     aTGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAg                                       >  2:1273957/1‑148 (MQ=255)
                                                                                                     aTGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTAGGTTTCACTGGCTGCGGCTGAg                                       >  1:2998724/1‑148 (MQ=255)
                                                                                                            gCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGc                                <  2:3479462/148‑1 (MQ=255)
                                                                                                                 tCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGt                           <  1:1025368/148‑1 (MQ=255)
                                                                                                                    aGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGa                        <  1:1776470/148‑1 (MQ=255)
                                                                                                                    aGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGa                        <  1:3538182/148‑1 (MQ=255)
                                                                                                                     gttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGAt                       <  2:916229/148‑1 (MQ=255)
                                                                                                                     gttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGAt                       <  2:3504031/148‑1 (MQ=255)
                                                                                                                      ttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATa                      >  2:742690/1‑148 (MQ=255)
                                                                                                                      ttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATa                      >  1:754517/1‑148 (MQ=255)
                                                                                                                      ttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATa                      >  1:221053/1‑148 (MQ=255)
                                                                                                                      ttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATa                      <  1:3274806/148‑1 (MQ=255)
                                                                                                                        gttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAAt                    <  1:446631/148‑1 (MQ=255)
                                                                                                                        gttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAAt                    <  1:2359329/148‑1 (MQ=255)
                                                                                                                        gttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAAt                    <  1:2359386/148‑1 (MQ=255)
                                                                                                                         ttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATc                   >  1:3070211/1‑148 (MQ=255)
                                                                                                                            cGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGc                <  1:3403967/148‑1 (MQ=255)
                                                                                                                               ggCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACg             >  2:1800773/1‑148 (MQ=255)
                                                                                                                                gCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGTGGCTGAGTGACTGCCGGGTCGATAATCCGCACGt            <  1:1856096/148‑1 (MQ=255)
                                                                                                                                 cTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGtt           <  1:2377360/148‑1 (MQ=255)
                                                                                                                                  ttCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTa          <  1:118390/148‑1 (MQ=255)
                                                                                                                                   tCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTAc         >  1:1603458/1‑148 (MQ=255)
                                                                                                                                       caccacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAAt     >  1:2902600/1‑148 (MQ=255)
                                                                                                                                       caccacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAAt     >  2:278792/1‑148 (MQ=255)
                                                                                                                                       caccacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAAt     >  2:1625410/1‑148 (MQ=255)
                                                                                                                                        accacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATg    >  1:2166330/1‑148 (MQ=255)
                                                                                                                                        accacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATg    <  1:530601/148‑1 (MQ=255)
                                                                                                                                          cacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCa               >  1:2675345/1‑135 (MQ=255)
                                                                                                                                          cacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGc  >  2:1991966/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                
TTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGC  >  CP009273/1038175‑1038460

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: