Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,223,153 G→T 6.4% intergenic (+46/+216) PP_3699 → / ← PP_3700 hypothetical protein/chromosome partitioning ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,223,1530GT 6.4% 542.6 / 17.3 164intergenic (+46/+216)PP_3699/PP_3700hypothetical protein/chromosome partitioning ATPase
Reads supporting (aligned to +/- strand):  ref base G (96/57);  new base T (3/8);  total (99/65)
Fisher's exact test for biased strand distribution p-value = 2.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

CACTTGCGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCCGCGGC  >  NC_002947/4223008‑4223302
                                                                                                                                                 |                                                                                                                                                     
  cTTGCGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGa                                                                                                                                                 <  1:4036506/150‑1 (MQ=255)
   ttGCGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGAt                                                                                                                                                >  2:5349836/1‑150 (MQ=255)
     gCGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGc                                                                                                                                              <  1:406088/150‑1 (MQ=255)
     gCGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGc                                                                                                                                              <  2:2847467/150‑1 (MQ=255)
      cGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGcc                                                                                                                                             <  1:5445316/150‑1 (MQ=255)
       gTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccc                                                                                                                                            <  2:4094755/150‑1 (MQ=255)
       gTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGAGATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGccc                                                                                                                                            >  1:2316528/1‑150 (MQ=255)
        tGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGcccc                                                                                                                                           <  1:2762955/150‑1 (MQ=255)
            aCGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccctga                                                                                                                                       <  2:2813407/150‑5 (MQ=255)
             cGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcct                                                                                                                                      >  1:812937/1‑150 (MQ=255)
               tCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctca                                                                                                                                    >  1:1805575/1‑150 (MQ=255)
                 gAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagc                                                                                                                                  <  1:3059464/150‑1 (MQ=255)
                   gTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCGAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccct                                                                                                                                         <  2:3735609/141‑3 (MQ=255)
                   gTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagcct                                                                                                                                >  2:494867/1‑150 (MQ=255)
                    ttGACGCTGTTCACGCTGAAGATCCTGAAGATTTGATCAGGGATGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagcctc                                                                                                                               >  2:3494464/1‑150 (MQ=255)
                         gCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagcctcagcct                                                                                                                          >  2:3824131/1‑150 (MQ=255)
                           tGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagcctcagcctca                                                                                                                        >  2:3810348/1‑150 (MQ=255)
                             ttCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc                                                                                                                      >  1:1722636/1‑150 (MQ=255)
                             ttCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc                                                                                                                      <  1:5349752/150‑1 (MQ=255)
                             ttCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc                                                                                                                      <  1:3722678/150‑1 (MQ=255)
                               cACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTcagcctcagcctcagcct                                                                                                                    >  1:3882085/1‑150 (MQ=255)
                                   cTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCaaa                                                                                                                <  1:4962860/150‑1 (MQ=255)
                                   cTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCaaa                                                                                                                <  1:4962891/150‑1 (MQ=255)
                                    tGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCaaaa                                                                                                               >  2:1277320/1‑150 (MQ=255)
                                      aaGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAAcc                                                                                                             >  1:2181252/1‑150 (MQ=255)
                                      aaGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAAcc                                                                                                             >  1:2072223/1‑150 (MQ=255)
                                        gATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcagc                                                                                                                                  >  1:2369631/1‑127 (MQ=255)
                                        gATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTa                                                                                                           <  2:1017046/150‑1 (MQ=255)
                                        gATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTa                                                                                                           <  1:5372611/150‑1 (MQ=255)
                                        gATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTa                                                                                                           <  1:4709478/150‑1 (MQ=255)
                                         aTCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAg                                                                                                          <  2:1397473/150‑1 (MQ=255)
                                             tGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTAt                                                                                                      >  1:3688304/1‑150 (MQ=255)
                                             tGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTAt                                                                                                      >  2:908698/1‑150 (MQ=255)
                                             tGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTAt                                                                                                      >  2:3132098/1‑150 (MQ=255)
                                              gAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCGGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATg                                                                                                     >  1:1922910/1‑150 (MQ=255)
                                               aaGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATgg                                                                                                    <  1:4562458/150‑1 (MQ=255)
                                               aaGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATgg                                                                                                    <  2:1722636/150‑1 (MQ=255)
                                                aGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAGGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGc                                                                                                   <  1:3710714/150‑1 (MQ=255)
                                                  cTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCg                                                                                                 >  2:2603732/1‑150 (MQ=255)
                                                   tttGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGa                                                                                                >  2:4311625/1‑150 (MQ=255)
                                                             gCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGaa                                                                                      <  2:706516/150‑1 (MQ=255)
                                                                  tatGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCcac                                                                                 >  1:1930656/1‑150 (MQ=255)
                                                                  tatGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCcac                                                                                 >  1:1689094/1‑150 (MQ=255)
                                                                      cTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCACAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAAtc                                                                             >  2:3358357/1‑150 (MQ=255)
                                                                       ttGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAAtct                                                                            >  2:3280708/1‑150 (MQ=255)
                                                                       ttGAGGAGTTTGAACAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAAtctc                                                                           >  2:2153664/1‑150 (MQ=255)
                                                                         gaggagTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCa                                                                          <  1:5349836/150‑1 (MQ=255)
                                                                         gaggagTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCa                                                                          >  2:3391575/1‑150 (MQ=255)
                                                                               tttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAAc                                                                    >  2:367975/1‑150 (MQ=255)
                                                                                ttGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAAcc                                                                   <  1:1101482/150‑1 (MQ=255)
                                                                                     cAAGACTGCATGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGc                                                              <  2:4051122/150‑1 (MQ=255)
                                                                                      aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGcc                                                             >  1:5143858/1‑150 (MQ=255)
                                                                                          cTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTAc                                                         <  2:2477816/150‑1 (MQ=255)
                                                                                           tGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCGGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTAc                                                         >  1:4156120/1‑149 (MQ=255)
                                                                                              aagtgaagtgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGt                                                     >  2:938140/1‑150 (MQ=255)
                                                                                               agtgaagtgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTc                                                    >  2:1234965/1‑150 (MQ=255)
                                                                                                gtgaagtgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCCCATGCAACCGCTGCCGTACGAGTCg                                                   <  2:4534956/150‑1 (MQ=255)
                                                                                                   aagtgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGca                                                >  2:2397400/1‑150 (MQ=255)
                                                                                                   aagtgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGca                                                >  2:2397457/1‑150 (MQ=255)
                                                                                                      tgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctc                                                                                                                                     >  1:2675710/1‑62 (MQ=255)
                                                                                                      tgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctc                                                                                                                                     <  2:2675710/62‑1 (MQ=255)
                                                                                                      tgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAg                                             >  1:4281016/1‑150 (MQ=255)
                                                                                                        gCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAATCCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTa                                           >  1:3753540/1‑150 (MQ=255)
                                                                                                          aaaTGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTAtt                                         >  1:2564361/1‑150 (MQ=255)
                                                                                                          aaaTGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTAtt                                         >  1:1671016/1‑150 (MQ=255)
                                                                                                          aaaTGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTAtt                                         >  1:4859075/1‑150 (MQ=255)
                                                                                                          aaaTGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTAtt                                         >  1:1115190/1‑150 (MQ=255)
                                                                                                            aTGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGc                                       <  1:4079678/150‑1 (MQ=255)
                                                                                                             tGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCa                                      <  1:2720121/150‑1 (MQ=255)
                                                                                                             tGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCa                                      >  2:3591434/1‑150 (MQ=255)
                                                                                                             tGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCa                                      >  2:3494502/1‑150 (MQ=255)
                                                                                                             tGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCa                                      <  1:1070302/150‑1 (MQ=255)
                                                                                                              gCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTa                                                                                                       >  1:4424979/1‑84 (MQ=255)
                                                                                                              gCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCAt                                     >  2:2037604/1‑150 (MQ=255)
                                                                                                               cTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCAtc                                    >  1:5155949/1‑150 (MQ=255)
                                                                                                                 gACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCAtctc                                  <  2:3793401/150‑1 (MQ=255)
                                                                                                                      tcgtGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCt                             >  2:854652/1‑150 (MQ=255)
                                                                                                                        gtGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCCGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTcc                           <  2:1646597/150‑1 (MQ=255)
                                                                                                                          gCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCCGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGa                         <  2:1986192/150‑1 (MQ=255)
                                                                                                                          gCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGa                         >  2:3052133/1‑150 (MQ=255)
                                                                                                                          gCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGa                         <  1:5014490/150‑1 (MQ=255)
                                                                                                                          gCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGa                         <  2:1805575/150‑1 (MQ=255)
                                                                                                                          gCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGa                         <  1:2904904/150‑1 (MQ=255)
                                                                                                                          gCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGa                         >  2:3973382/1‑150 (MQ=255)
                                                                                                                                   gcTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCt                <  2:2316528/150‑1 (MQ=255)
                                                                                                                                     tCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGa              >  1:1371345/1‑150 (MQ=255)
                                                                                                                                         ggggCATCTGTCCGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACAtt          >  1:963441/1‑150 (MQ=255)
                                                                                                                                         ggggCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACAtt          >  1:5220965/1‑150 (MQ=255)
                                                                                                                                          gggCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACACTCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTCTTTCATCTCTAGCTCCGAGCATGCCCTGACAttt         >  1:1013327/1‑150 (MQ=255)
                                                                                                                                          gggCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACAttt         <  1:1178369/150‑1 (MQ=255)
                                                                                                                                            gCATCTGTCCGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAAtctc                                                                           <  1:3413630/82‑1 (MQ=255)
                                                                                                                                            gCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTcc       <  2:3089209/150‑1 (MQ=255)
                                                                                                                                              aTCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCcgc     >  1:2902859/1‑150 (MQ=255)
                                                                                                                                              aTCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCcgc     >  2:4508400/1‑150 (MQ=255)
                                                                                                                                              aTCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCcgc     >  1:403316/1‑150 (MQ=255)
                                                                                                                                              aTCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCAGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCcgc     >  1:3494893/1‑150 (MQ=255)
                                                                                                                                               tCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGTCATTTCcgcg    >  1:4716991/1‑150 (MQ=255)
                                                                                                                                               tCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCcgcg    >  2:3678959/1‑150 (MQ=255)
                                                                                                                                               tCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCcgcg    >  2:262093/1‑150 (MQ=255)
                                                                                                                                                 ggACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCCGCGGc  >  2:4899316/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CACTTGCGTGCAACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACATTTCCGCGGC  >  NC_002947/4223008‑4223302

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: