Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,075,8770AC36.1% 75.9 / 12.7 61T109P (ACC→CCC) yihXalpha‑D‑glucose‑1‑phosphate phosphatase YihX
Reads supporting (aligned to +/- strand):  ref base A (18/21);  new base C (1/21);  total (19/42)
Fisher's exact test for biased strand distribution p-value = 5.52e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGCGTATTCCGCTGGCATCGCTTAAGAAGAGTTTTCATATGGGGGAGGCGTTTCATCAGCATGAGCGTGGGGAAATTAGCG‑ACGAAGCGTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGG  >  NC_000913/4075630‑4076126
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          
aaCGCGTATTCCGCTGGCATCGCTTAAGAAGAGTTTTCATATGGGGGAGGCGTTTCATCAGCATGAGCGTGGGGAAATTAGCG‑ACGAAGCGTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATAcc                                                                                                                                                                                                                                                          <  2:616454/250‑1 (MQ=255)
     gTATTCCGCTGGCATCGCTTAAGAAGAGTTTTCATATGGGGGAGGCGTTTCATCAGCATGAGCGTGGGGAAATTAGCG‑ACGAAGCGTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGc                                                                                                                                                                                                                                                   <  2:682280/251‑1 (MQ=255)
       aTTCCACTGGCATAGCTTACGCAGAGTTTTCATCTGGGGGAGGCGTTTCATCAGCATGAGCGTGGGGAACTTCGCG‑CCGAAGCGTTCACAGAGGCGCTGTGTCATGAGATGGCTCTACCGCTAAGATACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCACGCATAAACTGCGTGAGCAGGGGCATCGCGCGGGGGTGC‑TTTCCAATCCC‑ACCCGCCt                                                                                                                                                                                                                                                 >  2:672230/1‑251 (MQ=255)
                    cGCTTAAGAAGGGTTTTCATAGGGGGGAGGGGTTTCATCAGCAGGAG‑GGGGGGAAAATAGGGGACGAGGCTTTGGCGGAGGCGCTGGGTCATGAGGGGGTTCTACCGCTAAGTTACGGGCATTTCTCTCACGGCTGGCGGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACGGCGTGAGCAGGGGCATCGCGTGGGGGTGC‑TGTCCAAGCCA‑AACCGCCTGCATACCACCTTc                                                                                                                                                                                                                                    <  2:574095/251‑1 (MQ=255)
                        taaaaagaGTTTTCATATGGGGGAGGCGTTTCATCAGCATGAGCGTGGGGAAATTAGCG‑ACGAAGCTTTCGCAGGGGCGCTGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAATTG‑ATCGCCATCATGCAAAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTg                                                                                                                                                                                                                                  <  1:301864/245‑1 (MQ=255)
                        tAAGAAGAGTTTTCATATGGGGGAGGCGTTTCATCAGCATGAGCGTGGGGAAATTAGCG‑ACGAAGCTTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTCCCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGGG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGGGC‑TTTCCAATCCC‑AACCGCCTGCATACCACCTTCTgg                                                                                                                                                                                                                                 <  1:632856/250‑1 (MQ=255)
                           gaggagtttgcattgggggggggCGTTTCGTCAGCATGAGGGTGGGGAATTTGGCG‑ACGAAGCGTTCGCGGAGGGGCTGTGTCAGGAGAGGGCTCTCCCGCTAAGCTACGGGCAGTTCTCTCACGGCTGGCAGGCGTTGTT‑TGTGGCGCTGC‑GCCCGGAAGTG‑TTCGCCATCATGCATAAACTGCGTGGGCGGGGGCATCGGGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCgg                                                                                                                                                                                                                             <  2:66975/236‑1 (MQ=255)
                                    tCATATGGGGGAGGCGTTTCATCAGCATGAGCGTGGGGAAATTAGCG‑ACGAAGCGTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTTCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATAcc                                                                                                                                                                                                                    >  1:124185/1‑251 (MQ=255)
                                     cagttGGGGGAGGTGTTTCATCAGCATGGGCGTGGGAAAATTAGCG‑ACGAAGCGTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTCCCACTAAGCTACGAGCAGTTCTCGCCCGGCTGGCAGGCGGTGTT‑TGTTGGGCTGC‑GCCGGGAAGTA‑ATCGCCATCATCCATAAACTGCGTGAGCAGGGGCACCGCGGGGTGGTTC‑TTTCCAATCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATAcc                                                                                                                                                                                                                    <  2:577044/246‑1 (MQ=255)
                                                 gCGTTTCATCAGCATGAGCGTGGGGAAATTAGCG‑ACGAAGCGTCCGCAGAGGCGCTGTTTCACGAGATGGCTCTCCCGCTAAGCTACGAGCATTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCAAATCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGa                                                                                                                                                                                                       <  1:583595/251‑1 (MQ=255)
                                                       aataaGCATGGGCGTGGGGAAATTGGCG‑ACGAAGCTTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTCCCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCGGGCGGTTTT‑TGTGGCGCGGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAACCGGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTCCCAATCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGAtgctgc                                                                                                                                                                                                 <  2:200898/247‑1 (MQ=255)
                                                            gCATGAGCGTGGGGAATTTAGCG‑ACGAAGCGTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTTTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGAcc                                                                                                                                                                                            <  1:136424/251‑1 (MQ=255)
                                                               tGAGCGTGGGGCAATTAGCG‑ACGAAGCGTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGGTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCCTGCATAAACTGCGTGAGCCGGGGCATCGCGGGGGGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCatt                                                                                                                                                                                         >  2:28939/1‑250 (MQ=255)
                                                               gtggggTGGGTAAATTCGCG‑AAGAGGCGTGCGCAGGGGCGCTGGTTCAGGGAATGGCTGTACCCCCAAGCTACGAGATGTGCCCCCACGGCTGTAAGGCTGTGTT‑GGTGGCGCTGC‑GCCCGGAAGGG‑AGCGCCATCATGCATAAACTGGGTGAGGGGTGGCACCGCGTGGTGGGGC‑TTGCCACTCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCat                                                                                                                                                                                          <  2:600460/245‑1 (MQ=255)
                                                                  gCGTGGGGAAATGAGCG‑ACGAAGCGTTCGCAGGGGCGCGGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCATTTCTCTCACGGCTGGCACGCGGTTTT‑TGTGGCGCTGC‑GCCCGGAGGTG‑ATCGCCATCATGCAAAACCTGCGTGGGCGGGGGCAGCGCGTGGTGGTGC‑TTTCCAATCCC‑AACCGCCTGCATACCACCTGCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCAtatct                                                                                                                                                                                      <  1:33487/251‑1 (MQ=255)
                                                                                      gaAGCGTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTggg                                                                                                                                                                   >  1:1898/1‑251 (MQ=255)
                                                                                         gCGTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGGAGTG‑ATCGCCATCATGCATAAACTGCGTGCGCAGGGGCATCGCGTGGTGGGGC‑TTTCCCATACC‑ACACGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCAGCTGACCATATCTATCTGTCGCACGATCGGGGGAt                                                                                                                                                                >  2:447518/1‑251 (MQ=255)
                                                                                          cGTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATg                                                                                                                                                               >  1:548397/1‑251 (MQ=255)
                                                                                             ccGCAGAGGCGTTGTGTCAGGAGATGGCTCTCCCGCTAAGCTACGAGCATTTCTCTCCGGGCTGGCGGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAACCTGCGTGAGCAGGGGCACCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATgcgc                                                                                                                                                            <  1:62685/250‑1 (MQ=255)
                                                                                                      tgcTGTGTCATGAGAGGGCTCTACCGCTAAGCTACGAGCAGTTCTCTAACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑CCCCGGAAGTG‑ATCGCCCTCCTGCATACACTGCGTGAGCCGGGGCATCCCGTGGCGGCGC‑TTTCCAATACC‑AACCGCCTGCTTACCACCTTCTGGCCGGTAGAATTCCCGGGAATTCGGGCTGCTGCTGACCATATCTATCTGTTGCAAGAGCTGGGGATGTGCAAACCTGa                                                                                                                                                    >  2:173609/2‑250 (MQ=255)
                                                                                                       gcTGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGCCGTG‑ATCGCCATCATGCATAAAATGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCACCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGa                                                                                                                                                    >  1:69440/1‑249 (MQ=255)
                                                                                                              caggagagggctcgccctctaggcgaCGAGGAGTTCTCCCACGGCTGGCGGGGGGTTTT‑TGTTGCGCGGC‑GCCCGGAAGTG‑ATCGCCATCATGAAAAAACTGCGTCAGCAGGGGCATCGCGGGGTGGTGC‑TTTCAAATACA‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGGCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAAt                                                                                                                                           <  2:288911/226‑1 (MQ=255)
                                                                                                                    aTGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTAcc                                                                                                                                      >  1:58921/1‑250 (MQ=255)
                                                                                                                     tGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCCCCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCCTTTCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTAcc                                                                                                                                      >  2:664524/1‑249 (MQ=255)
                                                                                                                          ctACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGt                                                                                                                               >  2:255529/1‑251 (MQ=255)
                                                                                                                            cccGCGAAGCTACGCGCGGTCCTCTCCGGGGTGGCAGGCGGTGTT‑TGTGGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCTTGCATAACCGGCGTGAGCAGGGGCATCGCGTGGGGGTGT‑TTCCCAATCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGttt                                                                                                                             <  2:80868/250‑1 (MQ=255)
                                                                                                                            aCCGCTAAGCTACGAGCGGTTCTCACACGGCTGGCAGGCGGTGTT‑TGTGGCGCTGC‑GCCCGAAAGTG‑ATCCCCATCTTGCATAAACGGCGGGAGCAGGGGCATCGCGTGGTGGGGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAATAATACCCGGAAATTCGTGATGCTTCTTACCATATCTATCGGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATg                                                                                                                                <  1:277810/248‑1 (MQ=255)
                                                                                                                            aCCGCCAAGCAACGGGCAGTGCACTCACGGCGGGCGGGCGGTGTT‑GGTTGCG‑TGCCGCCGGGAAGGA‑ATCGCCACCAGGCAAAAACGGGGGGAGCAGGGGCAGCGCGTGGTGGGCC‑TTTCCAAGCCC‑AGCCGCCTGCAGACCACCTTCGGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTTAAGCACGAATTTACCAGCATGtt                                                                                                                              <  2:98513/250‑1 (MQ=255)
                                                                                                                             ccgtgaaGCGACGAGCAGTTCACTTTGGGCTGGCAGGCGGTGTT‑GGTGGCGCGGC‑GCCCGGAAGGC‑ATCGCCAACAGGCAAAAACTGCGAGAGCAGGGGAACCACGTGGGGGTGC‑TGTCAAACACC‑ATCCGCCTGCATACCACCGTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATGTACCAGCATGtttt                                                                                                                            <  1:477493/246‑1 (MQ=255)
                                                                                                                              cGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTg                                                                                                                           >  2:24508/1‑251 (MQ=255)
                                                                                                                                  aaGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAgg                                                                                                                       >  1:182510/1‑251 (MQ=255)
                                                                                                                                       aCGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGGGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATTTCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATTTTTTGCAGGCGGa                                                                                                                   >  1:166240/1‑250 (MQ=255)
                                                                                                                                            cattcctcgcACGGCTGGCGGGGGGGGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ACCGCCTCCATGCATAAGCTGCGTGAGGAGGGGCTGCGCG‑GGGGGGGCGTTGCCAATCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGAGCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGttt                                                                                                             <  2:37470/242‑1 (MQ=255)
                                                                                                                                                         gCTGGCAGGCGGTGTT‑TGTGGCGCTGC‑GCCCGGAAGGG‑ACCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGGGGTGC‑TTTCCAAACCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCATCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGAt                                                                                                <  2:645068/251‑1 (MQ=255)
                                                                                                                                                           gggcggggggtttt‑tgttGCGCTGC‑GCCCGGAAGTG‑ATGGCCATCAGGCAAAACCTGCGGGGGCGGGGGCACCGCGGGGTGGTGC‑TTCCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGAGCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATAc                                                                                              <  2:649398/239‑1 (MQ=255)
                                                                                                                                                              cAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑ACCCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGt                                                                                           >  1:445946/1‑251 (MQ=255)
                                                                                                                                                               ggggggTGTG‑TGTGGCGCGGC‑GCCCGGAAGTG‑TTCGCCAGCATGCAAAAACTGCGTGAGCAGGGGCATCGCGGGGGGGGCC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTc                                                                                          <  2:166240/247‑1 (MQ=255)
                                                                                                                                                                 gCGGTGTTGTGT‑GCGCTGC‑GCCCGGA‑GTGGATCGCCATCATGAAAAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TGTCCAATCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCtt                                                                                        <  2:561007/251‑1 (MQ=255)
                                                                                                                                                                 gCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACGGCGTGAGCAGGGGCATCGCGGGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCtt                                                                                        <  1:179154/251‑1 (MQ=255)
                                                                                                                                                                    ggtgtttggtgggtggc‑gcCCGTAAGTG‑CTCGCCATCATGCATAAACTGCGGGAGCGGGTGCATCGCGGGGTGGTTC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCtttt                                                                                      <  2:498489/236‑1 (MQ=255)
                                                                                                                                                                            ttgCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGAt                                                                              >  1:169248/1‑251 (MQ=255)
                                                                                                                                                                             ggcgcggc‑gcCCGGAAGTG‑ATCGCAATCATGCAAAAACGGCGTGAGCAGGGGCATCGGGTGGGGGTGT‑TTTCCAATCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATa                                                                             <  2:1898/250‑1 (MQ=255)
                                                                                                                                                                                gctgc‑gcCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACg                                                                          <  1:180230/251‑1 (MQ=255)
                                                                                                                                                                                cgtggcgccggaaaGTG‑ATCGCCACAATGAATAAACGGGGTGAGCGGGGGCATCGCGGGGGGGTGC‑TTTCCAATCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAAc                                                                           <  2:547753/239‑1 (MQ=255)
                                                                                                                                                                                       cccGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCCGGGGCATCGCGTGGTGGTGT‑TTCCCAATCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGata                                                                    <  1:28939/251‑1 (MQ=255)
                                                                                                                                                                                          gggaatggaccGCCATCATGCAAAAACGGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTCCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGataat                                                                  <  2:58921/241‑1 (MQ=255)
                                                                                                                                                                                                    cgccaccttgaaaaaactgggggggCAGGGGCTTCGCGTGGGGGGGC‑TTTCCAATACC‑ATCCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGa                                                         <  2:548397/226‑1 (MQ=255)
                                                                                                                                                                                                     gCCATCATGCATAACCTGGGTGGGCGGGGGCTGGGCGTGGTGGTTC‑TTTCAAATCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACTATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATTTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATGCGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGc                                                       <  1:654994/251‑1 (MQ=255)
                                                                                                                                                                                                        atcatGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCa                                                     >  1:221774/1‑250 (MQ=255)
                                                                                                                                                                                                             gcagaaacggggggggCAGGGGCACCGCGTGGGGGGTC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGc                                               <  2:574970/236‑1 (MQ=255)
                                                                                                                                                                                                                 aaacggcgggggCAGGGGCATCGCGGGGGGGTGC‑TTCCCAATCCA‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTGTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTggg                                           <  2:445946/240‑1 (MQ=255)
                                                                                                                                                                                                                      gCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCAtt                                       >  1:4125/1‑250 (MQ=255)
                                                                                                                                                                                                                       gggtggcgggggctccgcgGGGTGGTGC‑TTCCCAAGCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGAGAATATAGAAGGAGCCAATCAGCTGGTCATTAc                                     <  2:317805/236‑1 (MQ=255)
                                                                                                                                                                                                                          gAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATCCC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCGCCCAGCGATACGGTCTTTTTCGACTATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAg                                  <  1:120716/251‑1 (MQ=255)
                                                                                                                                                                                                                             cgggggCACCGCGTGGGGGGGC‑TTTCCAACACA‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGAGCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTAt                               <  2:221774/249‑1 (MQ=255)
                                                                                                                                                                                                                              gggggCACCGCGGGGGGGTGC‑TTTCCAAACCC‑AACCGCCTGCATACCACCTTCTTGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTAtt                              <  2:168048/250‑1 (MQ=255)
                                                                                                                                                                                                                              gaggggcctgcgcgtgtggttgtttCCAATACCCAACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGTAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAATGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTa                                <  2:105239/229‑1 (MQ=255)
                                                                                                                                                                                                                                gggCATCGCGTGGTGGTGC‑TTTCCAATCCC‑ATCCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTTGGCATTACCAGTATTCt                            <  1:170604/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gcgTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGaa                      >  1:598203/1‑250 (MQ=255)
                                                                                                                                                                                                                                        gtgggtggtgC‑TTTCCAACACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAg                    <  2:605748/248‑1 (MQ=255)
                                                                                                                                                                                                                                                         aaTACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATccc    <  2:625679/251‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCC‑AACAGCCTGCATAACACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTAACCCAGCGATACGGTCTTTTTCGGCGATAACGCCGCTAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCgg  >  1:491522/1‑250 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          
AACGCGTATTCCGCTGGCATCGCTTAAGAAGAGTTTTCATATGGGGGAGGCGTTTCATCAGCATGAGCGTGGGGAAATTAGCG‑ACGAAGCGTTCGCAGAGGCGCTGTGTCATGAGATGGCTCTACCGCTAAGCTACGAGCAGTTCTCTCACGGCTGGCAGGCGGTGTT‑TGTTGCGCTGC‑GCCCGGAAGTG‑ATCGCCATCATGCATAAACTGCGTGAGCAGGGGCATCGCGTGGTGGTGC‑TTTCCAATACC‑AACCGCCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGG  >  NC_000913/4075630‑4076126

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