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,3160GT99.6% 825.7 / ‑4.4 252R451S (CGT→AGT) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (123/128);  total (123/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 258 total aligned reads displayed.

GGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCAC  >  CP009273/1038169‑1038462
                                                                                                                                                   |                                                                                                                                                  
  atatTTTTAGTACGATGGCGCTGTTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGc                                                                                                                                                <  1:756359/150‑1 (MQ=255)
    atTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCaa                                                                                                                                              >  2:3833615/1‑150 (MQ=255)
    atTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCaa                                                                                                                                              >  1:4959859/1‑150 (MQ=255)
      ttttAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACa                                                                                                                                            <  1:1941096/150‑1 (MQ=255)
       tttAGTACGATNGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCAAGTTGTTCCGGGGCTTCTACACCACTACGCAACAt                                                                                                                                           >  1:2758289/1‑150 (MQ=255)
       tttAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACAt                                                                                                                                           <  2:3546452/150‑1 (MQ=255)
         tAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCg                                                                                                                                         >  2:3575554/1‑150 (MQ=255)
           gTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCa                                                                                                                                       >  2:4612196/1‑150 (MQ=255)
              cGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCAcgc                                                                                                                                    <  2:2650989/150‑1 (MQ=255)
                 tGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGcc                                                                                                                                 >  2:1720486/1‑150 (MQ=255)
                    cgcTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGc                                                                                                                              >  1:3935504/1‑150 (MQ=255)
                     gctgctGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCa                                                                                                                             <  2:2214715/150‑1 (MQ=255)
                       tgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACg                                                                                                                           <  1:5031597/150‑1 (MQ=255)
                          tgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCa                                                                                                                        >  1:4742288/1‑150 (MQ=255)
                           gATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCAc                                                                                                                       >  1:805252/1‑150 (MQ=255)
                           gATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCAc                                                                                                                       >  2:972990/1‑150 (MQ=255)
                             ttAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCAccc                                                                                                                     >  1:102949/1‑150 (MQ=255)
                             ttAGAAAATAAATTGTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCAccc                                                                                                                     <  2:4826638/150‑1 (MQ=255)
                              tAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCca                                                                                                                    <  1:4107031/150‑1 (MQ=255)
                                gAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCcaca                                                                                                                  <  2:3378181/150‑1 (MQ=255)
                                gAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCcaca                                                                                                                  >  2:1990927/1‑150 (MQ=255)
                                gAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCcaca                                                                                                                  >  1:1804664/1‑150 (MQ=255)
                                 aaaaTAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAg                                                                                                                 <  1:4504993/150‑1 (MQ=255)
                                   aataaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGaa                                                                                                               >  1:1248016/1‑150 (MQ=255)
                                    ataaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGaaa                                                                                                              >  1:1653979/1‑150 (MQ=255)
                                     taaatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGAGCCAGCACGGCACCCACAGAAAt                                                                                                             <  1:2917823/150‑1 (MQ=255)
                                       aatTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATaa                                                                                                           >  2:3301516/1‑150 (MQ=255)
                                        atTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATaaa                                                                                                          >  2:62523/1‑150 (MQ=255)
                                        atTTTTCTTAAGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATaaa                                                                                                          >  1:102092/1‑150 (MQ=255)
                                         tttttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAAc                                                                                                         <  2:1197897/150‑1 (MQ=255)
                                         tttttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAAc                                                                                                         <  2:1535960/150‑1 (MQ=255)
                                          ttttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACa                                                                                                        >  2:454905/1‑150 (MQ=255)
                                          ttttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACa                                                                                                        >  1:5268036/1‑150 (MQ=255)
                                          ttttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACa                                                                                                        <  2:1845545/150‑1 (MQ=255)
                                           tttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAg                                                                                                       <  2:3618697/150‑1 (MQ=255)
                                             tCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGc                                                                                                     <  1:4931061/150‑1 (MQ=255)
                                             tCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGc                                                                                                     >  2:2022159/1‑150 (MQ=255)
                                              cTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGcc                                                                                                    >  1:3516847/1‑150 (MQ=255)
                                              cTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGcc                                                                                                    >  1:650463/1‑150 (MQ=255)
                                              cTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGcc                                                                                                    <  1:2794903/150‑1 (MQ=255)
                                               ttACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCa                                                                                                   <  1:52228/150‑1 (MQ=255)
                                                     aGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATa                                                                                             <  1:3096310/150‑1 (MQ=255)
                                                     aGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATa                                                                                             <  2:403194/150‑1 (MQ=255)
                                                      gACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCTCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATaa                                                                                            >  1:1457300/1‑150 (MQ=255)
                                                      gACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCAGGGGCTTCTACACCTCTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATaa                                                                                            <  2:2980295/150‑1 (MQ=255)
                                                       aCGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATaaa                                                                                           <  2:3243711/150‑1 (MQ=255)
                                                        cGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATaaaa                                                                                          >  2:1332825/1‑150 (MQ=255)
                                                          ggTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAAcc                                                                                        >  1:4171643/1‑150 (MQ=255)
                                                             gCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAg                                                                                     <  2:2305568/150‑1 (MQ=255)
                                                             gCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAg                                                                                     <  1:4335087/150‑1 (MQ=255)
                                                              cGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGc                                                                                    >  2:5165751/1‑150 (MQ=255)
                                                              cGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGc                                                                                    >  1:1243202/1‑150 (MQ=255)
                                                               gTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGca                                                                                   >  2:4155214/1‑150 (MQ=255)
                                                               gTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGca                                                                                   >  2:4011731/1‑150 (MQ=255)
                                                                  tATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGcacca                                                                                >  2:4251217/1‑150 (MQ=255)
                                                                   aTCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGcaccac                                                                               <  2:1292037/150‑1 (MQ=255)
                                                                   aTCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGcaccac                                                                               <  1:3159043/150‑1 (MQ=255)
                                                                    tCCAGCCTCTCGGTCATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAAAATCGCACGCGCAAGCAAGTCACCAACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACg                                                                              <  1:3771285/150‑1 (MQ=255)
                                                                       aGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTc                                                                           >  1:286921/1‑150 (MQ=255)
                                                                       aGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTc                                                                           >  1:373395/1‑150 (MQ=255)
                                                                          cACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAAc                                                                        <  1:4579022/150‑1 (MQ=255)
                                                                          cACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAAc                                                                        <  2:3090745/150‑1 (MQ=255)
                                                                               ggACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCttt                                                                   <  1:305441/150‑1 (MQ=255)
                                                                               ggACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCttt                                                                   <  1:2560658/150‑1 (MQ=255)
                                                                               ggACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCatt                                                                   >  2:3033365/1‑150 (MQ=255)
                                                                                gACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtttt                                                                  <  1:2406642/150‑1 (MQ=255)
                                                                                  cATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtttttt                                                                >  2:2453784/1‑150 (MQ=255)
                                                                                      gggATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGt                                                            >  1:399342/1‑150 (MQ=255)
                                                                                        gATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGTGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGACAAGAATAAAACCAAGCACCACGATCAACGCTTTTTTCGGatt                                                          >  2:536272/1‑150 (MQ=255)
                                                                                          tAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCTCGTCGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCTACGCTTTTTTCGGTTTCa                                                        <  2:1588591/150‑1 (MQ=255)
                                                                                           aGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCAc                                                       <  1:3526137/150‑1 (MQ=255)
                                                                                              ggCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTgg                                                    <  1:2627320/150‑1 (MQ=255)
                                                                                                  tAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGc                                                <  1:686206/150‑1 (MQ=255)
                                                                                                   aaaCGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCg                                               <  1:1592480/150‑1 (MQ=255)
                                                                                                       gCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTg                                           >  1:4790066/1‑150 (MQ=255)
                                                                                                        cTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGa                                          <  2:5005428/150‑1 (MQ=255)
                                                                                                         tGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAg                                         >  1:1636778/1‑150 (MQ=255)
                                                                                                          gATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGt                                        <  1:3213518/150‑1 (MQ=255)
                                                                                                             gCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGAc                                     <  2:4959943/150‑1 (MQ=255)
                                                                                                              ccGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACt                                    >  2:3193166/1‑150 (MQ=255)
                                                                                                               cGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTg                                   >  1:4995511/1‑150 (MQ=255)
                                                                                                                 tGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGcc                                 >  2:5258954/1‑150 (MQ=255)
                                                                                                                      ttCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTc                            <  2:785776/150‑1 (MQ=255)
                                                                                                                         cAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGAt                         <  2:4466022/150‑1 (MQ=255)
                                                                                                                         cAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGAt                         <  2:1429601/150‑1 (MQ=255)
                                                                                                                           gttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATaa                       >  1:1404009/1‑150 (MQ=255)
                                                                                                                            ttgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAAt                      <  2:4441736/150‑1 (MQ=255)
                                                                                                                             tgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATc                     <  1:2163353/150‑1 (MQ=255)
                                                                                                                             tgttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATc                     <  2:4958942/150‑1 (MQ=255)
                                                                                                                              gttCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATcc                    <  2:886931/150‑1 (MQ=255)
                                                                                                                              gttCCGGGGATTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATACAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATcc                    <  2:1051792/150‑1 (MQ=255)
                                                                                                                                tCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGc                  >  2:593788/1‑150 (MQ=255)
                                                                                                                                     ggCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGtt             <  2:4205031/150‑1 (MQ=255)
                                                                                                                                      gCTTCTACACCACTACGCAACATCGCACGCTCCAGCACGGCACCCACAGAAATAAACATGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTTCCTGGTCGATAATCCGCACGTTa            >  2:1386553/1‑150 (MQ=255)
                                                                                                                                      gCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTa            <  1:3301516/150‑1 (MQ=255)
                                                                                                                                             caccacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATgg     <  2:2698223/150‑1 (MQ=255)
                                                                                                                                             caccacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATgg     <  1:4347979/150‑1 (MQ=255)
                                                                                                                                              accacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGc    >  1:2083104/1‑150 (MQ=255)
                                                                                                                                                cacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCAc  <  1:2852264/150‑1 (MQ=255)
                                                                                                                                                cacTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCAc  >  2:1871193/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCAC  >  CP009273/1038169‑1038462

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: