Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 732,548 A→T 9.7% E695D (GAA→GAT clpB → chaperone protein ClpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947732,5480AT 9.7% 818.9 / 59.5 267E695D (GAA→GATclpBchaperone protein ClpB
Reads supporting (aligned to +/- strand):  ref base A (140/100);  new base T (15/11);  total (156/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 273 total aligned reads displayed.

CTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTGGGGGC  >  NC_002947/732399‑732697
                                                                                                                                                     |                                                                                                                                                      
 tGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAg                                                                                                                                                       >  2:4894101/1‑150 (MQ=255)
    ttggTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACg                                                                                                                                                    <  1:7392211/148‑1 (MQ=255)
         gCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgccgc                                                                                                                                               <  2:1751409/150‑1 (MQ=255)
          caccacCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgccgcc                                                                                                                                              >  2:8708341/1‑150 (MQ=255)
           accaccAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCt                                                                                                                                             <  1:6263853/150‑1 (MQ=255)
              tcccggttATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                                                                                                                          >  2:1721377/8‑150 (MQ=255)
              tcccggttATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                                                                                                                          >  1:508177/8‑150 (MQ=255)
              tcccggttATGTCGGCTATGAAGAGGGTGGCGTTCTGAAGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                                                                                                                          >  2:7455129/8‑150 (MQ=255)
              accaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGAc                                                                                                                                          >  2:2872114/1‑150 (MQ=255)
               ccaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGAcc                                                                                                                                         <  2:2963930/150‑1 (MQ=255)
                  ggttATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGAc                                                                                                                                      >  1:8109158/4‑150 (MQ=255)
                  ggCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGAc                                                                                                                                      >  1:1260647/1‑150 (MQ=255)
                   gCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACa                                                                                                                                     >  1:2669523/1‑150 (MQ=255)
                    cTATGTAGGGTATGAAGAGGGCGGTGACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAg                                                                                                                                    <  1:8740818/150‑1 (MQ=255)
                           gggTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACggg                                                                                                                            >  2:4413221/1‑150 (MQ=255)
                           ggctATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                                                                                                                          <  1:7825273/134‑1 (MQ=255)
                            ggTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGgc                                                                                                                           <  1:5527061/150‑1 (MQ=255)
                             gTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACCGCCACGGgcg                                                                                                                          >  2:3900902/1‑150 (MQ=255)
                               aTGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCa                                                                                                                        >  1:8534463/1‑150 (MQ=255)
                                tGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCAc                                                                                                                       >  1:1728502/1‑150 (MQ=255)
                                  aaGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCg                                                                                                                     >  2:8734498/1‑150 (MQ=255)
                                  aaGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCg                                                                                                                     >  2:4671342/1‑150 (MQ=255)
                                    gagGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTg                                                                                                                   <  2:1122795/150‑1 (MQ=255)
                                     agGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGTGCGCCCGTggc                                                                                                                 <  1:1520500/150‑2 (MQ=255)
                                      gggCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGa                                                                                                                 <  2:360506/150‑1 (MQ=255)
                                       ggCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGTGCGCCCGTggct                                                                                                                <  2:5647422/150‑3 (MQ=255)
                                          ggTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTc                                                                                                             >  1:6805533/1‑150 (MQ=255)
                                           gcgttcTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTcc                                                                                                            <  2:3800141/145‑1 (MQ=255)
                                           gcgtacTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTcc                                                                                                            <  1:7571084/145‑1 (MQ=255)
                                           gTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTcc                                                                                                            >  2:6873872/1‑150 (MQ=255)
                                             tACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGc                                                                                                          <  1:4420156/150‑1 (MQ=255)
                                               tcTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTCCGCaa                                                                                                        >  1:5331113/2‑150 (MQ=255)
                                               ccTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCaa                                                                                                        <  2:357323/150‑1 (MQ=255)
                                                 tGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTCCGCAaca                                                                                                      <  2:4553908/150‑1 (MQ=255)
                                                   aCCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACAcc                                                                                                    >  1:2244064/1‑150 (MQ=255)
                                                     cGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgt                                                                                                  >  1:2535668/1‑150 (MQ=255)
                                                      gAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtg                                                                                                 >  2:8196699/1‑150 (MQ=255)
                                                            gTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtgatcgtg                                                                                           >  2:866860/1‑150 (MQ=255)
                                                             tGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtgatcgtga                                                                                          >  2:586315/1‑150 (MQ=255)
                                                             tGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtgatcgtga                                                                                          <  2:5481138/150‑1 (MQ=255)
                                                             tGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtgatcgtga                                                                                          <  2:4829005/150‑1 (MQ=255)
                                                                 gcgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGCCTTCCGCAACACCGTGATCGTGATGAc                                                                                      <  2:1465243/150‑1 (MQ=255)
                                                                 gcgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGAc                                                                                      <  1:4020200/150‑1 (MQ=255)
                                                                  cgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGAcc                                                                                     >  1:8068007/1‑150 (MQ=255)
                                                                   gcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCt                                                                                    >  1:5299426/1‑150 (MQ=255)
                                                                   gcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCt                                                                                    >  1:3378306/1‑150 (MQ=255)
                                                                     aaGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTcc                                                                                  <  1:3804996/150‑1 (MQ=255)
                                                                         catACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGATCAACCTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCAGCAACACCGTGATCGTGATGACCTCCAAcc                                                                              >  1:2884849/3‑150 (MQ=255)
                                                                         cTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAAcc                                                                              >  2:2760747/1‑150 (MQ=255)
                                                                         cTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAAcc                                                                              >  1:9105266/1‑150 (MQ=255)
                                                                          atatcaggtgatCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTCCGCAACACGGTCATCATCATGACCTCGAACCt                                                                             >  2:1809796/7‑150 (MQ=255)
                                                                           tatcaggtgatCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTCCGCAACACGGTCATCATCATGACCTCGAACCtc                                                                            >  2:7838171/6‑149 (MQ=255)
                                                                              tCGGTGGTGCTGCTGGACGTGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGACGACGGCCGCCTGACCGACAGCCACGGGCGCTCCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGc                                                                         >  2:1593344/1‑150 (MQ=255)
                                                                               aggtgatCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTCCGCAACACGGTCATCATCATGACCTCGAACCTCGGcg                                                                        >  1:8375626/2‑149 (MQ=255)
                                                                                ggtggtGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTc                                                                       >  2:7897116/1‑150 (MQ=255)
                                                                                 gtggtgCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACtt                                                                                                              <  1:363873/110‑1 (MQ=255)
                                                                                        tgctgGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGa                                                               <  1:5939715/150‑1 (MQ=255)
                                                                                        tgctgGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGa                                                               >  1:1934273/1‑150 (MQ=255)
                                                                                         gctgGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGAt                                                              >  1:2485306/1‑150 (MQ=255)
                                                                                            ggACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCa                                                           >  1:5871199/1‑150 (MQ=255)
                                                                                             gACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAg                                                          <  1:5825801/150‑1 (MQ=255)
                                                                                             gACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAg                                                          <  2:6708377/150‑1 (MQ=255)
                                                                                             gACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAg                                                          <  2:2432918/150‑1 (MQ=255)
                                                                                             gACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAg                                                          <  1:495342/150‑1 (MQ=255)
                                                                                             gACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAg                                                          <  1:6715528/150‑1 (MQ=255)
                                                                                                gAGGTGGAAAAAGTCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGaa                                                       >  2:6163824/1‑150 (MQ=255)
                                                                                                gAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGaa                                                       <  1:4673659/150‑1 (MQ=255)
                                                                                                   gTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtg                                                    >  1:1553005/1‑150 (MQ=255)
                                                                                                   gTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtg                                                    >  1:1432047/1‑150 (MQ=255)
                                                                                                    tGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtgg                                                   <  2:3499686/150‑1 (MQ=255)
                                                                                                     ggAAAAAGCCCACCCGGATGTGTTCAACGTGCTGGTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGAGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCAGGGCTCGGCGCAGATTCAGGAACtggt                                                  >  1:1167579/1‑150 (MQ=255)
                                                                                                      gAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtggtg                                                 <  1:4413209/150‑1 (MQ=255)
                                                                                                       aaaaaGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtggtgg                                                <  2:7212950/150‑1 (MQ=255)
                                                                                                       aaaaaGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtggtgg                                                <  2:1897390/150‑1 (MQ=255)
                                                                                                       aaaaaGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACtggtgg                                                <  1:7498262/150‑1 (MQ=255)
                                                                                                          aaGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTggtgggtg                                             >  1:296816/1‑150 (MQ=255)
                                                                                                          aaGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTggtgggtg                                             >  1:8843757/1‑150 (MQ=255)
                                                                                                           aGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGa                                            <  2:8873738/150‑1 (MQ=255)
                                                                                                            gCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGAc                                           >  1:6559361/1‑150 (MQ=255)
                                                                                                                aCCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACcgcg                                       >  2:3869627/1‑150 (MQ=255)
                                                                                                                aCCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACcgcg                                       >  2:3331149/1‑150 (MQ=255)
                                                                                                                 cccGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGa                                      <  1:5016700/150‑1 (MQ=255)
                                                                                                                  ccGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAg                                     >  2:5084290/1‑150 (MQ=255)
                                                                                                                  ccGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAg                                     <  2:419801/150‑1 (MQ=255)
                                                                                                                      aTGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGcac                                 >  2:7758604/1‑150 (MQ=255)
                                                                                                                           ttCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGACCTGGTGGGTGACCGCGAGGCACAGCGt                            >  2:9002037/1‑150 (MQ=255)
                                                                                                                           ttCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGt                            >  2:2519595/1‑150 (MQ=255)
                                                                                                                           ttCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGt                            >  2:6676640/1‑150 (MQ=255)
                                                                                                                             cAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAgg                                    >  2:6124058/1‑140 (MQ=255)
                                                                                                                             cAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGc                          <  1:2298003/150‑1 (MQ=255)
                                                                                                                              aaCGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCg                         >  2:1415776/1‑150 (MQ=255)
                                                                                                                                 gaggtGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCa                      >  2:2061545/5‑150 (MQ=255)
                                                                                                                                  tGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAg                     <  2:7397733/150‑1 (MQ=255)
                                                                                                                                      gTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTgat                 <  1:425405/150‑1 (MQ=255)
                                                                                                                                         gCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTgatgga              >  2:8449569/1‑150 (MQ=255)
                                                                                                                                         gCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTgatgga              <  2:212283/150‑1 (MQ=255)
                                                                                                                                             gTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCg          >  2:4315905/1‑150 (MQ=255)
                                                                                                                                              ttcTGGAAGACGGCCGACTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATAGTGATGACCTCCAACCTGGGATCGGCGCAGATTCAGGAACTGGAGGGTGACCGCGAGGCACAGCGTCCGGCAGTGATGGATGCgg         >  2:5343319/3‑150 (MQ=255)
                                                                                                                                               gCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGt        >  2:6776864/1‑150 (MQ=255)
                                                                                                                                                  ggAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTggg     >  1:8295936/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTGGGGGC  >  NC_002947/732399‑732697

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: