Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913288,3820AC31.8% 81.7 / 10.6 65T532P (ACC→CCC) yagHputative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base A (24/19);  new base C (1/19);  total (25/40)
Fisher's exact test for biased strand distribution p-value = 7.99e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGCCAGGGCG  >  NC_000913/288135‑288632
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                             
cGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCt                                                                                                                                                                                                                                                           <  1:375823/251‑1 (MQ=255)
cGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCt                                                                                                                                                                                                                                                           <  1:272100/251‑1 (MQ=255)
               ggtcgcGGACCATGCGGAGGGCGTCGTGCGGCGGGTGGACGGGGTTACGTGGGGCTCCCGCTACAGCTACTCGTTTGAGGGGGAGACGGGGCACACCGTGCCGGTGCCGTATGAGGCGTGGAAGCGGGCGGACGACTACACCGGCGGGCGCGGCTTCTTCCCCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGATTTCGACTACTTCACCTACGAGCCGGTCTAAc                                                                                                                                                                                                                                            <  2:27898/246‑1 (MQ=255)
               ggcggcGGAACGTGGGGAGTGCGTCTGGCGGCGGGGGGACGGGGACACGCTGGTCTACCGCTACAGCTACTCGTTTGAGGGGGAGACGTGCCCCACCGTGCCGGTGCCGTTGGAGGCGTGGAAGCTGTCGGACGACTCCACCGGCGGGCCCGGTTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGCCAGCACCGGCGACGGCTGCTCCGCGGCCTTCGACTACTTCCCCTCCGAGCCGGTCTAAc                                                                                                                                                                                                                                            <  2:507615/246‑1 (MQ=255)
               ggTGCCGGATCATGCGGACGGCGTCTGGCGGCGGGTGGACGTGGATCCGCTGGTCTCCCGCTACAGCTACTCGTTGGAGGGCGAGACGTGGCACACCGTGCCGGTGACGTAGGAGGCGGGGAAGCTTTCGGACGACTACACCGGCGGGCGCGGCTTCTTCCCCGGCGCGTTTGTGGGCCTGCACTGCGAGGACAGCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTCCGAGCCGGTCTAAc                                                                                                                                                                                                                                            <  2:240362/251‑1 (MQ=255)
                        aCATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGGGCGGCCTCTTTACCGGCGCGGTTTTTGGGCTGCACCGCGGGGGCACCCGCGGGGCCGGGTGCTACCCGGGCTTCGCATCCTTCACCtcccagacgggcttaaagccccagc                                                                                                                                                                                                                                   >  1:458319/1‑227 (MQ=255)
                         caggcggggggCGTCGGGCGGCGGGGGGACGTGGATACGCTGTTCACCCGCTACAGCCACCCGTTGGATGGCGAGAGGTGGCCCACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGGGGGCGCGGCTTCTTCACCGGCGCGTTTGGGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTGCGAGCCGGTCTAACGGCTCAGGcc                                                                                                                                                                                                                                  <  2:248555/241‑1 (MQ=255)
                                               ggggggACGTGGATACGCGGTTCTACCGCTACAGCTACTCGTTTGATGGCGGGACGTGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCCCCGGCGCGTTTTTGGGCCTGCACGGCGAGGACATCAGCGGCGTCGGTTGCTACGCGGACTTCGACTACTTCTCCTACGAGCCGGTCTAACGGCTCAGGCCTGGTAGCCCAGCGCGCGCgaga                                                                                                                                                                                                           <  1:429233/247‑1 (MQ=255)
                                                    gACGGGGATCCGCGGGTCTCCCGCTACAGCTCCTCGTTTGAGGGCGAGACGGGCCACCCCGGGCCGGTGCCGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGGGGACATCAGCGGCGACGGCTGCTACGCGGATTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGcgcg                                                                                                                                                                                                      <  2:340718/251‑1 (MQ=255)
                                                               gCTGGTCTCCCGTTACAGCTACTCGTTGGATGGCGAGACGTGGCCCCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGGGGGCGCGGCTTCTTCCCCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCctg                                                                                                                                                                                           <  1:61654/251‑1 (MQ=255)
                                                               gCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCctg                                                                                                                                                                                           <  1:227542/251‑1 (MQ=255)
                                                                  gttCTACCGCTACAGCTACCCGTTTGATGGCGAGACTTGGCACACCGTGCCGGTGTCGTAGGAGGCGTGGAAGCTGGCGGACAACTACATCGGCGGGCCCGGCTTCTTCACCGGCGCGTTGGTGTGCCTGCACTGCGAGTGCATCAGCGGCGGGGGCTGCTACGCGCACTGCGCCTTCTTCCCCGACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCctgc                                                                                                                                                                                          <  1:407723/247‑1 (MQ=255)
                                                                            tACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGGGCGGGCTCTTCTCCCGCGCGGTTTTGGGGCTGCCCTGCGAGGGAATCAACGGCGGCGGCCGCTACGCCGGCGTCGGCCACCTCACCTACCAACCGAGCCAACGGCCCAGGCCCGGGTGCCCCCCGCGCGCGAGAGGGCGGGCCCGGccccctgacgaaccccc                                                                                                                                                                              >  1:561154/1‑236 (MQ=255)
                                                                                     tCGTTTGATGGCGGGACGGGGCCCCCCGGGCGGGTGACGTATGAGGCGGGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCGTTTGGGGGCCTGCACAGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGAATACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAgct                                                                                                                                                                                   <  1:307487/237‑1 (MQ=255)
                                                                                      tgttggATGGCGAGACGGGGCACACGGCGCCGGTGACGTATGAGGCGGGGAAGTGGTCGGACTCCTACACGGGCGGTCGCGGCTTCTTCCCCGGCGGGTTTGTGGGCCTGCACTGCGAGTACTCCAGGGGCGATGGCTGCTACGCGGACTTCGACTCCTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCTGAAAtt                                                                                                                                                                      <  1:331082/244‑1 (MQ=255)
                                                                                       gttggatggggaggggaggCACGCCGTGCGGGTGCCGTATGGGGCGGGGAAGCTCTGGGCCGGCAAGTTGGGCGGGCGGGGCTTGTGCTCGGGCGCGTTTGGGGGCCTGCACTGCGGGCAGTTCGGCGGGGACGGCTGTGACGGGGACTTCTACTGCTTCACTTGCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGGGCGAGAGCGCGAGCCCGGCCTGCGGAAGCGGCTCGCGGAAATTAGc                                                                                                                                                                   <  1:575286/234‑1 (MQ=255)
                                                                                                    ccGTGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCGTTGGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTGCACCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGcgtcg                                                                                                                                                      <  1:303792/250‑1 (MQ=255)
                                                                                                       gggCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTcc                                                                                                                                                   <  1:447461/250‑1 (MQ=255)
                                                                                                                 ggCCGGTGACGTATGAGGCGGGGAAGCTGTCGGACGACTGCATCGGCGGGCGGGGCTTCTTCACCGGCGCGTTTGTGGGCGTGCGCGGGGAGGACATCAGCGGCGACGGCTGCTACGCGGATTTCGACTACTGCACCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCGTCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAgg                                                                                                                                         <  1:597131/250‑1 (MQ=255)
                                                                                                                          cGTAGGAGGCGTGGAAGCTGGCGGACGACTACATCGGCGGGCGCGGCTTCTTCCCCGGCGCGTGTGTGGGCCTGCACTGCGAGGACATCAGCGGCGCCGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCg                                                                                                                                <  1:385260/251‑1 (MQ=255)
                                                                                                                                         ggCTGTCGGGCGCCTCCATCGGCGGGCGCGGCTTCTTCCCCGGCGGGTGTGTGGGCCTGCACTGGGAGGACACCAGCGGCGAGGGCGGCTACGCGGACTTGGTCTACTTCACCAACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGggcgg                                                                                                                 <  1:129554/250‑1 (MQ=255)
                                                                                                                                               cGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGTCTTCGACTACTCCACCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTcc                                                                                                                                                   <  1:339152/211‑1 (MQ=255)
                                                                                                                                                ggACGCCTACAGCGGCGGGCGCGGCTTTTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTCCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGAc                                                                                                          <  1:512501/251‑1 (MQ=255)
                                                                                                                                                  acgacTACATCGGCGGGCGCGGCTTCTACACCGGCGCGTTGTTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTCCGCGGACTTCGACTACTTCCCCGACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGAcgc                                                                                                        <  1:26422/251‑1 (MQ=255)
                                                                                                                                                    gacTACACCGGCGGGCGCGGTTTCTCCCCCGGCGCGTTTGTGGGCCTGCCCTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTCCGACTACTCCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCgg                                                                                                      <  1:542743/251‑1 (MQ=255)
                                                                                                                                                     gttaccgcggcggGCCGGGCTTCTTCCCGGGCGCGTTTGTGGGCCTGCTCTGCGGGGGCATCAGCGGGTGCGGCTGGTGCGCGGACTGTGACTCCTTCCCCTACGAGCCGGTCTAACGGCGCAGTCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGTCCTTCTGAAGCTGCTCGCGGAAAGTAGCCGGCTCCGCGTCGTCCACGCGTGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCggg                                                                                                     <  1:623203/244‑1 (MQ=255)
                                                                                                                                                        ccATCGGCGGGCGCGGCTTCTTCCCCGGCGCGTTTGGGGGCCTGCACTGCGGGGACATCAGCGGCGGCGTCGGCTACGCGGACTTCGACTACTTCGCCTGCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGa                                                                                                    <  1:462814/248‑1 (MQ=255)
                                                                                                                                                          aTCGGCGGGCGCGGTTTCTTCACCGGCGCGTTTGGGGGCCTGCACTGCGAGGACATCAGCGGCGCCGGCTGCTACGCGGATTCCGACTCCTCCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCg                                                                                                <  1:1495/251‑1 (MQ=255)
                                                                                                                                                            cggcggGCGCGGCTTCTTCACCGGCGCGTTGGGGGGCCTGCACGGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACAGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTg                                                                                              <  1:529793/251‑1 (MQ=255)
                                                                                                                                                            cggcggGCGCGGCTTCTCCACCGGCGCGTTTGTGGGCCTGCATTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGATTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTg                                                                                              <  1:91471/251‑1 (MQ=255)
                                                                                                                                                              gggggcgCGGCTTCTTCCCCGCCTCGTTTGTGGGCCTCGCCTGCAAGCCCTCCAGCGGCGAGGGCTGCTCCGCGGACTTCAACTCCTTCCCCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTTAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGTATGTCATCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGt                                                                                            <  1:643686/249‑1 (MQ=255)
                                                                                                                                                                gggcgCGCCTTCGTCCCCGGCGGGTTTGTGGTACTGCTCCGCGAGTATAGCGGTGGCGACGGCTCCAACGCGGACTTCGACGACCTCACGCTCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGTGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTc                                                                                          <  1:356568/251‑1 (MQ=255)
                                                                                                                                                                    gcgGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCcac                                                                                       >  1:545416/1‑250 (MQ=255)
                                                                                                                                                                    gcgGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCGCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGACCTTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCCGCGCCGACAGGCTCCGGGCGGGGGTGACGCGGGACCCGTGGTTCcaac                                                                                      >  2:25554/1‑249 (MQ=255)
                                                                                                                                                                            ttcACCGGCGCTTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCTTGGTTCCACACCGGCAc                                                                              <  1:543266/250‑1 (MQ=255)
                                                                                                                                                                                   gcgcGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGCGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACCGGCCCAGTGCGGCGATGACGCGGGTCTCGGGGTTCCCCACCGGGCCCGCCACGc                                                                      >  2:4058/1‑251 (MQ=255)
                                                                                                                                                                                     gcgTTTGGGGGCCTCCACTGCGAGGACATCGGGGGCGCCGGCTGCTACGCGGACTTCGACTCCTCCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGGGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAg                                                                    <  1:68939/251‑1 (MQ=255)
                                                                                                                                                                                       gTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCCCACCGGCCCCGCCACGCcgcg                                                                  >  2:130260/1‑251 (MQ=255)
                                                                                                                                                                                            gtgGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCAGGTCTAACGGCTCAGGCCGGGTAGCCCCGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACg                                                              >  1:645609/1‑251 (MQ=255)
                                                                                                                                                                                                 ccTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGCGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTACCCACCGGCCCCGCCCCGCAGCGCACGCCCTg                                                         >  2:245405/1‑251 (MQ=255)
                                                                                                                                                                                                   tGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGCAGCTGCTCGCGGGAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGGGGTTCCCCACCGGCACCGCCACGCCGCGCACGCCCTGcc                                                       >  2:409574/1‑250 (MQ=255)
                                                                                                                                                                                                     gcccggCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTCCTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTc                                                      <  1:145761/246‑1 (MQ=255)
                                                                                                                                                                                                     gCACGGCGAGGACATCAGCGGCGCCGGCTGCTACGCGGACTTCGACTCCTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCTCCGCCACGCAGCGCACGCCCTGCTc                                                      <  1:162489/251‑1 (MQ=255)
                                                                                                                                                                                                      cccTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAg                                                                                                                                     <  1:277673/170‑1 (MQ=255)
                                                                                                                                                                                                      cACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAg                                                                                                                                     >  2:277673/1‑172 (MQ=255)
                                                                                                                                                                                                           ggAGGACACCAGGGGCGCCGGCTGCTACGCGGACTTCGATTACTTCCCCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTtctc                                                <  1:333319/250‑1 (MQ=255)
                                                                                                                                                                                                             aGGACACCGGCGGCGACGGCGGCTCCGCGGACTTCGACTCCTCCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCtt                                              <  1:95916/251‑1 (MQ=255)
                                                                                                                                                                                                              ggACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTc                                             >  1:415003/1‑251 (MQ=255)
                                                                                                                                                                                                                  accggcggcgACGCCTGCTCCGCGGCCGTCGACTACTCCACCGACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTc                                                                                                                                       <  2:433109/154‑1 (MQ=255)
                                                                                                                                                                                                                  aTCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTc                                                                                                                                       >  1:433109/1‑158 (MQ=255)
                                                                                                                                                                                                                     aGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGACGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGCGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGCAATTCGCCCGCTACGCGTCGTCCACGCGGGAGGTCAGCGTCGCCAGGCTCCGGGCGGCGATGACGCGGGACTCGTGGGTCCACCACGGGCCCGCCACGCCGCGCCAGCCCTGCTCGTTCTCTTGGCTGTcc                                      >  2:352883/1‑251 (MQ=255)
                                                                                                                                                                                                                        ggcgACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAgg                                   >  1:157064/1‑251 (MQ=255)
                                                                                                                                                                                                                         gcgACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAggg                                  >  1:251738/1‑251 (MQ=255)
                                                                                                                                                                                                                           gACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAAAGGCTCAGGCCGGGTAGCCCAGCGCGCGAGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCg                                >  1:41984/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggCTGCTACGCGGACATCGACTAATTCACCTACGAGCCGGTCTAACGCCTCAAGCCGGGGAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGACATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAACGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCCCACCGGCACCGCCACGCCGCGCACGCCCTGCCCGCGCTCTTCGCTGTCCCGGGCgg                               >  2:437449/1‑248 (MQ=255)
                                                                                                                                                                                                                                          gACTTCGAATACTTCACCTCCGAGACGGTCTCACGGCTCAGGCAGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTACTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACAACGGCACCGCCACGCCGCGCAAGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCgg                 >  2:661170/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aCTTCGACTACTTCACCTAAGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGGTCCacac                                                                                    >  2:94476/1‑183 (MQ=255)
                                                                                                                                                                                                                                           aCTTCGACTACTCCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCacac                                                                                    <  1:94476/183‑1 (MQ=255)
                                                                                                                                                                                                                                              tCGACTACTTCCCCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTc              <  1:206657/250‑1 (MQ=255)
                                                                                                                                                                                                                                               cGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTg            >  1:148744/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGCCTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTgc           >  2:154535/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCCCACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTgg           >  2:639255/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    actTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCACAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGccc       >  1:55488/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    actTCACCTACGAGCCGGTCTAACGGATCAGGCAGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGACCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGGGGTTCCACACCGGCTCCGCCACGCAGCGCACGCCCTGCGCGTTCTCTTCGCGGTCCAGGGCGTAGCCTTGCTCGCGGGTCCGCGcc        >  2:269601/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       tCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGGCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGCCTggg    >  2:215580/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         aCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGCCAGGgcg  >  1:118233/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                             
CGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGCCAGGGCG  >  NC_000913/288135‑288632

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: