Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,717,084 T→C 69.3% F116S (TTC→TCC)  PP_4173 → two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,717,0840TC69.3% 320.9 / 198.6 238F116S (TTC→TCC) PP_4173two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base T (45/28);  new base C (78/87);  total (123/115)
Fisher's exact test for biased strand distribution p-value = 4.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 238 total aligned reads displayed.

CAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGCCTTCAGTCAGCGCCAGCAG  >  NC_002947/4716938‑4717229
                                                                                                                                                  |                                                                                                                                                 
caacGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCa                                                                                                                                                >  2:1859126/1‑150 (MQ=255)
 aacGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAg                                                                                                                                               >  1:7758638/1‑150 (MQ=255)
 aacGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAg                                                                                                                                               >  2:3159792/1‑150 (MQ=255)
  acGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAgg                                                                                                                                              >  1:6824640/1‑150 (MQ=255)
       ctGCTTGAAGGGCTGATGCGTGCCACCCACAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGccc                                                                                                                                         <  2:1689560/149‑1 (MQ=255)
          cTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCAtt                                                                                                                                      >  2:6137191/1‑150 (MQ=255)
           ttGAAGGGCTGCTGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACCGGTCGAGGTGTCCCAGGCGCCCATTc                                                                                                                                     <  2:531132/150‑1 (MQ=255)
           ttGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTc                                                                                                                                     <  1:2733199/150‑1 (MQ=255)
            tGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGAAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCg                                                                                                                                    >  2:3794184/1‑150 (MQ=255)
             gAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGc                                                                                                                                   >  1:360612/1‑150 (MQ=255)
               aGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCt                                                                                                                                 <  1:6066782/150‑1 (MQ=255)
                gggCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTg                                                                                                                                <  2:7012699/150‑1 (MQ=255)
                  gCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTgca                                                                                                                              <  2:7696380/150‑1 (MQ=255)
                    tGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTgcagc                                                                                                                            >  1:7992452/1‑150 (MQ=255)
                     gATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTgcagca                                                                                                                           >  1:2852650/1‑150 (MQ=255)
                      aTGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCaa                                                                                                                          <  2:2711306/150‑1 (MQ=255)
                       tGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCaaa                                                                                                                         <  1:6319764/150‑1 (MQ=255)
                       tGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCaaa                                                                                                                         <  2:1030536/150‑1 (MQ=255)
                        gCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAAt                                                                                                                        >  2:3570859/1‑150 (MQ=255)
                          gTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGAAACAGGTCGAGGTGTCCCAGGCGCACATTCGCCTGCAGCAAATcc                                                                                                                      >  1:3193354/1‑150 (MQ=255)
                           tGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCg                                                                                                                     >  1:807425/1‑150 (MQ=255)
                            gCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGc                                                                                                                    <  1:4470529/150‑1 (MQ=255)
                            gCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGc                                                                                                                    <  2:6809492/150‑1 (MQ=255)
                            gCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGc                                                                                                                    <  2:6094867/150‑1 (MQ=255)
                              cACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCt                                                                                                                  <  1:7515932/150‑1 (MQ=255)
                              cACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCt                                                                                                                  <  2:5510799/150‑1 (MQ=255)
                                 ccTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGa                                                                                                               >  2:503479/1‑150 (MQ=255)
                                   tCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACa                                                                                                             >  1:3624761/1‑150 (MQ=255)
                                    cAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACaa                                                                                                            <  2:6414420/150‑1 (MQ=255)
                                     aGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAAc                                                                                                           <  1:2226074/150‑1 (MQ=255)
                                      gCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACg                                                                                                          <  2:8408567/150‑1 (MQ=255)
                                         tCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACt                                                                                                       >  1:3106601/1‑150 (MQ=255)
                                           ccGCACGTGCGCTTTCTGTAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTc                                                                                                     <  1:1033461/150‑1 (MQ=255)
                                           ccGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTc                                                                                                     >  2:6403165/1‑150 (MQ=255)
                                           ccGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTc                                                                                                     >  1:4743322/1‑150 (MQ=255)
                                             gCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCt                                                                                                   >  2:3687904/1‑150 (MQ=255)
                                             gCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCt                                                                                                   <  2:823061/150‑1 (MQ=255)
                                              cACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTg                                                                                                  <  1:1943203/150‑1 (MQ=255)
                                                 gTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCaa                                                                                               <  2:4260874/150‑1 (MQ=255)
                                                 gTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCaa                                                                                               <  2:4416071/150‑1 (MQ=255)
                                                   gcgcTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACa                                                                                             >  1:3436004/1‑150 (MQ=255)
                                                   gcgcTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACa                                                                                             <  1:2221834/150‑1 (MQ=255)
                                                    cgcTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAg                                                                                            <  2:807425/150‑1 (MQ=255)
                                                     gcTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAgg                                                                                           >  2:2008479/1‑150 (MQ=255)
                                                     gcTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAgg                                                                                           >  2:1399573/1‑150 (MQ=255)
                                                          cTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCaaa                                                                                      <  1:832644/150‑1 (MQ=255)
                                                           tGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATACGCCTGCAGAAAATCCGCCTGGACAACGACTTCCTGCAACAGTGCAAAg                                                                                     >  2:4612010/1‑150 (MQ=255)
                                                           tGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAg                                                                                     <  1:1663131/150‑1 (MQ=255)
                                                             gAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGcc                                                                                   >  2:1347763/1‑150 (MQ=255)
                                                               aGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGccac                                                                                 <  2:2343504/150‑1 (MQ=255)
                                                               aGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGccac                                                                                 <  1:6261458/150‑1 (MQ=255)
                                                                gTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGccacc                                                                                >  1:1168393/1‑150 (MQ=255)
                                                                 tGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGccacca                                                                               <  2:3034432/150‑1 (MQ=255)
                                                                      gACAGCCGCAACCACATCCTGGTGTATGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCg                                                                          >  2:8248784/1‑150 (MQ=255)
                                                                      gACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCg                                                                          >  2:6769465/1‑150 (MQ=255)
                                                                      gACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCg                                                                          <  2:6305146/150‑1 (MQ=255)
                                                                      gACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCg                                                                          <  2:2047226/150‑1 (MQ=255)
                                                                           ccGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCAt                                                                     >  2:778278/1‑150 (MQ=255)
                                                                             gCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCg                                                                   >  2:5359747/1‑150 (MQ=255)
                                                                              cAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCgg                                                                  >  1:8276338/1‑150 (MQ=255)
                                                                                 ccACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGGCTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTga                                                               <  1:8162366/150‑1 (MQ=255)
                                                                                  cacaTCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTgac                                                              >  1:7623617/1‑150 (MQ=255)
                                                                                     aTCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTgacgac                                                           <  1:5638253/150‑1 (MQ=255)
                                                                                           gtgtACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTg                                                     <  1:2012693/150‑1 (MQ=255)
                                                                                           gtgtACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTg                                                     <  1:6537299/150‑1 (MQ=255)
                                                                                                   cGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGt                                             <  1:8985863/150‑1 (MQ=255)
                                                                                                   cGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGt                                             <  2:8021367/150‑1 (MQ=255)
                                                                                                      gCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGAt                                          >  1:7128183/1‑150 (MQ=255)
                                                                                                      gCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGAt                                          >  1:5060340/1‑150 (MQ=255)
                                                                                                          tCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGAtcgtc                                      >  2:1231659/1‑150 (MQ=255)
                                                                                                          tCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGAtcgtc                                      >  2:4694840/1‑150 (MQ=255)
                                                                                                          tCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGAtcgtc                                      <  1:5732328/150‑1 (MQ=255)
                                                                                                          tCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCAGTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGAtcgtc                                      <  1:5591172/150‑1 (MQ=255)
                                                                                                            cgacgaAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCgg                                    >  2:4598454/1‑150 (MQ=255)
                                                                                                             gacggAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGGGATCGTCGGc                                   >  1:7884718/1‑150 (MQ=255)
                                                                                                             gacgaAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGc                                   >  1:4376821/1‑150 (MQ=255)
                                                                                                              acgaAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCa                                  >  2:5188716/1‑150 (MQ=255)
                                                                                                              acgaAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCa                                  >  2:2501531/1‑150 (MQ=255)
                                                                                                                  aaGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTc                              >  1:5465389/1‑150 (MQ=255)
                                                                                                                  aaGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTc                              <  1:1859126/150‑1 (MQ=255)
                                                                                                                  aaGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTc                              >  1:6916211/1‑150 (MQ=255)
                                                                                                                  aaGTGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTc                              >  1:5016806/1‑150 (MQ=255)
                                                                                                                     tGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGa                           >  2:3046066/1‑150 (MQ=255)
                                                                                                                     tGCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGa                           >  1:226031/1‑150 (MQ=255)
                                                                                                                      gCCAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGAt                          <  1:848295/150‑1 (MQ=255)
                                                                                                                       ccAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATg                         >  2:7774866/1‑150 (MQ=255)
                                                                                                                       ccAACCGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATg                         >  2:8880115/1‑150 (MQ=255)
                                                                                                                        cAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGa                        >  2:8590854/1‑150 (MQ=255)
                                                                                                                         aaCCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCGGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGAc                       <  2:6340540/150‑1 (MQ=255)
                                                                                                                           ccGCGCGCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGc                     <  1:921369/150‑1 (MQ=255)
                                                                                                                            cgcgcgCAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGcc                    >  2:7338731/1‑150 (MQ=255)
                                                                                                                              cgcgcAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGCCtt                  <  1:4995930/150‑1 (MQ=255)
                                                                                                                                cgcAACAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGCCTtca                <  1:5978198/150‑1 (MQ=255)
                                                                                                                                   aaCAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGCCTtcagtc             >  1:1369044/1‑150 (MQ=255)
                                                                                                                                   aaCAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGCCTtcagtc             >  2:5885790/1‑150 (MQ=255)
                                                                                                                                     cAGGTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGCCTtcagtcag           >  2:6885150/1‑150 (MQ=255)
                                                                                                                                        gTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGCCTTCAGTCAgcgc        >  2:4916933/1‑150 (MQ=255)
                                                                                                                                        gTCGAGGTGTCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGCCTTCAGTCAgcgc        <  1:2040695/150‑1 (MQ=255)
                                                                                                                                              gtgtCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGCCTTCAGTCAGCGCcagcag  >  1:2328099/1‑150 (MQ=255)
                                                                                                                                              gtgtCCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGCCTTCAGTCAGCGCcagcag  >  1:4042130/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGGCGCCCATTCGCCTGCAGCAAATCCGCCTGGACAACGACTTCCTGCAACAGGGCAAAGCCACCAGCCCGGCCATCGGTGACGACTACCTGGGCCGCGTGATCGTCGGCATGTCCGATGACGCCTTCAGTCAGCGCCAGCAG  >  NC_002947/4716938‑4717229

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: