Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 85,078:1 +T 8.5% intergenic (+56/+74) 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,0781.T 8.5% 1202.7 / 102.6 281intergenic (+56/+74)aroE/PP_0075shikimate 5‑dehydrogenase/choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base . (131/126);  new base T (10/14);  total (141/140)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-02

Only 100 of 291 total aligned reads displayed.

CAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGC  >  NC_002947/84930‑85226
                                                                                                                                                     |                                                                                                                                                     
caggcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAc                                                                                                                                                      <  1:3530760/150‑1 (MQ=255)
caggcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAc                                                                                                                                                      <  2:2617638/150‑1 (MQ=255)
  ggcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGTTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgt                                                                                                                                                    <  2:3191182/150‑1 (MQ=255)
   gcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtc                                                                                                                                                   <  1:5183768/150‑1 (MQ=255)
   gcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtc                                                                                                                                                   <  1:2527153/150‑1 (MQ=255)
        tGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCCATGCCGTTCCCACAAATTAAGTAGGACg                                                                                                                                              >  1:1306771/1‑150 (MQ=255)
        tGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACg                                                                                                                                              >  2:4918535/1‑150 (MQ=255)
        tGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACg                                                                                                                                              <  2:4767385/150‑1 (MQ=255)
         gAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACCAAGTACGTCGTACGc                                                                                                                                             >  2:1958159/1‑150 (MQ=255)
               ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGAAcc                                                                                                                                       >  1:6510890/1‑150 (MQ=255)
                tcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATcca                                                                                                                                      >  2:6667146/1‑149 (MQ=255)
                 cttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATccac                                                                                                                                     >  1:4394002/1‑148 (MQ=255)
                  tttATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCctct                                                                                                                                    >  1:4884457/1‑150 (MQ=255)
                    tATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGGAAGGGCGTTGGCACAAAGTACGTCGTTCGCCGGTGctcttt                                                                                                                                  >  1:4984986/1‑148 (MQ=255)
                    tATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                  <  2:1424006/150‑1 (MQ=255)
                      tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGTCTGCTTTGCAGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTgg                                                                                                                                <  2:7045477/150‑1 (MQ=255)
                      tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTgg                                                                                                                                <  1:1182465/150‑1 (MQ=255)
                       cTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTggg                                                                                                                               >  1:422504/1‑150 (MQ=255)
                        tggcgtggcgTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGa                                                                                                                              <  1:3026446/150‑1 (MQ=255)
                         ggcgtggcgtGTGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGTGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAg                                                                                                                             <  1:6967153/150‑1 (MQ=255)
                          gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGc                                                                                                                            >  2:542277/1‑150 (MQ=255)
                          gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGc                                                                                                                            >  2:100205/1‑150 (MQ=255)
                          gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCCGTCCCACAAAGTACGTCGTACGCCCATACACAGTGGGAGc                                                                                                                            >  1:3915096/1‑150 (MQ=39)
                           cgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACgg                                                                                                                          <  2:3318920/150‑1 (MQ=17)
                           cgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCg                                                                                                                           >  1:7247161/1‑150 (MQ=255)
                              ggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGcc                                                                                                                        <  1:6764586/150‑1 (MQ=255)
                              ggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCt                                                                                                                       <  1:1263192/150‑1 (MQ=18)
                                cgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCtt                                                                                                                      <  2:2410208/150‑1 (MQ=255)
                                  tgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGc                                                                                                                    <  2:5055909/150‑1 (MQ=255)
                                     ggCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTc                                                                                                                 <  1:3695177/150‑1 (MQ=255)
                                      gCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTcg                                                                                                                >  2:4282755/1‑150 (MQ=255)
                                      gCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTcg                                                                                                                >  2:6476408/1‑150 (MQ=255)
                                       ccGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTcgc                                                                                                               <  2:3917794/150‑1 (MQ=18)
                                       ccGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCAGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTcgc                                                                                                               <  1:1226579/150‑1 (MQ=255)
                                         gtaTAGGGGGCGGGTGTTGGCTGAACTGCGCCGGAAACTGGCTCGGGGGTGAGCTGGCGGGGGCCGCTTGGCGGCCATATCGCGACGAAAGGCGGTTCCCACAAAGTACGTCTTACGCCGATCATCTGTGGGAGCGGCCTTGCGTCGCGa                                                                                                             <  1:1596967/148‑1 (MQ=255)
                                         ggATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGa                                                                                                             <  2:5268636/150‑1 (MQ=255)
                                         ggATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGa                                                                                                             >  1:2988206/1‑150 (MQ=255)
                                          gATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGAt                                                                                                            >  2:2196518/1‑150 (MQ=255)
                                             aCGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATggg                                                                                                        <  2:600375/150‑1 (MQ=255)
                                              cGGCGCCGGTGTTGGCTGACCGGCGCCGGAAACTCGCTCGGGGCTGAGGTGGCCGGGGCCGCTTGGCGGCCCTAGCGCGACGCATGGCCGCTCCCACAGAGGATCGGCGACCGACGTACTTTGTGGTAACGGCCTTGCGTCGCGATaggg                                                                                                       <  2:6234825/150‑2 (MQ=9)
                                                    cGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCaaa                                                                                                 >  2:1538279/1‑150 (MQ=255)
                                                    cGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCaaa                                                                                                 <  2:514522/150‑1 (MQ=255)
                                                     ggTGTTGGCTGAACTGCGCCGGCAACTCGCTGGGGGCTGAGCTGGCGGGGGCCGCTTGGCGGCCCTAGCGCGCCGCACGGCGGTTCCCACAACGTACGTCGTAGGCCGATCCTCTGTGGGCGGGGCCTTGCGTCGCGATGGGCTGCAAAg                                                                                                <  1:2868488/150‑1 (MQ=255)
                                                      gtgtTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGCGGATCGGCGTACGCCGTACTTTGTGGGAACGGCCTTGCGTCGCGATAGGGCCGCCAAGc                                                                                               <  2:1756117/150‑1 (MQ=11)
                                                          tGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTGGGCGGCCCTATGGCGACGCAAGGCCGTTCCCACAAAGTCCGTCGAACGCCGATCCCCTGTGGGAGCGGCCTTGGGTCGCGATGGGCTGCAAAGCAGAc                                                                                           <  2:1575254/150‑1 (MQ=39)
                                                           ggCTGACCTGCGCCGGCACCGCGCTCGGGGCTGAGCTGGCGGGGGCCGCTTGGCGGCCCTATCGCGCCGCAAGGCGGTTCCCACAAAGTACGCCGTACGCCGAACCTCTGTGGGAGCGGCCTTGCGCCGCGATGGGCTGCAAAGCAGAcc                                                                                          <  2:2524305/150‑1 (MQ=255)
                                                                  cTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATc                                                                                   <  2:3442367/150‑1 (MQ=255)
                                                                  cTGCGCCGGCAACTCGCTCGGGGCTGAGCGGGCCGGGGCCGCTTGGCGGCGCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATc                                                                                   <  2:6488043/150‑1 (MQ=255)
                                                                              ctcgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGcc                                                                                                                                            >  2:2721437/1‑82 (MQ=255)
                                                                               tcgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTGGCACAATGTAAGTCGTACGCCGTTGCTCTGTGGGAGGGGGCTTGGGTGGCGATGGGGTGCTAAGCAGTCCCTGTATGCGAACCTCAttgg                                                                      >  2:5822900/1‑148 (MQ=255)
                                                                               tcgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACACAGTACGTCGTACGCCGATCCTCTGTGGGTGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTcc                                                                      >  2:4840940/1‑150 (MQ=255)
                                                                                          tGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGAt                                                           >  2:3817170/1‑150 (MQ=255)
                                                                                           gAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATc                                                          <  2:6522315/150‑1 (MQ=255)
                                                                                               tGGCGGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCGGTGGGAGCGGCCTTGCGACGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGc                                                      <  1:5252735/150‑1 (MQ=255)
                                                                                                ggCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                     <  2:7320920/150‑1 (MQ=255)
                                                                                                ggCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                     <  1:1814928/150‑1 (MQ=255)
                                                                                                 gCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGcac                                                    >  2:5556551/1‑150 (MQ=255)
                                                                                                        ccGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGCGCGGCCGTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCa                                             >  2:2231679/1‑150 (MQ=255)
                                                                                                           cTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGCTCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAccc                                          >  1:6179127/1‑150 (MQ=255)
                                                                                                            gtGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAcccc                                         <  2:2814861/149‑1 (MQ=255)
                                                                                                                   cccTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGc                                  <  1:7288326/150‑1 (MQ=255)
                                                                                                                   cccTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGc                                  <  1:5445273/150‑1 (MQ=255)
                                                                                                                     cTATGGCTACGCAAGGCCGGTCCCACAAAGTACGTCGTACGCCGTTCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCtt                                <  1:7419995/150‑1 (MQ=255)
                                                                                                                      caTCGCGACGCAAGGCCGCTCCCACAGAGGATCGGAGTACTACGTACTTTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCCCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTc                               >  2:1317561/2‑150 (MQ=18)
                                                                                                                        tCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTg                             <  2:3702751/150‑1 (MQ=255)
                                                                                                                         cgcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                            <  1:2953667/150‑1 (MQ=255)
                                                                                                                          gcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCa                           >  2:2146782/1‑150 (MQ=255)
                                                                                                                          gcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCa                           >  1:5911210/1‑150 (MQ=255)
                                                                                                                          gcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCa                           >  1:2349879/1‑150 (MQ=255)
                                                                                                                           cgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCaa                          >  2:6833513/1‑150 (MQ=255)
                                                                                                                            gACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAAc                         >  1:66103/1‑150 (MQ=255)
                                                                                                                            gACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAAc                         >  1:4418609/1‑150 (MQ=255)
                                                                                                                            gACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAAc                         >  2:4718024/1‑150 (MQ=255)
                                                                                                                               gCAAGGCCGTCCCCAAAAAGTACGTCTTACGCCGACCCTGTGGGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                      <  1:2700012/150‑1 (MQ=39)
                                                                                                                                 aaGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcctc                    >  1:4539219/1‑150 (MQ=255)
                                                                                                                                 aaGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcctc                    >  1:321144/1‑150 (MQ=255)
                                                                                                                                 aaGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcctc                    >  2:3363747/1‑150 (MQ=21)
                                                                                                                                  aGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCg                   >  1:3611570/1‑150 (MQ=255)
                                                                                                                                  aGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCg                   <  1:1152225/150‑1 (MQ=255)
                                                                                                                                  aGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCg                   >  2:890311/1‑150 (MQ=21)
                                                                                                                                       gTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAc              >  1:6295914/1‑150 (MQ=255)
                                                                                                                                       gTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAc              <  2:3618997/150‑1 (MQ=255)
                                                                                                                                        ttCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAcc             >  2:5899022/1‑150 (MQ=255)
                                                                                                                                         tCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGa                  <  1:2088114/144‑1 (MQ=25)
                                                                                                                                         tCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGa                  >  2:2088114/1‑144 (MQ=25)
                                                                                                                                           ccACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTcc                                                                      >  2:3543833/1‑90 (MQ=255)
                                                                                                                                           ccACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                            <  2:6304622/132‑1 (MQ=255)
                                                                                                                                           ccACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                            >  1:6304622/1‑132 (MQ=255)
                                                                                                                                           ccACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGc          >  2:1139257/1‑150 (MQ=255)
                                                                                                                                            cacaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg         >  2:6193118/1‑150 (MQ=255)
                                                                                                                                            cacaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg         <  2:2988206/150‑1 (MQ=255)
                                                                                                                                             acaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATc                <  1:3323789/142‑1 (MQ=255)
                                                                                                                                             acaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCgg        <  2:1306771/150‑1 (MQ=255)
                                                                                                                                             acaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCgg        >  2:4512355/1‑150 (MQ=255)
                                                                                                                                               aaaGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCaa                          <  2:2019331/130‑1 (MQ=255)
                                                                                                                                                 aGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTg    <  2:7165874/150‑1 (MQ=255)
                                                                                                                                                  gTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTgg   >  2:4157522/1‑150 (MQ=255)
                                                                                                                                                   tACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGc  >  2:4577939/1‑150 (MQ=255)
                                                                                                                                                   tACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGc  >  2:1856136/1‑150 (MQ=255)
                                                                                                                                                   tACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGc  >  1:4647059/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGC  >  NC_002947/84930‑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: