Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,051,268 G→T 6.0% intergenic (‑32/+28) 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,051,2680GT 6.0% 442.3 / 12.0 133intergenic (‑32/+28)PP_1822/folEA‑Ihypothetical protein/GTP cyclohydrolase I
Reads supporting (aligned to +/- strand):  ref base G (64/61);  new base T (4/4);  total (68/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

CCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACATGT  >  NC_002947/2051119‑2051417
                                                                                                                                                     |                                                                                                                                                     
 cGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGa                                                                                                                                                      >  1:820330/1‑150 (MQ=255)
    ttGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAgg                                                                                                                                                   <  1:4806170/150‑1 (MQ=255)
      gCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGt                                                                                                                                                 >  1:3881306/1‑150 (MQ=255)
      gCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGt                                                                                                                                                 <  2:3029709/150‑1 (MQ=255)
       cTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTc                                                                                                                                                >  2:4430401/1‑150 (MQ=255)
       cTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTc                                                                                                                                                <  2:3108672/150‑1 (MQ=255)
       cTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTc                                                                                                                                                <  1:2469279/150‑1 (MQ=255)
       cTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTc                                                                                                                                                <  1:2468653/150‑1 (MQ=255)
       cTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTc                                                                                                                                                <  1:5014420/150‑1 (MQ=255)
       cTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGACCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTc                                                                                                                                                <  2:52575/150‑1 (MQ=255)
       cTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGACCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTc                                                                                                                                                <  2:52554/150‑1 (MQ=255)
        tCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCg                                                                                                                                               >  2:59172/1‑150 (MQ=255)
        tCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCg                                                                                                                                               >  1:1856438/1‑150 (MQ=255)
        tCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCg                                                                                                                                               >  1:4015986/1‑150 (MQ=255)
        tCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCg                                                                                                                                               >  1:3643501/1‑150 (MQ=255)
        tCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCg                                                                                                                                               >  1:2668686/1‑150 (MQ=255)
         cGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCgg                                                                                                                                              <  1:2119740/150‑1 (MQ=255)
          gATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGc                                                                                                                                             <  1:59172/150‑1 (MQ=255)
            tGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTc                                                                                                                                           <  1:3213693/150‑1 (MQ=255)
            tGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTc                                                                                                                                           <  1:174359/150‑1 (MQ=255)
              cGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGcattt                                                                                                                                         <  2:1485490/150‑5 (MQ=255)
              cGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCgg                                                                                                                                         <  1:3407749/150‑1 (MQ=255)
              cGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCgg                                                                                                                                         <  1:3407745/150‑1 (MQ=255)
               gggCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCACAAGCCGACCCTGAAGGGTCGGCTCggg                                                                                                                                        <  2:4233181/150‑1 (MQ=255)
               gggCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTTCAGGGTCGGCttttg                                                                                                                                        <  2:2933099/150‑5 (MQ=255)
               gggCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCggg                                                                                                                                        >  1:424475/1‑150 (MQ=255)
               gggCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCggg                                                                                                                                        <  2:4088636/150‑1 (MQ=255)
               gggCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCggg                                                                                                                                        <  2:4232902/150‑1 (MQ=255)
                  ccTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGAGCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGcc                                                                                                                                     <  1:2155067/150‑1 (MQ=255)
                    tGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAg                                                                                                                                   >  1:2118494/1‑150 (MQ=255)
                           cggcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCaa                                                                                                                            <  1:2356766/150‑1 (MQ=255)
                                  gcgTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTg                                                                                                                     >  2:1529845/1‑150 (MQ=255)
                                     tttGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATc                                                                                                                  >  1:101313/1‑150 (MQ=255)
                                      ttGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCa                                                                                                                 <  2:820330/150‑1 (MQ=255)
                                          gcTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCt                                                                                                             <  2:3187996/150‑1 (MQ=255)
                                          gcTGGGCCGGGCGGCGGAGCCAGGCACCGAAGACGCGGAGCTGGCCGAGGCCAACCCGCTGGGCTTTGACCAGCACGGGCAACTCCTGAAAGGCATTAGCCGACCCGGATGGGTCGGCTCGGGTCCAGCAACCACTCACGTGTTCAGGCt                                                                                                             <  2:3188059/150‑1 (MQ=255)
                                            ttGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGa                                                                                                           <  2:1623982/150‑1 (MQ=255)
                                             tGGCCGGGCGGCTGAGCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAg                                                                                                          <  2:4014776/150‑1 (MQ=255)
                                                          gATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGTGGAACTCGCTGCgg                                                                                             <  2:1039361/150‑1 (MQ=255)
                                                           aTCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCggg                                                                                            <  2:4941179/150‑1 (MQ=255)
                                                            tCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGt                                                                                           >  2:2325616/1‑150 (MQ=255)
                                                             ccAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTg                                                                                          >  1:4214170/1‑150 (MQ=255)
                                                               aGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGAGCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTggc                                                                                        >  1:1850541/1‑150 (MQ=255)
                                                                  cACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTggcggc                                                                                     <  2:3713260/150‑1 (MQ=255)
                                                                     cGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGtt                                                                                  >  1:2911822/1‑150 (MQ=255)
                                                                     cGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGtt                                                                                  >  1:1747500/1‑150 (MQ=255)
                                                                       aaGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGACCTCGCTGCGGGTGGCGGCGTTGt                                                                                >  2:4653942/1‑150 (MQ=255)
                                                                       aaGACGATGAGCTGGCCGAGGCCAACCCGAAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGACCTCGCTGCGGGTGGCGGCGTTGt                                                                                >  2:36355/1‑150 (MQ=255)
                                                                            gATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCgg                                                                           >  1:1087275/1‑150 (MQ=255)
                                                                                 gCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGACCTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTc                                                                      >  2:4671378/1‑150 (MQ=255)
                                                                                     gCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCAccc                                                                  >  2:448201/1‑150 (MQ=255)
                                                                                           gCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATc                                                            >  1:1153218/1‑150 (MQ=255)
                                                                                             cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAc                                                          >  1:22828/1‑150 (MQ=255)
                                                                                             cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAc                                                          >  1:1206437/1‑150 (MQ=255)
                                                                                                cccGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGa                                                       <  1:1082874/150‑1 (MQ=255)
                                                                                                 ccGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAg                                                      <  1:4109096/150‑1 (MQ=255)
                                                                                                     aGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAACAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGa                                                  <  2:4111715/150‑1 (MQ=255)
                                                                                                      gggCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGAt                                                 >  2:874122/1‑150 (MQ=255)
                                                                                                           ttGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATgg                                            <  2:2582597/150‑1 (MQ=255)
                                                                                                           ttGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATgg                                            <  2:1336684/150‑1 (MQ=255)
                                                                                                             gACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTc                                          >  1:1816688/1‑150 (MQ=255)
                                                                                                                gATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGt                                           >  1:3112381/1‑146 (MQ=255)
                                                                                                                gATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGaa                                       >  1:4578936/1‑150 (MQ=255)
                                                                                                                gATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGaa                                       >  2:779050/1‑150 (MQ=255)
                                                                                                                gATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGCAACTCACCCAGCAACACCGAGGTGATCATGGTCGaa                                       >  2:664840/1‑150 (MQ=255)
                                                                                                                 aTCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAAt                                      <  1:3775697/150‑1 (MQ=255)
                                                                                                                   cATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTc                                    <  2:3302966/150‑1 (MQ=255)
                                                                                                                   cATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTc                                    <  2:2961390/150‑1 (MQ=255)
                                                                                                                     tGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTg                                  <  2:4193524/150‑1 (MQ=255)
                                                                                                                     tGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTg                                  <  2:694110/150‑1 (MQ=255)
                                                                                                                        gCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCtt                               >  2:1814365/1‑150 (MQ=255)
                                                                                                                          aaCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtct                             >  1:597498/1‑150 (MQ=255)
                                                                                                                           aCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                            <  2:1239925/150‑1 (MQ=255)
                                                                                                                           aCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                            <  2:3315788/150‑1 (MQ=255)
                                                                                                                           aCTCATGAAATCCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                            <  1:4700202/150‑1 (MQ=255)
                                                                                                                             tCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGa                          >  2:4368317/1‑150 (MQ=255)
                                                                                                                             tCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGa                          >  2:3206858/1‑150 (MQ=255)
                                                                                                                             tCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGa                          >  1:3637810/1‑150 (MQ=255)
                                                                                                                              ccTGAAATGCAAATGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGac                         >  2:2697164/1‑150 (MQ=255)
                                                                                                                               cTGAAATGCAAAAGCCGACACTGAAGGGTCGGCTCGGGGCCAGAAACCAATCACTTGATAAGGCTGAGCAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTCCTTCTCGaca                        >  2:3469969/1‑150 (MQ=255)
                                                                                                                                tGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGacac                       >  2:2625733/1‑150 (MQ=255)
                                                                                                                                    aTGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACcgcg                   >  1:227135/1‑150 (MQ=255)
                                                                                                                                        ccgaGCCGACCCTTCAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc               >  2:4080470/4‑150 (MQ=255)
                                                                                                                                        ccgaGCCGACCCTTCAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc               >  1:3825657/4‑150 (MQ=255)
                                                                                                                                        aaaaGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc               >  1:3605153/1‑150 (MQ=255)
                                                                                                                                        aaaaGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc               >  1:4120057/1‑150 (MQ=255)
                                                                                                                                        aaaaGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc               >  2:2960071/1‑150 (MQ=255)
                                                                                                                                           aGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatc            <  2:2118494/150‑1 (MQ=255)
                                                                                                                                           aGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatc            <  2:3643501/150‑1 (MQ=255)
                                                                                                                                            gCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatca           >  1:894817/1‑150 (MQ=255)
                                                                                                                                             ccGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatcat          >  1:281029/1‑150 (MQ=255)
                                                                                                                                               gACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGc        >  2:3103637/1‑150 (MQ=255)
                                                                                                                                                aCCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGca       <  1:3208954/150‑1 (MQ=255)
                                                                                                                                                aCCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGca       <  2:3297531/150‑1 (MQ=255)
                                                                                                                                                  ccTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGcaca     >  1:166842/1‑150 (MQ=255)
                                                                                                                                                  ccTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGcaca     >  1:3679919/1‑150 (MQ=255)
                                                                                                                                                   cTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGcacag    <  2:1362386/150‑2 (MQ=255)
                                                                                                                                                    tGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCgg                                                                                             >  2:194493/1‑60 (MQ=255)
                                                                                                                                                    tGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACAtg   >  1:1728030/1‑150 (MQ=255)
                                                                                                                                                     gAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATACTGCTTCTCGACACCGCGCATCATCATGCACAtgt  <  2:1856438/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACATGT  >  NC_002947/2051119‑2051417

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