Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,055,789 A→G 7.3% E549G (GAA→GGA)  spoT → bifunctional (p)ppGpp synthase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,055,7890AG 7.3% 611.3 / 7.6 212E549G (GAA→GGA) spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base A (41/141);  new base G (14/2);  total (69/143)
Fisher's exact test for biased strand distribution p-value = 3.21e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

AGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACT  >  NC_002947/6055641‑6055938
                                                                                                                                                    |                                                                                                                                                     
agcagcCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGaa                                                                                                                                                      <  1:4114652/150‑1 (MQ=255)
agcagcCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGaa                                                                                                                                                      <  2:2876675/150‑1 (MQ=255)
 gcagcCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAAc                                                                                                                                                     <  1:3124604/150‑1 (MQ=255)
  cagcCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACa                                                                                                                                                    <  1:3717850/150‑1 (MQ=255)
  cagcCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCAGCGGGGGGCCGAAGGCGAACa                                                                                                                                                    >  1:4566647/1‑150 (MQ=255)
    gcCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCCGCGGTCGGCCGGCGGCGGACacc                                                                                                                                                  >  1:1072103/1‑148 (MQ=255)
    gcCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCAGCAGTCGGCCGAAGGCCAACaac                                                                                                                                                  >  1:6137357/1‑148 (MQ=255)
       tGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCGCCAGCAGTCGGGCCGAGGCGcacaccagc                                                                                                                                               >  1:5986584/1‑145 (MQ=255)
       tGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCCGCAGTCGGCCGAAGGCGGAAagcagc                                                                                                                                               >  1:2127357/1‑150 (MQ=255)
         gagaAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCGGCTGTCGGCCGCAGGCGAACAGCGGCCg                                                                                                                                             >  2:3238561/1‑150 (MQ=255)
         gagaAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCGGCGGGCGGCCCGAGGCGGACAGCAGccc                                                                                                                                             >  2:2599957/1‑149 (MQ=255)
          agaAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCgg                                                                                                                                            <  1:1555575/150‑1 (MQ=255)
              aaTCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCAGCGGGCGGCCGCAGGCGGCCAGCGGCCGGCGcc                                                                                                                                        >  2:4081975/1‑150 (MQ=255)
                tCCCCCAGGAACGCATCCAGTCTATTCTCGCCTTTTACCGCGTGGAGCTCATATAAGACCTGCTCGAAGACATGGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCGGGCAGCTGTTCGCCTTCGGCCGACAGCAGGCGgcgcgcg                                                                                                                                      <  1:3658711/150‑4 (MQ=21)
                  ccccAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCCCCCGCGGGCGGCCGCAGGCCGAAACCAGCCGGCGCCCGac                                                                                                                                    >  1:1748874/1‑149 (MQ=255)
                  ccccAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCCCCGGTGGGCCGCACGCGCACAACAGCCCGCGccccac                                                                                                                                    >  1:4110108/1‑146 (MQ=255)
                  ccccAGGAACGCATCCAGTCTATTCTCGCCGAGTACAGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGaa                                                                                                                                    >  1:1105503/1‑150 (MQ=255)
                   cccAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCCCCCGCAGTCGGCCCCACGCGCACAGCGGCCGGCGccccaac                                                                                                                                   >  1:4695260/1‑145 (MQ=255)
                     cAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGc                                                                                                                                 <  2:3206494/150‑1 (MQ=255)
                      aGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCCCCAGCGGGGGGGCGAAGGCCAACAACAGCCGGCCCCCGAAGGcc                                                                                                                                >  2:2873026/1‑150 (MQ=255)
                      aGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGcc                                                                                                                                <  2:2469278/150‑1 (MQ=255)
                         aaCGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGTGGGccttc                                                                                                                             <  1:2224005/150‑4 (MQ=255)
                         aaCGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCAc                                                                                                                             <  1:6470259/150‑1 (MQ=255)
                          aCGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCGGCGGTCGGCCGGAGGCGCACAGCGGGCCGCGCCCGAAGGccccca                                                                                                                            >  2:2984723/1‑147 (MQ=255)
                            gCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTgg                                                                                                                          <  2:5997418/150‑1 (MQ=255)
                            gCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTgg                                                                                                                          <  1:5800414/150‑1 (MQ=255)
                            gCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTgg                                                                                                                          <  2:4913609/150‑1 (MQ=255)
                              aTCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCCGCGGGCGGCCGCAGGCGGCCACCGGCCCGCGCCCGCAGGCCCCCAGGCg                                                                                                                        >  2:2952215/1‑150 (MQ=255)
                                 cAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCAGCTGGCGGCCGCAGGCGGACAGCGGCCGGCGCCCCAAGGCCCCCAGGCGCTc                                                                                                                     >  2:5697877/1‑150 (MQ=255)
                                 cAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCAGCTGTCGGCCGCAGGCGAACAGCTGCCGGCGCCCGAAGGCCCCCAGGCGATc                                                                                                                     >  2:3988842/1‑150 (MQ=255)
                                  aGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCGCCAGCGGTGGGCCCCAGGCGCCCACCAGACGGCGCCCCAAGGCCCCCAGGCGAtgc                                                                                                                    >  2:2335030/1‑148 (MQ=255)
                                  aGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATcc                                                                                                                    <  2:838259/150‑1 (MQ=255)
                                   gTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCCCCGGCGGGGGGGCGCAGGCGACCCCCGGCCCGGGCCCGCAGGCCCCCCGGCgcggcg                                                                                                                   >  1:6451619/1‑145 (MQ=255)
                                      tATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCgg                                                                                                                <  2:4420788/150‑1 (MQ=255)
                                       aTTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCCGCCGTCGGCCGAAGGCGGACAGCAGGCCGCCCCCGAAGGCCCACCCGCGGTCCCCGGc                                                                                                               >  2:2649282/1‑150 (MQ=255)
                                                   tACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTg                                                                                                   <  1:1191414/150‑1 (MQ=255)
                                                        ccTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCt                                                                                              <  1:1832883/150‑1 (MQ=255)
                                                          tGGAGCTCATAGGAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCa                                                                                            <  1:5494903/150‑1 (MQ=255)
                                                               cTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTAc                                                                                       <  1:1482998/150‑1 (MQ=255)
                                                               cTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTAc                                                                                       <  2:252604/150‑1 (MQ=255)
                                                                 cATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTAcac                                                                                     <  1:5650552/150‑3 (MQ=255)
                                                                    agaagaCCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCaa                                                                                  <  2:2166815/150‑1 (MQ=255)
                                                                    agaagaCCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCaa                                                                                  <  2:3652236/150‑1 (MQ=255)
                                                                    agaagaCCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCGGCTGTGGGGCGCAGGCGACAAGCTGCCCGCGCCCCAAGGCCCACAGGCGATCCGCGGCACCCCAGGCCAGGTGCTCAGCGACGCCaa                                                                                  >  2:2050696/1‑150 (MQ=255)
                                                                       agaCCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGGTGTCGGCCGAAGGCGACAAGCTGCCCGCGCCAGAAGGCCCCCTGGCGCGCCGCGGCCCCGCCGGCCGGGGGCTCCGCGACGCCaatgg                                                                               >  2:6590017/1‑147 (MQ=255)
                                                                           cTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGtgctgc                                                                           <  1:2106187/150‑1 (MQ=255)
                                                                             gCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCAc                                                                         <  1:4483305/150‑1 (MQ=255)
                                                                               tCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCGCCGGCGGTCGGCCCCAGGCGCCCAGCGGCCGGCGCCCCAACGCCCACAGGCGCTCCGCCGCACCGCAGGCCGGGGGGGCCGCGACGCCAACagctgggccc                                                                       >  2:76791/1‑145 (MQ=255)
                                                                               tCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGc                                                                       <  2:5407382/150‑1 (MQ=255)
                                                                               tCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGc                                                                       <  1:1632174/150‑1 (MQ=255)
                                                                               tCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCGGCGGTCGGCCGCAGGCGGACCCCTGCCGGCGCCCCAAGGCCCACAGGCGCTCCGCGGCACCCAACGCCAGGGGGGCAGCGACGCCCACTGCTGCAccc                                                                       >  2:1539105/1‑148 (MQ=255)
                                                                                 gAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCGCCTGCTGTGGGCCGGAGGCGTTCAGCTGCCGGGGGCCGTAGGCCCCCAGGCGGTGCGCGGCACCCAAAGCCGGGTGGTCAGCGACGCCCAGTGGTGCACGCCg                                                                     >  2:3797793/1‑150 (MQ=255)
                                                                                 gAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCAGCTGTCGGCCGAACGCGAACACCGGCCCGCGCCCCAAGGCCCACAGGCGCTCCCCGGCACCCAAGGCCGGGTGCTCAGCAACCCCAACGGCGGCACGccc                                                                     >  2:6088231/1‑149 (MQ=255)
                                                                                       aTCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCCAATGCGAAAAGCTGCCGGCGCCACAAGGCCCACTGGCGCTACGCGGCACCGAAGGCCAGGTGCTCAGCAAAGCCAAGTGCTGGACGCCCATcccc                                                               >  2:6762344/1‑149 (MQ=255)
                                                                                           gCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTg                                                           <  1:2147606/150‑1 (MQ=255)
                                                                                            ccTGGGCAACCGCATGGCCTGCGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGa                                                          <  1:2798590/150‑1 (MQ=255)
                                                                                              tGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGAcc                                                        <  1:717550/150‑1 (MQ=255)
                                                                                                ggCAACCGCATGGCCTACGTGGTCGCGGGCCCCCCGCGGTCGGCCGCAGGCGGACAACTGCCGGCGCCCCAAGGCCCCCCGGCGGTCCGCCGCCCCCCAGGCCGGGTGCGCCGCGACGCCCACTGCGGTTCGCCGAGCCCCGGTGAcccc                                                      >  2:1325936/1‑149 (MQ=255)
                                                                                                   aaCCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGAtt                                                   <  1:4444501/150‑1 (MQ=255)
                                                                                                     ccGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGt                                                 <  1:4434221/150‑1 (MQ=255)
                                                                                                       gCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCg                                               <  1:5379418/150‑1 (MQ=255)
                                                                                                       gCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCg                                               <  1:1915310/150‑1 (MQ=255)
                                                                                                         cgggCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGc                                             <  2:5227178/148‑1 (MQ=255)
                                                                                                         aTGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGc                                             <  1:2555501/150‑1 (MQ=255)
                                                                                                         aTGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGc                                             <  1:5978933/150‑1 (MQ=255)
                                                                                                         aTGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGc                                             <  1:6622673/150‑1 (MQ=255)
                                                                                                         aTGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGc                                             <  1:2872655/150‑1 (MQ=255)
                                                                                                            gCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCAc                                          <  2:1158569/150‑1 (MQ=255)
                                                                                                            gCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCAc                                          <  2:5791769/150‑1 (MQ=255)
                                                                                                              cTACGTGGTCGCGGGCCCCCAGCTGTCGGCCGAACGCGAACAGCTTCCGGCGCCCGAAGGCCCACTGGCGGTCCGCGGCACCCAAGGCCTGGTGCTCAGCTACGCCAAATGGTGCACGCCCATAACGGGGGACACGATTGTAGGCCAccg                                        >  2:516888/1‑149 (MQ=255)
                                                                                                              cTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCt                                        <  1:3584429/150‑1 (MQ=255)
                                                                                                                    ggTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggc                                  <  2:413540/150‑1 (MQ=255)
                                                                                                                    ggTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggc                                  <  2:4682080/150‑1 (MQ=255)
                                                                                                                     gTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggcc                                 <  2:229058/150‑1 (MQ=255)
                                                                                                                     gTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggcc                                 <  2:5027696/150‑1 (MQ=255)
                                                                                                                      tCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggccg                                <  1:2383558/150‑1 (MQ=255)
                                                                                                                            gccgccTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAgg                          <  2:3987927/150‑1 (MQ=255)
                                                                                                                             ccgccTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAggg                         <  1:4235708/150‑1 (MQ=255)
                                                                                                                               gccTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCa                       >  2:1143291/1‑150 (MQ=255)
                                                                                                                                ccTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCAt                      >  2:2792948/1‑150 (MQ=255)
                                                                                                                                  tgctgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATgg                    >  1:1058616/1‑150 (MQ=255)
                                                                                                                                    ctgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGc                                                                       >  2:3403156/1‑97 (MQ=255)
                                                                                                                                    ctgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGc                                                                       <  1:3403156/97‑1 (MQ=255)
                                                                                                                                     tgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCg                 >  2:7205/1‑150 (MQ=255)
                                                                                                                                      gTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGGGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGt                <  1:3238825/150‑1 (MQ=255)
                                                                                                                                       tCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTg               <  2:2372750/150‑1 (MQ=255)
                                                                                                                                       tCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTg               <  1:4081975/150‑1 (MQ=255)
                                                                                                                                        cGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGc              <  1:891796/150‑1 (MQ=255)
                                                                                                                                        cGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGc              <  2:4566647/150‑1 (MQ=255)
                                                                                                                                        cGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGc              <  1:5112202/150‑1 (MQ=255)
                                                                                                                                         ggCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCa             >  1:4304932/1‑150 (MQ=255)
                                                                                                                                          gCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCCGCTCCGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGGCGTGCAc            >  1:2508632/1‑150 (MQ=255)
                                                                                                                                           ccGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCAcc           >  1:1395528/1‑150 (MQ=255)
                                                                                                                                            cGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCt          >  2:109896/1‑150 (MQ=255)
                                                                                                                                            cGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCt          >  1:2973274/1‑150 (MQ=255)
                                                                                                                                             gAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTg         <  2:5649080/150‑1 (MQ=255)
                                                                                                                                               aGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCATCTACGACAAGTGCTGCACGCCGATCCCGGGTGACCCGATTTTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGga       <  2:3832445/150‑1 (MQ=255)
                                                                                                                                               aGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGga       <  2:6137357/150‑1 (MQ=255)
                                                                                                                                               aGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGga       <  1:496008/150‑1 (MQ=255)
                                                                                                                                                  cGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGgagta    <  2:1748874/150‑3 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAACT  >  NC_002947/6055641‑6055938

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