Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6670GA33.3% 62.4 / 19.6 54pseudogene (937/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base G (26/10);  new base A (0/18);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.52e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGTTGCGGGAATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATG  >  NC_000913/3862430‑3862902
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑AGGTGCGGGa                                                                                                                                                                                                                                                  >  1:48787/1‑239 (MQ=255)
tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑AGGTAATTCCCACCAGCATGGCGGGTT‑AA‑GG‑TGCgggg                                                                                                                                                                                                                                                  >  1:926614/1‑238 (MQ=255)
tGGAGATTTGTGGCAGCTTCATATCCGCGTGGTTGCTGAACATCGGGTTCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCGCCGCCCAGGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGTTGCGGGa                                                                                                                                                                                                                                                  <  2:1259094/239‑1 (MQ=255)
     aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGCGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑CT‑CGGTAATTCCCACCAGCATGGC‑GGCT‑AA‑AGTTGCggggaacc                                                                                                                                                                                                                                              >  1:719433/1‑233 (MQ=255)
                          gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGGTGCGGGAAACAGC‑AAGCCCGCCA‑CCTTTACCCg                                                                                                                                                                                                                      >  1:1805086/1‑239 (MQ=255)
                          gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑AGTTGCGGGAAAAAGC‑AAGACCCCCA‑CCTTTAtccc                                                                                                                                                                                                                      >  2:2506911/1‑238 (MQ=255)
                           cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGCGGGTT‑AA‑GGTTGGGGGAACAGGC‑AAACCCGCCC‑CTTTTTccga                                                                                                                                                                                                                      >  1:2765707/1‑237 (MQ=255)
                               gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGTTACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGCT‑AA‑AGTTTCGGGAAACACG‑AAGCCCCCCACCTTTTAccgcatt                                                                                                                                                                                                                   >  1:89577/1‑233 (MQ=255)
                                          tCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTTAAA‑GGTTG‑GGGGATCAACAAACCCCCCCC‑CTTTTCCCCGTTTTTCCGGGGGGGc                                                                                                                                                                                                       >  1:692300/1‑239 (MQ=255)
                                           cggtgtCCAGTTTTGCAGTAAAA‑CAGAGCCAGACAGACCTCCGCCCATCGTC‑CCCCCCACCCAAAAGAAATAC‑TTACGGGTCGA‑AATTGGGGCCGCCAGCAGGGCGTGTACCCAAAAAAGAAACG‑GGGAAAAGAACA‑GAAAGGAAAATCCCGGCGT‑TT‑CGTTAATCCCCACCAGCAGGGG‑GTTT‑AG‑GTTTGGGGAAATCAAC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCa                                                                                                                                                                                                      <  2:1308690/233‑1 (MQ=255)
                                                   aGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑AGGTGGGGGGAACAGC‑AAAACCCCCC‑CTTTTTCCCGCTTTTTCCGGGCGGGAGGGaacaaa                                                                                                                                                                                             >  1:2434753/1‑235 (MQ=255)
                                                      tttGCAGTAAAACCCTGTTCAACAAGCCCCCCG‑CCATGGTGCCCCACCCCACAAAAGAGATAAATT‑C‑GGGCGACCATGGAGGCCGCAAGCAGGGGGTCTACCGAAAACAGAAAGG‑GAAAAATGAACA‑GGAAGGT‑AATTCCCGGGG‑TTGCGGTAACCCCCCCCACCAGGGG‑GGTA‑AG‑GGTGGCGGAATTCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTa                                                                                                                                                                                             <  2:1600021/236‑1 (MQ=255)
                                                      ggtGCGGTAAAA‑CCTACACAATCAGCCCACGGCCCGGGCAC‑CCCACCACACAGAAGGCTTCCACTCG‑GGCCGC‑CATGGGGCCCGCAAGCAGGGCGTGTCCCACAACCGGCAAGG‑GTAAAATGAAAA‑GGAAGGTAAATCCCAGGGT‑TT‑GGGTAATCCCCCCCAGCCTGGC‑GTTT‑AG‑GGTGGGGGAAATCAGC‑AAACCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTa                                                                                                                                                                                             <  2:692300/233‑1 (MQ=255)
                                                        tGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑GGGTGCGGGAATCAAC‑AAACCCCCCA‑CTTTTTCCCGATTTTTCCGTGCGGCAGTggaaaaaaac                                                                                                                                                                                          >  1:340525/1‑228 (MQ=255)
                                                              aaaa‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGCGGGTT‑AG‑GGTTCGGGGaacaaac                                                                                                                                                                                                                                            >  2:1168283/1‑179 (MQ=255)
                                                                            tCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑GGTTGCGGGGATCAGC‑AAGCCCCCCC‑CCTTTTTCCGCGTTTTCGGTGCGGCCGGGAAGGAAAATCGCCAAAGAAATAccccggg                                                                                                                                                                      >  2:135909/1‑234 (MQ=255)
                                                                             gcagccTCCGCCCATGTCG‑CCCCCCACGCACAAGAGATAGACCAC‑GGTCAG‑AATTGGGCCCCACGACCCGGGGGGCCCCAAAAACGCGAACG‑GAAAAAAGAACA‑GGAAGATAATCCCCAGCGG‑TA‑CGGAAATCCCCACCAGCAGGCC‑GGTT‑GA‑GGTTGCGGAAATCGGC‑ATGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACgg                                                                                                                                                                      <  2:2315545/232‑1 (MQ=255)
                                                                                                     ccacaccaAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑AGTTGCGGGGATACGC‑AAACCCCCCA‑CTTTTACACCACTTTTCGGTGCCGCAGTGAA‑GTACAACTCTCAAAAAAAAGCCACGGCA‑GCCAAAACTTTTGTGTTatcatgt                                                                                                                                           >  2:1829381/1‑233 (MQ=255)
                                                                                                                              tCGA‑CGTTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGGT‑AAGGGTTGGGGGAA‑AAGCAAACCCCCCCC‑CTTTTACCCCATTTTCCGG‑GGGGCCGGGAAGGAAAACCGCT‑AAGAAAAGCCCCCGGCACCCAAAGATTTTTAGGTTACGGGCCA‑GGGCAAAACCCCATCCCGGGaa                                                                                                                   >  2:1725538/1‑239 (MQ=255)
                                                                                                                                 a‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑AGGTGCGGGAATCCGC‑AAGCCCGCCC‑CTTTTACCCCGTTTTTCGGTGCGGGAGTGGA‑AGACATCGCTAAAGAAAAGCCCCCGGCC‑CCCAAGAATTTTTAGTTTACCTGCACGGGAAAAACCCCCTCCCGGGaacac                                                                                                                >  2:463228/1‑238 (MQ=255)
                                                                                                                                      gcagcccgccAGCGGGGCGGGTACCGCAAA‑AACAACGTGTAAGATAACAA‑AGATGATACCTCCCGGCGG‑TT‑GGTCAATCCCCACCAGCAGGGC‑GGTT‑AA‑GGTTGCGGGAAACAGC‑AACCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTGAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGaac                                                                                                                  <  2:2765707/227‑1 (MQ=255)
                                                                                                                                                     ggCTTGAACCGAAAACAAAAACG‑GT‑AAAAGAACAGGGATGGAAAATGCCAGCGG‑TT‑CAGTAATCCCCACCAGCATGGG‑GGTT‑AG‑GGTGGGGGAAAACAAA‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGgcgc                                                                                           <  2:463885/239‑1 (MQ=255)
                                                                                                                                                       cGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTG‑AA‑GGTTGGGGGGATCAAC‑AAACCCGCCC‑CTTTTTTCCGAATTTCCGGGGCGGCAGTGGA‑AAACATCCCTAAAGAAAAGCACACGGCA‑CCCAAGAATTTTTAGTTAACCGGGA‑GGGCGAAACAGCCCTACGCGGACAACGATATTCA‑CGGCTCCgcgcgc                                                                                         >  2:2500554/1‑237 (MQ=255)
                                                                                                                                                              cGCAAACAGAAACG‑GTAAAAGAAAGA‑GGAAGTAAAATCCCAGCGG‑TT‑CGGTAATCCCCCCCAGCATGGG‑GGTT‑AG‑GGTGGCGGAAATCAGC‑AACCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCa                                                                                                                                                                                                      <  2:544234/127‑1 (MQ=255)
                                                                                                                                                              cGCAAACAACAACG‑GTAAAAAGAACA‑GGAAGGAAATTCCCAGCGG‑TT‑CGGAAATCCCCCCCAGCAGGGC‑GGTT‑AG‑GGTTGCGGAAATCAGA‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGa                                                                                  <  1:463228/239‑1 (MQ=255)
                                                                                                                                                               aaaaaCAGCACGG‑GAAAAAGGAACG‑GGAAGTAAAATCCCGGGGT‑TT‑CGGAAATTCCCACCAGCTTGGG‑GTTG‑AA‑GGTGGGGGGAATCAAC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGaa                                                                                 <  1:2500554/237‑1 (MQ=255)
                                                                                                                                                                 aaaCAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑AGGTGGGGGGATCACG‑AAGCCCCCCA‑CTTTTACCCCATTTTTCGGTGCGGGAGTGAA‑GGACATCGCTAAAGAAAAGCCCACGCGA‑CCCAAGAATTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCCGGAACAAAGAT‑TTCTGCGCGGCGGCGCTCACACAGAAc                                                                                >  2:525475/1‑238 (MQ=255)
                                                                                                                                                                   aCAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGGTGCGGGAAACAGC‑AAGCCCGCCC‑CCTTTACCCGATTTTTCCGGGCGGCGGTGGA‑GGACATAGCTAAAAAAATACCCACGGCC‑CCCAAGATTTTTGAGTTACCGTGCA‑GGGCGAAAACGCCTTCCCGGAACAAACAT‑TTCCTCGGGTCGGCGCTCACACAGacaac                                                                              >  2:2502053/1‑234 (MQ=255)
                                                                                                                                                                        aaCG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑AGTTGCGGGGATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTTCGGGGCGGGAGTGAA‑GTACa                                                                                                                                                                                           >  1:544234/1‑127 (MQ=255)
                                                                                                                                                                        aaCG‑GAGAAAGGAACA‑GGAAGGAAAATCCCGGGGT‑TT‑CGATAATCCCCACCAGCATGGC‑GGTT‑AA‑GGTGGCGGGAATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAg                                                                        <  1:1090756/239‑1 (MQ=255)
                                                                                                                                                                             ggaaaaagaaaa‑gGAAGAAAAATTCCGGCGG‑TT‑GGTTAATCCCCACCAGCAGGGG‑GGTT‑AA‑GGTTGGGGGAATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAg                                                                        <  1:27383/232‑1 (MQ=255)
                                                                                                                                                                               aaaaTGAACG‑GGAAGGGAAATCCCAGCGG‑TT‑GGGTAATTCCCCCCAGCGGGGC‑GGTA‑AA‑GGTGGCGGGAATCAGA‑AGGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                   <  1:124248/237‑1 (MQ=255)
                                                                                                                                                                                   tGACCA‑GAAAGGAAAATCCCAGCGG‑TG‑CGGTAATCCCCACCAGCTTGGC‑GGTT‑AG‑GGTAGCGGAAAGCAGA‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  1:1061147/239‑1 (MQ=255)
                                                                                                                                                                                   tGAAGA‑GGAAGAAAATTCCCGGCGT‑TT‑CGGTAATCCCCACCAGCAGGGC‑GGTT‑AG‑GGTTGGGGAAATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  1:2242422/239‑1 (MQ=255)
                                                                                                                                                                                     aaaa‑GGAAGTAAAATTCCGGCGT‑TT‑CGGTAATCCCCCCCAGCAGGGC‑GGTT‑AG‑GGTTGGGGGAACCAGC‑AACCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  1:2040202/234‑1 (MQ=255)
                                                                                                                                                                                           gagggTAAATTCCAGGGG‑TA‑CGGTAATTCCCACCAGCTGGGG‑GGTT‑AG‑GGTGGGGGGAACAAGA‑AGGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  1:2738873/229‑1 (MQ=255)
                                                                                                                                                                                            gaGGTAAAATCCCGGGGTTT‑CGTTAATCCCCCCCAGCGGGGC‑GGTT‑AG‑GGTGGGGGAAACCAAC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGaa                                                                                 <  1:1626033/212‑1 (MQ=255)
                                                                                                                                                                                            aaGGTAATTCCCAGCGG‑TT‑CGGTAATCCCCCCCAGCGTGGC‑GGTT‑AG‑GGTGGCGGAAATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  1:1834308/231‑1 (MQ=255)
                                                                                                                                                                                            aaGGTAAATCCCGGGGT‑TT‑CGGTAACTCCCCCCAGCATGGG‑GTTG‑AA‑GGTGGGGGAAACCAGC‑ACGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                           <  1:2021630/234‑1 (MQ=255)
                                                                                                                                                                                                tAAATCCCAGCGG‑TC‑CGTTAATTCCCCCCAGCAGGGC‑GTTT‑AA‑GGTGGGGGAAATAAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                           <  1:2649972/230‑1 (MQ=255)
                                                                                                                                                                                                  aaTTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑GGGTTTGGGAATCAAA‑AAAATCCCCA‑ATTTTTACCGCTTTTCCg                                                                                                                                                                                                              >  1:920188/1‑85 (MQ=255)
                                                                                                                                                                                                     tCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGTTGCGGGGATCAGC‑AAGCCCGCCC‑CTTTTACCCGGTTTTTCGGGGCGGCAGTGAA‑GGACATCGCTAAAGAAATGcc                                                                                                                                                                           >  1:584809/1‑116 (MQ=255)
                                                                                                                                                                                                         gCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGTTGCGGGGATCAGC‑AAACCCGCCC‑CTTTTACCCGGTTTTTCGGTGCGGCAGGGAA‑GTACATCGCTAAAAAAAAGCCCACGGCA‑CCAAAGATTTTTGa                                                                                                                                                     >  1:895088/1‑133 (MQ=255)
                                                                                                                                                                                                                             ccgtaaTGGG‑GGTT‑AA‑GTTTGCGGAAATCGAA‑AACCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                       <  1:1866628/234‑1 (MQ=255)
                                                                                                                                                                                                                             ccaGGAGGGC‑GGTT‑AA‑GGTTGCGGAAACCAGC‑AACCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                       <  2:1742183/239‑1 (MQ=255)
                                                                                                                                                                                                                             ccaGCAGGGC‑GGTT‑AA‑GGTGGGGGGAACCAGA‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                       <  1:135909/239‑1 (MQ=255)
                                                                                                                                                                                                                                     ggtggTT‑AA‑GGTGGCGGAAAGCAGC‑AACCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATa                 <  1:2080204/235‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTT‑AATGGTGGCGGAAAGCAGA‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa               <  1:2614236/238‑1 (MQ=255)
                                                                                                                                                                                                                                            aaa‑GGTGGCGAAAATCACA‑AACCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc          <  1:2095503/238‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCGGGAATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGaa                                                                                                                                                                                >  1:367584/1‑69 (MQ=255)
                                                                                                                                                                                                                                                     cGGGTATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GCACACCGCTAAAGAAATGCCCAAGGCA‑CGAAAGATTTTTGAGTTAACGTGTA‑GGGCCAAACCGCCTACCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCCAACTGAAATCTTGCAGATGCTGCCTCCAGGACACTTGACTGCCGCCTTCCACAGCTGGCGGACCAAAGCTAACCTGTCCGtt      >  2:984561/1‑234 (MQ=255)
                                                                                                                                                                                                                                                     cGGGAATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAAct                                                                               <  1:848666/162‑1 (MQ=255)
                                                                                                                                                                                                                                                     cGGGAATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAAct                                                                               >  2:848666/1‑162 (MQ=255)
                                                                                                                                                                                                                                                     cGGGAATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg  >  1:2096498/1‑239 (MQ=255)
                                                                                                                                                                                                                                                     cGGGAATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg  >  1:784119/1‑239 (MQ=255)
                                                                                                                                                                                                                                                     cGGGAATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACACCGCTAAAGAAATACCCACGGCA‑CCAAAGATTTTTGAGGTACCGCGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACAACCGAAAGCCGCCTTCAAAAGCTCACAGGCACAAGCAAAACTGTACGTAGaat  >  2:709623/1‑237 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
TGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAA‑CCTAGTCAATCAGACCTCCGCCCATGTTG‑CCCACCACACCAAAGAGATACATTAC‑GGTCGA‑CATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACG‑GTGAAATGAACA‑GGAAGGTAAATTCCAGCGG‑TT‑CGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGTTGCGGGAATCAGC‑AAGCCCGCCA‑CTTTTACCCGATTTTCCGGTGCGGCAGTGAA‑GTACATCGCTAAAGAAATGCCCACGGCA‑CCAAAGATTTTTGAGTTACCGTGCA‑GGGCAAAACCGCCTTCCGGGAACAACGAT‑TTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATG  >  NC_000913/3862430‑3862902

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