Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,633,379 G→T 5.0% G17C (GGC→TGC)  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,633,3790GT 5.0% 875.1 / 24.4 258G17C (GGC→TGC) clpPATP‑dependent Clp protease proteolytic subunit
Reads supporting (aligned to +/- strand):  ref base G (115/129);  new base T (3/10);  total (119/139)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 278 total aligned reads displayed.

TGCTGAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAACTGCTGTTC  >  NC_002947/2633230‑2633528
                                                                                                                                                     |                                                                                                                                                     
tgctgAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAg                                                                                                                                                       <  2:3619944/150‑1 (MQ=255)
   tgAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAggcg                                                                                                                                                    <  2:1050588/150‑1 (MQ=255)
   tgAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAggcg                                                                                                                                                    >  2:758620/1‑150 (MQ=255)
   tgAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAggcg                                                                                                                                                    >  1:4332715/1‑150 (MQ=255)
        aaCGGTCGTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTg                                                                                                                                               <  2:5168066/150‑1 (MQ=255)
          cGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTgga                                                                                                                                             <  2:2238647/150‑2 (MQ=255)
           ggTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTggat                                                                                                                                            <  1:4356606/150‑3 (MQ=255)
           ggTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTggat                                                                                                                                            <  2:379844/150‑3 (MQ=255)
             tcttctCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTggatgt                                                                                                                                          <  1:527482/150‑5 (MQ=255)
             tcttctCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTccc                                                                                                                                          >  2:1640660/1‑150 (MQ=255)
                 ctctcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATg                                                                                                                                      <  2:1523698/150‑1 (MQ=255)
                       tAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATc                                                                                                                                <  2:3750876/150‑1 (MQ=255)
                        aGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCg                                                                                                                               >  2:2706702/1‑150 (MQ=255)
                          gAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAg                                                                                                                             <  2:2104466/150‑1 (MQ=255)
                            aaCCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCa                                                                                                                           <  2:3608800/150‑1 (MQ=255)
                            aaCCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCa                                                                                                                           <  2:4727350/150‑1 (MQ=255)
                              ccTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGt                                                                                                                         <  1:3277758/150‑1 (MQ=255)
                               cTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTc                                                                                                                        >  2:5514102/1‑150 (MQ=255)
                               cTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCAGGTCCCGATGGTTATCGAGCAGTc                                                                                                                        >  1:1457537/1‑150 (MQ=255)
                                   acaAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGcc                                                                                                                    <  2:1560270/150‑1 (MQ=255)
                                     aaGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCg                                                                                                                  >  2:3198924/1‑150 (MQ=255)
                                       gccagccTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTg                                                                                                                <  2:2268405/150‑1 (MQ=255)
                                        ccagccTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTgg                                                                                                               <  1:1435439/150‑1 (MQ=255)
                                         cagccTTCGATCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGc                                                                                                              <  2:1314510/150‑1 (MQ=255)
                                               tCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAAcgc                                                                                                        >  2:5356529/1‑150 (MQ=255)
                                               tCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAAcgc                                                                                                        >  2:647285/1‑150 (MQ=255)
                                                  aGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCt                                                                                                     >  1:3385791/1‑150 (MQ=255)
                                                    cTGGCTTGTGCGTATTCACGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTa                                                                                                   <  1:3850728/150‑1 (MQ=255)
                                                    cTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTa                                                                                                   <  1:3074550/150‑1 (MQ=255)
                                                     tGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTAc                                                                                                  >  1:2344782/1‑150 (MQ=255)
                                                     tGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTAc                                                                                                  <  2:2206479/150‑1 (MQ=255)
                                                      ggCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACg                                                                                                 >  2:609761/1‑150 (MQ=255)
                                                            tgCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCt                                                                                           >  2:2257052/1‑150 (MQ=255)
                                                             gCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTa                                                                                          >  2:507108/1‑150 (MQ=255)
                                                              cGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTAc                                                                                         <  1:2923828/150‑1 (MQ=255)
                                                               gTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACt                                                                                        >  2:3938298/1‑150 (MQ=255)
                                                                tATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTc                                                                                       <  1:3694450/150‑1 (MQ=255)
                                                                tATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTc                                                                                       <  1:684445/150‑1 (MQ=255)
                                                                     aaGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGcc                                                                                  >  2:4144261/1‑150 (MQ=255)
                                                                        ccATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGt                                                                               >  2:3422031/1‑150 (MQ=255)
                                                                         cATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGtt                                                                              <  2:399840/150‑1 (MQ=255)
                                                                            gACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCATAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGccc                                                                                                                   <  1:5355868/110‑1 (MQ=255)
                                                                            gACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGaa                                                                           >  2:2466484/1‑150 (MQ=255)
                                                                            gACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGaa                                                                           <  2:5412326/150‑1 (MQ=255)
                                                                            gACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGAGGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGaa                                                                           <  2:83686/150‑1 (MQ=255)
                                                                            gACTATTTGCGAGTCAGCCCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCACGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGaa                                                                           <  2:4472108/150‑1 (MQ=255)
                                                                             aCTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAg                                                                          <  1:1919865/150‑1 (MQ=255)
                                                                               tATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGa                                                                        >  1:4183600/1‑150 (MQ=255)
                                                                               tATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGa                                                                        >  2:311090/1‑150 (MQ=255)
                                                                               tATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGa                                                                        >  2:5624197/1‑150 (MQ=255)
                                                                               tATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGa                                                                        >  1:2642955/1‑150 (MQ=255)
                                                                               tATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGa                                                                        >  2:3966012/1‑150 (MQ=255)
                                                                               tATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGa                                                                        >  2:479563/1‑150 (MQ=255)
                                                                                 tttGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGc                                                                      <  1:495953/150‑1 (MQ=255)
                                                                                  ttGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGTAGGAGCg                                                                     >  1:3504843/1‑150 (MQ=255)
                                                                                   tGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACACGCGCCTGTTGAAGGAGCgt                                                                    <  2:4843705/150‑1 (MQ=255)
                                                                                   tGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCAGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCgt                                                                    >  2:832959/1‑150 (MQ=255)
                                                                                      gAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGtgtg                                                                 <  1:4504620/150‑1 (MQ=255)
                                                                                      gAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGtgtg                                                                 <  1:2135162/150‑1 (MQ=255)
                                                                                       aGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGa                                                                <  1:4606455/150‑1 (MQ=255)
                                                                                         tGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATc                                                              >  1:1010872/1‑150 (MQ=255)
                                                                                         tGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATc                                                              <  1:758620/150‑1 (MQ=255)
                                                                                         tGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATc                                                              >  1:1616751/1‑150 (MQ=255)
                                                                                            gcgcAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTc                                                           >  1:4002912/1‑150 (MQ=255)
                                                                                             cgcAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTcc                                                          <  1:391840/150‑1 (MQ=255)
                                                                                               cAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTg                                                        <  1:3524850/150‑1 (MQ=255)
                                                                                                aGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGAt                                                                                                                                       >  1:486782/1‑70 (MQ=255)
                                                                                                aGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGAt                                                                                                                                       <  2:486782/70‑1 (MQ=255)
                                                                                                aGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTgg                                                       >  1:64548/1‑150 (MQ=255)
                                                                                                aGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTgg                                                       >  2:4407628/1‑150 (MQ=255)
                                                                                                  gACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGtt                                                     <  2:2303881/150‑1 (MQ=255)
                                                                                                  gACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGtt                                                     <  1:1702167/150‑1 (MQ=255)
                                                                                                     aTGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGc                                                  <  2:1852912/150‑1 (MQ=255)
                                                                                                      tGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGaa                                                                           <  2:1152205/124‑1 (MQ=255)
                                                                                                      tGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGaa                                                                           >  1:1152158/1‑124 (MQ=255)
                                                                                                        tCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACgcgc                                                                                                      <  2:4304096/95‑1 (MQ=255)
                                                                                                        tCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGccat                                               <  2:2114767/150‑3 (MQ=255)
                                                                                                         cccGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCgg                                              <  1:692520/150‑1 (MQ=255)
                                                                                                           cGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTa                                            >  1:5564392/1‑150 (MQ=255)
                                                                                                            gCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTag                                           >  1:2719315/1‑150 (MQ=255)
                                                                                                            gCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTag                                           >  1:197879/1‑150 (MQ=255)
                                                                                                            gCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTag                                           >  1:1820139/1‑150 (MQ=255)
                                                                                                             cAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTaga                                          <  1:2220183/150‑1 (MQ=255)
                                                                                                              aaTTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTagag                                         <  1:2444626/150‑1 (MQ=255)
                                                                                                              aaTTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTagag                                         <  2:1061239/150‑1 (MQ=255)
                                                                                                                   ttATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTa                                    >  1:939606/1‑150 (MQ=255)
                                                                                                                   ttATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTa                                    >  2:1978692/1‑150 (MQ=255)
                                                                                                                   ttATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTa                                    >  2:2906973/1‑150 (MQ=255)
                                                                                                                       atTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATg                                <  2:1211545/150‑1 (MQ=255)
                                                                                                                                   tctGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCtcgtc                    >  1:523622/1‑150 (MQ=255)
                                                                                                                                   tctGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCATGTAGTCCTCTACCGGGCCAAccaggaa                    <  2:145257/150‑6 (MQ=38)
                                                                                                                                     tGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGtt                                                     >  2:292581/1‑115 (MQ=255)
                                                                                                                                     tGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGtt                                                     <  1:292581/115‑1 (MQ=255)
                                                                                                                                     tGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTcgtcgt                  >  1:1013728/1‑150 (MQ=255)
                                                                                                                                     tGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTcgtcgt                  >  1:3139609/1‑150 (MQ=255)
                                                                                                                                        cATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTgc               <  1:1057931/150‑1 (MQ=255)
                                                                                                                                            cAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCaa           >  1:2469709/1‑150 (MQ=255)
                                                                                                                                            cAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCaa           >  2:4955955/1‑150 (MQ=255)
                                                                                                                                            cAGGCCGCACTCGCCCTGTTCCCCATGTTTAAAGACCAGTCCGCCCGTGGCGAAAGACCTAACGACATCTACTCGCGCCTCTTTAACGAGCCTGTGATCTTCATGGTTGGCCAGTTAGAGGACAACATGGCCAACCTCGTCGTTGCGCaa           >  2:1037491/1‑150 (MQ=255)
                                                                                                                                             aGGCCGCAGGAGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGtt                 >  1:510053/1‑143 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCTGAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAACTGCTGTTC  >  NC_002947/2633230‑2633528

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: