Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,041 A→T 95.0% pseudogene (1393/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,461,0410AT95.0% 546.2 / 13.7 178pseudogene (1393/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (5/4);  new base T (104/65);  total (109/69)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 179 total aligned reads displayed.

CGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCA  >  CP009273/1460896‑1461188
                                                                                                                                                 |                                                                                                                                                   
cGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTcc                                                                                                                                                   >  2:530159/1‑148 (MQ=255)
      tgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGAc                                                                                                                                             >  1:3310691/1‑148 (MQ=255)
       gaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACa                                                                                                                                            <  2:2615747/148‑1 (MQ=255)
          gCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAgcg                                                                                                                                         <  1:2522986/148‑1 (MQ=255)
          gCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAgcg                                                                                                                                         <  2:1885679/148‑1 (MQ=255)
           ccACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAgcgc                                                                                                                                        >  2:1498920/1‑148 (MQ=255)
            cACGGCAAACAACATCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGATATTGCTGGAATTAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCa                                                                                                                                       <  1:2667217/148‑1 (MQ=14)
            cACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCa                                                                                                                                       >  1:1866634/1‑148 (MQ=255)
             aCGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCaa                                                                                                                                      <  1:2625861/148‑1 (MQ=255)
                gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACgg                                                                                                                                   <  2:508026/148‑1 (MQ=255)
                gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACgg                                                                                                                                   <  2:1866634/148‑1 (MQ=255)
                 cAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGt                                                                                                                                  <  2:481954/148‑1 (MQ=255)
                    acaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGa                                                                                                                               <  2:3481932/148‑1 (MQ=255)
                    acaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGa                                                                                                                               >  1:3481932/1‑148 (MQ=255)
                    acaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGa                                                                                                                               >  1:1306298/1‑148 (MQ=255)
                      aacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATa                                                                                                                             <  2:1562329/148‑1 (MQ=255)
                        caacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAAc                                                                                                                           >  2:2549847/1‑148 (MQ=255)
                         acccgCAACACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACa                                                                                                                          <  2:2719196/144‑1 (MQ=255)
                          acGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACaa                                                                                                                         >  2:3004141/1‑148 (MQ=255)
                           cGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAg                                                                                                                        >  2:316342/1‑148 (MQ=255)
                           cGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAg                                                                                                                        >  2:2864991/1‑148 (MQ=255)
                              cAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGCCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATGGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGc                                                                                                                     <  1:3435720/148‑1 (MQ=255)
                                aaaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCAc                                                                                                                   <  2:1255718/148‑1 (MQ=255)
                                 aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCAcc                                                                                                                  >  2:660157/1‑148 (MQ=255)
                                 aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCAcc                                                                                                                  <  1:2870678/148‑1 (MQ=255)
                                 aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCAcc                                                                                                                  <  1:277878/148‑1 (MQ=255)
                                 aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGAGAACAAGGGCAcc                                                                                                                  >  1:1622604/1‑148 (MQ=255)
                                  accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCa                                                                                                                 >  2:1715480/1‑148 (MQ=255)
                                   ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAt                                                                                                                >  2:2971804/1‑148 (MQ=255)
                                     accGTTGACGGCAAGGATTCCACCGGTACGGAACTTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGa                                                                                                              <  2:3230987/148‑1 (MQ=255)
                                         ttGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGt                                                                                                          >  2:3211754/1‑148 (MQ=255)
                                          tGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTc                                                                                                         >  2:2784573/1‑148 (MQ=255)
                                          tGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTc                                                                                                         >  2:2062102/1‑148 (MQ=255)
                                           gACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCa                                                                                                        <  1:292752/148‑1 (MQ=255)
                                           gACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCa                                                                                                        <  1:3155036/148‑1 (MQ=255)
                                           gACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCa                                                                                                        <  1:3155046/148‑1 (MQ=255)
                                            aCGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAc                                                                                                       >  2:3493131/1‑148 (MQ=255)
                                                cAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                                                                                   >  1:2701961/1‑148 (MQ=255)
                                                cAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                                                                                   >  1:1723220/1‑148 (MQ=255)
                                                  aGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATcc                                                                                                 >  1:241057/1‑148 (MQ=255)
                                                   ggATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                >  1:1877680/1‑148 (MQ=255)
                                                    gATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                               >  1:42236/1‑148 (MQ=255)
                                                    gATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                               >  1:1247149/1‑148 (MQ=255)
                                                    gATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                               >  1:1247142/1‑148 (MQ=255)
                                                    gATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                               >  1:1845584/1‑148 (MQ=255)
                                                    gATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                               <  1:2084198/148‑1 (MQ=255)
                                                     aTTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                                                                              <  1:2615325/148‑1 (MQ=255)
                                                       tCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                                                                            >  1:79413/1‑148 (MQ=255)
                                                          aCCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCa                                                                                         <  2:228991/148‑1 (MQ=255)
                                                               tACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGt                                                                                    >  2:3310018/1‑148 (MQ=255)
                                                               tACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGt                                                                                    >  2:3309962/1‑148 (MQ=255)
                                                                aCGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTa                                                                                   >  1:1680189/1‑148 (MQ=255)
                                                                 cGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAt                                                                                  <  2:421868/148‑1 (MQ=255)
                                                                 cGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAt                                                                                  >  1:3158359/1‑148 (MQ=255)
                                                                  ggAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATc                                                                                 >  2:1584517/1‑148 (MQ=255)
                                                                        tGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGAtt                                                                           >  2:2915847/1‑148 (MQ=255)
                                                                         gCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTg                                                                          <  2:1362999/148‑1 (MQ=255)
                                                                           tGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGAc                                                                        <  1:2784573/148‑1 (MQ=255)
                                                                             gCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACgg                                                                      >  2:483720/1‑148 (MQ=255)
                                                                                      ggAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGc                                                             >  1:2878913/1‑148 (MQ=38)
                                                                                      ggAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGc                                                             >  1:886090/1‑148 (MQ=38)
                                                                                         aGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCAt                                                          <  1:3004141/148‑1 (MQ=38)
                                                                                         aGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCAt                                                          <  1:1490989/148‑1 (MQ=38)
                                                                                          ggTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATc                                                         >  1:1912886/1‑148 (MQ=38)
                                                                                          ggTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATc                                                         >  1:1912747/1‑148 (MQ=38)
                                                                                           gTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCg                                                        <  2:2277902/148‑1 (MQ=38)
                                                                                           gTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAt                                                                                                                >  1:3391882/1‑92 (MQ=33)
                                                                                           gTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAt                                                                                                                <  2:3391882/92‑1 (MQ=33)
                                                                                             gATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTc                                                      >  1:1153160/1‑148 (MQ=38)
                                                                                             gATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTc                                                      >  1:1152048/1‑148 (MQ=38)
                                                                                             gATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTc                                                      >  1:10533/1‑148 (MQ=38)
                                                                                                tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAAt                                                   >  2:13534/1‑148 (MQ=255)
                                                                                                  aGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATaa                                                 >  1:1501373/1‑148 (MQ=255)
                                                                                                  aGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATaa                                                 >  1:3272713/1‑148 (MQ=255)
                                                                                                    gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACg                                               >  2:205884/1‑148 (MQ=255)
                                                                                                    gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACg                                               >  2:2006979/1‑148 (MQ=255)
                                                                                                    gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACg                                               >  1:1979533/1‑148 (MQ=255)
                                                                                                       ggCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAg                                            <  1:2062102/148‑1 (MQ=255)
                                                                                                         cGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGc                                          >  2:220425/1‑148 (MQ=255)
                                                                                                           aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACCGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCga                                        >  2:1311923/1‑148 (MQ=255)
                                                                                                           aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGaa                                             <  2:2517594/143‑1 (MQ=255)
                                                                                                           aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGaa                                             >  1:2517594/1‑143 (MQ=255)
                                                                                                                gATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCACTAACGAAGGCGAGAGCa                                   >  1:677100/1‑148 (MQ=255)
                                                                                                                gATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCa                                   >  2:2445240/1‑148 (MQ=255)
                                                                                                                gATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCa                                   >  2:2444354/1‑148 (MQ=255)
                                                                                                                 aTGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAgg                                                              <  2:1744631/120‑1 (MQ=17)
                                                                                                                 aTGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAc                                  >  1:2494238/1‑148 (MQ=255)
                                                                                                                 aTGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAc                                  >  1:668780/1‑148 (MQ=255)
                                                                                                                       gcggcggcgGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAc                            >  2:459833/1‑148 (MQ=255)
                                                                                                                           cggcggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAAt                        >  1:1025986/1‑148 (MQ=255)
                                                                                                                            ggcggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATg                       <  2:103275/148‑1 (MQ=255)
                                                                                                                                 tCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcggc                  >  2:3271285/1‑148 (MQ=255)
                                                                                                                                 tCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcggc                  >  2:1499628/1‑148 (MQ=255)
                                                                                                                                 tCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcggc                  >  1:1614699/1‑148 (MQ=255)
                                                                                                                                       tATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc            <  1:2553516/148‑1 (MQ=255)
                                                                                                                                        aTTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggca           <  2:668780/148‑1 (MQ=255)
                                                                                                                                          tGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACt         >  2:3117226/1‑148 (MQ=255)
                                                                                                                                           gATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTc        <  2:1247142/148‑1 (MQ=255)
                                                                                                                                              atCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa     >  1:2818510/1‑148 (MQ=255)
                                                                                                                                                 aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGCTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCa  <  2:149271/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCA  >  CP009273/1460896‑1461188

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: