Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,1640CG67.2% 83.1 / 63.0 62A39G (GCC→GGC) stfRputative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base C (13/7);  new base G (23/18);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGC  >  NC_000913/1428915‑1429410
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
ttATGAGGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGCAATGCCCAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACAAGCTGGcc                                                                                                                                                                                                                                                       >  2:167829/1‑251 (MQ=255)
  agtatgtcccgggCAGTGGCGACTGGCGCATTGACGGTTTCTGCGTGGGTTTCGGTTATAAGTTCTGATTTGCCTGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCt                                                                                                                                                                                                                                                      <  2:354781/239‑1 (MQ=255)
     aaGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTc                                                                                                                                                                                                                                                     >  2:124966/1‑248 (MQ=255)
         ctccGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGAACCACGGTGGTGGTGAACACGGTGGGCTCAGagact                                                                                                                                                                                                                                              >  2:126145/3‑247 (MQ=255)
           ccGGCAGTGGCGACTGGCGGACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTTAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATcc                                                                                                                                                                                                                                            <  1:292426/251‑1 (MQ=255)
               cAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGAGTAGCCAGGTAACACAGTTTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCGCCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGAt                                                                                                                                                                                                                                        <  1:286365/251‑1 (MQ=255)
                  tGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGa                                                                                                                                                                                                                                      >  1:114762/1‑250 (MQ=255)
                         tGGCGCACTGACGGTTTAATCGTGGGTGCCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGCCAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGGGTACTGAAAGACGGCCCAGGAAAACCGGTACCGAACTGCACAATCCAGCTGAAAGCAAAACGTAAAAACACCACGGTGGGGGTGAACACGCAGGCCTCAGAAAATCCGGATGAAGCAgggg                                                                                                                                                                                                                              >  2:366390/1‑250 (MQ=255)
                             gtaCTGACGGCTTCATCGTTGGCGTCGGATATAACTTCTGATTAGCCAGGTAACACAGCGTTATGACAACACGCCGGAACCGGTAGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCATGAAAGACGGAACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACAAGGTGGGCTCAGAGCATCAGGATCAAGCCGGGCGtt                                                                                                                                                                                                                           >  2:48353/3‑250 (MQ=255)
                                ttGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGGAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGACAGACGGCCCAGGAAAACCGGTACAGAACTGCCCCATTCAGCTGAAAGCAAGACGTAACAGCACCACGGTGGTGGTGAGCACCGGGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACa                                                                                                                                                                                                                        >  2:80997/2‑250 (MQ=255)
                                 tGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGc                                                                                                                                                                                                                      >  1:227407/1‑251 (MQ=255)
                                              gttggggtCGGTTATAAATTATGATTCGCCAGGCAACACAGTGTTATGACAGCCCGCCGGACCCGGTGGGCTTTTTTGTGGGGTGCATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCACAGAGAATCCGGATGAAGCCGGGCGTTACAACATggctgtggagg                                                                                                                                                                                                         >  2:298073/7‑242 (MQ=255)
                                                 ggggtCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTAc                                                                                                                                                                                                       <  1:426892/247‑1 (MQ=255)
                                                  tggtctgTTATAAATTCTGATTAGCCAGGTACCATAGTGTTATGACAGCCCGTCGGACCCGGTGGGCTTTGTTGTGGGGTGAATATGGTAGTAAAGATTTCAGGCGTCCTGAAAGATGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACGCGGTGGGCTCAGAGAAGCCGGTTGAAGCTGGGCGTTACAGCATGGATGTGGAGTACg                                                                                                                                                                                                      <  2:223864/244‑1 (MQ=255)
                                                   ggtCGGTTATAAGTTCTGAGTAGCCAGGTAATACAGGGTGATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTTAAAGCAAAACGTAACAGGACCACGGTGGTGGTGAGCACGCTGGCGTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGt                                                                                                                                                                                                    <  2:191792/250‑1 (MQ=255)
                                                   ggtCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGt                                                                                                                                                                                                    >  1:72875/2‑251 (MQ=255)
                                                             aaattCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtg                                                                                                                                                                                          <  1:350811/248‑3 (MQ=255)
                                                               attCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGTTTTTTTGTGGGGTGAATATGGCAGTATAGATTTCAGGAGGCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGTATGTGGAGTACGGTCAGTACAgtgtc                                                                                                                                                                                        <  2:192808/250‑5 (MQ=255)
                                                                    gATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcc                                                                                                                                                                                    <  1:14988/250‑9 (MQ=255)
                                                                    gATTAGCCAGGTAACACAGTGTTATGACAGCCCACCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAAAGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCt                                                                                                                                                                                   >  1:344460/1‑251 (MQ=255)
                                                                    gATTAGCAAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGAGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcct                                                                                                                                                                                   >  1:324221/1‑242 (MQ=255)
                                                                              gTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttgac                                                                                                                                                                         <  2:9245/251‑20 (MQ=255)
                                                                                 acacAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCACGACGTAACGGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTCCGGTCCGTACAgtgtcatcctgcaggatgacggt                                                                                                                                                                      >  1:356569/1‑229 (MQ=255)
                                                                                 acacAGTGTTATGACAGCCAGCCGGAACCGGTGGGCTTTTTTGCGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGCGAATCCGGATGCAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTAccgtgtcatcctgcaggttgacggt                                                                                                                                                                      >  2:54950/1‑227 (MQ=255)
                                                                                  cacaGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttgacggtt                                                                                                                                                                     >  1:391761/1‑228 (MQ=255)
                                                                                    caGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttgacggtttt                                                                                                                                                                   <  2:314857/251‑26 (MQ=255)
                                                                                                                cggacTTTTTTGTGGGGTGAATATGGCAATAAAGATTTCAGATGTACTGAACGACGGGACAGGAAAACCGGTACAGAACTGCACAATACAGCTGAACGCAAAAAGTAACAGCACCACGGTGGGGGTGAACACGCTGGCCTCAGAAACTCCGGCTGAAGCAGGGCGTTACAGCATGGACGTTGAGTACGGACAGTACAGCGTTATTCGGTTGGTGTAAGGATTCCCGCCGTCACATGCCGGGACCATTCCCg                                                                                                                                       >  2:456643/5‑251 (MQ=255)
                                                                                                                 gggtttttttGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGAGGGCATAGGAAAACCGTGACAGAACTGAACAACCCAGCTGAAAGCAAAACTTAGCAGCTCCGGGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCGGGGTGTTACATCTTGGACTTTGGGTACGGTAAGGGTAGGGTTATTCTGTTGTGGGAAGTGTTTCCGCCGTCACATGCCGGGAGCATTACCgt                                                                                                                                      <  2:251722/247‑1 (MQ=255)
                                                                                                                     ttttttGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTAt                                                                                                                                  >  1:12564/1‑251 (MQ=255)
                                                                                                                      ttgtgggggggTGAATAGGGCAGTAAAAATTTCAGGTGTACTAAAAGACGGCCCAGGAAAACCGGTACAAAACTGCACAATCCAGCTAAAAGCAAAACGTAACGGCACCACGGGGGTGGTGAACACGCGGGCCTCAGAAAATCCGGAGGAAGCCGGGTGTTACAGCATGGACGTTGAGTCCGGTCAGTACAGCGTTATTCTTTTGGTGTAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTa                                                                                                                                   <  2:377348/242‑1 (MQ=255)
                                                                                                                             ggggAGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGAAAAACCGGTACAGACCTGCACACTCCAACTGAAAGCAAAACGTAACAGCACCACGGTGGAGGTGAACACGCTGGCCTAAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGCCGTTGCGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATtc                                                                                                                          >  2:297954/1‑251 (MQ=255)
                                                                                                                                    atatGGCAGTAATGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGAGTATGAAGATTCACAAccg                                                                                                                   >  1:228729/1‑250 (MQ=255)
                                                                                                                                    atatGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAAccc                                                                                                                   >  1:200756/1‑251 (MQ=255)
                                                                                                                                    atatGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCAGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccg                                                                                                                   >  1:216607/1‑250 (MQ=255)
                                                                                                                                      atGGCAGTAAAGATTTCAGGAGTTCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccgg                                                                                                                  <  2:336621/250‑3 (MQ=255)
                                                                                                                                        ggCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccgggg                                                                                                                <  2:107834/250‑5 (MQ=255)
                                                                                                                                         gCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccggggac                                                                                                              <  2:271516/251‑7 (MQ=255)
                                                                                                                                            gTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCt                                                                                                           >  1:50414/1‑251 (MQ=255)
                                                                                                                                                aGATTTCAGGAGTTCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCTGGATTAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGGACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCTCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAAt                                                                                                       <  1:444898/251‑1 (MQ=255)
                                                                                                                                                aGATTTCAGGAGTCCTGAAAGACGGCACAGTAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGGGGAGAACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAAt                                                                                                       >  2:107083/1‑251 (MQ=255)
                                                                                                                                                   tttCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGACCTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGAGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGAt                                                                                                    >  1:194147/1‑251 (MQ=255)
                                                                                                                                                     tCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGAttt                                                                                                  >  1:204163/1‑251 (MQ=255)
                                                                                                                                                             aCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTg                                                                                          >  1:251064/1‑251 (MQ=255)
                                                                                                                                                                gAAAGAAGGCACAGGAAAACCGGTACAGAACAGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCGGGCCACAGAAAATCCGGATGAAGCCGGGCGTTACAGACTGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATACCCGCCGTCACATGCCGGGACCATTCCCGTGTATGAAGATTCTCAACCCGGGACACTGAATGATTTTCTCggagca                                                                                        >  2:356267/1‑246 (MQ=255)
                                                                                                                                                                     aCGGCACACGAAAACCGGTACAGAACTGCACAATCCAGCTGCAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACCCGCTGGCCTCAGAAAATCCGGATAAAGCCGGGCGTTACAGCCTGGACGTTGAGTACGGTCCGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGCTTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCCTGACg                                                                                  >  2:380228/1‑251 (MQ=255)
                                                                                                                                                                      cGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgg                                                                                 >  1:102459/1‑251 (MQ=255)
                                                                                                                                                                       ggCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgga                                                                                <  2:72875/251‑1 (MQ=255)
                                                                                                                                                                       ggCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgg                                                                                 >  1:73145/1‑250 (MQ=255)
                                                                                                                                                                       ggCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACgg                                                                                 >  1:399198/1‑250 (MQ=255)
                                                                                                                                                                                  aCCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGccc                                                                     <  1:90255/251‑1 (MQ=255)
                                                                                                                                                                                      gTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCggcc                                                                 <  1:348194/251‑4 (MQ=255)
                                                                                                                                                                                           gAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACAGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCTGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTGCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAgg                                                            <  1:54486/251‑1 (MQ=255)
                                                                                                                                                                                           gAACTGCACAATCCAGCTGAAAGCAAAACGTCACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAgg                                                            >  1:223011/1‑251 (MQ=255)
                                                                                                                                                                                                     acccAGTTGAAACCAAAACGTACAAGCCCCAGGGTGGTGGTGAACACGCGGGCCTCAGAAAATCCGGATGAGCCCGGGCGTTACAGCATGGACGTTGAGTTCGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCg                                                  <  2:317456/249‑1 (MQ=255)
                                                                                                                                                                                                       ccAGCTGAAAGCAAAACGTAACAGAACCAGGGGGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGTGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGtt                                                <  1:16027/251‑1 (MQ=255)
                                                                                                                                                                                                              aaagccagaCGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCACGGCCGGAGGTGCTGCGTCGTCTTGAACt                                         >  1:86258/9‑251 (MQ=255)
                                                                                                                                                                                                               aagccagaCGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTg                                        >  1:24147/8‑251 (MQ=255)
                                                                                                                                                                                                                          aaCAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAag                              <  1:141190/250‑1 (MQ=255)
                                                                                                                                                                                                                            cAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACTGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAgg                           >  1:48075/1‑251 (MQ=255)
                                                                                                                                                                                                                                accacGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGgc                       >  1:63478/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gtggtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAAcgc              <  2:12564/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ggtgAACACGCTGGACTCAGCAAATCCGGATGAAGCCGGGAGATACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCCTGACGGAGGATGATGCCCGTCCGGAGGCACTGAGCCGCTTTGAACTGAGGGGGGAAGAGGTGGCGCGTAACGCGt            >  2:274331/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      tGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGc  >  1:52594/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
TTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGC  >  NC_000913/1428915‑1429410

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