Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,206,908 A→C 5.5% Q119H (CAA→CAC xcpW → type II secretion pathway protein XcpW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,206,9080AC 5.5% 2222.4 / 13.1 710Q119H (CAA→CACxcpWtype II secretion pathway protein XcpW
Reads supporting (aligned to +/- strand):  ref base A (417/228);  new base C (2/63);  total (419/291)
Fisher's exact test for biased strand distribution p-value = 4.23e-24
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 733 total aligned reads displayed.

TACGACCAGCGTCTGCCGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGCTGGCCTGGCCTTCTGACCAACCGGTGCA  >  NC_002947/1206759‑1207057
                                                                                                                                                     |                                                                                                                                                     
tACGACCAGCGTCTGCCGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCaa                                                                                                                                                       >  1:8083564/1‑150 (MQ=255)
 aCGACCAGCGTCTGCCGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAAc                                                                                                                                                      <  1:7220593/150‑1 (MQ=255)
    aCCAGCGTCTGCCGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGcc                                                                                                                                                   >  2:6362152/1‑150 (MQ=255)
     ccAGCGTCTGCCGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCct                                                                                                                                                  >  1:5973767/1‑150 (MQ=255)
      cAGCGTCTGCCGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCctg                                                                                                                                                 >  2:7403778/1‑150 (MQ=255)
         cGTCTGCCGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCctgctg                                                                                                                                              <  2:2108294/150‑1 (MQ=255)
            cTGCCGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCctgctgctg                                                                                                                                           <  1:5175584/150‑1 (MQ=255)
            cTGCCGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCctgctgctg                                                                                                                                           <  1:6172602/150‑1 (MQ=255)
            cTGCCGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCctgctgctg                                                                                                                                           <  1:1241007/150‑1 (MQ=255)
               ccGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGaa                                                                                                                                        >  1:6492928/1‑150 (MQ=255)
                cGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGaaa                                                                                                                                       <  2:4363082/150‑1 (MQ=255)
                 gCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAg                                                                                                                                      >  1:1507300/1‑150 (MQ=255)
                   tGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCt                                                                                                                                    <  1:2202107/150‑1 (MQ=255)
                   agtccAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCt                                                                                                                                    >  2:7259926/4‑150 (MQ=255)
                      ccAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTcgg                                                                                                                                 <  1:7928300/150‑1 (MQ=255)
                       cAGCCGTTGGCTGGACTGGCGGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTcggc                                                                                                                                >  2:8338562/1‑150 (MQ=255)
                       cAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTcggc                                                                                                                                >  2:2832321/1‑150 (MQ=255)
                       cAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTcggc                                                                                                                                <  1:8779466/150‑1 (MQ=255)
                       cAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTcggc                                                                                                                                <  1:7644982/150‑1 (MQ=255)
                       cAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTcggc                                                                                                                                >  2:12540/1‑150 (MQ=255)
                       cAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTcggc                                                                                                                                >  1:2581992/1‑150 (MQ=255)
                        aGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTcggcg                                                                                                                               >  2:4901893/1‑150 (MQ=255)
                           cGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGAc                                                                                                                            >  2:7796115/1‑150 (MQ=255)
                           cGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGAc                                                                                                                            >  2:4577984/1‑150 (MQ=255)
                           cGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGAc                                                                                                                            <  2:3071450/150‑1 (MQ=255)
                             ttGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCt                                                                                                                          <  2:326896/150‑1 (MQ=255)
                                gCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTa                                                                                                                       <  1:141986/150‑1 (MQ=255)
                                 cTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTAc                                                                                                                      >  1:7721878/1‑150 (MQ=255)
                                  tGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTAcg                                                                                                                     >  1:4944878/1‑150 (MQ=255)
                                    gACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTAcgcc                                                                                                                   >  1:8495436/1‑150 (MQ=255)
                                    gACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTAcgcc                                                                                                                   >  1:3070413/1‑150 (MQ=255)
                                      cTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTAcgccgc                                                                                                                 >  1:1783473/1‑150 (MQ=255)
                                       tGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGTACAGGTccgccgag                                                                                                                <  2:6868521/150‑3 (MQ=255)
                                         gCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTgg                                                                                                              >  2:6228156/1‑150 (MQ=255)
                                          cAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGgc                                                                                                             <  2:3788175/150‑1 (MQ=255)
                                            gCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGTACAGGTccgccgagctttc                                                                                                           <  1:19286/150‑8 (MQ=255)
                                            gCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGgcgc                                                                                                           <  2:3884931/150‑1 (MQ=255)
                                              ccGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGccc                                                                                                         >  2:4646954/1‑150 (MQ=255)
                                              ccGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGccc                                                                                                         >  1:2800507/1‑150 (MQ=255)
                                               ccacAGCAACTGGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCTGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGcccc                                                                                                        >  2:564633/4‑150 (MQ=255)
                                               cGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGcccc                                                                                                        <  1:532684/150‑1 (MQ=255)
                                               cGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGcccc                                                                                                        >  2:2131123/1‑150 (MQ=255)
                                               cGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGcccc                                                                                                        >  2:80203/1‑150 (MQ=255)
                                                gCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCcg                                                                                                       >  1:4603022/1‑150 (MQ=255)
                                                   aGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCgcgc                                                                                                    <  1:5701673/150‑1 (MQ=255)
                                                   aGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCgcgc                                                                                                    >  1:4556436/1‑150 (MQ=255)
                                                    gCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCt                                                                                                   >  1:3189918/1‑150 (MQ=255)
                                                     cAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTg                                                                                                  >  2:8714992/1‑150 (MQ=255)
                                                        ccAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGt                                                                                               >  2:6827202/1‑150 (MQ=255)
                                                            ttGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAg                                                                                           <  2:4117517/150‑1 (MQ=255)
                                                             tGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAgc                                                                                          <  1:1609479/150‑1 (MQ=255)
                                                             tGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAgc                                                                                          <  1:3584256/150‑1 (MQ=255)
                                                                  tgGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGaa                                                                                     <  1:5879794/150‑1 (MQ=255)
                                                                  tgGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGaa                                                                                     <  1:6396141/150‑1 (MQ=255)
                                                                  tgGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGaa                                                                                     <  2:6065755/150‑1 (MQ=255)
                                                                  tgGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGaa                                                                                     <  2:5029594/150‑1 (MQ=255)
                                                                  tgGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGACTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGaa                                                                                     <  2:6618336/150‑1 (MQ=255)
                                                                   gggTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAAc                                                                                    <  2:560272/150‑1 (MQ=255)
                                                                   gggTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAAc                                                                                    <  1:2918910/150‑1 (MQ=255)
                                                                   gggTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCACCTGCAACGCCAGTGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAAc                                                                                    >  1:6956815/1‑150 (MQ=255)
                                                                     gTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGa                                                                                  >  1:2308880/1‑150 (MQ=255)
                                                                      tCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGAc                                                                                 >  1:8168221/1‑150 (MQ=255)
                                                                        gCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTc                                                                               >  2:5183011/1‑150 (MQ=255)
                                                                          cAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGa                                                                             >  2:338232/1‑150 (MQ=255)
                                                                           aGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGAc                                                                            >  1:8302740/1‑150 (MQ=255)
                                                                               ccGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAgg                                                                        >  1:3867101/1‑150 (MQ=255)
                                                                                cGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGt                                                                       >  1:418458/1‑150 (MQ=255)
                                                                                     tGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAgcg                                                                  >  1:7847204/1‑150 (MQ=255)
                                                                                         ggCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCAt                                                              >  2:473867/1‑150 (MQ=255)
                                                                                             gCCGACGGTGGCTCGCGCCTGCAACGCCAGCGCCTGGGTTTCGCGGCCCGCGAGCAACGCCTGCTGCTGTAACGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCGGGTTACCGGGCCAACGACTCGACCAGGTCAGCGCCATGAAc                                                          <  2:3649054/150‑1 (MQ=255)
                                                                                                   ggTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTCCAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACgttctg                                                    >  2:182991/1‑150 (MQ=255)
                                                                                                    ggggCCCGCGCGTGCACGCCCCGCGCGTGGGGTGCGGGCCCGGCGAGCACGGCCTGCGGGTGGAACGCTGGGCGGCCGTGTCCGCCGGTGGCGCCCGGCGGGGTTTACGGCGCCCACGACTCGACCAGGTCAGCGCCATGAACgttctgt                                                   <  2:293450/148‑1 (MQ=255)
                                                                                                     gggCTCGCGCCTCCACCGCCCGCGCCTGCGGTTGGTGGCCCGGGCGCAACGCCTGGTGCTGGAACGCTCGGCGGACGTGTACGCCGCTGGCGCCCCGCGGTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACgttctgtt                                                  <  2:6441281/149‑1 (MQ=255)
                                                                                                      ggCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACgttctgttc                                                 >  2:8555313/1‑150 (MQ=255)
                                                                                                         tCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTAc                                              >  1:5400884/1‑150 (MQ=255)
                                                                                                         gcgcgCCTGCAACGCCAGCGCCGGCGGTTCGCGGCCCGCGAGCAACGCCTGCGGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTAc                                              <  2:5721490/149‑1 (MQ=255)
                                                                                                           gcgcCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCa                                            >  1:3443675/1‑150 (MQ=255)
                                                                                                           gcgcCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCa                                            >  1:5055514/1‑150 (MQ=255)
                                                                                                             gcCTGCACCGCCCGCGCTGGGGTTTCGCGGCCCGCGACAACCGCGGGGGGGGGCAACGCGGGGCGCACGTGTCCGCCGGGGGCGCCCGGGGGTGGGTCCCGCGCCCACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAgg                                          <  2:2581992/150‑1 (MQ=255)
                                                                                                                   aaCGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGcc                                    >  2:7944489/1‑150 (MQ=255)
                                                                                                                   aaCGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGcc                                    >  2:8136213/1‑150 (MQ=255)
                                                                                                                   aaCGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGcc                                    >  2:7152257/1‑150 (MQ=255)
                                                                                                                        cAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCgctgg                               >  2:7592640/1‑150 (MQ=255)
                                                                                                                          gcgcGTGCGTTTCGCGGCCCGCGACCAACGCGTGCTGCTGTAAAGCGCGGGGCACCTTTACGCCGCGGGCGCCCCGCGCTGGTTACCGCGCCCACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCgctggct                             <  2:5345787/150‑1 (MQ=255)
                                                                                                                          gcgcCTGGGGTTCGCGGCCCGCGCGCACCGCGTGGTGGTGGAAGGGTCGGGGGACCTGGCCGCGGGGGGCGCCCCGCGGGGGGTCCCGCGCGCACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCgctggct                             <  2:449605/150‑1 (MQ=255)
                                                                                                                           cgcgTGGGGTTCGCGGCCCGCTAGCAACGCCTGCTGCTGGAATGCTGGGCGGACGTGTACGCCGCTGGCGCCCGGCGGTGGGTACCGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCgctggctg                            <  2:1922094/150‑1 (MQ=255)
                                                                                                                           cgcCTGCGGTGCGCGGCCCGGTACCAACCCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCGGGCGCCCCGCGGGGGGTCCTGCGCGCACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCgctggctg                            <  1:8668531/150‑1 (MQ=255)
                                                                                                                            gcCGGCGGTTCGCGGCCCGCGCGCACCGCCTGCTGGTGGAACGCTCGGCGGCCCTGCACGCCGCTGGCGCCCCGCGCTGGGTACCGCGCCAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGctggctgg                           <  2:5845480/150‑1 (MQ=255)
                                                                                                                               tGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGctggcct                        >  1:7354487/1‑150 (MQ=255)
                                                                                                                                ggggTTCGTGGCCCGCGAGCACCGCCTGGTGGTTGAAAGCGGGGGGGACCTTTACGCCGGGGGCGCCCCGGGGTGGTTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGctggcctg                       <  2:1578128/148‑1 (MQ=255)
                                                                                                                                ggggTTCGCGGCCCGCGAGCACCGCCTGCTGCGGCAACGCGCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCGGGGCACCGCGCCAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGctggcctg                       <  2:1353429/148‑1 (MQ=255)
                                                                                                                                       gcggCCCGCGCGCACGGCCTGCGGCGGGACCGCGCGGGGGCCCTGGACGCCGGGGGCGCCCCGCGGGGGGTCCCGCGCGCACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGCTGGCCTGGCCTTCt                <  2:7719827/148‑1 (MQ=255)
                                                                                                                                       gcgcccggcgACCACCGCCTGGTGGTGTAACGGGCGGGGGACCTGTCCGCGGGGGGCGCCCCGCGGGGGGTCCCGCGCCCACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGCTGGCCTGGCCTTCt                <  1:552455/143‑1 (MQ=255)
                                                                                                                                        aGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGCTGGCCTGGCCTTCTg               >  2:6969516/1‑150 (MQ=255)
                                                                                                                                         ggCCCGCGCCCACCGCCTGCTGGTGGAACGCGGGGGGGACCTGCCCGCCGGTGGCGCCCCGGGGTGGGTACCGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGCTGGCCTGGCCTTCTGa              <  1:6576435/150‑1 (MQ=255)
                                                                                                                                         ggCCCGCGAGCACCGCGTGGTGCTGGACCGCGGGGCGCACCTGTCGGCCGCGGGCGCCCGGGGGTGGGTACCGCGCCCACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGCTGGCCTGGCCTTCTGa              <  2:1044334/150‑1 (MQ=255)
                                                                                                                                         ggCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGCTGGCCTGGCCTTCTGa              >  1:1031828/1‑150 (MQ=255)
                                                                                                                                          gCCCGCGACCACCGCGTGCTGCTGGAAAGCGCGGCGCACCTGTCCGCCGCGGGCGCCCCGCGGGGGTCCCCGCGCGCACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGCTGGCCTGGCCTTCTGAc             <  1:6763759/150‑1 (MQ=255)
                                                                                                                                             cgcgAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGCTGGCCTGGCCTTCTGACCaa          >  1:5266326/1‑150 (MQ=255)
                                                                                                                                                     aCGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGCTGGCCTGGCCTTCTGACCAACCGGTGCa  >  1:2741616/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACGACCAGCGTCTGCCGCTGGCCAGCCGTTGGCTGGACTGGCAGCCCGCCAGCAACCAGTTGCTGTGGGTCGCCAGCACCGGCCTGCCGGCAGCCGAAGGTGGCTCGCGCCTGCAACGCCAGCGCCTGCGGTTCGAGGCCCGCGAGCAACGCCTGCTGCTGGAAAGCTCGGCGGACCTGTACGCCGCTGGCGCCCCGCGCTGGGTACAGCGCGAACGACTCGACCAGGTCAGCGCCATGAACGTTCTGTTCTACCAGGGCGGCCGCTGGCTGGCCTGGCCTTCTGACCAACCGGTGCA  >  NC_002947/1206759‑1207057

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: