Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,147,328 A→G 5.9% Q691R (CAG→CGG)  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,147,3280AG 5.9% 892.3 / 27.6 309Q691R (CAG→CGG) PP_1897DNA internalization‑related competence protein ComEC/Rec2
Reads supporting (aligned to +/- strand):  ref base A (162/108);  new base G (14/3);  total (176/133)
Fisher's exact test for biased strand distribution p-value = 7.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 316 total aligned reads displayed.

GGACACTGAAGTGCCGCGTATCGACTGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAG  >  NC_002947/2147179‑2147477
                                                                                                                                                     |                                                                                                                                                     
ggACACTGAAGTGCCGCGTATCGACTGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCa                                                                                                                                                       >  1:6585870/1‑150 (MQ=255)
 gACACTGAAGTGCCGCGTATCGACTGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAg                                                                                                                                                      <  2:4231780/150‑1 (MQ=255)
       gAAGTGCCGCGTATCGACTGGCTGCAGGCACCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTc                                                                                                                                                <  2:5887181/150‑1 (MQ=255)
        aaGTGCCGCGTATCGACTGGCTGCAGGCGCCGCATCATGCCAGCCGAAGTTCGTCCACCGAGGCGGTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCg                                                                                                                                               >  2:5510531/1‑150 (MQ=255)
              cgcgTATCGACTGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGt                                                                                                                                         <  2:2657293/150‑1 (MQ=255)
                cgTATCGACTGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTAt                                                                                                                                       >  1:8226398/1‑150 (MQ=255)
                cgTATCGACTGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTAt                                                                                                                                       <  1:8554808/150‑1 (MQ=255)
                 gTATCGACTGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATg                                                                                                                                      >  1:4324877/1‑150 (MQ=255)
                     cGACTGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGggcg                                                                                                                                  >  1:8342351/1‑150 (MQ=255)
                      gACTGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGggcgg                                                                                                                                 >  2:8127921/1‑150 (MQ=255)
                       aCTGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGggcggc                                                                                                                                >  1:4572830/1‑150 (MQ=255)
                       aCTGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGggcggc                                                                                                                                >  1:2958754/1‑150 (MQ=255)
                         tGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCAt                                                                                                                              >  2:2985828/1‑150 (MQ=255)
                           gCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTCCCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTCTGGGCGGCAggg                                                                                                                            >  2:7764168/1‑150 (MQ=255)
                           gCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAg                                                                                                                                                      <  2:7565517/124‑1 (MQ=255)
                            cTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATggg                                                                                                                           >  2:1932006/1‑150 (MQ=255)
                             tGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATgggg                                                                                                                          >  1:4829860/1‑150 (MQ=255)
                                 ggCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGggtggt                                                                                                                      <  1:4157950/150‑1 (MQ=255)
                                  gcgcCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGgtggtg                                                                                                                     <  1:7007416/150‑1 (MQ=255)
                                    gccgcATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGAt                                                                                                                   <  1:5383366/150‑1 (MQ=255)
                                    gccgcATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTGCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGAt                                                                                                                   <  2:7984891/150‑1 (MQ=255)
                                          tcatGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCTGCACGTGCGGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGa                                                                                                             <  1:5440668/150‑1 (MQ=255)
                                          tcatGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGa                                                                                                             >  1:6858113/1‑150 (MQ=255)
                                          tcatGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGa                                                                                                             >  1:1816638/1‑150 (MQ=255)
                                           catGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGCGGTGGTGATGCATGAt                                                                                                            <  2:2534668/150‑1 (MQ=255)
                                            atGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATa                                                                                                           >  2:839207/1‑150 (MQ=255)
                                             tGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATAc                                                                                                          >  1:1491729/1‑150 (MQ=255)
                                                cAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATAcggg                                                                                                       <  2:6046554/150‑2 (MQ=255)
                                                  gCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATAcggcgg                                                                                                     <  2:8601105/150‑1 (MQ=255)
                                                   ccGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCgga                                                                                                    <  2:3708195/150‑1 (MQ=255)
                                                   ccGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCgga                                                                                                    <  2:5067047/150‑1 (MQ=255)
                                                     gCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCggagg                                                                                                  >  1:5170985/1‑150 (MQ=255)
                                                     gCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCggagg                                                                                                  >  1:2583931/1‑150 (MQ=255)
                                                            gTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGgaggag                                                                                                >  2:1492889/1‑145 (MQ=255)
                                                             tCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCa                                                                                          >  1:7465586/1‑150 (MQ=255)
                                                                 ccGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGc                                                                                      <  1:3033454/150‑1 (MQ=255)
                                                                     ggCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGtt                                                                                  <  2:2319096/150‑1 (MQ=255)
                                                                     ggCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGtt                                                                                  <  1:1790263/150‑1 (MQ=255)
                                                                           cGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCt                                                                            <  2:4331520/150‑1 (MQ=255)
                                                                            gTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGTTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGTGGAGGGGGCATTGCGGTTGGTGCTg                                                                           <  1:3676315/150‑1 (MQ=255)
                                                                              ccGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTggg                                                                         >  2:7425530/1‑150 (MQ=255)
                                                                               cGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTgggg                                                                        >  2:7176574/1‑150 (MQ=255)
                                                                               cGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTgggg                                                                        >  1:4232167/1‑150 (MQ=255)
                                                                                gggCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGggc                                                                       <  1:8747641/150‑1 (MQ=255)
                                                                                  gCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGggcgg                                                                     >  2:4915542/1‑150 (MQ=255)
                                                                                    tACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCa                                                                   >  1:5055670/1‑150 (MQ=255)
                                                                                         caccgcGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATgggg                                                              >  1:979352/1‑150 (MQ=255)
                                                                                           ccgcGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTa                                                            <  1:2404174/150‑1 (MQ=255)
                                                                                               gcgGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGAAGCGTGATACGGCGCAGGGGGGGTCCTTGCGGGTGGTGCTGGGGAGGCGAGGTGTAGTCg                                                        >  2:5324731/1‑150 (MQ=255)
                                                                                                 ggggTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAg                                                      <  1:676194/150‑1 (MQ=255)
                                                                                                       ttGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCAGTGCGGTTGGTGCTGGGGCGGCATGGGGTCGTCGAGGGGGtt                                                >  2:6252268/1‑150 (MQ=255)
                                                                                                        tGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGTTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGTTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTa                                               <  2:3929946/150‑1 (MQ=255)
                                                                                                        tGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTCGTCGAggggggtc                                               >  2:7362542/1‑147 (MQ=255)
                                                                                                         gATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGTGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAg                                              <  2:525593/150‑1 (MQ=255)
                                                                                                         gATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAg                                              <  2:3171031/150‑1 (MQ=255)
                                                                                                          ttttCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCTGGTGGTCGAGCGGTATGGGCGGCTTGGGGTGGTGATGCTTGATACGGGGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAgg                                             <  1:4240746/149‑1 (MQ=255)
                                                                                                          aTTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAgg                                             <  1:1266921/150‑1 (MQ=255)
                                                                                                            ttCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAgggg                                           >  1:3857319/1‑150 (MQ=255)
                                                                                                               gcgGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCa                                        >  1:1161275/1‑150 (MQ=255)
                                                                                                               gcgGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCa                                        <  2:1814660/150‑1 (MQ=255)
                                                                                                                   gggCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAgcgg                                    >  1:1639813/1‑150 (MQ=255)
                                                                                                                     gcgcAGCAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATAATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAgcggcg                                  >  2:4046907/1‑150 (MQ=255)
                                                                                                                       gcAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGTATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGt                                >  1:3176169/1‑150 (MQ=255)
                                                                                                                       gcAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGt                                >  1:6084405/1‑150 (MQ=255)
                                                                                                                       gcAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGt                                >  2:7191755/1‑150 (MQ=255)
                                                                                                                       gcAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGAGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGt                                >  1:7112791/1‑150 (MQ=255)
                                                                                                                       gcAAAAACTGCTTGGGGCACCCGCACGTGCTGGTGGTCGAGGGGTTTGGGCGGCTTGGGGTGGTGATGCATGATCCGGGGGTGGTGGGGGCTTTGCGGTTGGTGGTGGGGCGGCATGGGGTAGGCGTGGGGGTTAGGGGGCAGCGGCGGt                                <  1:8065640/150‑1 (MQ=255)
                                                                                                                        caacaaCAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGtt                               >  1:3934947/1‑150 (MQ=255)
                                                                                                                         aacaacAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTc                              >  1:6396645/1‑150 (MQ=255)
                                                                                                                          acaacaGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCt                             <  2:1360814/150‑1 (MQ=255)
                                                                                                                              ccgCTGCGGGCACCCGCACGTGCTGGGGGTCGCGCGGTAGGGGCGGCTTGGGGTGGGGATGCAGGATCCGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTCGTCGCGGGGGTTAGGGGGCAGCGGCGGTTCTGGCg                         <  2:6361194/148‑1 (MQ=255)
                                                                                                                               cgCTTGGGGCACCGGCAGGTGCTGGGGGTGTTGGGGTTTGGGGGGATTGGGGTGGTGTTGATTGTTCGGGGGGTGGCGGGGGCTTTGGGTTTGGTGGTGGGGGGGATTGGGGTAGTCGCGGGGGTTTGGGGGCAGCGGCGGTTCTGGCgg                        <  1:6529926/149‑1 (MQ=255)
                                                                                                                                 cTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCCTGGGGTCGTCGAGGGGGGTAGGGGGCAGCGGCGGTTCTGGCGGGt                      >  2:6324582/1‑150 (MQ=255)
                                                                                                                                 cTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGt                      >  2:14588/1‑150 (MQ=255)
                                                                                                                                  ttCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGtt                     >  2:1988081/1‑150 (MQ=255)
                                                                                                                                  ttCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGtt                     >  1:4053855/1‑150 (MQ=255)
                                                                                                                                  ttCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGtt                     >  2:6452275/1‑150 (MQ=255)
                                                                                                                                     gggCACCCGCACGTGCAGGTGGTCGCGGGGTATGGGCGGCATGGGGTGGTGATGCATGATCCGGGGGAGGCGGGGGCTTGGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGCGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGa                  <  2:162554/150‑1 (MQ=255)
                                                                                                                                     gggCACCCGCACGTGCAGGTGGTCGAGCTGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGTAGCGGCGGTTCTGGCGGGTTCGa                  >  1:5655585/1‑150 (MQ=255)
                                                                                                                                     gggCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGa                  <  2:6554532/150‑1 (MQ=255)
                                                                                                                                     gaccACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTCGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGa                  >  2:6064833/4‑150 (MQ=255)
                                                                                                                                       gCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCTTTGGGGTTGGTGCTGGGGGGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAgg                <  2:5009627/150‑1 (MQ=255)
                                                                                                                                       gCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAgg                                                     <  2:4599983/113‑1 (MQ=255)
                                                                                                                                       gCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAgg                                                     >  1:4599983/1‑113 (MQ=255)
                                                                                                                                       gCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGAAGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAgg                >  2:2134410/1‑150 (MQ=255)
                                                                                                                                        cACCCGCAGGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAggg               <  2:6585870/150‑1 (MQ=255)
                                                                                                                                        cACCCGCACGTGCCGGTGGTCGCGCGGTAGGGGCGGCCTGGGGTGGTGATGCATGATCCGGGGGCGGCGGGGGCTTTGCGGTTGGGGCTGGGGCGGCTTGGGGTAGTCGAGGGGGTTCGGGGGCAGCGGCGGTTCTGGCGGGTTCGAggg               <  2:1163037/150‑1 (MQ=255)
                                                                                                                                          cctgCACGTGCGGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGa             >  2:7301119/4‑150 (MQ=255)
                                                                                                                                          cccGCTGGTGCTGGTGGTCGAGCGGTTTGGGGGGCATGGGGTGGTGATGCATGATACGGCGGAGGCGGGGGCATTGCGGTTGGTGGTGGGGCGGCTTGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGa             <  1:7308266/150‑1 (MQ=255)
                                                                                                                                          cccGCAGGTGCCGGTGGTCGCGCGGTATGGGCGGCCTGGGGGGGTGAGGCTTGATCCGGCGGCGGCGGGGGCATTGGGGTTGGTGGGGGGGCGGCTTGGGGTAGTCGTGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGa             <  2:4161877/150‑1 (MQ=255)
                                                                                                                                          cccGCAGGTGCAGGTGGTCGAGCGGTTTGGGCGGCATGGGGTGGTGAGGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCTTGGGTTAGTCGTGGGGGTTTGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGa             <  2:2592179/150‑1 (MQ=255)
                                                                                                                                          cccGCACGTGCTGGTGGTCGAGCGGTTGGGGCGGCATGGGGTGGTGATGCATGAGACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGa             <  2:4760451/150‑1 (MQ=255)
                                                                                                                                          cccGCACGTGCTGGTGGTCGAGCGGTAGGGGCGGCATGGGGTGGTGATGCATGATACGGCGGTGGCGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGa             <  2:7035538/150‑1 (MQ=255)
                                                                                                                                          cccGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGa             >  1:6042562/1‑150 (MQ=255)
                                                                                                                                           ctgCACGTGCGGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAg            >  1:572332/3‑150 (MQ=255)
                                                                                                                                           ccGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAg            >  2:6741795/1‑150 (MQ=255)
                                                                                                                                            cGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGt           >  2:2731209/1‑150 (MQ=255)
                                                                                                                                                 gTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTgaggc      >  1:2008119/1‑150 (MQ=255)
                                                                                                                                                   gCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTgaggcag    >  1:8138218/1‑150 (MQ=255)
                                                                                                                                                     aGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTgaggcagag  >  1:7677119/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGACACTGAAGTGCCGCGTATCGACTGGCTGCAGGCGCCGCATCATGGCAGCCGCAGTTCGTCCACCGAGGCGTTCGTCCGGGCTACCGCACCGCGCGGGGTATTGATTTCGCGGGGGCGCAACAACAGCTTCGGGCACCCGCACGTGCAGGTGGTCGAGCGGTATGGGCGGCATGGGGTGGTGATGCATGATACGGCGGAGGAGGGGGCATTGCGGTTGGTGCTGGGGCGGCATGGGGTAGTCGAGGGGGTTAGGGGGCAGCGGCGGTTCTGGCGGGTTCGAGGGGAGTGAGGCAGAG  >  NC_002947/2147179‑2147477

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: