Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,377,809:1 +A 7.7% intergenic (+267/+25) PP_2970 → / ← PP_2971 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,377,8091.A 7.7% 1102.5 / 87.1 272intergenic (+267/+25)PP_2970/PP_2971hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base . (145/106);  new base A (6/15);  total (151/121)
Fisher's exact test for biased strand distribution p-value = 1.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 294 total aligned reads displayed.

TGTTGTTTATGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACTTGCCGCCCA  >  NC_002947/3377660‑3377958
                                                                                                                                                      |                                                                                                                                                     
tgttgtTTATGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAg                                                                                                                                                        >  2:6860012/1‑150 (MQ=255)
 gttgttTATGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAgc                                                                                                                                                      <  1:3696492/150‑1 (MQ=255)
    gttTATGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAgcggc                                                                                                                                                   <  2:541337/150‑1 (MQ=255)
      ttATGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGa                                                                                                                                                 <  2:3453181/150‑1 (MQ=255)
         tGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCGCTGCAGAGGGCGATTc                                                                                                                                              <  2:1335114/150‑1 (MQ=255)
          gggCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCg                                                                                                                                             >  2:5770926/1‑150 (MQ=255)
           ggCACCTCGAATAACGCGAGGTTTTCAGTGCCTCCGACCGCCAAATGAAGGCCATGATCGACTCCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCgg                                                                                                                                            <  2:3582181/150‑1 (MQ=255)
           ggCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCgg                                                                                                                                            <  1:6623299/150‑1 (MQ=255)
            gCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCTGAATCGCCGCTGCAGAGGGCGATTCgga                                                                                                                                           <  1:1307681/150‑2 (MQ=255)
            gCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGc                                                                                                                                           <  1:1716877/150‑1 (MQ=255)
            gCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCGCTGCAGAGGGCGATTCgga                                                                                                                                           <  1:2702458/150‑2 (MQ=255)
             cACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCGCTGCAGAGGGCGATTCggaa                                                                                                                                          <  1:1740328/150‑3 (MQ=255)
             cACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCGCTGCAGAGGGCGATTCggaa                                                                                                                                          <  2:6819408/150‑3 (MQ=255)
             cACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCGCTGCAGAGGGCGATTCggaa                                                                                                                                          <  1:3071802/150‑3 (MQ=255)
             cACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTGCCGAATCGCCGCTGCAGAGGGCGATTCggaa                                                                                                                                          <  1:5225867/150‑3 (MQ=255)
                cTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCAc                                                                                                                                       >  2:2490418/1‑150 (MQ=255)
                cTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCAc                                                                                                                                       >  1:2492351/1‑150 (MQ=255)
                cTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCAc                                                                                                                                       >  2:5096920/1‑150 (MQ=255)
                cTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCAc                                                                                                                                       >  2:6003745/1‑150 (MQ=255)
                cTCGAATAACGCGAGGTTGTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGAGTCGGCCCAc                                                                                                                                       >  1:5863542/1‑150 (MQ=255)
                 tCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACt                                                                                                                                      >  1:2796822/1‑150 (MQ=255)
                 tCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACt                                                                                                                                      >  1:5095220/1‑150 (MQ=255)
                 tCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCACCCCGTTTTTCTTCCGAATCGCCGCTGCAGAGGGCGATTCggaa                                                                                                                                          >  1:1874083/1‑144 (MQ=255)
                 tCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCACCCCGTTTTTCTTCCGAATCGCCGCTGCAGAGGGCGATTCggaa                                                                                                                                          <  2:1874083/146‑3 (MQ=255)
                   gAATAACGCGTGGTTTTCAGTTAATCCGACCGCAAAATATTGGCCTTGATCGAATTGGGTTCTTTTTCTGGGGCTAATCACCGCTTTCTGGGGGGTGTCGCCCCTTTTTTCTTCCGAATCGCCCTCTGCTGCGGCGATTCGGCCCACTcc                                                                                                                                    <  2:3847319/150‑1 (MQ=255)
                            gAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAg                                                                                                                           >  1:2463358/1‑150 (MQ=255)
                            gAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAg                                                                                                                           >  2:2392980/1‑150 (MQ=255)
                               gTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGTTTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAg                                                                                                                        <  2:1118105/150‑1 (MQ=255)
                                          aTCCGCCCGCCAATTGATGGCCATGACCGAATTCGGTTCGTTTTCTGCGGCGAACCCCCGCTTTCTGGGCGGTGTCGCCCCTTTTTCCTCCCGAATCGCCCTGTGCAGGGGCTTTTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTcc                                                                                                             <  2:2309713/150‑1 (MQ=255)
                                           tCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTccc                                                                                                            >  1:5675997/1‑150 (MQ=255)
                                            ccGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCg                                                                                                           >  2:4049848/1‑150 (MQ=255)
                                             cGACCGCCAATTGTGGGCCTTGATCGATTGGGGTTCGTTTTGTGGGGCGAATCACGGGTTTGTGGGGGTTGTCGCCCGTTTTTTCTCCCGACTCGCCCTCGGCAGGGGGTTTTGGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGc                                                                                                          <  1:6181231/150‑1 (MQ=255)
                                             cGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGc                                                                                                          >  2:899316/1‑150 (MQ=255)
                                              gACCCCCAAATGTTGCCCTTTATCTATTTCTGGTCGTTTTCTGCGGCTAATCACCGCTTTCTGGGGGGTGTCGCCCCGTTTTTCTTCCCAATCGCCCCCCGCACCGGCCATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGcc                                                                                                         <  1:6689810/150‑1 (MQ=255)
                                                ccGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCTAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCtt                                                                                                       <  2:2813007/150‑1 (MQ=255)
                                                  ccaaaaTGTGGGCCATGAGCGATTGCGGTTCGATTTCGGGGGCGAATCACCGCTTTCTGGGCGGTGCCGCCCCTTTTTTCTTCCGAATCGCCCTCTGCAGGGGCTATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCt                                                                                                     <  2:6537059/147‑1 (MQ=255)
                                                   ccAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTg                                                                                                    >  1:6373645/1‑150 (MQ=255)
                                                     aaaTGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGa                                                                                                  >  2:4742892/1‑150 (MQ=255)
                                                         gatgGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGGGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTa                                                                                              <  2:1575850/150‑1 (MQ=255)
                                                            ggCCATGATCGAATTCGGTTCGTTTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCTTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTacac                                                                                           <  2:2236589/150‑1 (MQ=255)
                                                            ggCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTacac                                                                                           >  1:748345/1‑150 (MQ=255)
                                                             gCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACAcc                                                                                          >  1:3529002/1‑150 (MQ=255)
                                                             gCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACAcc                                                                                          >  2:6655447/1‑150 (MQ=255)
                                                              ccATTATCGAATTCGGTTCGTTTTCTGCGGCGAATCACGTCTTTGTGGGCGGTGTCGCCGCTTTTTTCTTCCTACTCGCCCTCGGCAGGGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCa                                                                                         <  2:5555413/150‑1 (MQ=255)
                                                              ccATGAGCGTAGCCGGTTCGATTGCGGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCCGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCa                                                                                         <  2:1540780/150‑1 (MQ=255)
                                                                aTGATCTAATTGGGTTCGTTTTCTGCGGCTAATCACCGGTTTGGGGGGGGTGTCGCCCCGTTTTTGTTCCGAATCGCCCTCTGCAGCGGCTATTGGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGc                                                                                       <  2:4188022/150‑1 (MQ=255)
                                                                   ctCGATTTCGGTTCGTTTTGTGCGGCGAACCCCCGCTTGCGGGGGGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCGGCTGGGGCGATTGGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGc                                                                                    <  2:3645249/149‑1 (MQ=255)
                                                                      gAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTc                                                                                 >  1:4298880/1‑150 (MQ=255)
                                                                        tttgggTTCTTTTTCTGGGGCTAAACCCCGCTTTCTGGGGGTTGGCGCCCCTTTTTTCTCCCCAATCGCCCTCTGCCGCGGCGTTTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCaa                                                                               <  1:983306/146‑1 (MQ=255)
                                                                            ggTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTg                                                                           >  2:3373259/1‑150 (MQ=255)
                                                                            ggTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTg                                                                           >  2:3199580/1‑150 (MQ=255)
                                                                             ggtctttttGTGCGGCGACCCCCCGCTGTGGGGGGGGTGGCGCCCGTTTTTCCTCCCGAACCGCCCTCGGCCGGGGCGTTCCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGt                                                                          <  1:3113279/144‑1 (MQ=255)
                                                                              ttCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTa                                                                         >  1:4671937/1‑150 (MQ=255)
                                                                              ttCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTa                                                                         >  2:2052227/1‑150 (MQ=255)
                                                                                 gATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCa                                                                      >  2:5654404/1‑150 (MQ=255)
                                                                                  aTTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAg                                                                     >  1:5949224/1‑150 (MQ=255)
                                                                                     tgggcggctAATCCCGGTTTTGTGGGGGGTGGCGCCCGTTTTTTCTCCCTAATCGCCCTCGGCCGGGGCTATCCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGACGCTTCAAGTTGTAGCAGgcc                                                                  <  1:6542182/147‑1 (MQ=255)
                                                                                           gcgACTCACCGCTTGCGGGGCGGGGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCcatc                                                            <  2:2669726/150‑1 (MQ=255)
                                                                                            cgAACCACCGCTTACTGGGCGGTGTCGCCCCGATTTTCTTCCGAATCGCCCCCTGCAGCGGCGATACGGCCCACTCCGGGCGTCAGAAGGCCTAGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCcatca                                                           <  2:2128171/150‑1 (MQ=255)
                                                                                                tCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgt                                                       >  1:4444065/1‑150 (MQ=255)
                                                                                                tCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgt                                                       >  1:2059396/1‑150 (MQ=255)
                                                                                                 cACCGCTTTCTGGGCGGTGTCTCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCATGATTCCCGCCTTCAGGAAGTACACCAGCCGCTACAAGTTGTAGCAGGCCGCCAtcatcgtc                                                      >  2:3666090/1‑150 (MQ=255)
                                                                                                   ccGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCcatc                                                            <  2:4088743/142‑1 (MQ=255)
                                                                                                   ccGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCcatc                                                            >  1:4088743/1‑142 (MQ=255)
                                                                                                   ccGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgtcat                                                    >  2:4619872/1‑150 (MQ=255)
                                                                                                   ccGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgtcat                                                    >  2:6760926/1‑150 (MQ=255)
                                                                                                   ccGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgtcat                                                    >  1:3923939/1‑150 (MQ=255)
                                                                                                   ccGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgtcat                                                    >  1:3023978/1‑150 (MQ=255)
                                                                                                   ccGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgtcat                                                    >  1:1402833/1‑150 (MQ=255)
                                                                                                    cGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgtcatc                                                   >  1:1113179/1‑150 (MQ=255)
                                                                                                    cGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCAtcatcgtcatc                                                   >  1:4030695/1‑150 (MQ=255)
                                                                                                         tggggggggTGTCGCCCCTTTTTTCTTCCGAATCGCCCTCTGCCGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCaaa                                              <  2:1764549/147‑1 (MQ=255)
                                                                                                         tggggggggTGGCGCCCTTTTTTTCTTCCCACTCGCCCCGTGCCGGGGCTTTTCGCCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCaaa                                              <  1:2392980/147‑1 (MQ=255)
                                                                                                           tGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGt                                            >  2:3329962/1‑150 (MQ=255)
                                                                                                             ggCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTc                                          >  2:498532/1‑150 (MQ=255)
                                                                                                               gggTGTCGCCCCTTTTTTCTTCCTAATCGCCCTCTGCCGCGGCTTTTGGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTcgc                                        <  1:2219125/149‑1 (MQ=255)
                                                                                                               gggTGTCGCCCCGTTTTTCGCCCCAACCGCCCTCTGCCGCGGCGTTTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTcgc                                        <  1:3244335/149‑1 (MQ=255)
                                                                                                                     cGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGccc                                  <  2:4936478/150‑1 (MQ=255)
                                                                                                                      gCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGcccc                                 >  2:4373732/1‑150 (MQ=255)
                                                                                                                       ccccGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCt                                >  2:1951866/1‑150 (MQ=255)
                                                                                                                        cccGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTg                               >  1:1052920/1‑150 (MQ=255)
                                                                                                                             tttaCTTCCGAAACGCCCTCTGCTGGGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGa                          <  1:4569074/150‑1 (MQ=255)
                                                                                                                                 cttcCGACTCGCCCTCTGCAGGGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGt                      <  2:1405926/150‑1 (MQ=255)
                                                                                                                                  ttcCTAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTa                     <  2:6230263/150‑1 (MQ=255)
                                                                                                                                  ttcCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTa                     >  2:2946416/1‑150 (MQ=255)
                                                                                                                                  ttcCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATAGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTa                     >  1:3337294/1‑150 (MQ=255)
                                                                                                                                    ccGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACg                   >  2:696152/1‑150 (MQ=255)
                                                                                                                                     cGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACgg                  >  2:5280371/1‑150 (MQ=255)
                                                                                                                                     cGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACgg                  >  2:6219996/1‑150 (MQ=255)
                                                                                                                                     cGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGGTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACgg                  >  2:5138198/1‑150 (MQ=255)
                                                                                                                                     cGAATCGCCGCTGCAGAGGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACgg                  >  2:2717172/1‑150 (MQ=255)
                                                                                                                                       aaTCGCCCTCTGCAGCGGGTATTCGCCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGAt                <  1:6003745/150‑1 (MQ=255)
                                                                                                                                          cGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCaa             <  1:498532/150‑1 (MQ=255)
                                                                                                                                            cccTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGAttcact           <  2:4850375/150‑6 (MQ=255)
                                                                                                                                            cccTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACt           >  1:5809625/1‑150 (MQ=255)
                                                                                                                                             ccTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACtt          <  2:5863542/150‑1 (MQ=255)
                                                                                                                                              ctctGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACTTg         >  2:3093298/1‑150 (MQ=255)
                                                                                                                                               tctGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACTTgc        >  2:4303905/1‑150 (MQ=255)
                                                                                                                                                 tGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACTTgccg      >  2:5309983/1‑150 (MQ=255)
                                                                                                                                                  gCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACTTgccgc     <  2:3529002/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGTTGTTTATGGGCACCTCGAATAACGCGAGGTTTTCAGTGAATCCGACCGCCAAATGATGGCCATGATCGAATTCGGTTCGATTTCTGCGGCGAATCACCGCTTTCTGGGCGGTGTCGCCCCGTTTTTCTTCCGAATCGCCCTCTGCAGCGGCGATTCGGCCCACTCCGGGCGTCAGAAGGCCTTGATTCCCGCCTTCTGGAAGTACACCAGCCGCTTCAAGTTGTAGCAGGCCGCCATCATCGTCATCGCAAAGTTCGCGCGCCCCTGTCCGATGGTACGGATCAACTTGCCGCCCA  >  NC_002947/3377660‑3377958

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