Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,594,526 A→C 5.5% R313R (CGT→CGG kgtP ← alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,594,5260AC 5.5% 270.4 / 6.7 91R313R (CGT→CGGkgtPalpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base A (39/47);  new base C (2/3);  total (41/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

ACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTGGTGTAGGTG  >  NC_002947/1594377‑1594673
                                                                                                                                                     |                                                                                                                                                   
aCGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCAGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGa                                                                                                                                                     >  1:2277134/1‑150 (MQ=255)
     gTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGAcggcg                                                                                                                                                >  2:58810/1‑150 (MQ=255)
     gTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGAcggcg                                                                                                                                                <  1:2384843/150‑1 (MQ=255)
      ttGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGAcggcgg                                                                                                                                               >  1:1288724/1‑150 (MQ=255)
       tGATCGAGGTATCGCCGCTTACGATGATCAGGGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACggcggc                                                                                                                                              <  1:306047/150‑1 (MQ=255)
       tGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACggcggc                                                                                                                                              <  2:1861412/150‑1 (MQ=255)
            gAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGCCGCCGTCCGATc                                                                                                                                         <  2:1476011/150‑1 (MQ=255)
             aGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCt                                                                                                                                        >  2:451416/1‑150 (MQ=255)
              ggTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCtt                                                                                                                                       <  1:58810/150‑1 (MQ=255)
                 ataGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTc                                                                                                                                    <  2:798459/150‑1 (MQ=255)
                    gCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGa                                                                                                                                 <  1:2120230/150‑1 (MQ=255)
                      cGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACa                                                                                                                               <  2:384995/150‑1 (MQ=255)
                       gCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAg                                                                                                                              >  1:1663718/1‑150 (MQ=255)
                        cTTACGATGATCAGTGCGGCCATGATCAGGAAGAAGGCACCCCACCAGGTCGGGAGGGTGTGAAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAGGGCGATAAGGGTCGGAGGGGGGCCGAGCCTGTCGGACAGa                                                                                                                             >  1:1917317/1‑150 (MQ=255)
                          tACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGacc                                                                                                                           <  1:451416/150‑1 (MQ=255)
                            cGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCAGCCTGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGaccac                                                                                                                         <  1:23301/150‑1 (MQ=255)
                             gatgatCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGaccacc                                                                                                                        >  1:192690/1‑150 (MQ=255)
                              atgatCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCg                                                                                                                       >  2:720583/1‑150 (MQ=255)
                              atgatCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCg                                                                                                                       >  2:2219310/1‑150 (MQ=255)
                                  tCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGatca                                                                                                                   <  1:1472736/150‑1 (MQ=255)
                                  tCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGatca                                                                                                                   <  1:207303/150‑1 (MQ=255)
                                   cAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGatcat                                                                                                                  <  1:1743300/150‑1 (MQ=255)
                                    aGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGatcatc                                                                                                                 <  2:1288724/150‑1 (MQ=255)
                                     gTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCg                                                                                                                >  1:2253466/1‑150 (MQ=255)
                                         ggCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCgtaag                                                                                                            <  2:1917317/150‑5 (MQ=255)
                                           ccATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACTGATCATCGGCTGCa                                                                                                          >  1:1789424/1‑150 (MQ=255)
                                           ccATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCa                                                                                                          <  2:192690/150‑1 (MQ=255)
                                           ccATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCa                                                                                                          <  2:2325778/150‑1 (MQ=255)
                                           ccATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCa                                                                                                          >  1:978054/1‑150 (MQ=255)
                                               gATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCTGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGca                                                                                                      >  2:756700/1‑150 (MQ=255)
                                               gATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCACGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACAGCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGct                                                                                                      <  2:1649779/150‑2 (MQ=255)
                                               gATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCACGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGct                                                                                                      <  2:1649628/150‑2 (MQ=255)
                                                aTCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGcac                                                                                                     <  2:730714/150‑1 (MQ=255)
                                                  cAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGCCGCCGTCCGAt                                                                                                                                          >  1:2345866/1‑111 (MQ=255)
                                                  cAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGCCGCCGTCCGAt                                                                                                                                          <  2:2345866/111‑1 (MQ=255)
                                                   aGGAAGAACGCACCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGAc                                                                                                                                                    <  1:630947/99‑1 (MQ=255)
                                                   aGGAAGAACGCACCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGAc                                                                                                                                                    >  2:630947/1‑99 (MQ=255)
                                                     gaagaaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCtt                                                                                                                                       >  1:145909/1‑111 (MQ=255)
                                                     gaagaaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCtt                                                                                                                                       <  2:145909/111‑1 (MQ=255)
                                                       agaaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATgaaca                                                                                              <  1:282465/150‑1 (MQ=255)
                                                        gaaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATgaacag                                                                                             <  2:581779/150‑1 (MQ=255)
                                                        gaaCGCACACCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGTACAGACCACCGATCATCGGCTGCAGGCACATgaacag                                                                                             <  2:181020/150‑1 (MQ=255)
                                                                 ccaccaGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGt                                                                                    >  2:59822/1‑150 (MQ=255)
                                                                 ccaccaGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGt                                                                                    >  1:2117761/1‑150 (MQ=255)
                                                                    ccaGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGCCGCCGTCCGATc                                                                                                                                         <  1:1537280/94‑1 (MQ=255)
                                                                    ccaGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGCCGCCGTCCGATc                                                                                                                                         >  2:1537280/1‑94 (MQ=255)
                                                                             gATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGa                                                                                                                                           >  1:1636260/1‑83 (MQ=255)
                                                                             gATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGa                                                                                                                                           <  2:1636260/83‑1 (MQ=255)
                                                                             gATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGatc                                                                                                                    <  1:2010907/106‑1 (MQ=255)
                                                                             gATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGatc                                                                                                                    >  2:2010907/1‑106 (MQ=255)
                                                                             gATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGAt                                                                        <  2:2277134/150‑1 (MQ=255)
                                                                                ggTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATAGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGAgggg                                                                     >  1:2051395/1‑149 (MQ=255)
                                                                                 gtgtGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTc                                                                                                                                    <  1:1372584/86‑1 (MQ=255)
                                                                                 gtgtGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTc                                                                                                                                    >  2:1372584/1‑86 (MQ=255)
                                                                                    tgCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGa                                                                                                                                 >  2:1811622/1‑86 (MQ=255)
                                                                                    tgCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGa                                                                                                                                 <  1:1811622/86‑1 (MQ=255)
                                                                                      cAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGa                                                               >  1:1440768/1‑150 (MQ=255)
                                                                                           tGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGa                                                               >  1:2182702/1‑145 (MQ=255)
                                                                                           tGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGa                                                               <  2:2182702/145‑1 (MQ=255)
                                                                                              tGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGa                                                                                                                             >  2:903455/1‑80 (MQ=255)
                                                                                              tGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGa                                                                                                                             <  1:903455/80‑1 (MQ=255)
                                                                                              tGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGa                                                       >  1:746785/1‑150 (MQ=255)
                                                                                                  gATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCt                                                   >  2:708943/1‑150 (MQ=255)
                                                                                                     ggggACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCAt                                                <  2:746785/150‑1 (MQ=255)
                                                                                                     ggggACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGAACTTGTCAGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGAGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCAt                                                >  1:2389185/1‑150 (MQ=255)
                                                                                                           ggTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCAc                                          <  1:708943/150‑1 (MQ=255)
                                                                                                                aCAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGgtgt                                     <  2:2253466/150‑1 (MQ=255)
                                                                                                                 cAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGGACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGtt                                    <  2:1369604/150‑1 (MQ=255)
                                                                                                                   gggTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGc                                                                                      >  2:333509/1‑98 (MQ=255)
                                                                                                                   gggTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGc                                                                                      <  1:333509/98‑1 (MQ=255)
                                                                                                                    ggTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGTAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGGGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGAc                                 >  2:368200/1‑150 (MQ=255)
                                                                                                                      tGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGActa                               >  1:1096982/1‑148 (MQ=255)
                                                                                                                      tGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCa                               <  1:2043077/150‑1 (MQ=255)
                                                                                                                       gCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACTAg                              <  2:1096982/150‑1 (MQ=255)
                                                                                                                       gCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAg                              <  2:1440768/150‑1 (MQ=255)
                                                                                                                            gCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCAc                                                          <  2:1784963/117‑1 (MQ=255)
                                                                                                                            gCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCAc                                                          >  1:1784963/1‑117 (MQ=255)
                                                                                                                                 ccGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGc                    >  1:802409/1‑150 (MQ=255)
                                                                                                                                  cGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGAtggt                                                                     <  1:1346815/100‑1 (MQ=255)
                                                                                                                                  cGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGAtggt                                                                     >  2:1346815/1‑100 (MQ=255)
                                                                                                                                  cGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCa                   >  2:2226915/1‑150 (MQ=255)
                                                                                                                                  cGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCa                   <  1:756700/150‑1 (MQ=255)
                                                                                                                                    aaGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGt                                                                                    >  2:1982014/1‑83 (MQ=255)
                                                                                                                                    aaGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGt                                                                                    <  1:1982014/83‑1 (MQ=255)
                                                                                                                                          aTCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCga                                                                           >  2:2166774/1‑86 (MQ=255)
                                                                                                                                          aTCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCga                                                                           <  1:2166774/86‑1 (MQ=255)
                                                                                                                                               gATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAg                                                              >  2:106388/1‑94 (MQ=255)
                                                                                                                                               gATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAg                                                              <  1:106388/94‑1 (MQ=255)
                                                                                                                                               gATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTggtgta      >  1:1808976/1‑150 (MQ=255)
                                                                                                                                               gATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTggtgta      >  1:12148/1‑150 (MQ=255)
                                                                                                                                                aTCGGACGGCGGCCGATCTTGTCGGACAGACCACTGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTggtgtag     <  2:1789424/150‑1 (MQ=255)
                                                                                                                                                   ggACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGTTGCAGGCACATGAACAGGAACAGCGTGGCGGTCGTGATGGTGGTTGAGTCACTGATGGTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTggtgtaggtg  <  2:119442/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTGGTGTAGGTG  >  NC_002947/1594377‑1594673

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