Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,146,043:1 +A 7.8% coding (787/2214 nt) PP_1897 → DNA internalization‑related competence protein ComEC/Rec2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,146,0431.A 7.8% 1512.1 / 121.0 372coding (787/2214 nt)PP_1897DNA internalization‑related competence protein ComEC/Rec2
Reads supporting (aligned to +/- strand):  ref base . (175/168);  new base A (16/13);  total (191/181)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 383 total aligned reads displayed.

GTTGCAGGCCACCGGCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGTTGCACTGGTC  >  NC_002947/2145895‑2146192
                                                                                                                                                     |                                                                                                                                                     
gTTGCAGGCCACCGGCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCt                                                                                                                                                      <  2:179290/150‑1 (MQ=255)
gTTGCAGGCCACCGGCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCt                                                                                                                                                      <  2:3317825/150‑1 (MQ=255)
gTTGCAGGCCACCGGCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCt                                                                                                                                                      >  2:5371414/1‑150 (MQ=255)
  tGCAGGCCACCGGCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTgg                                                                                                                                                    <  2:5578142/150‑1 (MQ=255)
       gCCACCGGCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATg                                                                                                                                               >  2:8225999/1‑150 (MQ=255)
       gCCACCGGCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATg                                                                                                                                               >  1:2152509/1‑150 (MQ=255)
       gCCACCGGCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATg                                                                                                                                               >  2:3853823/1‑150 (MQ=255)
         cACCGGCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGc                                                                                                                                             >  2:6954739/1‑150 (MQ=255)
         cACCGGCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGc                                                                                                                                             >  2:5721147/1‑150 (MQ=255)
         cACCGGCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCCGGGCCTGCGGCCTGGCCATGGc                                                                                                                                             >  2:2065612/1‑150 (MQ=255)
         cACCGGCACGGTGCACCTGTTGGTGAATTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCCGGGCCTGCGGCCTGGCCATGGc                                                                                                                                             >  2:2250865/1‑150 (MQ=255)
             ggCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCATGGCCAGGCCGCa                                                                                                                                         <  1:193685/150‑1 (MQ=255)
              gCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCATGGCCAGGCCGCAg                                                                                                                                        <  1:1898511/150‑1 (MQ=255)
              gCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAg                                                                                                                                        >  1:4589862/1‑150 (MQ=255)
              gCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGACTACCGTGGTTACCCTGGGCCTGCGGCCATGGCCAGGCCGCAg                                                                                                                                        <  1:8786352/150‑1 (MQ=255)
                     gCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGc                                                                                                                                 <  1:2936466/150‑1 (MQ=255)
                             tGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTAGGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGt                                                                                                                         <  1:3556042/150‑1 (MQ=255)
                              ggTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGtg                                                                                                                        >  1:6399350/1‑150 (MQ=255)
                                 gATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGtggct                                                                                                                     >  1:6801692/1‑150 (MQ=255)
                                  aTTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGtggctg                                                                                                                    >  1:5311060/1‑150 (MQ=255)
                                    ttCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGtggctggc                                                                                                                  <  2:5084999/150‑1 (MQ=255)
                                       cGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGtggctggctgg                                                                                                               >  2:6821961/1‑150 (MQ=255)
                                       cGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGtggctggctgg                                                                                                               >  2:3550801/1‑150 (MQ=255)
                                       cGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGtggctggctgg                                                                                                               >  1:1095849/1‑150 (MQ=255)
                                            aGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTg                                                                                                          <  2:1427905/150‑1 (MQ=255)
                                            aGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTg                                                                                                          <  1:5967050/150‑1 (MQ=255)
                                            aGCACAGCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTg                                                                                                          <  2:2678742/150‑1 (MQ=255)
                                             gCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTgg                                                                                                         >  2:3223383/1‑150 (MQ=255)
                                             gCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTgg                                                                                                         <  1:3946618/150‑1 (MQ=255)
                                              cacaTCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTggg                                                                                                        <  2:1295019/150‑1 (MQ=255)
                                                caTCGGCCTGGTCGCCGGTTTGCTCTAGGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCATGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGt                                                                                                      <  1:1018175/150‑1 (MQ=255)
                                                  tCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGc                                                                                                    <  1:1553511/150‑1 (MQ=255)
                                                  tCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGc                                                                                                    <  2:7875654/150‑1 (MQ=255)
                                                        tGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACg                                                                                                                            <  2:7705995/120‑1 (MQ=255)
                                                        tGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGc                                                                                              <  1:3433739/150‑1 (MQ=255)
                                                          gTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCATGGCCAGGCCGCa                                                                                                                                         >  2:8763410/1‑105 (MQ=255)
                                                          gTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAg                                                                                            >  1:8115350/1‑150 (MQ=255)
                                                           tCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGc                                                                                           >  1:8487379/1‑150 (MQ=255)
                                                           tCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGc                                                                                           >  1:4054694/1‑150 (MQ=255)
                                                           tCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGc                                                                                           >  2:1161371/1‑150 (MQ=255)
                                                            cGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCg                                                                                          <  1:4782439/150‑1 (MQ=255)
                                                              ccGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCggg                                                                                        <  1:2519136/150‑1 (MQ=255)
                                                                  tttGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGgcctg                                                                                    >  2:4443311/1‑150 (MQ=255)
                                                                  tttGCTCTATGGCCTGGTAGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGgcctg                                                                                    >  2:8593317/1‑150 (MQ=255)
                                                                      ctctATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGcctgcctg                                                                                >  2:6328311/1‑150 (MQ=255)
                                                                       tctATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGc                                                                                                    >  1:137554/1‑129 (MQ=255)
                                                                       tctATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGa                                                                               <  1:824902/150‑1 (MQ=255)
                                                                       tctATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGa                                                                               >  1:7454140/1‑150 (MQ=255)
                                                                       tctAGGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGa                                                                               <  1:4304329/150‑1 (MQ=255)
                                                                         tATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATg                                                                             >  2:5929593/1‑150 (MQ=255)
                                                                          aTGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTAGGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGc                                                                            <  1:2410391/150‑1 (MQ=255)
                                                                             gCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTgg                                                                         <  2:4327422/150‑1 (MQ=255)
                                                                                  gTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCgtg                                                                    >  1:5588142/1‑150 (MQ=255)
                                                                                   tCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCgtgg                                                                   >  2:64727/1‑150 (MQ=255)
                                                                                   tCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCgtgg                                                                   <  1:5251845/150‑1 (MQ=255)
                                                                                   tCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCgtgg                                                                   <  2:3777191/150‑1 (MQ=255)
                                                                                    cGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTAc                                                                                                                             <  2:8308343/91‑1 (MQ=255)
                                                                                    cGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCgtggt                                                                  <  2:6259312/150‑1 (MQ=255)
                                                                                    cGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCgtggt                                                                  <  2:4656565/150‑1 (MQ=255)
                                                                                     gCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCgtggtg                                                                 >  2:6596202/1‑150 (MQ=255)
                                                                                     gCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCgtggtg                                                                 >  1:4857913/1‑150 (MQ=255)
                                                                                     gCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGAGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGAAGCTGGCCgtggtg                                                                 >  2:7276840/1‑150 (MQ=255)
                                                                                           cTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGctc                                                           <  1:7095330/150‑1 (MQ=255)
                                                                                              gcgcGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTgg                                                        <  1:8622948/150‑1 (MQ=255)
                                                                                                gcgTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCg                                                      >  2:7740320/1‑150 (MQ=255)
                                                                                                   ttGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCt                                                   <  2:4637848/150‑1 (MQ=255)
                                                                                                   ttGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCt                                                   <  2:2497353/150‑1 (MQ=255)
                                                                                                    tGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTAGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTg                                                  <  2:1599015/150‑1 (MQ=255)
                                                                                                      ggggCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTgcg                                                >  1:7757710/1‑150 (MQ=255)
                                                                                                       gggCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTgcgc                                               >  2:1331715/1‑150 (MQ=255)
                                                                                                       gggCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTgcgc                                               >  2:4772228/1‑150 (MQ=255)
                                                                                                          cTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCttt                                            >  1:5421433/1‑150 (MQ=255)
                                                                                                           tATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTc                                           >  2:6440451/1‑150 (MQ=255)
                                                                                                           tATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTc                                           >  2:5186945/1‑150 (MQ=255)
                                                                                                               gCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCa                                       >  2:2172780/1‑150 (MQ=255)
                                                                                                                   gCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCtt                                   >  1:5097225/1‑150 (MQ=255)
                                                                                                                    cAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTg                                  <  1:8386515/150‑1 (MQ=255)
                                                                                                                      aGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGgc                                <  2:5785302/150‑1 (MQ=255)
                                                                                                                      aGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGgc                                <  2:4812194/150‑1 (MQ=255)
                                                                                                                          tACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTGGGCGCtt                            >  1:6625479/1‑150 (MQ=255)
                                                                                                                           aCCGTGGTTACCCTGGGCCTGCGGCGTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCttt                           <  1:3223383/150‑1 (MQ=255)
                                                                                                                            ccGTGGTTACCCTGGGCCTGCGGCCTGGGCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACGTTGGCGCtttt                          >  2:80069/1‑150 (MQ=255)
                                                                                                                              gTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCtttttt                        <  2:712363/150‑1 (MQ=255)
                                                                                                                                ggTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTc                      >  1:6271169/1‑150 (MQ=255)
                                                                                                                                 gTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTcc                     <  1:4408114/150‑1 (MQ=255)
                                                                                                                                      ccTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGc                <  2:2981528/150‑1 (MQ=255)
                                                                                                                                      ccTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGc                >  2:6033039/1‑150 (MQ=255)
                                                                                                                                      ccTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGc                >  2:7189456/1‑150 (MQ=255)
                                                                                                                                           gCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGtt           <  1:2065612/150‑1 (MQ=255)
                                                                                                                                           gCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGtt           <  1:7906771/150‑1 (MQ=255)
                                                                                                                                            gcTGCGGCCATGGCCAGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGTTg          >  2:1562799/2‑150 (MQ=255)
                                                                                                                                            gcTGCGGCCATGGCCAGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGTTg          >  2:3952809/2‑150 (MQ=255)
                                                                                                                                             cTGCGGCCATGGCCAGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGTTGc         >  1:6121369/1‑150 (MQ=255)
                                                                                                                                              tGCGGCCATGGCCAGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGTTGCa        >  2:6679040/1‑150 (MQ=255)
                                                                                                                                              tGCGGCCATGGCCAGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGTTGCa        >  2:4426578/1‑150 (MQ=255)
                                                                                                                                              tGCGGCCATGGCCAGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTACGTTGCTGGTTGCa        >  1:3714829/1‑150 (MQ=255)
                                                                                                                                               gCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGTTTCAc       >  2:4558678/1‑150 (MQ=255)
                                                                                                                                                cGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGTTGCACt      <  1:6921397/150‑1 (MQ=255)
                                                                                                                                                    cTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTTGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGTTGCACTGGTc  <  2:1378318/150‑1 (MQ=255)
                                                                                                                                                    cTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGTTGCACTGGTc  >  2:7470580/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTGCAGGCCACCGGCACGGTGCACCTGTTGGTGATTTCCGGCCAGCACATCGGCCTGGTCGCCGGTTTGCTCTATGGCCTGGTCGCCGGGCTGGCGCGTTGGGGGCTATGGCCGGCAAGGCTACCGTGGTTACCCTGGGCCTGCGGCCTGGCCATGGCCGCAGCGCTGGCCTACGGGTGGCTGGCTGGCGGTGGGGTGCCAGTGCAGCGGGCCTGCCTGATGCTGGCCGTGGTGTTGCTCTGGCGTCTGCGCTTTCGCCACCTTGGCGCTTTTTTTCCGTTGCTGGTTGCACTGGTC  >  NC_002947/2145895‑2146192

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: