Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,627,507 C→T 6.0% intergenic (+38/+28) PP_1421 → / ← PP_1422 sensor histidine kinase TctE/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,627,5070CT 6.0% 913.0 / 29.0 285intergenic (+38/+28)PP_1421/PP_1422sensor histidine kinase TctE/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (154/113);  new base T (6/11);  total (160/125)
Fisher's exact test for biased strand distribution p-value = 8.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

GTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGA  >  NC_002947/1627358‑1627656
                                                                                                                                                     |                                                                                                                                                     
gTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCcgc                                                                                                                                                       <  2:1865880/150‑1 (MQ=255)
 ttAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCcgcg                                                                                                                                                      >  1:3283990/1‑150 (MQ=255)
 ttAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCcgcg                                                                                                                                                      >  2:4583626/1‑150 (MQ=255)
    gggCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGaaa                                                                                                                                                   <  2:445852/150‑1 (MQ=255)
     ggCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAg                                                                                                                                                  >  1:859701/1‑150 (MQ=255)
      gCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAgg                                                                                                                                                 <  1:2655578/150‑1 (MQ=255)
       cTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAggg                                                                                                                                                <  1:1173379/150‑1 (MQ=255)
       cTGGCGATCGTGGGCGAGATCGGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAggg                                                                                                                                                <  2:5918713/150‑1 (MQ=255)
        tcgacaTCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAgggg                                                                                                                                               >  2:2177540/6‑150 (MQ=255)
        tGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAgggg                                                                                                                                               <  2:2353488/150‑1 (MQ=255)
         ggCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGc                                                                                                                                              <  2:6604890/150‑1 (MQ=255)
           cGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGa                                                                                                                                            <  2:1838776/150‑1 (MQ=38)
            gATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGac                                                                                                                                           >  1:6783137/1‑150 (MQ=255)
            gATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGac                                                                                                                                           >  2:6041481/1‑150 (MQ=255)
              tCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGacac                                                                                                                                         <  1:2856474/150‑1 (MQ=255)
                   ggCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGAGGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCt                                                                                                                                    <  1:279319/150‑1 (MQ=255)
                     cGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGACACAGGcaaa                                                                                                                                  <  2:3741632/150‑4 (MQ=37)
                     cGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTgg                                                                                                                                  <  2:6690722/150‑1 (MQ=255)
                       agaTCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTa                                                                                                                                <  2:6433849/150‑1 (MQ=255)
                            tGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGTTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGacac                                                                                                                                         <  2:5027465/136‑1 (MQ=37)
                            tGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGacac                                                                                                                                         >  1:5027465/1‑136 (MQ=37)
                                gggCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGACACAg                                                                                                                                       >  1:4863250/1‑134 (MQ=37)
                                    gcACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCaa                                                                                                                   <  1:6863844/150‑1 (MQ=255)
                                      aCCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAACGGGGCGACACAGGCTGGTACAGCTCAGCGCAACa                                                                                                                 >  1:5049269/1‑150 (MQ=255)
                                        cTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATg                                                                                                               >  1:5193380/1‑150 (MQ=255)
                                         tGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGc                                                                                                              >  2:1560453/1‑150 (MQ=255)
                                                        tGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcagca                                                                                               >  1:3763321/1‑150 (MQ=255)
                                                        tGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcagca                                                                                               <  2:1408698/150‑1 (MQ=255)
                                                        tGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcagca                                                                                               <  1:2328694/150‑1 (MQ=255)
                                                        tGCACGATGGCGAGAAGGGCGGGGTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCGGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcagca                                                                                               <  2:276860/150‑1 (MQ=255)
                                                             gagggCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCa                                                                                          <  2:6905704/147‑1 (MQ=255)
                                                               tGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATc                                                                                        >  2:265720/1‑150 (MQ=255)
                                                                   gagaAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGAcc                                                                                    >  1:550122/1‑150 (MQ=255)
                                                                   gagaAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGAcc                                                                                    >  2:4639691/1‑150 (MQ=255)
                                                                   gagaAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGAcc                                                                                    >  1:1091761/1‑150 (MQ=255)
                                                                         ggCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGc                                                                              >  2:3001677/1‑150 (MQ=255)
                                                                              gTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGt                                                                         <  1:3360027/150‑1 (MQ=255)
                                                                                tgcgggtgcgggTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGc                                                                       <  2:537619/150‑1 (MQ=255)
                                                                                    ggtgcgggtgAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCaaa                                                                   >  2:4827747/1‑150 (MQ=255)
                                                                                    ggtgagggtgagCTGTATCGCGGATTGATGGTTTGCATGTGTCGCACCTTTCATGGGCTCACCCGAGAAAGGGGCGAAACAGGCTGGTACAGCTCAGCGAAACATTATCCGCGCCTCAAGCAGTTCATCGACATCCAGCTCTGTGCaaaa                                                                   <  2:1110766/150‑1 (MQ=255)
                                                                                      tgcgggtgAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACg                                                                 <  2:178187/150‑1 (MQ=255)
                                                                                          ggtgAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCaa                                                             >  2:249854/1‑150 (MQ=255)
                                                                                            tgAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAAcc                                                           >  2:787354/1‑150 (MQ=255)
                                                                                            tgAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAAcc                                                           >  1:6047757/1‑150 (MQ=255)
                                                                                              aGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAAcccc                                                         <  2:550122/150‑1 (MQ=255)
                                                                                              aGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCAAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAAcccc                                                         >  1:2825204/1‑150 (MQ=255)
                                                                                                cTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAAAATGCTCCGCGCCTCCAGCAGTTCAACGACCTCCAGCTCGGTGCCAAACGGCAACCCCAg                                                       >  1:2553417/1‑150 (MQ=255)
                                                                                                 tttATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGGTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGa                                                      >  2:5356283/1‑150 (MQ=255)
                                                                                                 tttATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGa                                                      >  2:6359110/1‑150 (MQ=255)
                                                                                                  ttATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcagca                                                                                               <  1:1323447/108‑1 (MQ=255)
                                                                                                  ttATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCcagca                                                                                               >  2:1322883/1‑108 (MQ=255)
                                                                                                               aTGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACa                                                                                                                 <  2:7265702/77‑1 (MQ=255)
                                                                                                               aTGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACa                                                                                                                 >  1:3529077/1‑77 (MQ=255)
                                                                                                               aTGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCgg                                        <  1:6638554/150‑1 (MQ=255)
                                                                                                               aTGGTTTGCCTGTGTCGCCCCTTTCAAGGGCTCACCCGCGAAAGGGGCGACAAAGGCTGGTAAAGCTCAGCGCAACAAGCTACGCGCCTCCACCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCACAATGCGGTAAGcccg                                        >  1:1934165/1‑148 (MQ=255)
                                                                                                               aTGGTTTGCCTGGGTCGCCCCTTTCATGGGCTCACCCGCGACAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCgg                                        <  2:3870539/150‑1 (MQ=255)
                                                                                                                   tttGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGgcagc                                    >  2:423987/1‑150 (MQ=255)
                                                                                                                      gCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGgcagcgca                                 >  2:5275005/1‑150 (MQ=255)
                                                                                                                        cTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGgcagcgcagc                               >  2:2127856/1‑150 (MQ=255)
                                                                                                                        cTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGgcagcgcagc                               <  2:1598751/150‑1 (MQ=255)
                                                                                                                          gtgtCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCa                             >  1:2749412/1‑150 (MQ=255)
                                                                                                                          gtgtCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCa                             >  1:7414891/1‑150 (MQ=255)
                                                                                                                            gtCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAAt                           >  1:1786119/1‑150 (MQ=255)
                                                                                                                             tCGCCCCTTTCGCGGGTGAGCCCATGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGcc                                          >  1:3039601/1‑134 (MQ=37)
                                                                                                                             tCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATg                          <  1:6275702/150‑1 (MQ=255)
                                                                                                                             tCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATg                          >  1:7183009/1‑150 (MQ=255)
                                                                                                                             tCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATg                          >  2:4197154/1‑150 (MQ=255)
                                                                                                                               gCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCGCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCCGCCAATGGc                        >  1:149747/1‑150 (MQ=255)
                                                                                                                                ccccTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGcc                       <  1:7327736/150‑1 (MQ=255)
                                                                                                                                 cccTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCg                      >  1:5775470/1‑150 (MQ=255)
                                                                                                                                  ccTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCgg                     >  1:6404313/1‑150 (MQ=255)
                                                                                                                                  ccTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCgg                     <  1:2515932/150‑1 (MQ=255)
                                                                                                                                  ccTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCgg                     >  2:3108407/1‑150 (MQ=255)
                                                                                                                                    tttCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCt                                                                                                             >  1:6940035/1‑60 (MQ=255)
                                                                                                                                    tttCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTg                   >  1:4397822/1‑150 (MQ=255)
                                                                                                                                     ttCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGa                  <  2:5645589/150‑1 (MQ=255)
                                                                                                                                      tCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGAc                 >  2:2846746/1‑150 (MQ=255)
                                                                                                                                        aTGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGc               <  2:4067500/150‑1 (MQ=255)
                                                                                                                                        aTGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGc               <  2:5193380/150‑1 (MQ=255)
                                                                                                                                         tGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGcc              >  1:4629094/1‑150 (MQ=255)
                                                                                                                                          tggtgagcccACGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGccc             >  2:5133215/8‑150 (MQ=255)
                                                                                                                                          gggCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGccc             >  2:3361518/1‑150 (MQ=255)
                                                                                                                                           ggCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCcg            <  2:249046/150‑1 (MQ=255)
                                                                                                                                           ggCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCcg            >  1:4076328/1‑150 (MQ=255)
                                                                                                                                             cTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCcgcg          >  1:614169/1‑150 (MQ=255)
                                                                                                                                              tCACCCGCTATTGGGTCGACACAGGCTGGTACATCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCgcgc         <  1:6396914/150‑1 (MQ=255)
                                                                                                                                              tCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCgcgc         >  1:680767/1‑150 (MQ=255)
                                                                                                                                              tCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCgcgc         >  2:357148/1‑150 (MQ=255)
                                                                                                                                                aCCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGctg       >  1:2089066/1‑150 (MQ=255)
                                                                                                                                                aCCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGctg       >  2:4357483/1‑150 (MQ=255)
                                                                                                                                                aCCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGctg       <  1:1862215/150‑1 (MQ=255)
                                                                                                                                                aCCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGctg       >  1:1604098/1‑150 (MQ=255)
                                                                                                                                                 cccGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGctga      <  2:7347367/150‑1 (MQ=255)
                                                                                                                                                  ccGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTAc                                             >  1:2142849/1‑110 (MQ=255)
                                                                                                                                                  ccGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGctgac     >  1:5620056/1‑150 (MQ=255)
                                                                                                                                                  ccGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGctgac     >  2:5568905/1‑150 (MQ=255)
                                                                                                                                                  ccGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGctgac     >  1:3740650/1‑150 (MQ=255)
                                                                                                                                                  ccGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGctgac     >  1:2937864/1‑150 (MQ=255)
                                                                                                                                                   cgcgAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGctgact    >  2:6797520/1‑150 (MQ=255)
                                                                                                                                                    gcgAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGctgactg   >  1:4213317/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTAGGGCTGGCGATCGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTGCACGATGGCGAGAAGGGCGGGTTGCGGGTGCGGGTGAGCTTTATCGCGGATTGATGGTTTGCCTGTGTCGCCCCTTTCATGGGCTCACCCGCGAAAGGGGCGACACAGGCTGGTACAGCTCAGCGCAACATGCTCCGCGCCTCCAGCAGTTCATCGACCTCCAGCTCGGTGCCAAACGGCAACCCCAGATGGCGGTACGCCGGCAGCGCAGCCAATGGCCGGTGACGCCCGCGCTGACTGA  >  NC_002947/1627358‑1627656

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: