Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6700TA28.3% 46.3 / 12.8 63pseudogene (934/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 T (17/21);  new base A (10/5);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGTTGCGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGG  >  NC_000913/3862441‑3862908
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
gTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGAT‑AA‑GGGTTCGGGGAT‑caa‑c‑aaag                                                                                                                                                                                                                                              >  1:3426405/1‑236 (MQ=255)
     catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑TA‑AGTTGGGGGGAA‑CAA‑C‑AAACCCGcc                                                                                                                                                                                                                                         >  1:1924177/1‑238 (MQ=255)
                cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑AGGTTCGGGGAA‑CAG‑C‑AATCCCCCCACATTTTCCCGa                                                                                                                                                                                                                             >  1:800336/1‑239 (MQ=255)
                           aaCATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑AGGTGCGGGGAA‑CAA‑CAAAACCCCCCCTTTTTCCCCGCTTTTCCgtgg                                                                                                                                                                                                                   >  1:3400569/1‑236 (MQ=255)
                                          ttttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAA‑GGTGGGGGGACA‑AAC‑A‑AAGCCCCCCCCTTTTCCCCATTTTTCCCGTGGGGCAG‑GGAGTAC‑a                                                                                                                                                                                                   >  1:1114516/1‑238 (MQ=255)
                                          tgtggggacaaaaCCTAATAAACCAGCCCTCCGCGCGTGCTCCCCCCCACAAAAAAGAAATACACTACGCGCGAAAGGGAGCCCCCCAACACGGCGTGCACCGAAAACGGCAGCGGTAAAT‑GGAACAGGAGGTTAAATTCAAGAGTTTCGGTAATCCCCACTGGGGGGGT‑GTTT‑AA‑GGTGGGGGGAAT‑AAC‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑At                                                                                                                                                                                                  <  2:2468702/230‑1 (MQ=255)
                                             tGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑GAGTTTGGGAAG‑AAG‑C‑ACACCCCCCACTTTTAACCCCTTTTTCGGTGCGGCAGAGAAGGAC‑Atag                                                                                                                                                                                                >  2:1769382/1‑236 (MQ=255)
                                                ggTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTT‑AA‑GG‑TGCGGGGAA‑CCG‑C‑AAACCCCCCACcttttc                                                                                                                                                                                                                                 >  1:3425762/1‑201 (MQ=255)
                                                           gTCAATCACCCCCCCGCCTATGTCGCCCCCCCCACAAAAGATATACATCCGGCTCAACATTGGGCCCGCCAGCACGGCGTGTCCCAAAAACAGAAGCGGAGAAA‑TGAACGGGAAGGTAAATTCCGGGGGTTCGTTAACCCCCCCCAGCATGGC‑GGTA‑AA‑GGTGGCGGAAAA‑CAC‑C‑AACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCa                                                                                                                                                                                <  2:500802/239‑1 (MQ=255)
                                                                                         cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑GGTTTCGGGAAT‑CAG‑A‑CAACCCGCCCCCTTTATCCGCTTTTCCCGTGCGGGCGGGAAGGAA‑AACGCTTAAAAAAAG‑CCCCCCCCACCAAAGC‑ATTTTTAG‑TTA‑CCCTGCa                                                                                                                                               >  2:1174966/1‑239 (MQ=255)
                                                                                           caccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGGTGGGGGGAT‑CAG‑C‑AAAACCCCCCCCTTTACCCGATTTTTCCGAGCGGCAGGGAAGTACAATGGCCAAAAAAATC‑CCCACGCCCCACAAAA‑TTTTTAGT‑TTA‑CCG‑GCAgg                                                                                                                                             >  2:1227261/1‑239 (MQ=255)
                                                                                                      tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCCAACAGCAACGGCGAAA‑TGAACAGGAAGGCAAATTCCAGCGGTTCTGTAATTCCCACCAGCATGGC‑GGTT‑AA‑TGTTTCGGGGAT‑CTGAC‑AACCCGCCCCCTTTTACCC‑ATTTTCCGGGGCGGGCGGGAAGGAA‑CTCACTCAAAAAAAA‑ACCCCCGCGCCCATAA‑CTTTTTAG‑TTATCCGTGCTGGGCaaaagacg                                                                                                                                   >  2:1054440/1‑235 (MQ=255)
                                                                                                        cATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑AGGTGCGGGAAA‑CAA‑C‑AAGACCGCCCCCTTTACCCCGTTTTTCCGTGCGGCCGGGAAGAAC‑AACGCGAAAGAAAAA‑CCCACGGCACCAAAAA‑ATTTTTAGTTTA‑CCG‑GCACGGCAAAAcaaccc                                                                                                                                 >  2:1794272/1‑233 (MQ=255)
                                                                                                                tCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑AGTTTCGGGGAA‑CAG‑C‑AAAACCCCCCCCTTTAACCGCTTTTTCGGGGCGGCAGTGGAgga                                                                                                                                                                                                      >  1:2519771/1‑164 (MQ=255)
                                                                                                                   aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑GGTTGGGGGGAT‑TAG‑C‑CAGCCCGGCACtttttc                                                                                                                                                                                                                                 >  1:2443745/1‑134 (MQ=255)
                                                                                                                   aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTT‑AA‑GG‑TGCGGGGAA‑CAA‑CAAAACCCCCCCCTTTTACCCGGTTTTCCGG‑GCGGGCGGGAAGGAA‑ACACCCAAAAAAAAG‑CCCCCCGCCCCAAAAATTTTTTGGG‑TTC‑CCGCG‑AGGGCAAAAACG‑CCTTCCGGGgaccac                                                                                                                    >  2:1757724/1‑239 (MQ=255)
                                                                                                                         aGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑AGTTGCGGGGAT‑CAA‑C‑CAACCCCCCCCCTTTTCCCCCTTTTTCGGTGGGGCCGGGAAAGAA‑AACGCTAAAAAAAAA‑CCCCCCGCCCCCAAAA‑TTTTTTAA‑GTT‑TCGTGGCGGGCAAAAACG‑CCCTCCTGGgacaaccaattt                                                                                                              >  2:2846435/1‑233 (MQ=255)
                                                                                                                                       gtgtctACCGCAAACACCCAGGGTAAAA‑TAAAAGGGAAGGTAAATCCCGGGGTTTCGGTAATTCCCACCGCCAGGGC‑GGTT‑GG‑GGTGGGGGGAAT‑CAG‑A‑AAACCCGCCACTATTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCAc                                                                                                                                                                               <  2:368282/159‑1 (MQ=255)
                                                                                                                                             aCCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑TA‑GGGTGCGGGGAA‑CAA‑CAAACCCCCCCCCTTTTACCCCTTTTTCCGGGGGGGCGGGGAGAAAA‑ACCCCAAAAAAAATC‑CCCCCGCACCCAAAAA‑TTTTTGGG‑TTA‑‑CGCGCAGGGGAAAAACC‑CCCTCCCGGGACAAAGCTTTCTCCGGCTCGgcgccccaaa                                                                                           >  2:1133718/1‑232 (MQ=255)
                                                                                                                                                gaaaaCAGCAA‑GGTGAAAGTGAAAGGGAGGGTAACTCCCGGGGGTGCGGAAACCCCCCCAGGCTGGGG‑GGTT‑AG‑GGGTGGGGGAAT‑CAC‑C‑ACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                                       <  2:1089777/237‑1 (MQ=255)
                                                                                                                                                 cAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑AGGTGCGGGAAT‑TAG‑C‑CAACCCGCCACTTTTACCCGATTTTCCCGTGCGGCAGTGAAAGAA‑ATCGCTAAAGAAAAGCCCCCCGGCCCCaaa                                                                                                                                                                      >  1:368282/1‑164 (MQ=255)
                                                                                                                                                    aCAGAAACGGTAAAA‑TGAACAGGAAGGAAAATCCCAGGGGTTCGGTAATCCCCACCAGCAGGGG‑GTTT‑AG‑GGTTGGGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTg                                                                                                                                                                                                                   <  1:2331164/117‑1 (MQ=255)
                                                                                                                                                      aaaaaCGGTAAAA‑TGAACGGGAAGTTAATTCCCGGCGGTTCGGTATTCCCCACCAGCTGGGC‑GTTT‑AG‑GGTTGGGAAAAA‑AAC‑A‑AACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCg                                                                                                          <  1:437563/211‑1 (MQ=255)
                                                                                                                                                         aaCGGGGAAA‑TGAACAGGAAGGAAATTCCCGGCGTTTCGTTAATCCCCACCAGCATGGC‑GGTT‑AG‑GGTGGGGGGAAC‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCg                                                                                                     <  1:2453067/216‑1 (MQ=255)
                                                                                                                                                             ggaaaa‑aGAACAGAAAGTAAAATCCCAGGGGTCCGGTAATCCCCACCACCGGGGG‑GGTA‑AG‑GGTGGGGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                                           <  1:1047547/235‑1 (MQ=255)
                                                                                                                                                             gTGAAA‑TGAACAGGAGAGTAAATGCCAGGGTTTCGGTAACCCCCACGGGCAGGGG‑GTTT‑AA‑GGTGGGGGAAAT‑CAG‑C‑AGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                                           <  2:279067/238‑1 (MQ=255)
                                                                                                                                                              tGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑AA‑AGGTTCGGGGAT‑CAG‑C‑AAACCCCCCACCTTTTCCCCATTTTTCCGTGCGGGAGGGaa                                                                                                                                                                                                         >  2:2331164/1‑117 (MQ=255)
                                                                                                                                                                aaa‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑AGTTGCGGGGAT‑CAG‑C‑AAGCCCCCCCCTTTTTCCCCATTTTTCGGTGCGGCAGTGAAGTAC‑ATCGCTTAAGAAATG‑CCCACGGCACCAAAGA‑tttttta                                                                                                                                                            >  1:640393/1‑155 (MQ=255)
                                                                                                                                                                aaa‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGT‑TA‑GGGTGGGGGGAT‑CAg                                                                                                                                                                                                                                                     >  2:2384981/1‑73 (MQ=255)
                                                                                                                                                                aaa‑TGAACAGGAAGGAAAATCCCGGGGTTTGGTTAATCCCCCCCAGCAGGGC‑GTTT‑AG‑GGTGGGGGAAAC‑CGG‑C‑AGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtg                                                                           <  1:2194497/235‑1 (MQ=255)
                                                                                                                                                                 aa‑TGAACAGGAAGGAAATTCCCGGGGGTTCGGTATTCCCCCCCAGCATGGG‑GGTT‑AG‑GGTGGGGGAAAT‑CAG‑C‑AACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                      <  1:153194/239‑1 (MQ=255)
                                                                                                                                                                    tGAAAAGGAAGGAAAATTCCGGCGTTTCGTTAATCCCCCCCAGCAGGGC‑GGTT‑AG‑GGTGGGGGAAAT‑AAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                    <  1:1133718/239‑1 (MQ=255)
                                                                                                                                                                     aaaagggAAGTTAAATCCCAGCGGTTCGTAAATTCCCCCCATGAGGGC‑GGTT‑AA‑GGTTGCGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGccc                                                                   <  1:2376252/234‑1 (MQ=255)
                                                                                                                                                                      aaCAGGAAGTTAAATTCCGGCGGTTGGGTAATCCCCCCCAGCAGGGC‑GTTT‑AG‑GGTGGGGGAAAC‑CAC‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGAt                                                                                                                 <  1:1279858/192‑1 (MQ=255)
                                                                                                                                                                              ggaaaaTCCCAGGGTTTCGTTAATCCCCCCCAGCAGGGG‑GGTT‑AA‑GGTGGGGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGAttt                                                                                                                                                                                                                          <  1:3561767/82‑1 (MQ=255)
                                                                                                                                                                               gaaaaTCCCGGCGTTTCGTTATTCCCCCCCAGCTTGGC‑GTTT‑AA‑GGTTGCGGAAAA‑CAG‑A‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                    <  1:1769382/226‑1 (MQ=255)
                                                                                                                                                                                tAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑TA‑GGTTGCGGGGGA‑CAG‑G‑AAGGCCGCCACTTTTTCCCCATTTTTCGGTGCGGCAGGGGAGGAC‑ATCCCCAAAGAAATG‑CCCACGGCACCCAAGA‑ATTTTTAG‑TTA‑CCGTGCAGGGCAAAACCGCCCTTCCGGGAACAA‑GATTTTAGCgggc                                                                                                       >  2:1279858/1‑190 (MQ=255)
                                                                                                                                                                                aaaaCTCCGGCGGTTGAGTAATCCCCGCAAGCGTGGG‑GTTT‑GA‑GTTGGGGGAAAC‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                                    <  2:1975488/226‑1 (MQ=255)
                                                                                                                                                                                             ttCGGTAATCCCCACCAGCGTGGC‑GGTT‑AA‑GGTGGCGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAAct                                                                                     <  1:566194/197‑1 (MQ=255)
                                                                                                                                                                                                  aactcccccccAGCATGGC‑GGTT‑AA‑GGTTGCGGAAAT‑CAG‑C‑ACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAg                                                                                                                                                           <  2:2023933/120‑1 (MQ=255)
                                                                                                                                                                                                    atccccACCAGCATGGC‑GGTT‑AA‑GGTGGCGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGc                                    <  1:1174966/236‑1 (MQ=255)
                                                                                                                                                                                                     tcccccccGGCAGGGC‑GTTT‑AA‑GGTGGGGGGAAC‑CAG‑A‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAAt                                                                                                                                                                                      <  1:1717048/95‑1 (MQ=255)
                                                                                                                                                                                                         caccaGCGGGGC‑GGTA‑AA‑TGTGGCGGAAAT‑CAG‑C‑AGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                               <  1:1237659/239‑1 (MQ=255)
                                                                                                                                                                                                         caccaGCAGGGG‑GGTT‑AA‑GGGGGGGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc                               <  1:884707/239‑1 (MQ=255)
                                                                                                                                                                                                          cccaGCTGGGC‑GGTT‑AA‑GTTTGGGGAAAA‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTgcg                                                                      <  1:2738076/198‑1 (MQ=255)
                                                                                                                                                                                                          cccaGCA‑GGG‑GGTT‑AAGGGTTGCGGAAAC‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCt                                                                                                <  1:138546/172‑1 (MQ=255)
                                                                                                                                                                                                           ccaGCTGGGC‑GGTT‑AG‑GTTGGCGGGAAT‑AAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                             <  1:3468314/239‑1 (MQ=255)
                                                                                                                                                                                                           ccaGAAGGGC‑GGTG‑AG‑GGTGGCGGAAAC‑CAA‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                             <  1:2199944/239‑1 (MQ=255)
                                                                                                                                                                                                            caGCATGGC‑GGTT‑AA‑GGGTGCGGGGAT‑CAG‑C‑AAACCCCCCCCTTTTTCCCGCTTTTCCCGTGCGGCCGGGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTTAG‑GTT‑CCGCGCAGGGCAAAACCC‑CCTTCCGGGGACAACGATTTTTGCCGCTCCGCCGTCCACCCGCACTCTTGCTGAAGACGCGCCCAGTACCTATCACTGCCCCCCTCTACAGCTccccggcaca                            >  2:2199944/1‑232 (MQ=255)
                                                                                                                                                                                                               cATGGC‑GGTT‑AA‑GGTTGCGGGGAT‑CAG‑C‑AAGCCCCCCCCTTTTTCCCGATTTTCCGGGGCGGGAGGGGAGGAC‑a                                                                                                                                                                                                   >  1:3029603/1‑74 (MQ=255)
                                                                                                                                                                                                               cATGGC‑GGTT‑AA‑GGTGGGGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCgg                                                                                                                                                                                                                     <  1:448399/57‑1 (MQ=255)
                                                                                                                                                                                                                    c‑GGTT‑TA‑GGGTGCGGGAAT‑CAG‑C‑AAGCCCGCCCCTTTTTCCCCATTTTTCCGGGCGGCAGTGGAGGAC‑ATCGCTTAAGAAATG‑CCCCCGGCACCCAAAA‑TTTTTTAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTTAGCGGCTCCGCGGTCCAAACGAAc                                                                                      >  2:138546/1‑173 (MQ=255)
                                                                                                                                                                                                                     ggttt‑AA‑GGTTGCGGAAAC‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                     <  1:2142843/235‑1 (MQ=255)
                                                                                                                                                                                                                             gagtggCGGGAAACCAG‑C‑AACACCGACACTTTTACCCGATTTTCCGGTGCGACAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa                     <  2:2922638/228‑1 (MQ=255)
                                                                                                                                                                                                                               gtgggggaaaT‑AAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt             <  1:1496621/230‑1 (MQ=255)
                                                                                                                                                                                                                               gtggCGGAAAC‑CAG‑C‑ACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt             <  1:3410793/235‑1 (MQ=255)
                                                                                                                                                                                                                                ttGCGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCa                                                                                                                                                                        >  1:2078901/1‑86 (MQ=255)
                                                                                                                                                                                                                                tggCGGAAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGt             <  1:2089341/235‑1 (MQ=255)
                                                                                                                                                                                                                                   cGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccgcc                                                <  1:2274838/199‑1 (MQ=255)
                                                                                                                                                                                                                                   cGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgccgcc                                                >  2:2274838/1‑199 (MQ=255)
                                                                                                                                                                                                                                    ggaaaT‑CAG‑C‑AACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg        <  1:2799748/235‑1 (MQ=255)
                                                                                                                                                                                                                                        aT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTg                                                                                                                                                 >  1:1649968/1‑98 (MQ=255)
                                                                                                                                                                                                                                        aT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAg           >  1:2922638/1‑231 (MQ=255)
                                                                                                                                                                                                                                         t‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAgtgg  >  1:2983729/1‑239 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
GTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAA‑TGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTT‑AA‑GGTTGCGGGAAT‑CAG‑C‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAC‑ATCGCTAAAGAAATG‑CCCACGGCACCAAAGA‑TTTTTGAG‑TTA‑CCGTGCAGGGCAAAACCG‑CCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATGAAGTGG  >  NC_000913/3862441‑3862908

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