Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 700,407 C→A 5.9% N336K (AAC→AAA nagE → fused N‑acetyl glucosamine specific PTS enzyme: IIC, IIB, and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273700,4070CA 5.9% 572.9 / 9.6 160N336K (AAC→AAAnagEfused N‑acetyl glucosamine specific PTS enzyme: IIC, IIB, and IIA components
Reads supporting (aligned to +/- strand):  ref base C (75/75);  new base A (7/3);  total (82/78)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 162 total aligned reads displayed.

TCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACT  >  CP009273/700263‑700554
                                                                                                                                                |                                                                                                                                                   
tCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTc                                                                                                                                                  >  1:1815266/1‑148 (MQ=255)
     tgctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAt                                                                                                                                             >  1:3454983/1‑148 (MQ=255)
      gctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtg                                                                                                                                            >  1:1351622/1‑148 (MQ=255)
      gctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtg                                                                                                                                            >  1:1351892/1‑148 (MQ=255)
      gctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtg                                                                                                                                            <  2:2712025/148‑1 (MQ=255)
      gctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtg                                                                                                                                            <  2:1329451/148‑1 (MQ=255)
       ctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtgc                                                                                                                                           >  1:1653005/1‑148 (MQ=255)
            cctcctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGctgctg                                                                                                                                      <  2:1046736/148‑1 (MQ=255)
             ctcctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGTTgg                                                                                                                                     <  2:3065558/148‑1 (MQ=255)
                  gcacgcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATg                                                                                                                                <  1:3025557/148‑1 (MQ=255)
                   cacgcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATgg                                                                                                                               >  2:1149899/1‑148 (MQ=255)
                   cacgcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAGGCTGCTGGTGATgg                                                                                                                               >  2:1149878/1‑148 (MQ=255)
                      gcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCg                                                                                                                            >  2:2036152/1‑148 (MQ=255)
                        acTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGtt                                                                                                                          <  2:2677038/148‑1 (MQ=255)
                         ctgctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTa                                                                                                                         >  1:968398/1‑148 (MQ=255)
                           gctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtc                                                                                                                       <  2:1476024/148‑1 (MQ=255)
                            ctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtct                                                                                                                      >  2:3122015/1‑148 (MQ=255)
                            ctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtct                                                                                                                      >  2:2006152/1‑148 (MQ=255)
                            ctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtct                                                                                                                      <  2:793533/148‑1 (MQ=255)
                            ctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtct                                                                                                                      <  2:1496236/148‑1 (MQ=255)
                              gACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtcttc                                                                                                                    >  2:2204160/1‑148 (MQ=255)
                               aCCGGTATCATCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtcttct                                                                                                                   >  1:1236165/1‑148 (MQ=255)
                               aCCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtcttct                                                                                                                   >  1:3183101/1‑148 (MQ=255)
                                    tATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCt                                                                                                              <  2:735335/148‑1 (MQ=255)
                                    tATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCt                                                                                                              <  2:735105/148‑1 (MQ=255)
                                    tATCAGCCTGTTTGTGGCAACGCTGATGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCt                                                                                                              >  2:3455377/1‑148 (MQ=255)
                                      tCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTAt                                                                                                            >  2:3235289/1‑148 (MQ=255)
                                      tCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTAt                                                                                                            <  1:1957674/148‑1 (MQ=255)
                                       cAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATc                                                                                                           <  2:250966/148‑1 (MQ=255)
                                       cAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATc                                                                                                           <  2:250895/148‑1 (MQ=255)
                                        aGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCt                                                                                                          >  2:434429/1‑148 (MQ=255)
                                        aGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCt                                                                                                          <  1:2896299/148‑1 (MQ=255)
                                         gCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTa                                                                                                         >  1:3248814/1‑148 (MQ=255)
                                          ccTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTAc                                                                                                        >  1:1345011/1‑148 (MQ=255)
                                           cTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACt                                                                                                       <  1:3469130/148‑1 (MQ=255)
                                           cTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACt                                                                                                       <  1:790728/148‑1 (MQ=255)
                                            tGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACtt                                                                                                      <  1:64602/148‑1 (MQ=255)
                                             gTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTc                                                                                                     >  2:1014189/1‑148 (MQ=255)
                                              tttGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCg                                                                                                    >  2:1701883/1‑148 (MQ=255)
                                                    gCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGgtgt                                                                                              <  2:3392217/148‑1 (MQ=255)
                                                     cAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGtt                                                                                             >  1:641478/1‑148 (MQ=255)
                                                     cAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGtt                                                                                             <  1:1009466/148‑1 (MQ=255)
                                                       aCGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCa                                                                                           >  1:3310733/1‑148 (MQ=255)
                                                       aCGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCa                                                                                           >  2:2729233/1‑148 (MQ=255)
                                                       aCGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCa                                                                                           >  2:2855471/1‑148 (MQ=255)
                                                       aCGCAGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCa                                                                                           >  2:2143171/1‑148 (MQ=255)
                                                        cGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAg                                                                                          <  2:281645/148‑1 (MQ=255)
                                                        cGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAg                                                                                          <  2:283877/148‑1 (MQ=255)
                                                        cGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAg                                                                                          <  2:3388693/148‑1 (MQ=255)
                                                           tgctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTc                                                                                       >  2:3553645/1‑148 (MQ=255)
                                                            gctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCg                                                                                      <  1:107699/148‑1 (MQ=255)
                                                             ctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCgg                                                                                     >  1:1595468/1‑148 (MQ=255)
                                                             ctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGGTCAGTTCgg                                                                                     >  1:1595490/1‑148 (MQ=255)
                                                                ggTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTa                                                                                  >  1:2432358/1‑148 (MQ=255)
                                                                ggTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTa                                                                                  >  2:1978166/1‑148 (MQ=255)
                                                                  tATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATc                                                                                <  1:2855331/148‑1 (MQ=255)
                                                                   aTCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATcc                                                                               >  1:516443/1‑148 (MQ=255)
                                                                        cgcgGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATg                                                                          >  1:341974/1‑148 (MQ=255)
                                                                         gcgGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGt                                                                         <  1:209230/148‑1 (MQ=255)
                                                                         gcgGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGt                                                                         >  2:2278777/1‑148 (MQ=255)
                                                                          cgGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGtt                                                                        <  2:2232154/148‑1 (MQ=255)
                                                                          cgGGCTTCTCGTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGtt                                                                        <  2:2232036/148‑1 (MQ=255)
                                                                              cTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAAc                                                                    >  2:141869/1‑148 (MQ=255)
                                                                               ttctctttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAAcc                                                                   <  2:2135622/148‑1 (MQ=255)
                                                                                    tttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaaa                                                              <  1:3445419/148‑1 (MQ=255)
                                                                                     ttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaaaa                                                             >  2:559091/1‑148 (MQ=255)
                                                                                      tctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAAc                                                            >  1:902555/1‑148 (MQ=255)
                                                                                       ctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACg                                                           <  1:1865488/148‑1 (MQ=255)
                                                                                         ctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGcc                                                         >  2:92587/1‑148 (MQ=255)
                                                                                         ctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGcc                                                         >  2:1443394/1‑148 (MQ=255)
                                                                                          tGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCg                                                        <  2:2739136/148‑1 (MQ=255)
                                                                                             gggggCTATCGACTACGCTTTGTTTTATAACCTGCCGGCCGCCAGCCATAACGTCTGGATGCTGGTGGTTATGGGCGTTATCTTCTTCGCTATCTACTTCGTGTTGTTCTGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGt                                                     <  2:1548863/148‑1 (MQ=255)
                                                                                              ggggCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTc                                                    >  1:2705171/1‑148 (MQ=255)
                                                                                                       gACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATa                                           >  1:299425/1‑148 (MQ=255)
                                                                                                        aCTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGTTGGGCTTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCTTGTTCAACCTGAAAACGCCGGGTCGTGAAGATaa                                          <  1:880985/148‑1 (MQ=255)
                                                                                                         cTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATaaa                                         <  2:3183101/148‑1 (MQ=255)
                                                                                                         cTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATaaa                                         >  2:356914/1‑148 (MQ=255)
                                                                                                           aCGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaaga                                       >  1:1285413/1‑148 (MQ=255)
                                                                                                           aCGCGTAGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaaga                                       <  2:3040615/148‑1 (MQ=255)
                                                                                                            cgcgTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCTTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaagaa                                      <  1:540568/148‑1 (MQ=255)
                                                                                                             gcgTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaagaag                                     <  2:1309728/148‑1 (MQ=255)
                                                                                                                 tGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACga                                 >  2:3005022/1‑148 (MQ=255)
                                                                                                                   aTGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCGGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACgaga                               <  2:188254/148‑1 (MQ=255)
                                                                                                                       ataACCTGCCGGCCGCCAGCCATAAAGTCTGTATGCTGCTGGTGATGGGCTTTATCTTCTTCGATATCTACTTCGTGGTGTTCTGTTCGTTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGt                           <  1:1978166/148‑1 (MQ=255)
                                                                                                                         aaCCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTa                         >  1:3146870/1‑148 (MQ=255)
                                                                                                                              gccggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaa                    >  1:225258/1‑148 (MQ=255)
                                                                                                                                cggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaaga                  <  2:272081/148‑1 (MQ=255)
                                                                                                                                  gccgccAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCACCCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGaagaag                <  2:1550374/148‑1 (MQ=255)
                                                                                                                                  gccgccAGCCAGAACGTCTGGATGCTGCGGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGaagaag                <  2:1815266/148‑1 (MQ=255)
                                                                                                                                      ccagccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCaa            <  2:588251/148‑1 (MQ=255)
                                                                                                                                      ccagccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCaa            >  2:289155/1‑148 (MQ=255)
                                                                                                                                      ccagccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCaa            <  1:3553645/148‑1 (MQ=255)
                                                                                                                                          ccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGc        >  2:1871197/1‑148 (MQ=255)
                                                                                                                                          ccagAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGc        >  2:1844957/1‑148 (MQ=255)
                                                                                                                                           cagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCa       <  2:3248814/148‑1 (MQ=255)
                                                                                                                                            agAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCaa      <  1:2487464/148‑1 (MQ=255)
                                                                                                                                            agAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCaa      >  2:911391/1‑148 (MQ=255)
                                                                                                                                            agAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCaa      >  2:911404/1‑148 (MQ=255)
                                                                                                                                                cGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACt  >  1:2136727/1‑148 (MQ=255)
                                                                                                                                                cGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAAAAGCAACACt  >  1:2136998/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACT  >  CP009273/700263‑700554

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: