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.

CGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGCCAGGGCG  >  NC_000913/288135‑288632
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                             
cGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑Ct                                                                                                                                                                                                                                                           <  1:375823/251‑1 (MQ=255)
cGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑Ct                                                                                                                                                                                                                                                           <  1:272100/251‑1 (MQ=255)
               ggtcgcGGACCATGCGGAGGGCGTCGTGCGGCGGGTGGACGGGGTTACGTGGGGCTCCCGCTACAGCTACTCGTTTGAGGGGGAGACGGGGCACACCGTGCCGGTGCCGTATGAGGCGTGGAAGCGGGCGGACGACTACACCGGCGGGCGCGGCTTCTTCCCCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGATTTCGACTACTTCAC‑CTACGAGCCGGTCTAAc                                                                                                                                                                                                                                            <  2:27898/246‑1 (MQ=255)
               ggcggcGGAACGTGGGGAGTGCGTCTGGCGGCGGGGGGACGGGGACACGCTGGTCTACCGCTACAGCTACTCGTTTGAGGGGGAGACGTGCCCCACCGTGCCGGTGCCGTTGGAGGCGTGGAAGCTGTCGGACGACTCCACCGGCGGGCCCGGTTTCTTCACCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGCCAGCACCGGCGACGGCTGCTCCGCGGCCTTCGACTACTTCCC‑CTCCGAGCCGGTCTAAc                                                                                                                                                                                                                                            <  2:507615/246‑1 (MQ=255)
               ggTGCCGGATCATGCGGACGGCGTCTGGCGGCGGGTGGACGTGGATCCGCTGGTCTCCCGCTACAGCTACTCGTTGGAGGGCGAGACGTGGCACACCGTGCCGGTGACGTAGGAGGCGGGGAAGCTTTCGGACGACTACACCGGCGGGCGCGGCTTCTTCCCCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACAGCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTCCGAGCCGGTCTAAc                                                                                                                                                                                                                                            <  2:240362/251‑1 (MQ=255)
                        aCATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGGGCGGCCTCTTTACCGGCGCGGTTTTTGGGC‑T‑GCACCGCGGGGGCACCCGCGGGGCCGGGTGCTACCCGGGCTTCGCATCCTTCAC‑Ctcccagacgggcttaaagccccagc                                                                                                                                                                                                                                   >  1:458319/1‑227 (MQ=255)
                         caggcggggggCGTCGGGCGGCGGGGGGACGTGGATACGCTGTTCACCCGCTACAGCCACCCGTTGGATGGCGAGAGGTGGCCCACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGGGGGCGCGGCTTCTTCACCGGCGCG‑TTTGGGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTGCGAGCCGGTCTAACGGCTCAGGcc                                                                                                                                                                                                                                  <  2:248555/241‑1 (MQ=255)
                                               ggggggACGTGGATACGCGGTTCTACCGCTACAGCTACTCGTTTGATGGCGGGACGTGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCCCCGGCGCG‑TTTTTGGGCCT‑GCACGGCGAGGACATCAGCGGCGTCGGTTGCTACGCGGACTTCGACTACTTCTC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑TGGTAGCCCAGCGCGCGCgaga                                                                                                                                                                                                           <  1:429233/247‑1 (MQ=255)
                                                    gACGGGGATCCGCGGGTCTCCCGCTACAGCTCCTCGTTTGAGGGCGAGACGGGCCACCCCGGGCCGGTGCCGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCG‑TTTGTGGGCCT‑GCACTGCGGGGACATCAGCGGCGACGGCTGCTACGCGGATTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGcgcg                                                                                                                                                                                                      <  2:340718/251‑1 (MQ=255)
                                                               gCTGGTCTCCCGTTACAGCTACTCGTTGGATGGCGAGACGTGGCCCCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGGGGGCGCGGCTTCTTCCCCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCctg                                                                                                                                                                                           <  1:61654/251‑1 (MQ=255)
                                                               gCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCctg                                                                                                                                                                                           <  1:227542/251‑1 (MQ=255)
                                                                  gttCTACCGCTACAGCTACCCGTTTGATGGCGAGACTTGGCACACCGTGCCGGTGTCGTAGGAGGCGTGGAAGCTGGCGGACAACTACATCGGCGGGCCCGGCTTCTTCACCGGCGCG‑TTGGTGTGCCT‑GCACTGCGAGTGCATCAGCGGCGGGGGCTGCTACGCGCACTGCGCCTTCTTCCC‑CGACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCctgc                                                                                                                                                                                          <  1:407723/247‑1 (MQ=255)
                                                                            tACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGGGCGGGCTCTTCTCCCGCGCG‑GTTTTGGGGCT‑GCCCTGCGAGGGAATCAACGGCGGCGGCCGCTACGCCGGCGTCGGCCACCTCAC‑CTACCAACCGAGCCAACGGCCCAGGCCCGGGT‑GCCCCCCGCGCGCGAGAGGGCGGGCCCGGccccctgacgaaccccc                                                                                                                                                                              >  1:561154/1‑236 (MQ=255)
                                                                                     tCGTTTGATGGCGGGACGGGGCCCCCCGGGCGGGTGACGTATGAGGCGGGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCG‑TTTGGGGGCCT‑GCACAGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGAATACTTCCC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAgct                                                                                                                                                                                   <  1:307487/237‑1 (MQ=255)
                                                                                      tgttggATGGCGAGACGGGGCACACGGCGCCGGTGACGTATGAGGCGGGGAAGTGGTCGGACTCCTACACGGGCGGTCGCGGCTTCTTCCCCGGCGGG‑TTTGTGGGCCT‑GCACTGCGAGTACTCCAGGGGCGATGGCTGCTACGCGGACTTCGACTCCTTCCC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCTGAAAtt                                                                                                                                                                      <  1:331082/244‑1 (MQ=255)
                                                                                       gttggatggggaggggaggCACGCCGTGCGGGTGCCGTATGGGGCGGGGAAGCTCTGGGCCGGCAAGTTGGGCGGGCGGGGCTTGTGCTCGGGCGCG‑TTTGGGGGCCT‑GCACTGCGGGCAGTTCGGCGGGGACGGCTGTGACGGGGACTTCTACTGCTTCAC‑TTGCGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGGGCGAGAGCGCGAGCCCGGCCTGCGGAAGCGGCTCGCGGAAATTAGc                                                                                                                                                                   <  1:575286/234‑1 (MQ=255)
                                                                                                    ccGTGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCG‑TTGGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTGCAC‑CTCCGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGcgtcg                                                                                                                                                      <  1:303792/250‑1 (MQ=255)
                                                                                                       gggCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTcc                                                                                                                                                   <  1:447461/250‑1 (MQ=255)
                                                                                                                 ggCCGGTGACGTATGAGGCGGGGAAGCTGTCGGACGACTGCATCGGCGGGCGGGGCTTCTTCACCGGCGCG‑TTTGTGGGCGT‑GCGCGGGGAGGACATCAGCGGCGACGGCTGCTACGCGGATTTCGACTACTGCAC‑CTCCGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCGTCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAgg                                                                                                                                         <  1:597131/250‑1 (MQ=255)
                                                                                                                          cGTAGGAGGCGTGGAAGCTGGCGGACGACTACATCGGCGGGCGCGGCTTCTTCCCCGGCGCG‑TGTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGCCGGCTGCTACGCGGACTTCGACTACTTCCC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCg                                                                                                                                <  1:385260/251‑1 (MQ=255)
                                                                                                                                         ggCTGTCGGGCGCCTCCATCGGCGGGCGCGGCTTCTTCCCCGGCGGG‑TGTGTGGGCCT‑GCACTGGGAGGACACCAGCGGCGAGGGCGGCTACGCGGACTTGGTCTACTTCAC‑CAACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGggcgg                                                                                                                 <  1:129554/250‑1 (MQ=255)
                                                                                                                                               cGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGTCTTCGACTACTCCAC‑CTCCGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTcc                                                                                                                                                   <  1:339152/211‑1 (MQ=255)
                                                                                                                                                ggACGCCTACAGCGGCGGGCGCGGCTTTTTCACCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTCCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGAc                                                                                                          <  1:512501/251‑1 (MQ=255)
                                                                                                                                                  acgacTACATCGGCGGGCGCGGCTTCTACACCGGCGCG‑TTGTTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTCCGCGGACTTCGACTACTTCCC‑CGACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGAcgc                                                                                                        <  1:26422/251‑1 (MQ=255)
                                                                                                                                                    gacTACACCGGCGGGCGCGGTTTCTCCCCCGGCGCG‑TTTGTGGGCCT‑GCCCTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTCCGACTACTCCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCgg                                                                                                      <  1:542743/251‑1 (MQ=255)
                                                                                                                                                     gttaccgcggcggGCCGGGCTTCTTCCCGGGCGCG‑TTTGTGGGCCT‑GCTCTGCGGGGGCATCAGCGGGTGCGGCTGGTGCGCGGACTGTGACTCCTTCCC‑CTACGAGCCGGTCTAACGGCGCAGTCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGTCCTTCTGAAGCTGCTCGCGGAAAGTAGCCGGCTCCGCGTCGTCCACGCGTGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCggg                                                                                                     <  1:623203/244‑1 (MQ=255)
                                                                                                                                                        ccATCGGCGGGCGCGGCTTCTTCCCCGGCGCG‑TTTGGGGGCCT‑GCACTGCGGGGACATCAGCGGCGGCGTCGGCTACGCGGACTTCGACTACTTCGC‑CTGCGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGa                                                                                                    <  1:462814/248‑1 (MQ=255)
                                                                                                                                                          aTCGGCGGGCGCGGTTTCTTCACCGGCGCG‑TTTGGGGGCCT‑GCACTGCGAGGACATCAGCGGCGCCGGCTGCTACGCGGATTCCGACTCCTCCCC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCg                                                                                                <  1:1495/251‑1 (MQ=255)
                                                                                                                                                            cggcggGCGCGGCTTCTTCACCGGCGCG‑TTGGGGGGCCT‑GCACGGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCC‑CTACGAGCCGGTCTAACAGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTg                                                                                              <  1:529793/251‑1 (MQ=255)
                                                                                                                                                            cggcggGCGCGGCTTCTCCACCGGCGCG‑TTTGTGGGCCT‑GCATTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGATTTCGACTACTTCCC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTg                                                                                              <  1:91471/251‑1 (MQ=255)
                                                                                                                                                              gggggcgCGGCTTCTTCCCCGCCTCG‑TTTGTGGGCCTCGC‑CTGCAAGCCCTCCAGCGGCGAGGGCTGCTCCGCGGACTTCAACTCCTTCCC‑CTCCGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTTAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGTATGTCATCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGt                                                                                            <  1:643686/249‑1 (MQ=255)
                                                                                                                                                                gggcgCGCCTTCGTCCCCGGCGGG‑TTTGTGGTACT‑GCTCCGCGAGTATAGCGGTGGCGACGGCTCCAACGCGGACTTCGACGACCTCACGCT‑CGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGTGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTc                                                                                          <  1:356568/251‑1 (MQ=255)
                                                                                                                                                                    gcgGCTTCTTCACCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCcac                                                                                       >  1:545416/1‑250 (MQ=255)
                                                                                                                                                                    gcgGCTTCTTCACCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCGCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGACCTTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCCGCGCCGACAGGCTCCGGGCGGGGGTGACGCGGGACCCGTGGTTCcaac                                                                                      >  2:25554/1‑249 (MQ=255)
                                                                                                                                                                            ttcACCGGCGCT‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCTTGGTTCCAC‑ACCGGCAc                                                                              <  1:543266/250‑1 (MQ=255)
                                                                                                                                                                                   gcgcG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGCGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACCGGCCCAGTGCGGCGATGACGCGGGTCTCGGGGTTCCCC‑ACCGGGCCCGCCACGc                                                                      >  2:4058/1‑251 (MQ=255)
                                                                                                                                                                                     gcg‑TTTGGGGGCCT‑CCACTGCGAGGACATCGGGGGCGCCGGCTGCTACGCGGACTTCGACTCCTCCCC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGGGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAg                                                                    <  1:68939/251‑1 (MQ=255)
                                                                                                                                                                                       g‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCCC‑ACCGGCCCCGCCACGCcgcg                                                                  >  2:130260/1‑251 (MQ=255)
                                                                                                                                                                                            gtgGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCAGGTCTAACGGCTCAGGCC‑GGGTAGCCCCGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACg                                                              >  1:645609/1‑251 (MQ=255)
                                                                                                                                                                                                 ccT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGCGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTACCC‑ACCGGCCCCGCCCCGCAGCGCACGCCCTg                                                         >  2:245405/1‑251 (MQ=255)
                                                                                                                                                                                                   t‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGCAGCTGCTCGCGGGAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGGGGTTCCCC‑ACCGGCACCGCCACGCCGCGCACGCCCTGcc                                                       >  2:409574/1‑250 (MQ=255)
                                                                                                                                                                                                     gcccggCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTCCTTCCC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTc                                                      <  1:145761/246‑1 (MQ=255)
                                                                                                                                                                                                     gCACGGCGAGGACATCAGCGGCGCCGGCTGCTACGCGGACTTCGACTCCTTCCC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCTCCGCCACGCAGCGCACGCCCTGCTc                                                      <  1:162489/251‑1 (MQ=255)
                                                                                                                                                                                                      cccTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAg                                                                                                                                     <  1:277673/170‑1 (MQ=255)
                                                                                                                                                                                                      cACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAg                                                                                                                                     >  2:277673/1‑172 (MQ=255)
                                                                                                                                                                                                           ggAGGACACCAGGGGCGCCGGCTGCTACGCGGACTTCGATTACTTCCC‑CTCCGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTtctc                                                <  1:333319/250‑1 (MQ=255)
                                                                                                                                                                                                             aGGACACCGGCGGCGACGGCGGCTCCGCGGACTTCGACTCCTCCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCtt                                              <  1:95916/251‑1 (MQ=255)
                                                                                                                                                                                                              ggACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTc                                             >  1:415003/1‑251 (MQ=255)
                                                                                                                                                                                                                  accggcggcgACGCCTGCTCCGCGGCCGTCGACTACTCCAC‑CGACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTc                                                                                                                                       <  2:433109/154‑1 (MQ=255)
                                                                                                                                                                                                                  aTCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTc                                                                                                                                       >  1:433109/1‑158 (MQ=255)
                                                                                                                                                                                                                     aGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGACGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGCGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGCAATTCGCCCGCTACGCGTCGTCCACGCGGGAGGTCAGCGTCGCCAGGCTCCGGGCGGCGATGACGCGGGACTCGTGGGTCCACCACGGGC‑CCGCCACGCCGCGCCAGCCCTGCTCGTTCTCTTGGCTGTcc                                      >  2:352883/1‑251 (MQ=255)
                                                                                                                                                                                                                        ggcgACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAgg                                   >  1:157064/1‑251 (MQ=255)
                                                                                                                                                                                                                         gcgACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAggg                                  >  1:251738/1‑251 (MQ=255)
                                                                                                                                                                                                                           gACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAAAGGCTCAGGCC‑GGGTAGCCCAGCGCGCGAGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCg                                >  1:41984/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggCTGCTACGCGGACATCGACTAATTCAC‑CTACGAGCCGGTCTAACGCCTCAAGCC‑GGGGAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGACATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAACGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCCC‑ACCGGCACCGCCACGCCGCGCACGCCCTGCCCGCGCTCTTCGCTGTCCCGGGCgg                               >  2:437449/1‑248 (MQ=255)
                                                                                                                                                                                                                                          gACTTCGAATACTTCAC‑CTCCGAGACGGTCTCACGGCTCAGGCA‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTACTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑AACGGCACCGCCACGCCGCGCAAGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCgg                 >  2:661170/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aCTTCGACTACTTCAC‑CTAAGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGGTCCac‑ac                                                                                    >  2:94476/1‑183 (MQ=255)
                                                                                                                                                                                                                                           aCTTCGACTACTCCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCac‑ac                                                                                    <  1:94476/183‑1 (MQ=255)
                                                                                                                                                                                                                                              tCGACTACTTCCC‑CTCCGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTc              <  1:206657/250‑1 (MQ=255)
                                                                                                                                                                                                                                               cGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTg            >  1:148744/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGCCTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTgc           >  2:154535/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCCC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTgg           >  2:639255/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    actTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCACAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGccc       >  1:55488/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    actTCAC‑CTACGAGCCGGTCTAACGGATCAGGCA‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGACCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGGGGTTCCAC‑ACCGGCTCCGCCACGCAGCGCACGCCCTGCGCGTTCTCTTCGCGGTCCAGGGCGTAGCCTTGCTCGCGGGTCCGCGcc        >  2:269601/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       tCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGGCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGCCTggg    >  2:215580/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         aC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGCCAGGgcg  >  1:118233/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                             
CGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGCGGCTTCTTCACCGGCGCG‑TTTGTGGGCCT‑GCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCAC‑CTACGAGCCGGTCTAACGGCTCAGGCC‑GGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCAC‑ACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGCGGGTCTGCGCCAGGGCG  >  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: ‑