Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,038,316 G→T 6.3% 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,3160GT 6.3% 906.6 / 24.8 236R451S (CGT→AGT) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base G (98/123);  new base T (9/6);  total (107/129)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 238 total aligned reads displayed.

TATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTG  >  CP009273/1038172‑1038464
                                                                                                                                                |                                                                                                                                                    
tatTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCa                                                                                                                                                 >  2:1600165/1‑150 (MQ=255)
    tttAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACAt                                                                                                                                             <  2:2175106/150‑1 (MQ=255)
           cGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCAcgc                                                                                                                                      >  2:1797235/1‑150 (MQ=255)
            gATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCAcgcg                                                                                                                                     >  2:2534469/1‑150 (MQ=255)
             aTGGCGCTGCAGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCAGAAACGCTGATGCCGTGCTCTTCCAGTTGTTACGGGGCTTCTACACCACTACGCAACATCGCACgcgc                                                                                                                                    >  1:3660781/1‑150 (MQ=255)
              tGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGcc                                                                                                                                   >  1:1072537/1‑150 (MQ=255)
              tGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGcc                                                                                                                                   >  1:352266/1‑150 (MQ=255)
                   ctgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCAc                                                                                                                              <  2:650299/150‑1 (MQ=255)
                   ctgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCAc                                                                                                                              <  1:3516238/150‑1 (MQ=255)
                    tgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACg                                                                                                                             >  1:2488149/1‑150 (MQ=255)
                    tgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACg                                                                                                                             >  2:4071581/1‑150 (MQ=255)
                    tgctgATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACg                                                                                                                             >  2:3269903/1‑150 (MQ=255)
                        gATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCAc                                                                                                                         >  2:4634435/1‑150 (MQ=255)
                          ttAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCAccc                                                                                                                       >  1:918115/1‑150 (MQ=255)
                            aGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCcac                                                                                                                     >  1:25246/1‑150 (MQ=255)
                              aaaaTAAATTTTTCTTACGCAGTCGGGTGCGTTTATCGAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCGTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAg                                                                                                                   <  2:4620771/150‑1 (MQ=255)
                                       ttttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACa                                                                                                          >  2:4543098/1‑150 (MQ=255)
                                       ttttCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACa                                                                                                          >  1:68115/1‑150 (MQ=255)
                                            ttACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCa                                                                                                     <  2:2042912/150‑1 (MQ=255)
                                             tACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCaa                                                                                                    >  1:1045275/1‑150 (MQ=255)
                                              aCGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAg                                                                                                   >  1:852498/1‑150 (MQ=255)
                                               cGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGa                                                                                                  >  1:2926552/1‑150 (MQ=255)
                                               cGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGa                                                                                                  >  2:539453/1‑150 (MQ=255)
                                                gCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGaa                                                                                                 <  2:674835/150‑1 (MQ=255)
                                                 cAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAAt                                                                                                <  2:4583845/150‑1 (MQ=255)
                                                 cAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAAt                                                                                                <  1:1198806/150‑1 (MQ=255)
                                                 cAGACGGGTGCGTTTATCCAGCCACTCGGACAGTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAAt                                                                                                >  2:1749978/1‑150 (MQ=255)
                                                  aGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCGTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATa                                                                                               >  2:1160124/1‑150 (MQ=255)
                                                      gggTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAAc                                                                                           >  1:2991778/1‑150 (MQ=255)
                                                       ggTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAAcc                                                                                          <  1:2716370/150‑1 (MQ=255)
                                                        gTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCa                                                                                         <  1:4380283/150‑1 (MQ=255)
                                                        gTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCa                                                                                         >  1:2667406/1‑150 (MQ=255)
                                                        gTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCa                                                                                         >  1:3240543/1‑150 (MQ=255)
                                                        gTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCa                                                                                         <  1:1880935/150‑1 (MQ=255)
                                                         tGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCaa                                                                                        <  1:2840517/150‑1 (MQ=255)
                                                         tGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCaa                                                                                        <  2:228218/150‑1 (MQ=255)
                                                         tGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAAAAAAACCaa                                                                                        >  2:4199291/1‑150 (MQ=255)
                                                          gCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAg                                                                                       >  1:1143499/1‑150 (MQ=255)
                                                          gCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAg                                                                                       <  2:1369465/150‑1 (MQ=255)
                                                                    aGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTc                                                                             <  1:1843538/150‑1 (MQ=255)
                                                                     gCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCa                                                                            <  1:1197841/150‑1 (MQ=255)
                                                                     gCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCa                                                                            <  1:2198667/150‑1 (MQ=255)
                                                                       cACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAAc                                                                          <  2:4299884/150‑1 (MQ=255)
                                                                       cACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAAc                                                                          <  1:1681297/150‑1 (MQ=255)
                                                                        aCTCGGACATTGGTATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGATTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACg                                                                         <  1:524376/150‑1 (MQ=255)
                                                                         cTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGc                                                                        >  1:2019129/1‑150 (MQ=255)
                                                                         cTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGc                                                                        >  1:821516/1‑150 (MQ=255)
                                                                          tCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCATGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCt                                                                       <  1:3757566/150‑1 (MQ=255)
                                                                          tCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCt                                                                       <  1:3934343/150‑1 (MQ=255)
                                                                          tCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCt                                                                       <  2:2676377/150‑1 (MQ=255)
                                                                          tCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCt                                                                       <  2:3363060/150‑1 (MQ=255)
                                                                           cGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtt                                                                      <  1:3627790/150‑1 (MQ=255)
                                                                           cGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtt                                                                      <  1:4596742/150‑1 (MQ=255)
                                                                              aCATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCttttt                                                                   >  1:569900/1‑150 (MQ=255)
                                                                               cATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCtttttt                                                                  <  1:3831145/150‑1 (MQ=255)
                                                                                   gggATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGt                                                              >  2:4011718/1‑150 (MQ=255)
                                                                                      aTAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTc                                                           <  2:1430309/150‑1 (MQ=255)
                                                                                       tAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCa                                                          <  1:2703975/150‑1 (MQ=255)
                                                                                        aGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCAc                                                         >  2:3810811/1‑150 (MQ=255)
                                                                                        aGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCAc                                                         >  2:2890520/1‑150 (MQ=255)
                                                                                        aGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCAc                                                         >  2:3524417/1‑150 (MQ=255)
                                                                                             cATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCt                                                    <  2:3742861/150‑1 (MQ=255)
                                                                                             cATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGACACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCt                                                    <  1:3410797/150‑1 (MQ=255)
                                                                                              aTAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTg                                                   >  2:3570043/1‑150 (MQ=255)
                                                                                                aaaCGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCg                                                 <  2:1977464/150‑1 (MQ=255)
                                                                                                 aaCGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCgg                                                <  2:3731637/150‑1 (MQ=255)
                                                                                                  aCGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGc                                               <  2:1935752/150‑1 (MQ=255)
                                                                                                  aCGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGc                                               <  2:4002831/150‑1 (MQ=255)
                                                                                                   cGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCt                                              <  2:2194571/150‑1 (MQ=255)
                                                                                                   cGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCt                                              >  1:4418040/1‑150 (MQ=255)
                                                                                                      tGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCATTTTTCGGTTTCACTGGCTGCGGCTGAg                                           >  1:1965727/1‑150 (MQ=255)
                                                                                                           ccGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACt                                      <  1:4212677/150‑1 (MQ=255)
                                                                                                           ccGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACt                                      <  1:2082665/150‑1 (MQ=255)
                                                                                                              tGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGcc                                   <  1:1518243/150‑1 (MQ=255)
                                                                                                               gCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCg                                  >  1:4534098/1‑150 (MQ=255)
                                                                                                               gCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCg                                  >  1:1413461/1‑150 (MQ=255)
                                                                                                                ctctTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCgg                                 >  1:2614297/1‑150 (MQ=255)
                                                                                                                ctctTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCgg                                 <  1:2162890/150‑1 (MQ=255)
                                                                                                                 tctTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCggg                                <  1:3960798/150‑1 (MQ=255)
                                                                                                                 tctTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCggg                                <  1:1874459/150‑1 (MQ=255)
                                                                                                                 tctTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCggg                                <  1:3046874/150‑1 (MQ=255)
                                                                                                                  ctTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGt                               <  1:2372948/150‑1 (MQ=255)
                                                                                                                  ctTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGt                               <  1:1307043/150‑1 (MQ=255)
                                                                                                                   ttCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTc                              <  2:1579566/150‑1 (MQ=255)
                                                                                                                   ttCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTc                              <  2:3734189/150‑1 (MQ=255)
                                                                                                                   ttCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTc                              <  2:1579267/150‑1 (MQ=255)
                                                                                                                   ttCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTc                              >  1:751840/1‑150 (MQ=255)
                                                                                                                    tCCAGTTGTTCCGGGGCTTCTACACCACTACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCg                             <  2:4074044/150‑1 (MQ=255)
                                                                                                                     ccAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGa                            >  2:3773115/1‑150 (MQ=255)
                                                                                                                        gttgttCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATaa                         >  1:2713688/1‑150 (MQ=255)
                                                                                                                         ttgttCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAAt                        <  1:1608269/150‑1 (MQ=255)
                                                                                                                            ttCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCg                     <  1:636292/150‑1 (MQ=255)
                                                                                                                                      tCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCa           <  1:3342284/150‑1 (MQ=255)
                                                                                                                                        tACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAAt         <  2:3172456/150‑1 (MQ=255)
                                                                                                                                        tACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAAt         <  2:314225/150‑1 (MQ=255)
                                                                                                                                        tACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAAt         <  2:2837091/150‑1 (MQ=255)
                                                                                                                                        tACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAAt         <  2:1172791/150‑1 (MQ=255)
                                                                                                                                        tACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAAt         <  1:1371197/150‑1 (MQ=255)
                                                                                                                                         acacCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATg        <  2:1309189/150‑1 (MQ=255)
                                                                                                                                         acacCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATg        <  1:3625338/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTG  >  CP009273/1038172‑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: