Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,628,386 A→C 7.5% A16A (GCA→GCC)  clpP → ATP‑dependent Clp protease proteolytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,628,3860AC 7.5% 220.7 / 7.2 76A16A (GCA→GCC) clpPATP‑dependent Clp protease proteolytic subunit
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base C (2/4);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

GTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAACTGCT  >  NC_002947/2628244‑2628532
                                                                                                                                              |                                                                                                                                                    
gTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAggcggc                                                                                                                                                <  1:143248/149‑1 (MQ=255)
 tAACGGTTTTCTCTCGTAGGAAACCTCCACATGCCAGCCTTCGAGCTGTCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAggcggc                                                                                                                                                <  1:131581/148‑1 (MQ=255)
   aCGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCt                                                                                                                                              <  8:170320/148‑1 (MQ=255)
   aCGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCt                                                                                                                                              <  1:394877/148‑1 (MQ=255)
      gTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTggat                                                                                                                                          <  7:149777/149‑3 (MQ=255)
      gTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTggat                                                                                                                                          <  7:119327/149‑3 (MQ=255)
       tcttctCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCCTGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTcc                                                                                                                                         <  4:106801/149‑1 (MQ=255)
       tcttctCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTcc                                                                                                                                         <  2:108644/149‑1 (MQ=255)
              tcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGt                                                                                                                                  >  5:168640/1‑149 (MQ=255)
              tcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGt                                                                                                                                  >  5:19615/1‑149 (MQ=255)
              tcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGt                                                                                                                                  >  3:138008/1‑149 (MQ=255)
               cGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGt                                                                                                                                  >  7:127391/1‑148 (MQ=255)
                 tAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTAt                                                                                                                               >  6:358209/1‑149 (MQ=255)
                  aGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTccc                                                                                                                                        >  5:358891/1‑139 (MQ=255)
                  aGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTccc                                                                                                                                        <  6:358891/139‑1 (MQ=255)
                   ggAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATc                                                                                                                              <  5:355900/148‑1 (MQ=255)
                   ggAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATc                                                                                                                              >  4:284313/1‑148 (MQ=255)
                    gAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGa                                                                                                                            <  5:358209/149‑1 (MQ=255)
                     aaaCCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAg                                                                                                                           >  2:99260/1‑149 (MQ=255)
                       aCCTACAGAAGCAATCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTTGGAGTGAGCTCATGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCa                                                                                                                         <  7:357460/149‑1 (MQ=255)
                         cTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGt                                                                                                                       >  3:205597/1‑149 (MQ=255)
                                agccagccTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTACCGAGGGTTATCGAGCAGTCCG‑CCCg                                                                                                               >  6:363637/1‑149 (MQ=255)
                                 gccagccTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGt                                                                                                              >  6:312202/1‑149 (MQ=255)
                                   cagccTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTgg                                                                                                            <  8:35099/149‑1 (MQ=255)
                                       cTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGaa                                                                                                        >  3:338552/1‑149 (MQ=255)
                                        ttCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAAc                                                                                                       >  2:325094/1‑149 (MQ=255)
                                                   ttGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGAAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGa                                                                                             >  2:351947/1‑148 (MQ=255)
                                                     gtgCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCAGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACa                                                                                           >  1:73415/1‑148 (MQ=255)
                                                               aaGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGa                                                                                                                                      <  6:8435/96‑1 (MQ=255)
                                                               aaGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGa                                                                                                                                      >  5:8435/1‑96 (MQ=255)
                                                               aaGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCgcgc                                                                                >  3:50349/1‑149 (MQ=255)
                                                               aaGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCgcgc                                                                                >  7:50451/1‑149 (MQ=255)
                                                               aaGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCgcgc                                                                                >  6:273254/1‑149 (MQ=255)
                                                                    aTGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGtt                                                                           >  6:282745/1‑149 (MQ=255)
                                                                     tGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTg                                                                          <  3:284313/149‑1 (MQ=255)
                                                                       aCTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTg                                                                          >  1:84210/1‑147 (MQ=255)
                                                                       aCTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGaa                                                                        <  4:102726/149‑1 (MQ=255)
                                                                        cTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAg                                                                       <  5:303506/149‑1 (MQ=255)
                                                                           tttGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCGGCGGCCGGGTCCCGATGGGTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGCCATCTACTCGCGGCTGTTGAAggcg                                                                    >  8:127317/1‑147 (MQ=255)
                                                                           tttGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑Ag                                                                   <  1:95240/149‑1 (MQ=255)
                                                                               ggAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCgtg                                                               <  5:282745/149‑1 (MQ=255)
                                                                               ggAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCgtg                                                               <  8:127391/149‑1 (MQ=255)
                                                                                 aGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATc                                                                                         <  7:362017/122‑1 (MQ=255)
                                                                                 aGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATc                                                                                         <  7:367224/122‑1 (MQ=255)
                                                                                 aGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATc                                                                                         >  8:362017/1‑122 (MQ=255)
                                                                                 aGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATc                                                                                         >  8:367224/1‑122 (MQ=255)
                                                                                       cgcAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTc                                                       >  2:229491/1‑149 (MQ=255)
                                                                                       cgcAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTc                                                       <  6:19615/149‑1 (MQ=255)
                                                                                           ggACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTgg                                                   <  8:50451/149‑1 (MQ=255)
                                                                                                  tCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCg                                                                                                          <  2:265405/88‑1 (MQ=255)
                                                                                                  tCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCg                                                                                                          >  1:265405/1‑88 (MQ=255)
                                                                                                    ccGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCgg                                          >  8:111466/1‑149 (MQ=255)
                                                                                                    ccGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCgg                                          <  1:355066/149‑1 (MQ=255)
                                                                                                     cGCAATTCTTATATTCAGCAGAGCTCTGTCATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTTTCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGCAGCGTGTGATCTTCCTGGTTGGCCCg                                           <  1:205789/148‑1 (MQ=255)
                                                                                                     cGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCgg                                          >  2:338669/1‑148 (MQ=255)
                                                                                                     cGCAAGTCGTATATTCAGCAGAGCTCTGACATCCATGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCAG‑CCCGTGGAGAACGCGCTTACGACATCTACACGCGCCTCTAGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCg                                           >  2:205789/1‑147 (MQ=255)
                                                                                                      gCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCgg                                          >  2:312852/1‑147 (MQ=255)
                                                                                                        aaTTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTaga                                      <  5:273254/149‑1 (MQ=255)
                                                                                                        aaTTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATAGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTa                                        >  7:128479/1‑147 (MQ=255)
                                                                                                         aTTCTTATATTGAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACTCTATTACGACATCTACTCGCGCCTGTTGAAGT‑AGCGTGTGATCTTCCTGGTTGGCCCGGTagag                                     <  1:71345/149‑1 (MQ=255)
                                                                                                         aTTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTagag                                     <  4:138008/149‑1 (MQ=255)
                                                                                                         aTTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCACGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTagag                                     <  4:205597/149‑1 (MQ=255)
                                                                                                               atatTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTa                                <  1:338669/148‑1 (MQ=255)
                                                                                                                  ttCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACAt                             <  2:84210/148‑1 (MQ=255)
                                                                                                                     agcagAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGc                          >  8:220586/1‑148 (MQ=255)
                                                                                                                      gcagAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGcc                         >  6:119823/1‑148 (MQ=255)
                                                                                                                              ctGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCtcggc                >  5:341508/1‑147 (MQ=255)
                                                                                                                                 aCATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATgg                           >  2:232515/1‑135 (MQ=255)
                                                                                                                                 aCATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATgg                           <  1:232515/135‑1 (MQ=255)
                                                                                                                                 aCATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGtt             >  3:102896/1‑149 (MQ=255)
                                                                                                                                 aCATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGtt             <  5:119823/149‑1 (MQ=255)
                                                                                                                                   aTCCACGCCGCAGGCGGGCTGGTCCCGATTGTTATCGATCATTCCGCCCCGCAATGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTg            <  1:325094/149‑1 (MQ=255)
                                                                                                                                     ccAGGCCGCCTGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTgcg          >  6:1053/1‑148 (MQ=255)
                                                                                                                                      cAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTgcgc         <  4:102896/148‑1 (MQ=255)
                                                                                                                                      cAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTgcgc         <  6:168640/148‑1 (MQ=255)
                                                                                                                                            gCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAActgct  >  3:317569/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCG‑CCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGG‑AGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAACTGCT  >  NC_002947/2628244‑2628532

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑