Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,526,645 A→C 7.5% D152A (GAC→GCC)  PP_4015 → HflD‑like high frequency lysogenization protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,526,6450AC 7.5% 217.5 / 3.0 72D152A (GAC→GCC) PP_4015HflD‑like high frequency lysogenization protein
Reads supporting (aligned to +/- strand):  ref base A (39/27);  new base C (2/4);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCTGGACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTGAACGAACTGTAC  >  NC_002947/4526497‑4526793
                                                                                                                                                    |                                                                                                                                                    
ctgctgGACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGa                                                                                                                                                      <  3:402960/149‑1 (MQ=255)
    tgGACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATg                                                                                                                                                  >  6:60346/1‑149 (MQ=255)
      gACACCATCGGCAAACGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCATTGTATAAGCACACCCTGAGCCCCTTGCGCCAGCGCATCCCGGGGCATGGCCACATgcg                                                                                                                                                >  4:65265/1‑149 (MQ=255)
        cacCATCGGCAAGCGCATGCCGCAAATCCAGTCGCAGGCCGATCATTTCGGCCTTGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGTTTCATGTCGCCATgcacc                                                                                                                                              <  8:193138/149‑4 (MQ=255)
        cacCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTTGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCt                                                                                                                                              <  1:124529/149‑1 (MQ=255)
           cATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTcc                                                                                                                                           >  5:295722/1‑149 (MQ=255)
                gCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTGCACGAAAACGTCATCGCTTCCAGTGGAGCCGGGTATCAGGCCACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCAACCTGCGCTGCCTgcag                                                                                                                                      >  7:344913/1‑149 (MQ=255)
                  aaCCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGCGCATGGCGACATGCGCTTCCTgcagca                                                                                                                                    >  2:203928/1‑149 (MQ=255)
                    ccGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTgc                                                                                                                                        <  1:87334/143‑1 (MQ=255)
                    ccGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTgc                                                                                                                                        >  2:87334/1‑143 (MQ=255)
                    ccGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCagcagc                                                                                                                                  <  4:296736/149‑1 (MQ=255)
                        cTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCcaa                                                                                                                              <  1:221736/149‑1 (MQ=255)
                        cTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCcaa                                                                                                                              <  7:10109/149‑1 (MQ=255)
                            cGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGCGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAAt                                                                                                                          <  1:203928/149‑1 (MQ=255)
                              cAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATg                                                                                                                         >  1:346430/1‑148 (MQ=255)
                               aaaTCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGAc                                                                                                                                                     <  3:44851/119‑1 (MQ=255)
                               aaaTCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGAc                                                                                                                                                     >  4:44851/1‑119 (MQ=255)
                                aaTCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCGCCATgca                                                                                                                                                <  7:21519/123‑2 (MQ=255)
                                aaTCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCGCCATgca                                                                                                                                                >  8:21519/1‑122 (MQ=255)
                                aaTCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCt                                                                                                                      <  8:292799/149‑1 (MQ=255)
                                     aGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCg                                                                                                                    >  5:232623/1‑146 (MQ=255)
                                          cAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCcg                                                                                                            >  6:169989/1‑149 (MQ=255)
                                             gCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCgcgc                                                                                                         >  4:195856/1‑149 (MQ=255)
                                               cGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGccc                                                                                                       >  4:159936/1‑149 (MQ=255)
                                               cGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGAGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGccc                                                                                                       >  1:154033/1‑149 (MQ=255)
                                                gAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCct                                                                                                      <  5:60346/149‑1 (MQ=255)
                                                gAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCct                                                                                                      >  3:287884/1‑149 (MQ=255)
                                                gAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCct                                                                                                      >  4:297414/1‑149 (MQ=255)
                                                 aGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCGCCATgcac                                                                                                                                               <  6:24104/107‑3 (MQ=255)
                                                 aGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCGCCATgcac                                                                                                                                               >  5:24104/1‑105 (MQ=255)
                                                           ccTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCg                                                                                           >  7:124525/1‑149 (MQ=255)
                                                               gTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGAGCATCCAGGTGCATGGCTACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCa                                                                                        <  6:295722/148‑1 (MQ=255)
                                                                 tCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATcc                                                                                     >  1:327464/1‑149 (MQ=255)
                                                                   aCGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATccgc                                                                                   >  5:318597/1‑149 (MQ=255)
                                                                   aCGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATc                                                                                      >  5:335910/1‑146 (MQ=255)
                                                                   aCGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATc                                                                                      <  6:335910/146‑1 (MQ=255)
                                                                        aaCGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATccgcgccgc                                                                              >  8:168229/1‑149 (MQ=255)
                                                                           gTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCcgcgccgcgc                                                                            >  6:132165/1‑148 (MQ=255)
                                                                           gTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCcgcgccgcgc                                                                            >  3:120729/1‑148 (MQ=255)
                                                                              aTCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATccgcgccgccacagt                                                                        >  1:18180/1‑143 (MQ=255)
                                                                              aTCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACt                                                                        <  3:297414/149‑1 (MQ=255)
                                                                              aTCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACt                                                                        <  7:168229/149‑1 (MQ=255)
                                                                               tCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACtg                                                                       >  5:203745/1‑149 (MQ=255)
                                                                                cGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGa                                                                          >  4:29938/1‑145 (MQ=255)
                                                                                cGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGa                                                                          <  3:29938/145‑1 (MQ=255)
                                                                                   ttCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTggc                                                                   >  3:351277/1‑149 (MQ=255)
                                                                                   ttCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTggc                                                                   >  5:316958/1‑149 (MQ=255)
                                                                                     ccAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTggc                                                                   <  4:287884/147‑1 (MQ=255)
                                                                                        gTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAg                                                              <  8:124525/149‑1 (MQ=255)
                                                                                           gAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCt                                                            <  3:159936/148‑1 (MQ=255)
                                                                                              ccTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTgggc                                                        >  5:178531/1‑149 (MQ=255)
                                                                                                    aTCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCAc                                                  >  8:237115/1‑149 (MQ=255)
                                                                                                       aGGACCCCCGGAGCTCCTGGCGCCAGCGCATCCAGGGTCAGCGCACCATGCGCGTCCTGCAGCAGCCCCACAAGGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGt                                               <  8:344913/149‑1 (MQ=255)
                                                                                                        ggACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGtt                                              <  2:346430/149‑1 (MQ=255)
                                                                                                          acacCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTgg                                            <  6:178531/149‑1 (MQ=255)
                                                                                                          acacCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTgg                                            <  6:318597/149‑1 (MQ=255)
                                                                                                              ccTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGc                                        >  3:85528/1‑149 (MQ=255)
                                                                                                              ccTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGc                                        >  3:91037/1‑149 (MQ=255)
                                                                                                                  aGCACCTTGCGCCAGCGCATCCAGGTGCATGGCTACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCGGCGCCCTGCTGCTGTCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCTTTGGCAGCTg                                      <  2:154033/147‑1 (MQ=255)
                                                                                                                  aGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTg                                      <  2:327464/147‑1 (MQ=255)
                                                                                                                      cattGCGCCAGCGCATCCAGGTGCATGGGGACATGCGCTTCCTGCAGCAGCCCAACAAAGCCTCGAAGATCCGCGCCCTGATGCTGGCCGGCATCCGCGCAGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGtt                                 >  1:9463/3‑148 (MQ=255)
                                                                                                                       cTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGtt                                              >  4:48970/1‑134 (MQ=255)
                                                                                                                       cTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGtt                                              <  3:48970/134‑1 (MQ=255)
                                                                                                                       cTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGtt                                 >  5:215652/1‑147 (MQ=255)
                                                                                                                          gcgcCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGcc                            >  7:292497/1‑149 (MQ=255)
                                                                                                                            gcCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCgg                          <  6:19592/149‑1 (MQ=255)
                                                                                                                              cAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGgcg                        >  7:19753/1‑149 (MQ=255)
                                                                                                                               aGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGgcgc                       <  4:120729/149‑1 (MQ=255)
                                                                                                                                       cAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGc                <  6:232623/148‑1 (MQ=255)
                                                                                                                                          gTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTgaa            >  3:258299/1‑149 (MQ=255)
                                                                                                                                                ggCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTGAACGAACt      >  3:158158/1‑149 (MQ=255)
                                                                                                                                                 gCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTGAACGAACt      <  6:203745/148‑1 (MQ=255)
                                                                                                                                                 gCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCAGCGCCCTGCTGCTGGCCGGCATCAGCGAAGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTGAACGAACTg     >  7:242534/1‑149 (MQ=255)
                                                                                                                                                    aCATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTGAACGAACTGTAc  <  4:158158/149‑1 (MQ=255)
                                                                                                                                                    aCATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTGAACGAACTGTAc  <  5:169989/149‑1 (MQ=255)
                                                                                                                                                    aCATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTGAACGAACTGTAc  <  8:19753/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGCTGGACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTGAACGAACTGTAC  >  NC_002947/4526497‑4526793

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑