Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 85,077 T→G 7.0% intergenic (+55/+75) aroE → / ← PP_0075 shikimate 5‑dehydrogenase/choline sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294785,0770TG 7.0% 750.2 / 33.6 256intergenic (+55/+75)aroE/PP_0075shikimate 5‑dehydrogenase/choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base T (124/110);  new base G (8/10);  total (133/123)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

Only 100 of 261 total aligned reads displayed.

AGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGC  >  NC_002947/84928‑85226
                                                                                                                                                     |                                                                                                                                                      
aGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGt                                                                                                                                                        <  1:969730/150‑1 (MQ=255)
  caggcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAc                                                                                                                                                     <  1:8298246/150‑1 (MQ=255)
     gcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtc                                                                                                                                                  <  1:1750298/150‑1 (MQ=255)
     gcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtc                                                                                                                                                  <  1:4813673/150‑1 (MQ=255)
     gcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtc                                                                                                                                                  >  1:3732684/1‑150 (MQ=255)
       agCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACCAAGTAcgtcgt                                                                                                                                                >  2:5756569/1‑150 (MQ=255)
       agCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtcgt                                                                                                                                                <  1:3835379/150‑1 (MQ=255)
          tGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAACTAAGTCGTACg                                                                                                                                             >  1:447405/1‑150 (MQ=255)
           gAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGc                                                                                                                                            >  2:4689159/1‑150 (MQ=255)
              gCCTTCTTTATCTGGCGTGGCGTGGGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGa                                                                                                                                         >  2:3417620/1‑150 (MQ=255)
              gCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGGTCCCACAAAGTACGTCGTACGCCGa                                                                                                                                         >  1:7082948/1‑150 (MQ=255)
               ccTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGaa                                                                                                                                        >  2:6861436/1‑149 (MQ=255)
                 ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAATGCCGGTCCCACAAAGTACGTCGTACGCCGATcc                                                                                                                                      >  1:3276864/1‑150 (MQ=255)
                 ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATcc                                                                                                                                      >  2:6111866/1‑150 (MQ=255)
                 ttcttTATCTGGCGTGGCATGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATcc                                                                                                                                      >  2:1796412/1‑150 (MQ=255)
                   cttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGTCTGCTTTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTacgacgtacttt                                                                                                                                    <  2:8668003/150‑7 (MQ=18)
                    tttATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGTGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCctct                                                                                                                                   <  2:4003339/150‑1 (MQ=255)
                    tttATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCctct                                                                                                                                   >  2:1465257/1‑150 (MQ=255)
                      tATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                 <  1:4044953/150‑1 (MQ=255)
                       aTCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgtg                                                                                                                                <  1:8027878/150‑1 (MQ=255)
                       aTCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgtg                                                                                                                                <  1:7415935/150‑1 (MQ=255)
                        tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTgg                                                                                                                               <  2:2420324/150‑1 (MQ=255)
                        tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTgg                                                                                                                               <  2:7227150/150‑1 (MQ=255)
                        tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTgg                                                                                                                               >  2:5381346/1‑150 (MQ=255)
                          tggcgtggcgTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGa                                                                                                                             <  1:7041075/150‑1 (MQ=255)
                          tggcgtggcgTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGa                                                                                                                             <  2:4774422/150‑1 (MQ=255)
                           ggcgtggcgtGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGaac                                                                                                                           <  2:7505720/150‑3 (MQ=18)
                           ggcgtggcgtGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGaac                                                                                                                           <  2:4728729/150‑3 (MQ=18)
                           ggcgtggcgtGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAg                                                                                                                            <  2:2789130/150‑1 (MQ=255)
                            gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGc                                                                                                                           >  2:6551557/1‑150 (MQ=255)
                            gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGaac                                                                                                                           <  1:2666426/150‑3 (MQ=25)
                              gtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGTCGGCTTTGCAGCCCAGCGCGACGCAAGGCCACTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGc                                                                                                                        <  2:6539766/150‑1 (MQ=12)
                              gtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCgg                                                                                                                         >  1:668737/1‑150 (MQ=255)
                              gtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATACTCTGTGGGAGCgg                                                                                                                         >  1:5580265/1‑150 (MQ=255)
                                        gCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTcg                                                                                                               >  2:3248618/1‑150 (MQ=255)
                                        gCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCGTGCGTcg                                                                                                               >  1:2957867/1‑150 (MQ=255)
                                            gATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGAt                                                                                                           <  2:5433144/150‑1 (MQ=255)
                                              tACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCACGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATgg                                                                                                         <  2:3800295/150‑1 (MQ=255)
                                              tACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCGGGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATgg                                                                                                         <  1:1306999/150‑1 (MQ=255)
                                               aCGGCGCCGGTGTGGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATggg                                                                                                        <  2:266981/150‑1 (MQ=255)
                                                 ggCGCCGGTGTGGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGGTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCACGGCCGTCCCCACAACGCACGTCGTACGCCGAGCCTCTGTGGGAGCGGCCTGGCGCCGCGATGGGCt                                                                                                      <  2:3088265/150‑1 (MQ=255)
                                                      cGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCaaa                                                                                                 <  2:7204876/150‑1 (MQ=255)
                                                        gtgtTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGc                                                                                               >  2:4556750/1‑150 (MQ=255)
                                                         tgtTGGCTGAACTGCGCCGGCAACTCGCTGGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAACTTACGTCGTCGGCCGATCCGCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCa                                                                                              <  2:259883/150‑1 (MQ=39)
                                                           ttGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCGGGCCGGGGCCGCGGGGCGGCCCTATCGCGACGCACGGCCGTTCCCCCAACGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGa                                                                                            <  2:5659967/150‑1 (MQ=39)
                                                            tGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAc                                                                                           >  2:8294836/1‑150 (MQ=255)
                                                              gCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAccc                                                                                         >  2:679359/1‑150 (MQ=255)
                                                                tGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGGCGCGGTGGGCTGCAAAGCAGACCCAg                                                                                       >  2:8428196/1‑150 (MQ=255)
                                                                   aCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTCCGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAAt                                                                                    <  2:90578/150‑1 (MQ=255)
                                                                   aCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGACCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAAt                                                                                    <  2:1399989/150‑1 (MQ=255)
                                                                   aCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTGGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAAt                                                                                    <  2:174681/150‑1 (MQ=255)
                                                                       cgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCa                                                                                <  2:5249350/150‑1 (MQ=255)
                                                                            gCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGGCGCAAGGCCGGTCCCCCAAAGTAACTCGTACGCCGCTCCTCGGTGGGGGCGGCCGTGCGTCGCCATGGGCTGCAAAGCAGACCCACAAaaccaaccgc                                                                           >  1:6390125/1‑148 (MQ=255)
                                                                              aaCTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCAt                                                                         >  1:2992586/1‑150 (MQ=255)
                                                                               aCTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCAtt                                                                        <  2:1166906/150‑1 (MQ=255)
                                                                                ctcgctggGGGCTGCGGTGGCCGGGGCCGCTTGGCGGCCCTATCGCGCCGCACGGCCGTCCCCACAAATTACGTCGTCCGCCGACCCTCTGTGGGAGCGGCCTTGCGCCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTc                                                                       <  1:5037450/150‑1 (MQ=255)
                                                                                ctcgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTc                                                                       >  2:5306317/1‑150 (MQ=255)
                                                                                cgcgcggggggCTGCGGGGGCGGGGGCCGCTTGGCGGCCCTAGCGCGACGAACGCCCGTCCCCACACCGTACGTCGTCCGCCGACCCACTGTGGGCGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTc                                                                       <  1:2596797/143‑1 (MQ=25)
                                                                                 tcgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGTCGCAAGGCCGTTGCCCCAAACTAAGGCGTACGCCCCTCCTCGGTGGGGGCGGCCGGGCGGCGCGGTGGGGGGCAAAGCCGACCCCGGAGCCCCACCACATTcc                                                                      >  1:4798001/1‑150 (MQ=255)
                                                                                  cgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCt                                                                     <  2:8805617/150‑1 (MQ=255)
                                                                                   gctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCCTTGCCACAAAGTACGTCGTACGCCCATCCACTGTGGGGGCGGCCTTGCGGCGCGATGGGCTGCAAAGCAGACCCAGAAACCCCACCACATTccac                                                                    >  1:5566841/1‑148 (MQ=255)
                                                                                     tcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTACCACAAAGTAAGTCGGACGCCGATCCACTGTGGGAGCGGCCGTTCCTCGCGATGGGCTGCAAAGCACAACCAGAATACCAAACACATTACACaa                                                                  >  1:8766548/1‑150 (MQ=37)
                                                                                                cgggCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCggg                                                       <  2:8060329/148‑1 (MQ=255)
                                                                                                cTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCggg                                                       <  2:8042658/150‑1 (MQ=255)
                                                                                                  ggCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCCCAAAGTACGTCGTACGCCGATCCTCTGTGGGGGCGGCCGTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAAACCCAACCACAATCCTCAAACCGGATCGGGca                                                     >  1:5247360/1‑150 (MQ=39)
                                                                                                  ggCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                     >  1:5249844/1‑150 (MQ=255)
                                                                                                  cgatagggCCGCCAAGCGGCCCTATCGCGACGCAAGGCCATTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                     >  1:9059909/6‑150 (MQ=39)
                                                                                                    ccGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAAAACCAACCTCATTCCTCAAACCGGATCGGGcaca                                                   >  1:6681306/1‑150 (MQ=255)
                                                                                                         gCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTcc                                              <  2:3707598/150‑1 (MQ=255)
                                                                                                         gCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTAAGCCGATCCTCTGTGGGAGCGGCCGTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCACAATCCTCAAACCGGATAGGGCACAACTcc                                              >  1:5911595/1‑150 (MQ=39)
                                                                                                         gCCGCTTGGCGGCCCTATCGAGACGCAAGGCCGTTCCCACAAACTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCATGCGTCCCGATGGGCTGCAAAGCCGACCCAGAATCACCACCCCATTCATCAAACACGGTCGGGCCCAACGcc                                              >  2:3333531/1‑150 (MQ=255)
                                                                                                           gggttGGGGGCCCTTTGGCGACGCAAGGCCGTTCCCACAAAGTACGTCTTACGCCGATCCGGTGTGGGAGGGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAc                                            <  1:8689720/147‑1 (MQ=38)
                                                                                                           cGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAc                                            <  2:2161977/150‑1 (MQ=255)
                                                                                                                  cggcCCTATCGCGACGCACGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTcc                                     <  1:6036535/150‑1 (MQ=255)
                                                                                                                     cccTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGc                                  <  2:5098685/150‑1 (MQ=255)
                                                                                                                         accgcgACGCACGGCCGTCCCCAAAAAGTACGTCGTACGCCTATCCTGTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCt                              <  1:7735384/148‑1 (MQ=39)
                                                                                                                         aTCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGGGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCt                              <  1:2012743/150‑1 (MQ=255)
                                                                                                                          tCGCGCCGCACGGCCGTTCCCACAAAGTAGGTCGTACGCCGAGCCTCTGTGGGTGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTg                             <  1:5724016/150‑1 (MQ=39)
                                                                                                                            gcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCa                           >  2:7466009/1‑150 (MQ=255)
                                                                                                                            gcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCATCCAGCTTCTGCa                           >  1:2425516/1‑150 (MQ=255)
                                                                                                                                 gCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                      >  2:7855383/1‑150 (MQ=255)
                                                                                                                                 gCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                      >  1:6240727/1‑150 (MQ=255)
                                                                                                                                 gCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                      >  2:687992/1‑150 (MQ=255)
                                                                                                                                      gCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGag                 >  1:4412029/1‑149 (MQ=255)
                                                                                                                                      gCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGAt                 >  2:4392812/1‑150 (MQ=21)
                                                                                                                                       cctttaccAAAATGTACGTCTTACGCCTATACTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATc                <  1:4767689/144‑1 (MQ=38)
                                                                                                                                       ccGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATc                >  2:6483139/1‑150 (MQ=21)
                                                                                                                                        cGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCa               >  2:4308928/1‑150 (MQ=255)
                                                                                                                                         gTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAc              >  2:3125841/1‑150 (MQ=255)
                                                                                                                                           tCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCt            >  2:579309/1‑150 (MQ=25)
                                                                                                                                             ccACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGc          >  1:4571970/1‑150 (MQ=255)
                                                                                                                                              cacaACGTCCGTGGTACGCCGATCCTCGGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg         <  1:1066657/150‑1 (MQ=255)
                                                                                                                                              cacaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg         <  1:687992/150‑1 (MQ=255)
                                                                                                                                              cacaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg         >  2:3751334/1‑150 (MQ=255)
                                                                                                                                              caaaaaGTACGTCGTACGCCTATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg         <  1:6960331/147‑1 (MQ=255)
                                                                                                                                               acaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCgg        >  2:589890/1‑150 (MQ=255)
                                                                                                                                                caAAGTACGTCGTACGCCGATCCTCTGTGGGAGTGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGc       >  1:6980441/1‑150 (MQ=255)
                                                                                                                                                caAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGc       <  2:6264000/150‑1 (MQ=255)
                                                                                                                                                    gTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTgg   >  2:8190733/1‑150 (MQ=255)
                                                                                                                                                     tACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGc  >  1:5310368/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AGCAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGC  >  NC_002947/84928‑85226

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: