Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 284,587 T→A 13.4% intergenic (+33/+30) betT‑II → / ← atsK choline/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947284,5870TA13.4% 562.3 / 71.4 223intergenic (+33/+30)betT‑II/atsKcholine/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (109/80);  new base A (20/10);  total (132/91)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-01

Only 100 of 230 total aligned reads displayed.

ACCTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCTG  >  NC_002947/284439‑284736
                                                                                                                                                    |                                                                                                                                                      
   tGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGcc                                                                                                                                                    >  1:55517/1‑150 (MQ=255)
    gAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGccc                                                                                                                                                   <  2:825100/150‑1 (MQ=255)
    gAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGccc                                                                                                                                                   <  2:3093644/150‑1 (MQ=255)
     aGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCcg                                                                                                                                                  <  1:2517324/150‑1 (MQ=255)
     aGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCcg                                                                                                                                                  <  1:888325/150‑1 (MQ=255)
        gtacccAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGa                                                                                                                                               <  2:1740635/145‑1 (MQ=255)
        gAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGa                                                                                                                                               >  1:1527741/1‑150 (MQ=255)
        gAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGa                                                                                                                                               <  1:1955917/150‑1 (MQ=255)
        gAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                               <  2:3777845/150‑1 (MQ=255)
         aaGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                              <  1:3745326/150‑1 (MQ=255)
          aGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAag                                                                                                                                             <  2:4181204/150‑1 (MQ=255)
          aGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGAag                                                                                                                                             <  2:239435/150‑1 (MQ=255)
           ggCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                            <  2:2881579/150‑1 (MQ=255)
           ggCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                            <  2:266600/150‑1 (MQ=255)
             cAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACTAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                                                 >  1:3056570/1‑143 (MQ=255)
             cAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACTAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                                                 <  2:3056570/143‑1 (MQ=255)
              aGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                         >  2:109136/1‑150 (MQ=255)
               gCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGcc                                                                                                                                        >  1:4180147/1‑150 (MQ=255)
                ccAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCa                                                                                                                                       >  2:1463172/1‑150 (MQ=255)
                ccAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCa                                                                                                                                       <  1:3644462/150‑1 (MQ=255)
                      cTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAgg                                                                                                                                >  1:4483467/1‑150 (MQ=255)
                       tATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTAGCGGGCTTGCCCGCGAAGAGGCCAGACAGGc                                                                                                                               >  1:443430/1‑150 (MQ=255)
                         tGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCaa                                                                                                                             >  1:182376/1‑150 (MQ=255)
                          gACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAAc                                                                                                                            >  1:2123717/1‑150 (MQ=255)
                              tggtggGCTACAGCATGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACTTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCt                                                                                                                        <  1:2631603/150‑1 (MQ=255)
                                gtggGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCCCGAAGAGGCCAGACAGGCAACCGCTAg                                                                                                                      >  1:788405/1‑150 (MQ=255)
                                      tACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGt                                                                                                                <  2:4479706/150‑1 (MQ=255)
                                       aCAGCAAGGAGCAGATAATCAACGACATCCTCGACCAGTACTAGCGTCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTTTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTg                                                                                                               <  2:1262055/150‑1 (MQ=255)
                                           cAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATc                                                                                                           >  2:1286344/1‑150 (MQ=255)
                                            aaGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGGCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATcc                                                                                                          >  1:1699100/1‑150 (MQ=255)
                                              ggAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTAGCGGGCGTGCCCCCGAAGCGGCCCGACAGGCACACGCGAGAAACGGGCAAccac                                                                                                        >  1:1918148/1‑148 (MQ=255)
                                              ggAGAAGATCATCAACGACATCCTCCACCAGTACTCGGGGCACTTGCTGTACCTGCACCTGGTCCGCTGACAACGCTCTCCTGTCCTGGCCTGTTGGGGGTCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATcccc                                                                                                        <  1:1021638/150‑1 (MQ=255)
                                               gAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCa                                                                                                       >  2:4430182/1‑150 (MQ=255)
                                                 gCAGATCATCACCGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAg                                                                                                     <  2:1152110/150‑1 (MQ=255)
                                                       catcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGGAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGt                                                                                               >  1:4222668/1‑150 (MQ=255)
                                                        atcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTc                                                                                              >  1:3169789/1‑150 (MQ=255)
                                                         tcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCg                                                                                             >  1:1083856/1‑150 (MQ=255)
                                                         tcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCCCGCAGAGGCCAGAAAGGCAACCGCAAGTAAGGTGCATCCCCAACGATTGTCg                                                                                             >  1:33744/1‑150 (MQ=255)
                                                         tcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCAGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCg                                                                                             >  1:2187969/1‑150 (MQ=255)
                                                           aaCGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGcc                                                                                           >  2:757839/1‑150 (MQ=255)
                                                           aaCGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTAGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGcc                                                                                           >  1:1353222/1‑150 (MQ=255)
                                                              gACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCgg                                                                                        >  2:1474110/1‑150 (MQ=255)
                                                              gACACCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCCCCTGTCCTGGCCTGTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCgg                                                                                        <  1:699830/150‑1 (MQ=255)
                                                                 aTCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGc                                                                                     <  2:2108647/150‑1 (MQ=255)
                                                                     tCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCa                                                                                 >  1:732957/1‑150 (MQ=255)
                                                                     tCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCa                                                                                 >  2:1931522/1‑150 (MQ=255)
                                                                        aCCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGCCCAAGCTCTCCTGTCCTGGCCTCTGCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCAccc                                                                              <  2:2531217/150‑1 (MQ=255)
                                                                        aCCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCAccc                                                                              <  1:3360148/150‑1 (MQ=255)
                                                                          cAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgc                                                                            >  2:339401/1‑150 (MQ=255)
                                                                            ttACGTGGGGCACTGGCCGTCCCTGCACCGGGCCCGCTCACCACGCTCGCCTGCCGGGGCCTCTTCGGGGGCTGCCCCGCTAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgccg                                                                          <  1:1300611/149‑1 (MQ=255)
                                                                            gtccgcgcgcCACAGGCAGTACCTCCACCGGGCCCGCTCACCACGCACCCCGGCCCGGGCCCCTGCGGGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgccg                                                                          <  1:36247/144‑1 (MQ=255)
                                                                             tACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCgg                                                                         >  2:337226/1‑150 (MQ=255)
                                                                              aCGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTCGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCggg                                                                        >  1:4087767/1‑149 (MQ=255)
                                                                               cgcgCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTGCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGtt                                                                       <  2:1893041/147‑1 (MQ=255)
                                                                                gAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtc                                                                      >  2:3973122/1‑150 (MQ=255)
                                                                                  gCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcgg                                                                    >  2:3219412/1‑150 (MQ=255)
                                                                                  gCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcgg                                                                    >  2:1691012/1‑150 (MQ=255)
                                                                                     gCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATACCCAAGGATGGTCGCCCGGTGCATCACCCCCCGGTtcgggcg                                                                 >  1:4219274/1‑150 (MQ=255)
                                                                                       acaTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTAGCGGGCTTGCCCCCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCAACCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGGTcggtcggt                                                               >  1:971277/1‑150 (MQ=255)
                                                                                             aGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                            <  2:2996564/68‑1 (MQ=255)
                                                                                             aGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                            >  1:2996564/1‑68 (MQ=255)
                                                                                              gTACCTGCACCTGGTCCGCTGACCATGCTCTCCTGGCCTGGCCTCTGCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatcatc                                                        <  2:3853070/150‑1 (MQ=255)
                                                                                              gTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatcatc                                                        >  2:726336/1‑150 (MQ=255)
                                                                                              gTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCCCGAAGAGGCCCGTCAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATTGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTGataaac                                                        >  1:2029099/1‑146 (MQ=255)
                                                                                               tACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGGGCATCACCCGCCGGGTCGGTCGGTAATCATcc                                                       >  1:3680621/1‑150 (MQ=255)
                                                                                                  cTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCAcc                                                    >  2:4271951/1‑150 (MQ=255)
                                                                                                      aCCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCaa                                                >  2:411019/1‑150 (MQ=255)
                                                                                                        cTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCaaaa                                              <  2:3022706/150‑1 (MQ=255)
                                                                                                              cGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGCCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGCTGCATCACCCGCCGGTTCGGACGGTAATCATCCACCGCAAAATGCTgg                                        >  2:1820270/1‑150 (MQ=255)
                                                                                                               gCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGGTCGGTCGGTAATCATCCACCGCAAAATGCTggg                                       >  2:381032/1‑150 (MQ=255)
                                                                                                                   aCCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCAc                                   >  2:2329455/1‑150 (MQ=255)
                                                                                                                            tcctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTc                          <  1:2329455/150‑1 (MQ=255)
                                                                                                                            tcctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTc                          >  2:2707747/1‑150 (MQ=255)
                                                                                                                            tcctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGGTCGGTCGGTAAAAATCCACCGCAAAATGCTGGGTCACGCGGTTGTc                          >  1:2753242/1‑150 (MQ=255)
                                                                                                                                 tcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGa                     >  2:1394414/1‑150 (MQ=255)
                                                                                                                                    tgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACg                  >  2:1621217/1‑150 (MQ=255)
                                                                                                                                     ggCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGc                 >  1:2683613/1‑150 (MQ=255)
                                                                                                                                     ggCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCAACCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGc                 >  1:3981990/1‑150 (MQ=255)
                                                                                                                                       ccTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCa               <  2:55517/150‑1 (MQ=255)
                                                                                                                                        ctctTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCAc              >  2:2594902/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACg             >  1:2301730/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACg             >  1:805802/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACg             >  1:2129558/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACg             >  2:2963160/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACg             >  2:3551519/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGAGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACg             >  1:4300143/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACg             >  1:392103/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACg             >  2:3309476/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACg             >  2:1570779/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTc           >  2:4509240/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCAGGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGt            >  2:2825680/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTc           >  2:2389996/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTc           >  2:1359545/1‑150 (MQ=255)
                                                                                                                                            tCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCa          <  2:1700970/150‑1 (MQ=255)
                                                                                                                                            tCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCa          >  2:163921/1‑150 (MQ=255)
                                                                                                                                            tCGCGGGCAAGCCCGCGAAGAGGCCAGGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTc           >  2:2431489/1‑150 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATtc       <  2:4376316/150‑1 (MQ=255)
                                                                                                                                                 ggCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATtctc     <  2:4059762/150‑1 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCTg  >  1:1739939/1‑150 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCTg  >  1:1646666/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
ACCTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCTG  >  NC_002947/284439‑284736

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: