Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,526,642 G→T 5.1% 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 5.1% 258.0 / 3.5 79G151V (GGC→GTC) PP_4015HflD‑like high frequency lysogenization protein
Reads supporting (aligned to +/- strand):  ref base G (41/34);  new base T (2/2);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

TGCTGGACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTGAACGAA  >  NC_002947/4526498‑4526787
                                                                                                                                                |                                                                                                                                                 
tgctgGACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGa                                                                                                                                                >  4:367360/1‑148 (MQ=255)
     gACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATgcg                                                                                                                                          >  7:213598/1‑149 (MQ=255)
     gACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATgcg                                                                                                                                          >  1:326221/1‑149 (MQ=255)
         ccATCCGCAACCGCCGGCCGCAAATCCAGTCGCAGGCGGAGCATTTCGGCCTGGTTGACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATTGCGACATGCGCTTc                                                                                                                                      <  4:313029/149‑1 (MQ=255)
              ggCAACCGCCTGCCTCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTgc                                                                                                                                  <  2:337265/148‑1 (MQ=255)
               gCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTgcag                                                                                                                                >  5:290975/1‑149 (MQ=255)
                 aaCCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCGCCATgcacc                                                                                                                                        <  3:299784/139‑4 (MQ=255)
                 aaCCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCGCCATgcacc                                                                                                                                        >  4:299784/1‑136 (MQ=255)
                 aaCCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCGCCATgcac                                                                                                                                         <  3:192545/138‑3 (MQ=255)
                 aaCCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGTCGCCATgcac                                                                                                                                         >  4:192545/1‑136 (MQ=255)
                     gCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGcc                                                                                                                           >  2:156592/1‑148 (MQ=255)
                       cTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCcaa                                                                                                                        >  5:98353/1‑149 (MQ=255)
                           cGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGTAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAAt                                                                                                                    <  6:417836/149‑1 (MQ=255)
                            gCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATg                                                                                                                   >  6:329227/1‑149 (MQ=255)
                            gCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATg                                                                                                                   >  8:229897/1‑149 (MQ=255)
                             cAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGc                                                                                                                  >  4:153832/1‑149 (MQ=255)
                               aaTCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCt                                                                                                                <  6:98353/149‑1 (MQ=255)
                                   cAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGaa                                                                                                            <  3:153832/149‑1 (MQ=255)
                                    aGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCg                                                                                                              >  1:282450/1‑146 (MQ=255)
                                    aGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCg                                                                                                              >  2:197782/1‑146 (MQ=255)
                                          aGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCcgc                                                                                                     <  7:375439/149‑1 (MQ=255)
                                          aGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCcgc                                                                                                     >  8:375439/1‑149 (MQ=255)
                                               gaaatgctCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCct                                                                                                >  4:190613/8‑149 (MQ=255)
                                                    tttCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCcaa                                                                                                                        >  5:154580/1‑120 (MQ=255)
                                                    tttCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCcaa                                                                                                                        <  6:154580/120‑1 (MQ=255)
                                                    tttCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCctgctgc                                                                                           >  7:116928/1‑149 (MQ=255)
                                                     ctCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCctgctgc                                                                                           <  6:15008/147‑1 (MQ=255)
                                                       cGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCctgctgctg                                                                                         <  5:329227/148‑1 (MQ=255)
                                                       cGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCctgctgctg                                                                                         <  2:136245/148‑1 (MQ=255)
                                                               ttCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATc                                                                                >  1:269053/1‑149 (MQ=255)
                                                                  aCGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATccgc                                                                             <  2:326221/149‑1 (MQ=255)
                                                                    gAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATccgcgc                                                                           >  6:261031/1‑149 (MQ=255)
                                                                    gAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATccgcgc                                                                           <  6:141144/149‑1 (MQ=255)
                                                                         cGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCcgcgccgcgc                                                                      <  8:213598/149‑1 (MQ=255)
                                                                             aTCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACt                                                                  >  8:243773/1‑149 (MQ=255)
                                                                                gCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACtgt                                                                <  3:190613/148‑1 (MQ=255)
                                                                                gCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACtgt                                                                <  2:282450/148‑1 (MQ=255)
                                                                                  ttCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTggc                                                             >  3:140162/1‑149 (MQ=255)
                                                                                  ttCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTggc                                                             >  1:50079/1‑149 (MQ=255)
                                                                                       gTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAg                                                        >  5:34233/1‑149 (MQ=255)
                                                                                             ccTTGTATCAGGACACACTGAGCACCTTGCGACAGCGCATCCAGTTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCTAATATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCTCGACTGTGGCGGCAGCTgggc                                                  <  2:365085/149‑1 (MQ=255)
                                                                                              cTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTgggcg                                                 <  5:261031/149‑1 (MQ=255)
                                                                                               ttGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTgggcgg                                                >  6:316674/1‑149 (MQ=255)
                                                                                               ttGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTgggcgg                                                >  1:279091/1‑149 (MQ=255)
                                                                                                 gTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTgggcgggc                                              >  8:73222/1‑149 (MQ=255)
                                                                                                  tATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCa                                             >  3:241540/1‑149 (MQ=255)
                                                                                                   aTCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCcgcgccgcgc                                                                      >  2:317233/1‑123 (MQ=255)
                                                                                                   aTCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCcgcgccgcgc                                                                      <  1:317233/123‑1 (MQ=255)
                                                                                                     cAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCg                                          >  4:164762/1‑149 (MQ=255)
                                                                                                      aGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGt                                         <  2:50079/149‑1 (MQ=255)
                                                                                                        gACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTg                                       >  5:161687/1‑149 (MQ=255)
                                                                                                        gACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTg                                       <  7:243773/149‑1 (MQ=255)
                                                                                                         acacCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTgg                                      <  6:99009/149‑1 (MQ=255)
                                                                                                          cacCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTgg                                      <  6:290975/148‑1 (MQ=255)
                                                                                                             ccTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGc                                  >  4:94970/1‑149 (MQ=255)
                                                                                                              cTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTg                                       >  3:316582/1‑143 (MQ=255)
                                                                                                              cTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTg                                       <  4:316582/143‑1 (MQ=255)
                                                                                                               tGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTg                                >  3:105375/1‑149 (MQ=255)
                                                                                                               tGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTg                                >  4:73165/1‑149 (MQ=255)
                                                                                                                 aGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTg                                <  4:140162/147‑1 (MQ=255)
                                                                                                                 aGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTg                                <  2:279091/147‑1 (MQ=255)
                                                                                                                    aCCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGtt                           <  5:316674/149‑1 (MQ=255)
                                                                                                                      cTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGt                                         >  7:130346/1‑133 (MQ=255)
                                                                                                                      cTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGt                                         <  8:130346/133‑1 (MQ=255)
                                                                                                                      cTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTg                                       >  8:67902/1‑135 (MQ=255)
                                                                                                                      cTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTg                                       <  7:67902/135‑1 (MQ=255)
                                                                                                                       ttGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGgt                              >  7:205155/1‑143 (MQ=255)
                                                                                                                       ttGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGgt                              <  8:205155/143‑1 (MQ=255)
                                                                                                                         gcgcCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAg                        <  1:156592/147‑1 (MQ=255)
                                                                                                                            ccAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGgc                   <  6:34233/149‑1 (MQ=255)
                                                                                                                                  cATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCAACGTTGGCAGCTGGTGTTCAGCCGGCGCAAGt             >  5:325819/1‑149 (MQ=255)
                                                                                                                                   aTCCAGGTGCATGGCGACCTGCGCATCCTGCAGCAGCCCAAAAATGCCACGAAGCTCCGCCCCCTGCAGATGGCCGGAATCCGCCCCGCGCGACAGTGGCGGCGGCGGGGCGGGCCCCGTTGGCAGCTGGTGTTCAGGCGGCGCAAGtt            >  8:259679/1‑149 (MQ=255)
                                                                                                                                      cAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGc          <  2:269053/148‑1 (MQ=255)
                                                                                                                                       aGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCt         >  7:290634/1‑148 (MQ=255)
                                                                                                                                        ggTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTg        >  2:165714/1‑148 (MQ=255)
                                                                                                                                           gCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGAGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTgaac     >  8:280814/1‑148 (MQ=255)
                                                                                                                                            cATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGgt                              >  2:382914/1‑122 (MQ=255)
                                                                                                                                            cATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGgt                              <  1:382914/122‑1 (MQ=255)
                                                                                                                                             aTGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTgaacgaa  <  3:164762/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TGCTGGACACCATCGGCAACCGCCTGCCGCAAATCCAGTCGCAGGCCGAGCATTTCGGCCTGGTTCACGAAAACGTCATCGCTTCCAGTGGAGCCTTGTATCAGGACACCCTGAGCACCTTGCGCCAGCGCATCCAGGTGCATGGCGACATGCGCTTCCTGCAGCAGCCCAACAATGCCTCGAAGATCCGCGCCCTGCTGCTGGCCGGCATCCGCGCCGCGCGACTGTGGCGGCAGCTGGGCGGGCACCGTTGGCAGCTGGTGTTCAGCCGGCGCAAGTTGCTGAACGAA  >  NC_002947/4526498‑4526787

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: ‑