Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 85,081 T→G 8.6% intergenic (+59/+71) 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,0810TG 8.6% 667.9 / 36.5 248intergenic (+59/+71)aroE/PP_0075shikimate 5‑dehydrogenase/choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base T (117/108);  new base G (14/9);  total (131/117)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 264 total aligned reads displayed.

GGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGG  >  NC_002947/84932‑85230
                                                                                                                                                       |                                                                                                                                                     
 gcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtc                                                                                                                                                      <  1:4813673/150‑1 (MQ=255)
 gcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtc                                                                                                                                                      <  1:1750298/150‑1 (MQ=255)
 gcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtc                                                                                                                                                      <  1:6939424/150‑1 (MQ=255)
   agCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtcgt                                                                                                                                                    <  1:3835379/150‑1 (MQ=255)
     cTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCCTTACCACAAACTAACTCGTAc                                                                                                                                                  >  1:3071815/1‑150 (MQ=255)
       gAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGc                                                                                                                                                >  2:4688914/1‑150 (MQ=255)
        aaGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACACAGTACGTCGTACGcc                                                                                                                                               >  1:4457411/1‑150 (MQ=255)
          gCCTTCTTTATCTGGCGTGGCGTGGGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGa                                                                                                                                             >  2:3417620/1‑150 (MQ=255)
          gCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGGTCCCACAAAGTACGTCGTACGCCGa                                                                                                                                             >  1:7082948/1‑150 (MQ=255)
           ccTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGag                                                                                                                                            >  2:7228067/1‑149 (MQ=255)
           ccTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGAt                                                                                                                                            >  2:6861898/1‑150 (MQ=255)
             ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCGAGGCCGTTCCCACAAAGTACGTCGTACGCCGATcc                                                                                                                                          >  2:5951832/1‑150 (MQ=255)
             ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATcc                                                                                                                                          >  2:7982216/1‑150 (MQ=255)
             ttcttTATCTGGCGTGGCATGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATcc                                                                                                                                          >  2:1796412/1‑150 (MQ=255)
                tttATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGTGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCctct                                                                                                                                       <  2:4003339/150‑1 (MQ=255)
                tttATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCctct                                                                                                                                       >  1:6227018/1‑150 (MQ=255)
                  tATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                     <  2:4384051/150‑1 (MQ=255)
                   aTCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgtg                                                                                                                                    <  1:8027878/150‑1 (MQ=255)
                   aTCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgtg                                                                                                                                    <  1:7415935/150‑1 (MQ=255)
                   aTCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGGTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCGTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgtg                                                                                                                                    <  2:5784884/150‑1 (MQ=255)
                    tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTgg                                                                                                                                   >  2:5381346/1‑150 (MQ=255)
                      tggcgtggcgTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGa                                                                                                                                 <  2:4774422/150‑1 (MQ=255)
                      tggcgtggcgTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGa                                                                                                                                 <  1:7041075/150‑1 (MQ=255)
                      tggcgtggcgTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGa                                                                                                                                 <  1:659959/150‑1 (MQ=255)
                       ggcgtggcgtGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACCAAATAACTCGTACGCCGATCCACAGTGGGAg                                                                                                                                >  1:4979113/1‑150 (MQ=39)
                       ggcgtggcgtGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAg                                                                                                                                <  2:2789130/150‑1 (MQ=255)
                        gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCGGAGGGAGc                                                                                                                               >  2:6551268/1‑150 (MQ=255)
                          gtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCgg                                                                                                                             <  2:6641731/150‑1 (MQ=255)
                          gtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCgg                                                                                                                             >  1:668737/1‑150 (MQ=255)
                          gtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATACTCTGTGGGAGCgg                                                                                                                             >  1:5580265/1‑150 (MQ=255)
                          gtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCCGTCCCCCAAAGTAAGGCGGACGCCGAGCCTCTGTGGGGGCgg                                                                                                                             >  1:2377830/1‑150 (MQ=255)
                              cgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCtt                                                                                                                         <  2:142328/150‑1 (MQ=255)
                              cgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCtt                                                                                                                         <  1:5263897/150‑1 (MQ=21)
                                tgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGc                                                                                                                       <  2:653298/150‑1 (MQ=255)
                                tgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGc                                                                                                                       <  1:2919430/150‑1 (MQ=255)
                                   ggCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCCACGCCGTTCCCACAAAGTACGTCGGACCCCGATCCTCTGTGGGGGCGGCCTGGCggc                                                                                                                    >  1:4735540/1‑148 (MQ=255)
                                    gCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTcg                                                                                                                   >  2:3248618/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)
                                             ggCGCCGGTGTTGGCTGAACGGCGCCGGCAACTCGCTGGGGGGTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATGGCGACGCTCGGCCTTTCCCACAAAGTACGTGGCACGCCGATCATCTGTGGGAGGGGCATTGCGACGCGATGGGCt                                                                                                          <  2:6895213/150‑1 (MQ=255)
                                             ggCGCCGGTGTGGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGGTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCACGGCCGTCCCCACAACGCACGTCGTACGCCGAGCCTCTGTGGGAGCGGCCTGGCGCCGCGATGGGCt                                                                                                          <  2:3088265/150‑1 (MQ=255)
                                                  cGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCaaa                                                                                                     <  2:7204876/150‑1 (MQ=255)
                                                    gtgtTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGc                                                                                                   >  2:4727246/1‑150 (MQ=255)
                                                    gtgtTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGc                                                                                                   <  1:8357858/150‑1 (MQ=255)
                                                    gtgtTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACCTCGTAAGCCCATCCGCAGTGGGAGCGGCCATGCGTCGCGATGGGCTGCAAAGc                                                                                                   >  1:2661244/1‑150 (MQ=39)
                                                     tgtTGGCTGAACTGCGCCGGCAACTCGCTGGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAACTTACGTCGTCGGCCGATCCGCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCa                                                                                                  <  2:259883/150‑1 (MQ=39)
                                                     tgtTGGCTGAACGGCGCCGGCAACTCGCTCGGGGCTGCGCTGGCCGGGGCCGCTTGGCGGCCCTAGCGCGACGCAAGGCCGTCCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCa                                                                                                  <  2:259996/150‑1 (MQ=255)
                                                        tGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAc                                                                                               >  2:8295618/1‑150 (MQ=255)
                                                          gCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAccc                                                                                             >  2:679359/1‑150 (MQ=255)
                                                          gCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCATGCGTCGCGATGGGCTGCAAAGCAGAccc                                                                                             >  1:566706/1‑150 (MQ=255)
                                                              aaCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGaa                                                                                         >  1:9091206/1‑150 (MQ=255)
                                                               aCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTCCGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAAt                                                                                        <  2:90578/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)
                                                                            ctcgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGAGCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCTACCTCATTc                                                                           <  1:5237199/150‑1 (MQ=255)
                                                                             tcgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGGCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCAATcc                                                                          >  1:7600802/1‑150 (MQ=255)
                                                                             tcgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGGTACCACAAAGTACGTCGTACGCCCATCCGCGGTGGGAGCGGCCGTGGGGCGCGATGGGGTGCAAAGCAGACCCACAAACCCAACcacaatac                                                                          >  1:4797985/1‑146 (MQ=37)
                                                                              cgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCt                                                                         <  2:8805617/150‑1 (MQ=255)
                                                                              cgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCt                                                                         <  2:4200383/150‑1 (MQ=255)
                                                                                  gggggCTGTGCTGGCGGGGGCGGCTTGGCGGCCCTATGGCGACGCAAGGCGTTCCCAACAAAGTAGGTCTTCGGCCGTGCCTGTGGGGGTGCGGCCTGGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCaaa                                                                     <  1:248448/149‑1 (MQ=255)
                                                                                       cTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGa                                                                <  2:1349190/150‑1 (MQ=255)
                                                                                            cTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCggg                                                           <  2:8042658/150‑1 (MQ=255)
                                                                                              ggCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                         >  1:5249844/1‑150 (MQ=255)
                                                                                               gCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGcac                                                        >  2:1136388/1‑150 (MQ=255)
                                                                                                 ggcggCCGCGTGGCGGCCCTAACGCGCCGCAAGGCCGTTCCCCCAACCCCCGTCGTACGCCGCTCCTCCGGGGGAGCGGCCTCGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCCAACCGGCTCGGGcacac                                                      <  2:8250409/147‑2 (MQ=255)
                                                                                                     gCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTcc                                                  <  2:1637959/150‑1 (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:567712/1‑150 (MQ=255)
                                                                                                       cGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAc                                                <  2:2161977/150‑1 (MQ=255)
                                                                                                         cTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACCAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAccc                                              >  2:7573765/1‑150 (MQ=255)
                                                                                                          ttGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCGGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCACCCTCATGCCGCAAACCGGAGCGGGCACAACTCCAcccc                                             >  2:5844715/1‑150 (MQ=39)
                                                                                                              cggcCCTATCGCGACGCAAGGCCGTTCCCACAATGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCAGTCCTCAAACCGGATCGGGCACAACTCCACCCCTTcc                                         <  2:6338012/150‑1 (MQ=255)
                                                                                                              cggcCCTATCGCGACGAAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTcc                                         <  1:8923919/150‑1 (MQ=255)
                                                                                                               ggcCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTGTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCa                                        <  1:6033632/150‑1 (MQ=255)
                                                                                                                 cccTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGc                                      <  2:5098685/150‑1 (MQ=255)
                                                                                                                  ccTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCt                                     >  2:5326342/1‑150 (MQ=255)
                                                                                                                      ccgctACGAACGGCCGTCCCCAAAAAGTACGTCGTACGCCGATCCTCTGTGGGTGGGGCCTTGCGGCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTg                                 <  1:3895118/149‑1 (MQ=255)
                                                                                                                        gcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCa                               >  2:7466009/1‑150 (MQ=255)
                                                                                                                             gCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                          >  2:8628476/1‑150 (MQ=255)
                                                                                                                               aaGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcctc                        <  2:2198448/150‑1 (MQ=255)
                                                                                                                                 ggCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGa                      <  2:8104801/150‑1 (MQ=255)
                                                                                                                                  gCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGag                     >  1:4412029/1‑149 (MQ=255)
                                                                                                                                   ccGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATc                    >  2:6143587/1‑150 (MQ=21)
                                                                                                                                   ccGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATc                    >  1:3611126/1‑150 (MQ=21)
                                                                                                                                        cccACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGAAGACCCATAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTg               >  1:5554170/1‑150 (MQ=255)
                                                                                                                                         ccACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGc              >  2:4523642/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)
                                                                                                                                             caAAGTACGTCGTACGCCGATCCTCTGTGGGAGTGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGc           >  1:6980441/1‑150 (MQ=255)
                                                                                                                                              aaaGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGcc          >  2:3107954/1‑150 (MQ=255)
                                                                                                                                                   aCGTCGTCCGCCTATCCTCTGTGGGTGGGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGcc     <  1:1841830/150‑1 (MQ=255)
                                                                                                                                                      gtcgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCg   <  2:8775611/150‑1 (MQ=255)
                                                                                                                                                       tcgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCgg  <  1:4493468/150‑1 (MQ=255)
                                                                                                                                                       tcgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCgg  >  1:4657193/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GGCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGG  >  NC_002947/84932‑85230

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: