Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913288,3820AC34.5% 67.6 / 12.0 56T532P (ACC→CCC) yagHputative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base A (15/21);  new base C (2/17);  total (17/39)
Fisher's exact test for biased strand distribution p-value = 2.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGAGCAGCCCATTCCGGTGCCGGAAC‑ATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGCC  >  NC_000913/288134‑288626
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
aCGAGCAGCCCATTCCGGTGCCGGAAC‑ATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATCCGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCAcc                                                                                                                                                                                                                                                    <  2:414030/251‑1 (MQ=255)
aCGAGCAGCCCAGTCCGGTGCCGGAAC‑ATGCGGGGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTGTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACGACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCAcc                                                                                                                                                                                                                                                    <  1:143709/251‑1 (MQ=255)
  gAGCAGCCCAGTCCGGTGCCGGAAC‑ATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCCCCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTa                                                                                                                                                                                                                                                  <  2:345820/251‑1 (MQ=255)
    gcagcCCATTCCGGTGCCGGAAC‑ATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACg                                                                                                                                                                                                                                                <  1:21584/251‑1 (MQ=255)
    gcagcCCATTCCGGTGCCGGAAC‑ATGCGGAGAGCGTCTGGCGGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACg                                                                                                                                                                                                                                                <  2:233875/251‑1 (MQ=255)
               cgggtgcGGAACCATGCGGAGAGCGTCTGGCTGCGGGTGGAGGTGGATACGGGGGTCTACCGCTACAGCTACTTGTGTGATGGCGAGACGTGGAACACCGTGCCGGGGACGTATGAGGCGTGGAAGCGGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTGGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCCCCTCCGAGCCGGTCTa                                                                                                                                                                                                                                      <  2:163266/245‑1 (MQ=255)
                 gggcgggAAC‑TTGGGGAGAGAGTCTGGCGGGGGGGGGAGGTGGATACGCTGGTCTCCCGCTACAGCGACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGGGTGGAAGCTGTCGGACGACTACAGCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGGGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCCCCTACGAGCCGGTCTAACg                                                                                                                                                                                                                                   <  2:125676/246‑1 (MQ=255)
                      ggAAC‑ATGCGGAGAGCGTCTGGCTGCGGGGGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCCCCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCa                                                                                                                                                                                                                              <  1:328644/251‑1 (MQ=255)
                                           ggcggggggTGGACGGGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGGGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGGCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACGTC‑GACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAgcgcg                                                                                                                                                                                                           <  1:421951/244‑1 (MQ=255)
                                                gggggggccgggtcgaCGCTGGTCTACGGCCACTGCTTCCGGTTGGAGGGCGAGGCGTGGCCCACCGTGCGGGGGGCGTTTGAGGTGGGGAAGCGGTCGGACGATTCCATCGGCGGGCGCGGCTTCTTCCCCGG‑CGGGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGGGGCGAGGGCTGCG‑ACGCGGACTTC‑GACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCgag                                                                                                                                                                                                     <  2:357135/236‑1 (MQ=255)
                                                cGGGTGGACGGGGATACGCTGGTCTACCGCTACAGCTACTCGTTGGATGGCGAGACGTGGCACCCCGTGCGGGTGACGTATGAGGGGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGt                                                                                                                                                                                                                      <  2:355327/234‑1 (MQ=255)
                                                   gTGGTCGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCCCCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAgc                                                                                                                                                                                                  <  1:51563/251‑1 (MQ=255)
                                                                cGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCct                                                                                                                                                                                     <  1:133802/251‑1 (MQ=255)
                                                                   ggtgcTACCGCTACAGCTACTCGTTTGATGGCGGGATGTGGCCCCCGGTGCCGGTGAGGTATGAGGTGTGGAGGCGGTCGGACGACTACTTGGGCGGGGGCGGCTTCGTCACGGG‑CGGGTTGGTGGGCCTGCAC‑TGCGAGAACATCGGCGGCGACGGCTGCT‑ACGCGGACGTC‑GACTACTTCCCCGGCGAGCCGGTCTAACGGCTCAGGCCTGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCctgct                                                                                                                                                                                  <  2:486108/247‑1 (MQ=255)
                                                                    ggTCTTCCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCGCACCTTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGATTCCAGCGGCGGGCGCGGCTTCTTCACGGG‑CGCGTTGGTGGGACGGCAC‑GGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGTCTTC‑GCCTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCctgc                                                                                                                                                                                   <  1:249826/249‑1 (MQ=255)
                                                                         cccGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACTTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGaa                                                                                                                                                                               <  1:253274/247‑1 (MQ=255)
                                                                               aCAGCTCCTCGTTTGAGGGCGAGACGTGGCACACCGGGCCGGTGACGTATGGGGCGTGGAAGCTGTCGGACGCTTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCGGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGTACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTcg                                                                                                                                                                       <  1:464622/250‑1 (MQ=255)
                                                                                cAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTcgcg                                                                                                                                                                     <  1:386216/251‑1 (MQ=255)
                                                                                           ttgagggggataCGTGGCACTCCGGGGGGGTGAGGTAGGAGGCGGGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCCTCTTCACCGG‑CGCGTTGGTGGGCCTGCAC‑TGCGAGGCCATCGGCGGCGACGGCTGCG‑CCGCGGACTTC‑GCCTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGcc                                                                                                                                                           <  2:129418/239‑1 (MQ=255)
                                                                                                ggcgagctgTGGCACTCCGTTGCGGTGACGTATGAGGCGTGGAAGTTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCCCCGG‑CGCGTTGGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGTCGGCTGCT‑ACGCGGCCTTC‑GACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTcc                                                                                                                                                     <  1:225917/243‑1 (MQ=255)
                                                                                                 gcgggaCGTGGCACACCGGGCTGGTGACGTATGAGGCGGGGATGCTGTCTGCCGACTTCATCGGGGGGCGCGGCTTCTTCACGGG‑CGCGTTGTTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGAGGGGTGCT‑ACGCGGACTGC‑GACTCCTTCATCTAGGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCTCGCGATATCGCGAGCCCGGGCGGTTGAAGCTGCTCGCGGAAATTAGCCAGCTCcg                                                                                                                                                    <  1:409674/247‑1 (MQ=255)
                                                                                                                                          aaGCTGTCGGACGACTACAGCGGGGGGCGCGGCTTCTTCCCCGG‑CGCGTTTGGGGGCCTGCAC‑GGCGAGGACATCAGCGGCGACGGCTGCT‑GCGGGGCCTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGggcg                                                                                                           <  2:133233/251‑1 (MQ=255)
                                                                                                                                                  ggtcgacTACCTCGGCGGGCGCGGCTTCTTCACCGT‑CGCGTTGGTGGGCCGGCAC‑GGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTCC‑GTCTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCTGCTCCGCTTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGTCGATGAc                                                                                                   <  1:298118/248‑1 (MQ=255)
                                                                                                                                                  ggccgattACACCGGCGGGCTCGGCTTCTTCCCCGG‑CGCGTTTGTGGGCCGGCAC‑TGCGAGGACATCAGCGGCGACGTCTGCT‑ACGCGGACTTC‑GACTACTGCACCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGAc                                                                                                   <  1:125806/244‑1 (MQ=255)
                                                                                                                                                    ccgacTACTCCGGCGGGCGCGGCTTCTTCCCCGG‑CGCGTTTTTGGGCCTGCAC‑TGCGAGGACACCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GCCTACTTCCCCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGAcgc                                                                                                 <  1:46042/250‑1 (MQ=255)
                                                                                                                                                     cgacTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGAcgcg                                                                                                <  1:362473/251‑1 (MQ=255)
                                                                                                                                                      cgctcccgcggcggGCGCG‑CTTCTGCCCCGGGCGCTTTTGTGGTCC‑GCACCTGGGAGGACCTCGGCGGGGACGGCTGCT‑ACGCGGGCTTC‑GCCTCCTCCCCCCCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCgg                                                                                               <  2:319533/243‑1 (MQ=255)
                                                                                                                                                          aCATCGGCGGGGGCGGCTGTTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑GGCGAGGACAGCGGCGGCGACGGCTGCT‑ACGCGGTCTTC‑GACTTCTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGGCGACAGGCTCAGGGCGGCGATGACGTGGGAc                                                                                            <  1:411840/250‑1 (MQ=255)
                                                                                                                                                              cgggggggggggCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTg                                                                                       <  1:375176/241‑1 (MQ=255)
                                                                                                                                                                    gcgcgtcttttttgccGG‑CGCGTTTGTGGGCCTGCAC‑TGCGCGGCTGGCAGCGGCGCCGGCT‑CTGCCGCGGACGTC‑GACTCCTCCCCCTTCGGGCCGGTCTAACGGCGCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGGTCCGCGTCGTCCACGCGGGAGGTCATCTTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGGGGTTCca                                                                                 <  2:448771/237‑1 (MQ=255)
                                                                                                                                                                     cgcgGCTTCTTCCCCGG‑CGCGTTTGTGGGCCGGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑CCGCGGACTTC‑GACTACTCCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCcac                                                                                <  1:377996/251‑1 (MQ=255)
                                                                                                                                                                          gttccttccccgggccgTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑CCGCGGACTTC‑GCCTACGTCCCCGACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGGGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCg                                                                            <  1:279599/236‑1 (MQ=255)
                                                                                                                                                                          cttcttCACCGG‑AGCGCTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCCCCTACCAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGCAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCCGCGTCGACCGGCTCAGGGCGGCGATGACGCGGGCCTCGTGGTTcccccccgg                                                                           >  2:239667/1‑244 (MQ=255)
                                                                                                                                                                            tcttcACCGG‑CGCGTTTGCGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACCTC‑GACTACTTACCCTAAGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGCGAGCGCGAGCACGGCCTGCTGAAGCTGCTCGCGGAAATTAGCAAGCTCCGAGTCGTCCACGCGGGGGGCCAGAGTCGACAGGCTCCGGGCGGGGAGGACGCGGGGCTCGTGGTGCCAAACCGGcc                                                                         >  2:453255/1‑250 (MQ=255)
                                                                                                                                                                             cttcACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACGGGCTCAGGGCGGCGATGACGCGGGACCCGTGGGTCCACACCGGCAc                                                                        >  2:98307/1‑251 (MQ=255)
                                                                                                                                                                                  cccggcggttgtttggcccTGCAC‑TGCTAGGATGTCGGTGGCGCCGGCTGCT‑ACGCGG‑CTTCCGCCTCCTTCTCCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGGGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGcc                                                                    <  2:121063/234‑1 (MQ=255)
                                                                                                                                                                                    gg‑CGCGCTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCAcc                                                                 >  1:131286/1‑250 (MQ=255)
                                                                                                                                                                                       ggtgTTTGGGGGCCTGCCC‑TGGGAGGACACCAGCGGCGCGGGCTGCT‑ACGCGGCCTTC‑GACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCa                                                               <  1:371787/248‑1 (MQ=255)
                                                                                                                                                                                         cgTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCAGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGACGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCCCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACCCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAgc                                                             >  1:150953/1‑251 (MQ=255)
                                                                                                                                                                                           tttGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAgcgc                                                           >  1:259342/1‑251 (MQ=255)
                                                                                                                                                                                            ttGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACCGCCGCT‑ACGCAGACTTC‑CACTACTTCACCTACGAGCCGGTCTAACGGCTCCGGCCGGGTAGCCCAGCCCGCCAGAGAGCGCGAGCCCGGCCTGCTGAAGCTTCTCGCGGAAATTAGCCAGCTCCGCGTCGCCCCCGCGGGAGGTCATCGTCGACCGTCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACACAgcgc                                                           >  1:68415/1‑250 (MQ=255)
                                                                                                                                                                                              gtgGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGAATTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACg                                                        >  1:159851/1‑251 (MQ=255)
                                                                                                                                                                                                  gCCTGCAC‑GGCGAGAACATCGGCGGCGCCGGCGGCT‑ACGCGGCCTTC‑GACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGTTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCt                                                    <  1:284222/251‑1 (MQ=255)
                                                                                                                                                                                                     ggCAC‑TGCAAGGACATCGGCGCCGACGGCTGCT‑CCGCGGCCTCC‑GACTCCTCCCCCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTg                                                   <  1:99080/248‑1 (MQ=255)
                                                                                                                                                                                                        aC‑TGCGAGGACACCAGCGGCGACGGCTGCT‑CCGCGGACTTC‑GACTACTCCACCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGt                                              <  1:462159/251‑1 (MQ=255)
                                                                                                                                                                                                                  aCATCAGAGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGc                                     >  1:355359/1‑251 (MQ=255)
                                                                                                                                                                                                                    aTCAGCGGCGACGGCTGCT‑ACGCGGCCTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGAGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCCGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAAACCGGCACCGCAACGCCGAGCAGGCAATGCTCGTTCTCTACGCTg                                   >  2:471962/1‑251 (MQ=255)
                                                                                                                                                                                                                             gACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCg                          >  1:357608/1‑251 (MQ=255)
                                                                                                                                                                                                                             gACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGAGCg                          >  1:116426/1‑251 (MQ=255)
                                                                                                                                                                                                                              cggtctgct‑CCGCGGATTCC‑GCCTACTCCCCCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCATGGCGt                         <  1:145843/247‑1 (MQ=255)
                                                                                                                                                                                                                                        ccgcgGTCTTC‑GTCTTCTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCtt                   <  1:60478/247‑1 (MQ=255)
                                                                                                                                                                                                                                                  c‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCCGGGCGTAGCCTTGCTCGGGGGTCTg      >  2:386889/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    gACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCCCGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGCAGCCTTGCTCGCGGGTCTgg     >  2:343721/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    gACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTgc     >  2:92855/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      ctactTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTgcgc   >  2:221551/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       tactTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGGTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGcc  >  2:302031/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
ACGAGCAGCCCATTCCGGTGCCGGAAC‑ATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGG‑CGCGTTTGTGGGCCTGCAC‑TGCGAGGACATCAGCGGCGACGGCTGCT‑ACGCGGACTTC‑GACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGCC  >  NC_000913/288134‑288626

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