Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,054,296 A→C 5.8% intergenic (‑33/+27) PP_1822 ← / ← folEA‑I hypothetical protein/GTP cyclohydrolase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,054,2960AC 5.8% 1032.4 / 33.9 306intergenic (‑33/+27)PP_1822/folEA‑Ihypothetical protein/GTP cyclohydrolase I
Reads supporting (aligned to +/- strand):  ref base A (128/159);  new base C (3/15);  total (131/175)
Fisher's exact test for biased strand distribution p-value = 2.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 309 total aligned reads displayed.

GCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACATGTG  >  NC_002947/2054148‑2054445
                                                                                                                                                    |                                                                                                                                                     
gCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGaa                                                                                                                                                      <  1:7452972/150‑1 (MQ=255)
gCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGaa                                                                                                                                                      <  1:4770787/150‑1 (MQ=255)
  ttGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAgg                                                                                                                                                    <  1:1217351/150‑1 (MQ=255)
  ttGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAgg                                                                                                                                                    >  2:7904546/1‑150 (MQ=255)
  ttGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAgg                                                                                                                                                    <  1:2719600/150‑1 (MQ=255)
  ttGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAgg                                                                                                                                                    >  2:751709/1‑150 (MQ=255)
    gCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGt                                                                                                                                                  <  2:3483326/150‑1 (MQ=255)
    gCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGt                                                                                                                                                  <  1:1957550/150‑1 (MQ=255)
     cTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTc                                                                                                                                                 <  2:5690323/150‑1 (MQ=255)
     cTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTc                                                                                                                                                 <  2:7634588/150‑1 (MQ=255)
        gATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGc                                                                                                                                              <  1:7027611/150‑1 (MQ=255)
        gATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGc                                                                                                                                              <  2:2608303/150‑1 (MQ=255)
          tGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTc                                                                                                                                            <  1:5471993/150‑1 (MQ=255)
           gCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCttt                                                                                                                                           <  2:7343578/150‑3 (MQ=255)
            cGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGCtttt                                                                                                                                          <  1:4208129/150‑4 (MQ=255)
            cGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCgg                                                                                                                                          <  2:4285195/150‑1 (MQ=255)
            cGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCgg                                                                                                                                          <  2:6084173/150‑1 (MQ=255)
             gggCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGCttttg                                                                                                                                         <  2:6532010/150‑5 (MQ=255)
             gggCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCCGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCggg                                                                                                                                         <  2:8005151/150‑1 (MQ=255)
              ggCGTGGGCGGCGGCGTCGCGTTTGCGGTTGGGCGGGCGGCTGATCCAGTCACCGAAGACGTTGAGCTGGCCGAGGCCAACCCGAAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGCttttgc                                                                                                                                        <  1:8308426/150‑6 (MQ=255)
              ggCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCgggg                                                                                                                                        >  1:3836515/1‑150 (MQ=255)
                     gcggcggcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCaa                                                                                                                                 >  1:632956/1‑150 (MQ=255)
                      cggcggcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGc                                                                                                                                       <  2:7755388/143‑1 (MQ=255)
                          ggcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAAt                                                                                                                            >  1:1826390/1‑150 (MQ=255)
                          ggcgTCGCGTTTGCGCTTGGCCGGGCGACTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAAt                                                                                                                            <  2:2822396/150‑1 (MQ=255)
                             gTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCAc                                                                                                                         <  2:4521825/150‑1 (MQ=255)
                              tCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACt                                                                                                                        >  1:2430048/1‑150 (MQ=255)
                               cgcgTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACtt                                                                                                                       <  2:4647902/150‑1 (MQ=255)
                               cgcgTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACtt                                                                                                                       <  1:671636/150‑1 (MQ=255)
                                 cgTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGa                                                                                                                     <  1:1408464/150‑1 (MQ=255)
                                      gggcTTGGCCGGGCGGCTGATCCAGTAACCGAAGACGTTGAGCTGGCCGTGGCCAACCCGCTTGGCTTTGACGTTCTTGGGCATCTCCTGATATGCAAAAGCCGACCGTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAgg                                                                                                                <  1:4066286/148‑1 (MQ=255)
                                          ttGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGa                                                                                                            >  2:1355359/1‑150 (MQ=255)
                                          ttGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGa                                                                                                            <  2:249099/150‑1 (MQ=255)
                                           tGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAg                                                                                                           <  1:6465750/150‑1 (MQ=255)
                                                     gcTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTg                                                                                                 >  1:7246962/1‑150 (MQ=255)
                                                         aTCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCggg                                                                                             <  1:5059201/150‑1 (MQ=255)
                                                          tCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGt                                                                                            <  1:7843644/150‑1 (MQ=255)
                                                             aGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTggc                                                                                         >  1:1345052/1‑150 (MQ=255)
                                                                cACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCACCGCCTGAACTGCACAAGCCGACCCTGACGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTggcggc                                                                                      <  2:377502/150‑1 (MQ=255)
                                                                cACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTggcggc                                                                                      <  2:3806213/150‑1 (MQ=255)
                                                                 aCCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGgcggcg                                                                                     >  1:7154641/1‑150 (MQ=255)
                                                                   cGAAGACGATGAGCTGGCCGAGGCCAACCCACAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGtt                                                                                   <  1:52164/150‑1 (MQ=255)
                                                                    gAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTg                                                                                  >  1:1020609/1‑150 (MQ=255)
                                                                    gAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTg                                                                                  >  1:3316612/1‑150 (MQ=255)
                                                                     aaGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAATTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGt                                                                                 >  2:1824826/1‑150 (MQ=255)
                                                                      aGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTc                                                                                >  2:3891919/1‑150 (MQ=255)
                                                                        ccGATGAGCTGGCCGAGGCCAACCCGCAGGGCTGGGACGCTCCTGGGCAACTCCTGAACGGCAACAGCCGACCCTGACGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGACCTCGCTGCGGGTGGCGGCGTTGTcgc                                                                              <  2:2944408/149‑1 (MQ=255)
                                                                        aCGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTcgc                                                                              <  1:3911994/150‑1 (MQ=255)
                                                                          gATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCgg                                                                            <  2:3719249/150‑1 (MQ=255)
                                                                            tGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGaa                                                                          <  2:265514/150‑1 (MQ=255)
                                                                               gCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTc                                                                       >  1:8596219/1‑150 (MQ=255)
                                                                                cTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCa                                                                      <  1:5026830/150‑1 (MQ=255)
                                                                                 tGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCAc                                                                     >  2:23401/1‑150 (MQ=255)
                                                                                    ccGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCa                                                                  >  1:2352646/1‑150 (MQ=255)
                                                                                     cGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAg                                                                 <  1:8671745/150‑1 (MQ=255)
                                                                                      gAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGc                                                                <  2:7861511/150‑1 (MQ=255)
                                                                                        ggCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCAt                                                              >  1:8736804/1‑150 (MQ=255)
                                                                                           cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAc                                                           >  1:495587/1‑150 (MQ=255)
                                                                                           cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAc                                                           >  1:6447635/1‑150 (MQ=255)
                                                                                            aaCCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAcc                                                          >  1:2878246/1‑150 (MQ=255)
                                                                                            aaCCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAcc                                                          <  1:1513184/150‑1 (MQ=255)
                                                                                               ccGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAg                                                       <  2:2128586/150‑1 (MQ=255)
                                                                                               ccGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAg                                                       <  1:2856742/150‑1 (MQ=255)
                                                                                               ccGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAg                                                       <  1:5961916/150‑1 (MQ=255)
                                                                                                cGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAgg                                                      >  2:5222587/1‑150 (MQ=255)
                                                                                                cGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAgg                                                      <  1:1727186/150‑1 (MQ=255)
                                                                                                 gCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGt                                                     >  1:4430337/1‑150 (MQ=255)
                                                                                                   aGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATTAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGa                                                   >  1:637252/1‑150 (MQ=255)
                                                                                                     ggCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATc                                                 <  1:3019866/150‑1 (MQ=255)
                                                                                                         ttGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATgg                                             <  2:8380626/150‑1 (MQ=255)
                                                                                                          tGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGt                                            <  1:816151/150‑1 (MQ=255)
                                                                                                           gACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTc                                           >  2:1620622/1‑150 (MQ=255)
                                                                                                              gATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGaa                                        >  1:5509275/1‑150 (MQ=255)
                                                                                                              gATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGaa                                        >  1:5756322/1‑150 (MQ=255)
                                                                                                                tCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGactg                                      >  2:8436400/1‑147 (MQ=255)
                                                                                                                   tGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTg                                   >  1:3076345/1‑150 (MQ=255)
                                                                                                                    gggCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGc                                  <  2:8200673/150‑1 (MQ=255)
                                                                                                                        aaCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtct                              <  2:2999945/150‑1 (MQ=255)
                                                                                                                           tCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGa                           <  2:2576012/150‑1 (MQ=255)
                                                                                                                           tCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGa                           >  2:6878519/1‑150 (MQ=255)
                                                                                                                            ccTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGac                          <  2:1428158/150‑1 (MQ=255)
                                                                                                                              tGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGacac                        >  2:1080035/1‑150 (MQ=255)
                                                                                                                              tGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGacac                        >  2:640943/1‑150 (MQ=255)
                                                                                                                                   tGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCgcgc                   <  2:75599/150‑1 (MQ=255)
                                                                                                                                      ccgaGCCGACCCTTCAGGGTCGGCTTGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc                >  2:269446/4‑150 (MQ=255)
                                                                                                                                      aaaaGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc                >  2:8626295/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc                >  1:4872906/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc                >  1:7744953/1‑150 (MQ=255)
                                                                                                                                         aGCCGACCCTTCAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatc             >  2:8096382/1‑150 (MQ=255)
                                                                                                                                         aGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatc             <  2:6356610/150‑1 (MQ=255)
                                                                                                                                         aGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatc             <  1:2159273/150‑1 (MQ=255)
                                                                                                                                          gCCGACCCTGTAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatca            <  1:8038941/150‑1 (MQ=255)
                                                                                                                                          gCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatca            <  1:3683769/150‑1 (MQ=255)
                                                                                                                                           ccGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatcat           <  2:4829593/150‑1 (MQ=255)
                                                                                                                                            cGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCAGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATg          >  2:6709864/1‑150 (MQ=255)
                                                                                                                                             gACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGACCTCACCCAGCATCACCGAGGTGATCAGGGTCGAATGCTGCTTCTCGACACCGCGCATCATCATGc         >  2:4070129/1‑150 (MQ=255)
                                                                                                                                             gACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGc         >  2:4124159/1‑150 (MQ=255)
                                                                                                                                                ccTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGcaca      >  2:7035174/1‑150 (MQ=255)
                                                                                                                                                ccTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGAGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGcaca      >  1:6633517/1‑150 (MQ=255)
                                                                                                                                                   gAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACAtgt   >  2:1977668/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACATGTG  >  NC_002947/2054148‑2054445

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