Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,054,296 A→C 6.6% 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 6.6% 694.5 / 27.3 208intergenic (‑33/+27)PP_1822/folEA‑Ihypothetical protein/GTP cyclohydrolase I
Reads supporting (aligned to +/- strand):  ref base A (98/95);  new base C (4/10);  total (102/106)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 213 total aligned reads displayed.

GCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACATGTG  >  NC_002947/2054148‑2054445
                                                                                                                                                    |                                                                                                                                                     
gCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGaa                                                                                                                                                      <  1:5139527/150‑1 (MQ=255)
gCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGaa                                                                                                                                                      <  1:79011/150‑1 (MQ=255)
 cTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAg                                                                                                                                                     <  2:4006124/150‑1 (MQ=255)
    gCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGt                                                                                                                                                  >  1:1693369/1‑150 (MQ=255)
      tCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCg                                                                                                                                                <  1:1338780/150‑1 (MQ=255)
      tCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCg                                                                                                                                                >  2:1522390/1‑150 (MQ=255)
      tCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCg                                                                                                                                                >  2:2286237/1‑150 (MQ=255)
      tCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCg                                                                                                                                                >  2:2904820/1‑150 (MQ=255)
           gCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGCttt                                                                                                                                           <  2:146050/150‑3 (MQ=255)
           gCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGCttt                                                                                                                                           <  1:4727062/150‑3 (MQ=255)
            cGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCTGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGCtttt                                                                                                                                          <  2:3128122/150‑4 (MQ=255)
            cGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCgg                                                                                                                                          <  1:3680217/150‑1 (MQ=255)
            cGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCgg                                                                                                                                          <  2:5181861/150‑1 (MQ=255)
            cGGGCCCGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGCtttt                                                                                                                                          <  1:4136958/150‑4 (MQ=255)
            cGGGCCCGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGCtttt                                                                                                                                          <  1:3178794/150‑4 (MQ=255)
              ggCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTcttgc                                                                                                                                        <  2:674358/150‑5 (MQ=255)
               gCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCTGTCACCGATGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGc                                                                                                                                       <  2:710239/150‑1 (MQ=255)
                 cTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCa                                                                                                                                     <  1:4447636/150‑1 (MQ=255)
                 cTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCa                                                                                                                                     <  2:2440941/150‑1 (MQ=255)
                     gcggcggcgTCGCGTTTGCGCTGGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGCTCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCaa                                                                                                                                 <  2:3736927/150‑1 (MQ=255)
                      cggcggcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAAc                                                                                                                                <  1:2662736/150‑1 (MQ=255)
                      cggcggcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAAc                                                                                                                                <  1:5264607/150‑1 (MQ=255)
                       ggcggcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAAcc                                                                                                                               <  2:5530772/150‑1 (MQ=255)
                               cgcgTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACtt                                                                                                                       <  1:4982152/150‑1 (MQ=255)
                               cgcgTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACtt                                                                                                                       <  1:3731708/150‑1 (MQ=255)
                                   tttGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATc                                                                                                                   >  2:1685581/1‑150 (MQ=255)
                                    ttGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCa                                                                                                                  <  2:4668157/150‑1 (MQ=255)
                                     tGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAg                                                                                                                 <  2:2579751/150‑1 (MQ=255)
                                        gcTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCt                                                                                                              >  2:3460461/1‑150 (MQ=255)
                                          ttGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGa                                                                                                            <  1:5254304/150‑1 (MQ=255)
                                          ttGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGa                                                                                                            <  1:3995700/150‑1 (MQ=255)
                                          ttGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGa                                                                                                            <  1:940737/150‑1 (MQ=255)
                                           tGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAg                                                                                                           <  2:3979544/150‑1 (MQ=255)
                                              ccGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGaa                                                                                                        >  1:4350156/1‑150 (MQ=255)
                                              ccGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGaa                                                                                                        <  2:1336660/150‑1 (MQ=255)
                                               cGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAAc                                                                                                       >  2:2604395/1‑150 (MQ=255)
                                                  gcggcTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCg                                                                                                    <  1:986812/150‑1 (MQ=255)
                                                     gcTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTg                                                                                                 <  2:1230827/150‑1 (MQ=255)
                                                      cTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGc                                                                                                >  1:5204061/1‑150 (MQ=255)
                                                           ccAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTg                                                                                           <  1:159175/150‑1 (MQ=255)
                                                             aGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTggc                                                                                         >  1:2292017/1‑150 (MQ=255)
                                                             aGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTggc                                                                                         >  1:3244251/1‑150 (MQ=255)
                                                              gTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTggcg                                                                                        <  1:1441279/150‑1 (MQ=255)
                                                                cACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTggcggc                                                                                      <  2:1280561/150‑1 (MQ=255)
                                                                   cGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGtt                                                                                   >  2:961792/1‑150 (MQ=255)
                                                                    gAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTg                                                                                  >  2:5081051/1‑150 (MQ=255)
                                                                     aaGACGATGAGCGGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGt                                                                                 <  1:2130119/150‑1 (MQ=255)
                                                                        aCGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTcgc                                                                              <  2:2334426/150‑1 (MQ=255)
                                                                          gATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCgg                                                                            >  2:1371012/1‑150 (MQ=255)
                                                                             gAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAAc                                                                         <  1:2133305/150‑1 (MQ=255)
                                                                                  ggCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCAcc                                                                    <  2:5429816/150‑1 (MQ=255)
                                                                                       aGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCa                                                               <  2:4532922/150‑1 (MQ=255)
                                                                                         gCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCAta                                                             >  1:2439479/1‑149 (MQ=255)
                                                                                          ccAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCa                                                            <  1:3895899/150‑1 (MQ=255)
                                                                                          ccAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCa                                                            >  1:2115267/1‑150 (MQ=255)
                                                                                           cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTg                                                                                                 <  2:877888/112‑1 (MQ=255)
                                                                                           cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTg                                                                                                 >  1:877888/1‑112 (MQ=255)
                                                                                           cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAc                                                           >  2:1890234/1‑150 (MQ=255)
                                                                                           cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAc                                                           >  2:5617405/1‑150 (MQ=255)
                                                                                            aaCCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAcc                                                          >  1:1929536/1‑150 (MQ=255)
                                                                                              cccGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGa                                                        <  1:3336167/150‑1 (MQ=255)
                                                                                               ccGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAg                                                       <  1:1776633/150‑1 (MQ=255)
                                                                                               ccGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAg                                                       >  2:4070458/1‑150 (MQ=255)
                                                                                               ccGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAg                                                       >  2:501772/1‑150 (MQ=255)
                                                                                                cGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAgg                                                      <  1:4847555/150‑1 (MQ=255)
                                                                                                cGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAgg                                                      <  1:1896671/150‑1 (MQ=255)
                                                                                                 gCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGt                                                     >  2:1738193/1‑150 (MQ=255)
                                                                                                  cAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTg                                                    <  1:3001946/150‑1 (MQ=255)
                                                                                                   aGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGa                                                   >  1:5397486/1‑150 (MQ=255)
                                                                                                        tttGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATg                                              >  1:1183247/1‑150 (MQ=255)
                                                                                                          tGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGt                                            >  2:3951129/1‑150 (MQ=255)
                                                                                                              gATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGaa                                        >  1:3309939/1‑150 (MQ=255)
                                                                                                              gATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGaa                                        <  1:5077020/150‑1 (MQ=255)
                                                                                                                 cATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTc                                     >  1:2554570/1‑150 (MQ=255)
                                                                                                                  aTGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAGGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCt                                    <  1:3579389/150‑1 (MQ=255)
                                                                                                                    gggCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGc                                  <  1:4634260/150‑1 (MQ=255)
                                                                                                                      gCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCtt                                >  2:5258612/1‑150 (MQ=255)
                                                                                                                        aaCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtct                              <  1:3716076/150‑1 (MQ=255)
                                                                                                                           tCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGa                           >  2:3306545/1‑150 (MQ=255)
                                                                                                                              tGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGacac                        >  2:3942534/1‑150 (MQ=255)
                                                                                                                              tGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGAGCCTGGTCGACTTCTGCTGCTCGacac                        >  2:3353883/1‑150 (MQ=255)
                                                                                                                               gAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACAcc                       >  1:4116495/1‑150 (MQ=255)
                                                                                                                                aaaTGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACcg                      >  2:531629/1‑150 (MQ=255)
                                                                                                                                      ccgaGCCGACCCTTCAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc                >  1:5334119/4‑150 (MQ=255)
                                                                                                                                      aattgCCGAGCCTGAAGGGTCGGCTCGGGGCCAGCAAGCAATCACTTGATGAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACGGAGGTGATCATGGTCGAATTCTGCTTCTCGACAGCGCGcatc                >  2:3559709/5‑150 (MQ=255)
                                                                                                                                      aaaaGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc                >  1:1293383/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc                >  2:3686297/1‑150 (MQ=255)
                                                                                                                                        gaGCCGACCCTTCAGGGTCGGCTTGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcat              >  1:812657/2‑150 (MQ=255)
                                                                                                                                          gCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatca            <  2:3015790/150‑1 (MQ=255)
                                                                                                                                           ccGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatcat           <  2:2163083/150‑1 (MQ=255)
                                                                                                                                           ccGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatcat           <  2:1787913/150‑1 (MQ=255)
                                                                                                                                           ccGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatcat           >  1:1530228/1‑150 (MQ=255)
                                                                                                                                            cGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATg          >  1:3950119/1‑150 (MQ=255)
                                                                                                                                            cGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATg          >  1:125445/1‑150 (MQ=255)
                                                                                                                                            cGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATg          >  1:1211482/1‑150 (MQ=255)
                                                                                                                                                ccTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGcaca      >  1:710052/1‑150 (MQ=255)
                                                                                                                                                  tGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACAtg    >  1:637688/1‑150 (MQ=255)
                                                                                                                                                  tGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACAtg    >  1:4586823/1‑150 (MQ=255)
                                                                                                                                                   gAAGGGTCGGCTCGGGGCCAGCAACCACTCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACAtgt   <  2:3935465/150‑1 (MQ=255)
                                                                                                                                                   gAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACAtgt   <  2:5433938/150‑1 (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: