Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 732,557 C→T 8.3% R698R (CGC→CGT clpB → chaperone protein ClpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947732,5570CT 8.3% 792.0 / 41.8 266R698R (CGC→CGTclpBchaperone protein ClpB
Reads supporting (aligned to +/- strand):  ref base C (146/97);  new base T (15/7);  total (162/104)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

GCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTGGGGGCGCACTTCCG  >  NC_002947/732408‑732706
                                                                                                                                                     |                                                                                                                                                      
  accaccAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCt                                                                                                                                                      >  2:5016700/1‑150 (MQ=255)
  accaccAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCt                                                                                                                                                      >  2:182672/1‑150 (MQ=255)
     tcccggttATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                                                                                                                                   >  2:1721377/8‑150 (MQ=255)
        aGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGa                                                                                                                                                <  1:8373731/150‑1 (MQ=255)
         ggCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGAc                                                                                                                                               >  1:1260647/1‑150 (MQ=255)
              tGTAGGGTATGAAGAGGCCGGTTACCTGACCTATGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGACGTGGAAAAAGCCCACCCGGATGAGTTCAACGTGCTGTTGCAGGTGCTGTAAGACGGCAGCCTGACCGACAGccc                                                                                                                                          >  2:5355103/1‑149 (MQ=255)
                tAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACg                                                                                                                                        >  2:8823311/1‑150 (MQ=255)
                 aGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACgg                                                                                                                                       >  2:5740745/1‑150 (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)
                       tGAAGAGGTCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCAc                                                                                                                                >  1:8197970/1‑150 (MQ=255)
                       tGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCAc                                                                                                                                >  1:1728502/1‑150 (MQ=255)
                        gAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCAcc                                                                                                                               >  1:2910125/1‑150 (MQ=255)
                        gAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCAcc                                                                                                                               >  1:3048433/1‑150 (MQ=255)
                         aaGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCg                                                                                                                              >  2:4671342/1‑150 (MQ=255)
                           gagGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTg                                                                                                                            >  1:3684047/1‑150 (MQ=255)
                           gagGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTg                                                                                                                            <  2:1122795/150‑1 (MQ=255)
                           gagGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTg                                                                                                                            >  1:4707677/1‑150 (MQ=255)
                           gagGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTg                                                                                                                            <  1:6434251/150‑1 (MQ=255)
                            agGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGTGCGCCCGTggc                                                                                                                          <  1:1520500/150‑2 (MQ=255)
                             gggCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGa                                                                                                                          <  2:360506/150‑1 (MQ=255)
                              ggCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGACAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGTGCGCCCGTggct                                                                                                                         <  2:5646240/150‑3 (MQ=255)
                              ggCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGAc                                                                                                                         >  1:5861312/1‑150 (MQ=255)
                                 ggTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTc                                                                                                                      >  2:1664258/1‑150 (MQ=255)
                                  gcgttcTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTcc                                                                                                                     <  2:3800141/145‑1 (MQ=255)
                                  gTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTcc                                                                                                                     <  1:5915628/150‑1 (MQ=255)
                                    tACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGc                                                                                                                   <  1:4420156/150‑1 (MQ=255)
                                      ccTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCaa                                                                                                                 <  2:357323/150‑1 (MQ=255)
                                        tGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTCCGCAaca                                                                                                               <  2:4553908/150‑1 (MQ=255)
                                        tGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAaca                                                                                                               <  2:6724865/150‑1 (MQ=255)
                                          aCCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACAcc                                                                                                             >  1:2244064/1‑150 (MQ=255)
                                            cGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgt                                                                                                           >  2:712285/1‑150 (MQ=255)
                                             gAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtg                                                                                                          >  1:2508804/1‑150 (MQ=255)
                                                   gTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtgatcgtg                                                                                                    <  1:4135697/150‑1 (MQ=255)
                                                    tGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCTGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtgatcgtga                                                                                                   <  1:4978725/150‑1 (MQ=255)
                                                    tGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtgatcgtga                                                                                                   >  2:586315/1‑150 (MQ=255)
                                                    tGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtgatcgtga                                                                                                   <  2:4829005/150‑1 (MQ=255)
                                                       ggcgCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGtgatga                                                                                                <  2:8671617/150‑1 (MQ=255)
                                                        gcgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGAc                                                                                               <  1:4020200/150‑1 (MQ=255)
                                                          gcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCt                                                                                             >  1:5299426/1‑150 (MQ=255)
                                                          gcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCt                                                                                             >  1:3378306/1‑150 (MQ=255)
                                                                cTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAAcc                                                                                       >  2:2760747/1‑150 (MQ=255)
                                                                cTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAAcc                                                                                       >  1:9105266/1‑150 (MQ=255)
                                                                  tatcaggtgatCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTCCGCAACACGGTCATCATCATGACCTCGAACCtc                                                                                     >  2:7838171/6‑149 (MQ=255)
                                                                  tatcaggtgatCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTCCGCAACACGGTCATCATCATGACCTCGAACCtc                                                                                     >  1:7210281/6‑149 (MQ=255)
                                                                   aCTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTgg                                                                                    <  2:2910125/150‑1 (MQ=255)
                                                                       ggtggtGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTc                                                                                >  2:5087768/1‑150 (MQ=255)
                                                                       ggtggtGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTc                                                                                >  2:8322058/1‑150 (MQ=255)
                                                                       ggtggtGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTc                                                                                >  2:7897116/1‑150 (MQ=255)
                                                                       ggtggtGCTGCTGGACGACCTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCCCCGTGGACTTCCGCACCACCGTGATCGTGATGACCTCCAACCTGGGcgc                                                                                >  2:6136072/1‑148 (MQ=255)
                                                                             gctgctGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCa                                                                          <  1:7466987/150‑1 (MQ=255)
                                                                               tgctgGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGgcg                                                                            <  1:2648319/146‑1 (MQ=255)
                                                                               tgctgGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGa                                                                        >  1:1934273/1‑150 (MQ=255)
                                                                                gctgGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGAt                                                                       >  1:2485306/1‑150 (MQ=255)
                                                                                   ggACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCa                                                                    <  2:6304960/150‑1 (MQ=255)
                                                                                   ggACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCa                                                                    <  2:4209830/150‑1 (MQ=255)
                                                                                   ggACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCa                                                                    <  1:7518381/150‑1 (MQ=255)
                                                                                   ggACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAACCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCa                                                                    <  1:1668960/150‑1 (MQ=255)
                                                                                    gACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAg                                                                   <  1:6715528/150‑1 (MQ=255)
                                                                                       gAGGTGGAAAAAGTCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGaa                                                                >  2:1560794/1‑150 (MQ=255)
                                                                                       gAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGaa                                                                <  1:4673659/150‑1 (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)
                                                                                                ataGCCTACGCGGATGTGTTCATCGTGCTGTTGCAGGTGCTGTTAGACGGCCGCCTTTAGGACTGCCTCGGGCGCTGCTTGGTGTTCGGCTACACCGTGATCGTGATGACCTCCTACCTGGGCTCGGCGCAGATTCAGGTAGTggtgggt                                                       <  2:8722537/148‑1 (MQ=255)
                                                                                                aaaGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTggtgggt                                                       <  1:8596529/150‑1 (MQ=255)
                                                                                                 aaGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTggtgggtg                                                      >  1:8843757/1‑150 (MQ=255)
                                                                                                  aGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGa                                                     <  2:8873738/150‑1 (MQ=255)
                                                                                                   gCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGAc                                                    >  1:6559361/1‑150 (MQ=255)
                                                                                                   gCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGAc                                                    >  1:8431857/1‑150 (MQ=255)
                                                                                                     ccACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACcg                                                  >  2:6702649/1‑150 (MQ=255)
                                                                                                        cccGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGa                                               <  2:8596929/150‑1 (MQ=255)
                                                                                                         ccGGATGTGTTCATCGTGCTGTTGCAGGTGCTGTAAGACGGCCGCCTGACCGACAGCCACGGGCTCAACGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCTCAGATTCAGGAACTGGTGGGTGACCGCGAg                                              <  1:5355103/150‑1 (MQ=255)
                                                                                                         ccGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAg                                              >  2:2083116/1‑150 (MQ=255)
                                                                                                             aTGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGcac                                          >  2:7758604/1‑150 (MQ=255)
                                                                                                                  ttCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGACCTGGTGGGTGACCGCGAGGCACAGCGt                                     >  2:9002037/1‑150 (MQ=255)
                                                                                                                  ttCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGt                                     >  2:7302848/1‑150 (MQ=255)
                                                                                                                   tCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTg                                    >  2:2013095/1‑150 (MQ=255)
                                                                                                                    cAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAgg                                             <  1:6124058/140‑1 (MQ=255)
                                                                                                                     aaCGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCg                                  >  2:3102400/1‑150 (MQ=255)
                                                                                                                     aaCGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCg                                  >  2:1415776/1‑150 (MQ=255)
                                                                                                                       cGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGACCTGGTGGGTGACCGCGAGGCACAGCGTGCGGc                                >  1:8078138/1‑150 (MQ=255)
                                                                                                                       cGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGc                                <  2:2508804/150‑1 (MQ=255)
                                                                                                                         tGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAg                              <  2:7397733/150‑1 (MQ=255)
                                                                                                                                 cAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTgatggat                      >  1:5068832/1‑150 (MQ=255)
                                                                                                                                      gCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGt                 >  2:6777372/1‑150 (MQ=255)
                                                                                                                                          gAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTgggg             <  2:1786265/150‑1 (MQ=255)
                                                                                                                                                ggCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTGGGGGCGCAc       >  2:7708199/1‑150 (MQ=255)
                                                                                                                                                ggCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTGGGGGCGCAc       >  1:7931950/1‑150 (MQ=255)
                                                                                                                                                ggCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTGGGGGCGCAc       >  1:7325546/1‑150 (MQ=255)
                                                                                                                                                ggCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTGGGGGCGCAc       >  2:2667882/1‑150 (MQ=255)
                                                                                                                                                ggCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTGGGGGCGCAc       >  2:3401123/1‑150 (MQ=255)
                                                                                                                                                ggCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTGGGGGCGCAc       >  1:8759315/1‑150 (MQ=255)
                                                                                                                                                ggCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTACGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTGGGGGCGCAc       >  2:220750/1‑150 (MQ=255)
                                                                                                                                                 gccgccTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGaa                                                                <  2:3850674/92‑1 (MQ=255)
                                                                                                                                                 gccgccTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGaa                                                                >  1:3850674/1‑92 (MQ=255)
                                                                                                                                                 gccgccTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTGGGGGCGCACt      >  2:4403519/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGCGCAGATTCAGGAACTGGTGGGTGACCGCGAGGCACAGCGTGCGGCAGTGATGGATGCGGTGGGGGCGCACTTCCG  >  NC_002947/732408‑732706

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: