Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,526,642 G→T 8.4% G151V (GGC→GTC)  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,6420GT 8.4% 220.4 / 8.9 71G151V (GGC→GTC) PP_4015HflD‑like high frequency lysogenization protein
Reads supporting (aligned to +/- strand):  ref base G (37/28);  new base T (3/3);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

CTGCTGGACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTGAACGAACTG  >  NC_002947/4526497‑4526790
                                                                                                                                                 |                                                                                                                                                    
ctgctgGACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGa                                                                                                                                                   <  3:402960/149‑1 (MQ=255)
 tgctgGACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCgc                                                                                                                                                   >  7:193138/1‑147 (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                                                                                                                           <  7:10109/149‑1 (MQ=255)
                        cTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCcaa                                                                                                                           <  1:221736/149‑1 (MQ=255)
                            cGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGCGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAAt                                                                                                                       <  1:203928/149‑1 (MQ=255)
                              cAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATg                                                                                                                      >  1:346430/1‑148 (MQ=255)
                               aaaTCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGAc                                                                                                                                                  >  4:44851/1‑119 (MQ=255)
                               aaaTCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGAc                                                                                                                                                  <  3:44851/119‑1 (MQ=255)
                                aaTCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCGCCATgca                                                                                                                                             >  8:21519/1‑122 (MQ=255)
                                aaTCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCGCCATgca                                                                                                                                             <  7:21519/123‑2 (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                                                                                                   >  3:287884/1‑149 (MQ=255)
                                                gAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCct                                                                                                   >  4:297414/1‑149 (MQ=255)
                                                gAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCct                                                                                                   <  5:60346/149‑1 (MQ=255)
                                                 aGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCGCCATgcac                                                                                                                                            >  5:24104/1‑105 (MQ=255)
                                                 aGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCGCCATgcac                                                                                                                                            <  6:24104/107‑3 (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                                                                                   <  6:335910/146‑1 (MQ=255)
                                                                   aCGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATc                                                                                   >  5:335910/1‑146 (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                                                                     <  7:168229/149‑1 (MQ=255)
                                                                              aTCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACt                                                                     <  3:297414/149‑1 (MQ=255)
                                                                               tCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACtg                                                                    >  5:203745/1‑149 (MQ=255)
                                                                                cGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGa                                                                       <  3:29938/145‑1 (MQ=255)
                                                                                cGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGa                                                                       >  4:29938/1‑145 (MQ=255)
                                                                                   ttCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTggc                                                                >  5:316958/1‑149 (MQ=255)
                                                                                   ttCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTggc                                                                >  3:351277/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)
                                                                                                      cAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATgg                                                                                                                                                      >  5:184734/1‑44 (MQ=255)
                                                                                                      cAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATgg                                                                                                                                                      <  6:184734/44‑1 (MQ=255)
                                                                                                       aGGACCCCCGGAGCTCCTGGCGCCAGCGCATCCAGGGTCAGCGCACCATGCGCGTCCTGCAGCAGCCCCACAAGGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGt                                            <  8:344913/149‑1 (MQ=255)
                                                                                                        ggACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGtt                                           <  2:346430/149‑1 (MQ=255)
                                                                                                          acacCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTgg                                         <  6:318597/149‑1 (MQ=255)
                                                                                                          acacCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTgg                                         <  6:178531/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)
                                                                                                                                                 |                                                                                                                                                    
CTGCTGGACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTGAACGAACTG  >  NC_002947/4526497‑4526790

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: