Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,038,316 G→T 85.7% 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,3160GT85.7% 677.2 / 123.2 286R451S (CGT→AGT) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base G (22/19);  new base T (143/102);  total (165/121)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 295 total aligned reads displayed.

GATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTG  >  CP009273/1038170‑1038464
                                                                                                                                                  |                                                                                                                                                    
gagatTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCactac                                                                                                                                                    <  2:4582807/146‑1 (MQ=255)
    tttttAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCaa                                                                                                                                                >  1:2217690/1‑149 (MQ=255)
       ttAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACAt                                                                                                                                             <  2:4391385/149‑1 (MQ=255)
        tAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATc                                                                                                                                            <  2:508429/149‑1 (MQ=255)
         aGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCg                                                                                                                                           >  1:408105/1‑149 (MQ=255)
         aGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCg                                                                                                                                           >  1:1277772/1‑149 (MQ=255)
           tACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCa                                                                                                                                         <  2:1584765/149‑1 (MQ=255)
               aTGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCAcgcg                                                                                                                                     >  1:3420898/1‑149 (MQ=255)
                 ggCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGcc                                                                                                                                   <  1:166772/149‑1 (MQ=255)
                    gctgctGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGc                                                                                                                                >  2:1980697/1‑149 (MQ=255)
                    gctgctGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGc                                                                                                                                >  1:1858079/1‑149 (MQ=255)
                    gctgctGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGc                                                                                                                                <  2:1372154/149‑1 (MQ=255)
                    gctgctGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGc                                                                                                                                <  2:3760575/149‑1 (MQ=255)
                     ctgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCa                                                                                                                               <  2:256845/149‑1 (MQ=255)
                     ctgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCa                                                                                                                               <  2:926000/149‑1 (MQ=255)
                      tgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCAc                                                                                                                              >  1:2788992/1‑149 (MQ=255)
                       gctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACg                                                                                                                             >  2:3898415/1‑149 (MQ=255)
                           aTTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCAc                                                                                                                         >  2:4116104/1‑149 (MQ=255)
                                   ataaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGaa                                                                                                                 >  2:1033735/1‑149 (MQ=255)
                                    taaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGaaa                                                                                                                >  1:3835523/1‑149 (MQ=255)
                                     aaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAAt                                                                                                               >  2:83131/1‑149 (MQ=255)
                                     aaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAAt                                                                                                               >  2:3337969/1‑149 (MQ=255)
                                        tttttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATaaa                                                                                                            <  2:2370920/149‑1 (MQ=255)
                                          tttCTTAGGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACa                                                                                                          >  1:2758462/1‑149 (MQ=255)
                                          tttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACa                                                                                                          <  1:1420172/149‑1 (MQ=255)
                                          tttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACa                                                                                                          >  1:1750820/1‑149 (MQ=255)
                                           ttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAg                                                                                                         <  1:2183403/149‑1 (MQ=255)
                                             cTTACGCCGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGAGAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCGTCTACACCACTCCGCAACATCGCACGCGCCAGCACGGCACCCACAGACATAATCAGGc                                                                                                       >  2:3041153/1‑149 (MQ=255)
                                              ttACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGcc                                                                                                      >  2:2420877/1‑149 (MQ=255)
                                               tACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCa                                                                                                     <  2:1784375/149‑1 (MQ=255)
                                                aCGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCaa                                                                                                    <  1:886601/149‑1 (MQ=255)
                                                  gCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGa                                                                                                  <  2:1669336/149‑1 (MQ=255)
                                                    aGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAAt                                                                                                <  1:2033636/149‑1 (MQ=255)
                                                     gACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATa                                                                                               >  1:1573157/1‑149 (MQ=255)
                                                     gACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATa                                                                                               >  2:3938204/1‑149 (MQ=255)
                                                     gACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCACGAATa                                                                                               >  2:3266196/1‑149 (MQ=255)
                                                       cGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTATACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATaaa                                                                                             <  2:4084776/149‑1 (MQ=255)
                                                         ggTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAAc                                                                                           <  1:2951308/149‑1 (MQ=255)
                                                         ggTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAAc                                                                                           >  2:4327431/1‑149 (MQ=255)
                                                          gTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAAcc                                                                                          >  1:3074118/1‑149 (MQ=255)
                                                           tGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCa                                                                                         >  2:1541579/1‑149 (MQ=255)
                                                            gCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCaa                                                                                        <  2:2323023/149‑1 (MQ=255)
                                                            gCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCaa                                                                                        >  2:2682360/1‑149 (MQ=255)
                                                              gTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGc                                                                                      >  2:231344/1‑149 (MQ=255)
                                                              gTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGc                                                                                      <  1:1088149/149‑1 (MQ=255)
                                                              gTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGc                                                                                      <  2:3551891/149‑1 (MQ=255)
                                                               tttATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGca                                                                                     >  2:4274457/1‑149 (MQ=255)
                                                                ttATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGcac                                                                                    <  1:4773266/149‑1 (MQ=255)
                                                                   tCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGcaccac                                                                                 >  2:1140062/1‑149 (MQ=255)
                                                                      aGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGtt                                                                              >  1:1392483/1‑149 (MQ=255)
                                                                      aGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGtt                                                                              <  2:2010698/149‑1 (MQ=255)
                                                                       gCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTc                                                                             >  1:3701991/1‑149 (MQ=255)
                                                                       gCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTc                                                                             <  2:741822/149‑1 (MQ=255)
                                                                       gCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTc                                                                             >  2:4521879/1‑149 (MQ=255)
                                                                       gCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTc                                                                             <  2:2129309/149‑1 (MQ=255)
                                                                          aCTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAAc                                                                          >  1:2725564/1‑149 (MQ=255)
                                                                          aCTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAAc                                                                          >  1:1825749/1‑149 (MQ=255)
                                                                           cTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACg                                                                         >  2:118043/1‑149 (MQ=255)
                                                                            tCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGc                                                                        <  1:4126541/149‑1 (MQ=255)
                                                                            tCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGc                                                                        <  1:1946968/149‑1 (MQ=255)
                                                                              ggACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtt                                                                      >  2:1116808/1‑149 (MQ=255)
                                                                              ggACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtt                                                                      <  2:653852/149‑1 (MQ=255)
                                                                                 cATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCttttt                                                                   >  1:756656/1‑149 (MQ=255)
                                                                                   ttGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTc                                                                 >  1:3784668/1‑149 (MQ=255)
                                                                                    tGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCg                                                                <  1:986401/149‑1 (MQ=255)
                                                                                     ggaaTAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCgg                                                               <  1:2415836/146‑1 (MQ=255)
                                                                                     ggaaTAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCgg                                                               <  1:1726414/146‑1 (MQ=255)
                                                                                                aTAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCt                                                    <  1:4229761/149‑1 (MQ=255)
                                                                                                 tAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTg                                                   >  1:1345371/1‑149 (MQ=255)
                                                                                                  aaaCGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGc                                                  <  1:3106813/149‑1 (MQ=255)
                                                                                                      gCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCt                                              >  2:4697640/1‑149 (MQ=255)
                                                                                                        tGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGa                                            <  1:3027673/149‑1 (MQ=255)
                                                                                                          aTGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGt                                          >  1:2690536/1‑149 (MQ=255)
                                                                                                           tGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTg                                         >  2:590128/1‑149 (MQ=255)
                                                                                                               gTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTg                                     >  1:2303573/1‑149 (MQ=255)
                                                                                                               gTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGATCAACGCTTTTATCGGTTTCACTGGCTGCGGCTGAGTGACTg                                     >  1:3333828/1‑149 (MQ=255)
                                                                                                                tGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGc                                    >  2:2023027/1‑149 (MQ=255)
                                                                                                                    ctTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCggg                                >  1:1138237/1‑149 (MQ=255)
                                                                                                                     ttCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGt                               >  1:3472081/1‑149 (MQ=255)
                                                                                                                        cAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGa                            <  1:449506/149‑1 (MQ=255)
                                                                                                                         aGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGAt                           >  1:1017121/1‑149 (MQ=255)
                                                                                                                          gttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATa                          >  1:598634/1‑149 (MQ=255)
                                                                                                                              ttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATcc                      <  2:779082/149‑1 (MQ=255)
                                                                                                                              ttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATcc                      >  1:1199761/1‑149 (MQ=255)
                                                                                                                                 cGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCa                   <  1:4338773/149‑1 (MQ=255)
                                                                                                                                  ggggCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCAc                  >  1:2532004/1‑149 (MQ=255)
                                                                                                                                   gggCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACg                 >  1:4717163/1‑149 (MQ=255)
                                                                                                                                     gCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGtt               <  2:1277772/149‑1 (MQ=255)
                                                                                                                                     gCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGtt               <  1:711913/149‑1 (MQ=255)
                                                                                                                                     gCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGtt               <  2:4410024/149‑1 (MQ=255)
                                                                                                                                     gCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGtt               >  1:1492583/1‑149 (MQ=255)
                                                                                                                                      cTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTa              <  1:3198200/149‑1 (MQ=255)
                                                                                                                                       ttCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTAc             >  2:3156428/1‑149 (MQ=255)
                                                                                                                                       ttCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTAc             >  2:4462537/1‑149 (MQ=255)
                                                                                                                                         cTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCa           >  1:3255488/1‑149 (MQ=255)
                                                                                                                                          tACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCaa          <  2:3502988/149‑1 (MQ=255)
                                                                                                                                          tACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCaa          >  1:1037866/1‑149 (MQ=255)
                                                                                                                                           acacCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAAt         <  1:2525739/149‑1 (MQ=255)
                                                                                                                                              ccacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGc      >  1:82518/1‑149 (MQ=255)
                                                                                                                                                  gacgCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTg  >  2:491636/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTG  >  CP009273/1038170‑1038464

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: