Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,038,158 C→A 13.4% A503A (GCG→GCT)  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,1580CA13.4% 678.3 / 53.8 209A503A (GCG→GCT) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base C (82/99);  new base A (6/22);  total (88/121)
Fisher's exact test for biased strand distribution p-value = 2.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 216 total aligned reads displayed.

AGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCT  >  CP009273/1038014‑1038305
                                                                                                                                                |                                                                                                                                                   
agTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGcc                                                                                                                                                  <  2:1648851/148‑1 (MQ=255)
agTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGcc                                                                                                                                                  <  2:691394/148‑1 (MQ=255)
 gTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCa                                                                                                                                                 <  2:1060588/148‑1 (MQ=255)
 gTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCa                                                                                                                                                 <  1:802157/148‑1 (MQ=255)
     aaCTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGaa                                                                                                                                             <  2:334652/148‑1 (MQ=255)
      aCTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGacg                                                                                                                                            >  1:2698627/1‑146 (MQ=255)
      aCTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAg                                                                                                                                            <  2:1428863/148‑1 (MQ=255)
      aCTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAg                                                                                                                                            <  1:3424185/148‑1 (MQ=255)
      aCTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAg                                                                                                                                            <  1:917871/148‑1 (MQ=255)
            aaaCGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGata                                                                                                                                      <  2:521720/148‑1 (MQ=255)
             aaCGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGatat                                                                                                                                     >  2:2531104/1‑148 (MQ=255)
             aaCGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGatat                                                                                                                                     >  1:2811437/1‑148 (MQ=255)
             aaCGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGatat                                                                                                                                     >  1:2607088/1‑148 (MQ=255)
              aCGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATAtt                                                                                                                                    >  1:916184/1‑148 (MQ=255)
              aCGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATAtt                                                                                                                                    >  2:1145065/1‑148 (MQ=255)
               cGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATAttt                                                                                                                                   >  1:3498697/1‑148 (MQ=255)
                  tttACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTa                                                                                                                                <  1:1784248/148‑1 (MQ=255)
                  tttACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATAATAGCGAAATGCAGACTGGTTCGTAGCGCACTTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTa                                                                                                                                <  1:832695/148‑1 (MQ=255)
                       cACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACg                                                                                                                           <  2:2798627/148‑1 (MQ=255)
                        aCTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGa                                                                                                                          <  1:427131/148‑1 (MQ=255)
                        aCTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGa                                                                                                                          <  1:3148630/148‑1 (MQ=255)
                             cTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGgcg                                                                                                                     >  1:894678/1‑148 (MQ=255)
                             cTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGgcg                                                                                                                     >  1:3443039/1‑148 (MQ=255)
                                 cGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCgctgc                                                                                                                 <  1:2783075/148‑1 (MQ=255)
                                 cGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCgctgc                                                                                                                 <  1:1471839/148‑1 (MQ=255)
                                   tCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGctgctg                                                                                                               <  2:2548632/148‑1 (MQ=255)
                                    cGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTGCGATGGCGCTGCTGa                                                                                                              <  1:3076162/148‑1 (MQ=255)
                                    cGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGa                                                                                                              <  1:2876652/148‑1 (MQ=255)
                                    cGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGa                                                                                                              >  2:1901965/1‑148 (MQ=255)
                                    cGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGa                                                                                                              <  1:2560048/148‑1 (MQ=255)
                                      cACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGAtt                                                                                                            >  2:2970506/1‑148 (MQ=255)
                                        ccGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAg                                                                                                          <  1:463877/148‑1 (MQ=255)
                                         cGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATGGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGa                                                                                                         <  1:3173910/148‑1 (MQ=255)
                                         cGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGa                                                                                                         >  2:3436092/1‑148 (MQ=255)
                                            tGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaaaa                                                                                                      <  1:2269126/148‑1 (MQ=255)
                                            tGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaaaa                                                                                                      <  2:616556/148‑1 (MQ=255)
                                            tGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaaaa                                                                                                      <  2:3562037/148‑1 (MQ=255)
                                             gATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaat                                                                                                     >  2:1217225/1‑148 (MQ=255)
                                                cATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataaa                                                                                                  <  2:915853/148‑1 (MQ=255)
                                                cATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataaa                                                                                                  <  2:904111/148‑1 (MQ=255)
                                                 aTTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataaat                                                                                                 >  1:311161/1‑148 (MQ=255)
                                                   tAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttt                                                                                               <  2:1499415/148‑1 (MQ=255)
                                                   tAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttt                                                                                               <  1:1454557/148‑1 (MQ=255)
                                                   tAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttt                                                                                               <  1:622152/148‑1 (MQ=255)
                                                    aGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAtttt                                                                                              <  2:847163/148‑1 (MQ=255)
                                                     gAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttttt                                                                                             <  1:1901965/148‑1 (MQ=255)
                                                     gAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttttt                                                                                             <  1:3059096/148‑1 (MQ=255)
                                                        tGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCtt                                                                                          <  2:2986791/148‑1 (MQ=255)
                                                        tGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCtt                                                                                          <  1:2004804/148‑1 (MQ=255)
                                                          ttATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCTCAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTAc                                                                                        <  1:549219/148‑1 (MQ=255)
                                                          ttATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTAc                                                                                        <  1:3023985/148‑1 (MQ=255)
                                                            aTTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGc                                                                                      >  1:3079469/1‑148 (MQ=255)
                                                            aTTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGc                                                                                      >  1:408901/1‑148 (MQ=255)
                                                             ttCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCa                                                                                     >  1:1202022/1‑148 (MQ=255)
                                                              tctcCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAg                                                                                    <  2:2306694/148‑1 (MQ=255)
                                                               ctcCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGa                                                                                   <  1:2429788/148‑1 (MQ=255)
                                                                tcCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGAc                                                                                  >  1:3313956/1‑148 (MQ=255)
                                                                  cGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACgg                                                                                >  1:3549102/1‑148 (MQ=255)
                                                                      tcCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGc                                                                            >  1:734098/1‑148 (MQ=255)
                                                                       ccATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCg                                                                           <  1:902494/148‑1 (MQ=255)
                                                                       ccATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCg                                                                           <  1:100176/148‑1 (MQ=255)
                                                                        catcatAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGt                                                                          <  2:154464/148‑1 (MQ=255)
                                                                        catcatAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGt                                                                          >  1:697428/1‑148 (MQ=255)
                                                                         atcatAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGtt                                                                         <  2:2749110/148‑1 (MQ=255)
                                                                         atcatAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGtt                                                                         <  2:1906660/148‑1 (MQ=255)
                                                                            atAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTAt                                                                      >  2:2879892/1‑148 (MQ=255)
                                                                             tAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATc                                                                     <  1:2468305/148‑1 (MQ=255)
                                                                               gCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCa                                                                   >  1:780659/1‑148 (MQ=255)
                                                                                    aTGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCAc                                                              >  1:2949561/1‑148 (MQ=255)
                                                                                      gCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGACGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTc                                                            >  2:1393638/1‑148 (MQ=255)
                                                                                       cAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCg                                                           >  2:1025368/1‑148 (MQ=255)
                                                                                         gACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGa                                                         <  1:1478156/148‑1 (MQ=255)
                                                                                          aCTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGAc                                                        <  1:1905009/148‑1 (MQ=255)
                                                                                          aCTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGAc                                                        <  1:3271736/148‑1 (MQ=255)
                                                                                             ggTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACAtt                                                     <  2:212/148‑1 (MQ=255)
                                                                                             ggTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACAtt                                                     >  1:1232048/1‑148 (MQ=255)
                                                                                                  gTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGAt                                                <  2:787275/148‑1 (MQ=255)
                                                                                                     gcgcACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGt                                             <  2:2676640/148‑1 (MQ=255)
                                                                                                      cgcACGTACGGCTTCCACAGCAGAATCCGCCGGTTTATCCACCGCCAGGAAGTTGATATTTTTAGTACGATGGCTCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCAAGCCACTCGGACATTGGGATAGTg                                            <  2:310344/148‑1 (MQ=255)
                                                                                                      cgcACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTg                                            <  2:311161/148‑1 (MQ=255)
                                                                                                       gcACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTgg                                           <  1:1190023/148‑1 (MQ=255)
                                                                                                       gcACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTgg                                           <  1:2927494/148‑1 (MQ=255)
                                                                                                        cACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGc                                          <  2:3313630/148‑1 (MQ=255)
                                                                                                        cACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGc                                          <  2:813333/148‑1 (MQ=255)
                                                                                                            tACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATaa                                      <  2:2949561/148‑1 (MQ=255)
                                                                                                            tACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATaa                                      >  2:407943/1‑148 (MQ=255)
                                                                                                             aCGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATaaa                                     <  1:1963663/148‑1 (MQ=255)
                                                                                                             aCGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATaaa                                     >  2:3403967/1‑148 (MQ=255)
                                                                                                                gCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGc                                  >  1:2811944/1‑148 (MQ=255)
                                                                                                                 cTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCt                                 >  2:2438149/1‑148 (MQ=255)
                                                                                                                     cacaGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATg                             <  2:115618/148‑1 (MQ=255)
                                                                                                                            gAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGct                      <  2:780659/148‑1 (MQ=255)
                                                                                                                             aaTCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGctc                     <  2:2641030/148‑1 (MQ=255)
                                                                                                                              aTCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACTCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGctct                    >  2:2735944/1‑148 (MQ=255)
                                                                                                                                         tATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTc         >  1:2093245/1‑148 (MQ=255)
                                                                                                                                           tCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACAATGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCg       >  2:3017128/1‑148 (MQ=255)
                                                                                                                                             cACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCggg     >  2:446631/1‑148 (MQ=255)
                                                                                                                                             cACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCggg     >  2:2382449/1‑148 (MQ=255)
                                                                                                                                               ccGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGc   >  1:434613/1‑148 (MQ=255)
                                                                                                                                                cGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCt  >  2:3274806/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCT  >  CP009273/1038014‑1038305

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: ‑