Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,465,7870GA39.2% 54.9 / 24.6 53R132C (CGC→TGC) treBtrehalose‑specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base G (22/9);  new base A (0/20);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 1.68e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGG‑CGTGCCGCCCATTTTTTT‑CACC‑GCTGACCA‑GCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGT  >  NC_000913/4465542‑4466031
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
cGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGT‑CGTGCCGCCCATTTTTTT‑CACC‑GCTGACCA‑GCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAAc                                                                                                                                                                                                                                                    <  1:429227/251‑1 (MQ=255)
 ggCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGG‑CGTGCCGCCCATTTTTTT‑CACC‑GCTGACCA‑GCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAAcc                                                                                                                                                                                                                                                   >  2:15708/1‑251 (MQ=255)
 ggCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGG‑CGTGCCGCCCATTTTTTT‑CACC‑GCTGACCA‑GCAAATACCGACCGGCAGGGAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACCGAAAATCGTAAATCGGTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTTACCGTT‑‑GCTCATGGGCCAA‑T‑CGCCGATCACATTGCGAAAAcc                                                                                                                                                                                                                                                   >  2:45954/1‑251 (MQ=255)
    aGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGG‑CGTGCCGCCCATTTTTTT‑CACC‑GCTGACCA‑GCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑Ag                                                                                                                                                                                                                                               <  1:291278/251‑1 (MQ=255)
       tgctgcCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGG‑CGTGCCGCCCATTTTTTT‑CACC‑GCTGACCA‑GCAAATACCGACCGGCAGGGAGAAGAAGATTGCTTC‑ACCGATCAA‑CCACAGAAAATCGGAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCCGACCGTT‑‑GATCATGGGCAAA‑T‑CCCCGATCACATTGCGAAAACCG‑AGG‑At                                                                                                                                                                                                                                           >  2:571019/1‑251 (MQ=255)
                 gcataTAAGCGTTCACCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGACCGG‑CGTGCCGCCCATTTTTTT‑CACC‑GCTGACCA‑GCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACAGAAAATCGTAGATCGTTTGCGGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAACCCG‑AGG‑ATAAAA‑ccgccg                                                                                                                                                                                                                                <  2:131594/247‑1 (MQ=255)
                                          aaaaaCCAGGGTCACGCAAACCACGATACCAAGGATGGG‑GGTGCCCCCCTTTTTTTC‑CCCC‑GCGGACAA‑CAAAATCCCGCCCGGCAGGAAGAAGAAG‑TCCTTTCCACCGATAAA‑CCAAAAAAAATCGTAGATCTTTGGCAGGGAAGGTTACATGGGCGCCAGCGTCTGCCCGTT‑‑GCTCAGGGGCAAA‑T‑CGCCGACCACATTGGAAAAACCG‑AGG‑ATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGga                                                                                                                                                                                                       <  1:288208/249‑1 (MQ=255)
                                           aaaaCCAGTGCCCCGCCAAGCACGAAACCAAGAAGCGG‑GGTGCCCCCAATTTTTTC‑CCCC‑GCTAACCA‑GAAAACCCCGGCCGGCAGAAAGAAAAACACCGTTTC‑ACAGATAAA‑CCAAAAAAAATCGAAGATGGTTAGCAGGGAAGGGTACATGCGCGCAAGCTTCTGCCCGTG‑‑CCTGAGGGGAAAA‑C‑CGCCAACCACATGGAAAAAACCA‑AGA‑ACCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGgaa                                                                                                                                                                                                      <  2:79887/250‑1 (MQ=255)
                                            aaaCCAGTGCCCCCCCAACCACGA‑ACCCAGGGTCGGGCGTCCCGCCCATTTTTTC‑CCCC‑CCTGCCCA‑GAAAATACCGCCGGGCGGGAAAAAGAAGATCGCTTC‑CCCAACCAA‑CCACAAAAAATCTGAGTTGGTTTGGAGGAAAGGGTACTTTTGCGCAAGGGTCGGCCCGGT‑‑GCCCATGGGCAAA‑T‑CGCCGATCCCATTGC‑AAAACCG‑GGGGATCAAC‑CCGCCGCTAATCAACGCGGGCAGCAACGGgaag                                                                                                                                                                                                     <  2:446552/251‑1 (MQ=255)
                                               ccAGTGCCCCGCAAAGCACAATCCCAAGGATCGG‑CGTCCCGCC‑ATTTTTTT‑CCCCCGCTGACCA‑GAAAATCCCGCCCGGCGGGTAGAAGAAGACCGCTCC‑CCCGATAAC‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGTTACATTGGCGCCAGCTTCTGACCGTT‑‑GCTCATGGGAAAA‑T‑CGCCGATCACATTGCGAAACCCG‑AGG‑ATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGaagaag                                                                                                                                                                                                  <  1:485205/251‑1 (MQ=255)
                                                 ggtgtCACGCCAAGCACGATCCCAAGGATCGG‑CGTGCCGCCCTTTTTTTC‑CCCC‑GCTGACCA‑GCAAATCCCGACGGGCAGGTAGAAGAAGATCGCTTC‑ACCAATCAA‑CCACAGAAAATCGTAGATCTTTTGCGGGGAGGGTTACATTTGCGCCAGCGTCTGACCTTT‑‑GCTCAGGGGCAAA‑T‑CGCCGATCACATGGCAAAACCCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGAt                                                                                                                                                                                                <  1:526482/250‑1 (MQ=255)
                                                 aggtacACGCAAAGCACGATCCAAAGGATGGG‑CGGCCCGCCTATTTTTGT‑CCCC‑GCT‑ACCAAGCAAAACCCGCCCGGGAGGAAGAAAAAAACCGCTAC‑ACCGAACAA‑CCACAGAAAATCGAAGACCGTTTGCGGGGAAGGGTACATTGGCCCCAGCGGCTGACCGTT‑‑GCACATGGGAAAA‑T‑CGCCAAGCACATGGCAAAAACCG‑AGG‑AAAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGAt                                                                                                                                                                                                <  1:483884/246‑1 (MQ=255)
                                                    gtccccccAAGCAGGATCCCAAGGAGCGG‑CGTGCCGCCAATTTTTTC‑CACC‑GCGGCCCA‑ACAAATCCCGACCGGCAGGAAAAAGAAGATCTTTCC‑CCCGATCAC‑CCCAAGAAAAGGGTAGATCGTTGGCAGGGAAGGTTAAATTGGCGCCAGCGTGTCAGCTTG‑‑GCTCAGGGGAAAA‑G‑CGCCGATCACATTGCAAAAGCCG‑AGG‑AACAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCAc                                                                                                                                                                                             <  1:583854/245‑1 (MQ=255)
                                                           cacacccGATACCAAGGATGGG‑CGTCCCTCCCTTTTTTTC‑CCCC‑GCTGACCC‑TCAAATCCCGACCGGCAGGTAAAAGAGGACCTCTAC‑CCCAACCAC‑CCACAAAAAATCGTATACCGTTTGGGGTTAAGGTTCCATTTGCCCCAGCGTTTGCCCTTT‑‑TCCCTTGGGCAAG‑T‑CGCCTATCACATGGCAAAAACGG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGaaa                                                                                                                                                                                      <  1:390033/245‑1 (MQ=255)
                                                           cAAGCACGATACCAAGGATCGG‑CGTGCCGCCCATTTTTTT‑CACC‑GCTGACCA‑GCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGaa                                                                                                                                                                                       >  2:416988/1‑250 (MQ=255)
                                                                   aTACCAAGGATCGG‑CGTGCCGCCCATTTTTTT‑CACC‑GCTGACCA‑GCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑ACACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTTCATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgag                                                                                                                                                                                  >  2:396286/1‑247 (MQ=255)
                                                                      aaaaGGAGCGG‑CGTGCCGACCATTTTTTA‑CACC‑GGGAAACA‑ACAAAAACCGACGGGAAGAAAGAAAAAGAGCGCTTA‑ACC‑ATCAAACCAGAAAAAATCGTAGAGTTTTGGCGGGGAGGGGAACAGTTGCGCCAGCGTCTGACCGTT‑‑GGCGAGGGGAAAA‑T‑CGCGGATAACATGGAAAAAACAG‑AGA‑ACAAAA‑CCGCCGATAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCaa                                                                                                                                                                           <  2:187056/249‑1 (MQ=255)
                                                                       aagggATCGG‑CGCGCC‑CCCATTTTTTTTCCCC‑GCTAACCA‑GAAAATCCCGCCGGGCAGGTAGAAGAAGATCGCTCC‑ACGAATCAA‑CCAAAGAAAAGCGTAGATCGTTGGCAGGGAAGGTTAAATTGGCGCCAGCGTCTGCCCGTG‑‑GCCAATGGGAAAA‑G‑CCCCGATCACATTGCAAAAACCG‑AGT‑ACAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAAc                                                                                                                                                                          <  1:359221/248‑1 (MQ=255)
                                                                                   gTGCCGCCCATTTTTTT‑CACC‑GCTGACCA‑GCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGcc                                                                                                                                                                >  1:148033/1‑250 (MQ=255)
                                                                                    tGCCGCCCATTTTTTT‑CACC‑GCTGACCA‑GCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCGGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCAt                                                                                                                                                              >  1:30497/1‑251 (MQ=255)
                                                                                         ccctttttttt‑CCCC‑GCTCACCA‑GCAAACCCCGCCCGGCAGGAAGAAGAAGATCGCTTC‑ACCGAACAA‑CCACAAAAAACCGTAGATCGTTTGCGGGGAGGGTTACATTTGCGCCAGCTTCTGACCGTT‑‑GCTCAGGGGAAAA‑T‑CGCCGATCACATTGGGAAACCCG‑GGA‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCat                                                                                                                                                         <  1:491895/247‑1 (MQ=255)
                                                                                            attttttcccACC‑GCTGACCA‑GCAAAGCCCCACCGGAAGTAAAAAAAAGACCGCTCC‑CCCGAACAA‑CAACAAAATATCGTAGTTCGTTGGCAGGGAGGGGTACATTGGCGCCA‑‑GTTTCAGCGTTTCGCTCAGGGAGAAG‑T‑CGCCGTTCACATTGC‑AAAACCG‑GAG‑TTCAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCTCGAAATGAGAGATCAACTGCTCATGCCATTTCatat                                                                                                                                                       <  1:36225/242‑1 (MQ=255)
                                                                                                       cc‑GCTGACCA‑GCAAACCCCACCCGGCAGGTAGAAGAAGATCGCTTA‑CCCGACCAA‑CCACAAAAAAGCGTAGATCGTTTGCAGGGAAGGTTACTTTTGCGCCAGCGTCTGCCCGTT‑‑GCTCATGGGCAAA‑T‑CGCCAATCACATGGCGAAACCCG‑AGA‑AACAAA‑CCGCCGCTAATCAACGCGGGCGGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGggcg                                                                                                                                            <  1:447854/251‑1 (MQ=255)
                                                                                                             gACCA‑GCAAATCCCGACCGGCAGGTAGAAGAAGAGCCCTCC‑CCCAATCAA‑CCAAAAAAAATCGTAGATGGTTGGGAGGGAGGGTTACATTTGCGCCAGCGTCTGCCCTTT‑‑GCTCA‑GGGCAAAAC‑CGCCGATCCCAT‑GGGAAAACCGGAGG‑AAAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCttt                                                                                                                                       <  1:519576/251‑1 (MQ=255)
                                                                                                             gACCA‑GAAAATCCCGCCCGGCGGGTAGAAGAAGACCGCTTC‑ACCGATCAA‑CCAAAAAAAATCGTAGATGTTTTGCGGGGAGGGGTACTTTGGCGCCGGCGTCTGACCGTT‑‑GCTCAGGGGAAAA‑C‑CGCCGACCAAATGGAAAAACCCG‑AGA‑AAAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCtt                                                                                                                                        <  1:307261/250‑1 (MQ=255)
                                                                                                                  aaaaaaaccccacgggCAGGTAAAAGAAGATCCCTTC‑CCCGATCAC‑CCAAAAAAACCCGTAGTTCGTTGGCAGGAAGGGGTACATTCGCGCCAGCGTCTGCCGGTT‑‑TCGAAGGGGCAAA‑C‑CGCCAATCAAATTGCAAAAACCT‑AGG‑AAAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCttttttt                                                                                                                                   <  1:263416/237‑1 (MQ=255)
                                                                                                                        tACCGACCGGCAGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGc                                                                                                                             >  2:249732/1‑251 (MQ=255)
                                                                                                                                  gaGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACAAAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGCCCGTT‑‑GCTCAGGGAAAAA‑T‑CGCCGATCACATTGCAAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCaa                                                                                                                   <  1:416988/250‑1 (MQ=255)
                                                                                                                                   aGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCACTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAAc                                                                                                                  >  1:449816/1‑251 (MQ=255)
                                                                                                                                        gaagaagaTGGCTTC‑ACCGAAAAA‑ACAAAAAAAAGCGTAGAGTGTTTGAAGGGAAGGGTACATTTGCGCGAGCTTCTGACGGTT‑‑GCTCAGGGGAAAA‑T‑CGCCAATAAAATGGCAAAAACCA‑AGA‑AAAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTgcg                                                                                                             <  1:389082/251‑1 (MQ=255)
                                                                                                                                         aaaaagtgcGCTCC‑ACCGACCAA‑CCACAAAAAATCGTAGGTCGTTTGCGGGGAAGGATACATTGGCGCCAGCGTCTGGCGTTT‑‑GCTCA‑GGGAAAA‑TTCCCCGATCACATTGCGAAAACCG‑AGG‑ATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTgcgc                                                                                                            <  1:441920/243‑1 (MQ=255)
                                                                                                                                               aTCGCTTC‑ACCGATCAA‑CCACAGAAAATCGTAGATCGATTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GATCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACATGCTCTTTATCATCCTGCGCCTGTcc                                                                                                      >  1:556816/1‑251 (MQ=255)
                                                                                                                                                     cc‑ACCGAGCAA‑CCACAGAAAATCGTAGTTCTTTTGGGGGGAGGGGTACATGTGCGCCAGCGTCTGCCCGTT‑‑GCTCATGGGCAAA‑C‑CGCCGATCACATTGCAAAAACCG‑GGG‑ATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGa                                                                                                <  1:459384/250‑1 (MQ=255)
                                                                                                                                                              cAA‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATc                                                                                         >  2:25364/1‑250 (MQ=255)
                                                                                                                                                                a‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGt                                                                                      >  1:556137/1‑251 (MQ=255)
                                                                                                                                                                     caGAAAATCGTAGATTGTTTGCACGGAAAGGACCATTTGCGGCAGCGTAAGAACGTC‑‑GCACATGGCCAAA‑T‑CCCCGAACACATTGCGAAAACCG‑AAG‑ATCAAA‑CCGCCGCTAATCACCGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATTAGAGCTCAACTGCTCATGCCATTTCATATTAAGCCGCGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTgttt                                                                                  >  2:339987/1‑248 (MQ=255)
                                                                                                                                                                      aaaaattcttagatctttggaggggAGGGGTACATTGGCGCCGGCGTCTGCCCGTT‑‑GCCAAGGGGAAAA‑C‑CGCCGACCACATTGCAAAAACCG‑AGG‑ACCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTg                                                                                 <  1:249732/229‑1 (MQ=255)
                                                                                                                                                                          aaTCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGAt                                                                               >  2:455213/1‑249 (MQ=255)
                                                                                                                                                                          aaTCGTAGATCGCTTGAAGGGAAGGGTACATTTGCGCAAGAGTCTGACCGTC‑‑GCTACTGGGAAAA‑T‑CGCAGATAACATTGTGAAAACCC‑AGG‑ATCAAA‑CCGCCGATAATCAACGCGGGCAGAAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGACATTTCATATTATGTCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGCCGCCATAAGTGCTTGAt                                                                               >  1:290044/1‑249 (MQ=255)
                                                                                                                                                                             cGTAGATCTTTTGCGGGGAAGGTTACATTGGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAC‑C‑CGCGGATCAAATGGCAAAACCCG‑AGA‑ATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTa                                                                           <  1:403070/250‑1 (MQ=255)
                                                                                                                                                                              gTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAAt                                                                         >  1:568006/1‑251 (MQ=255)
                                                                                                                                                                                            gggAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTg                                                           >  2:557181/1‑251 (MQ=255)
                                                                                                                                                                                               aggggtccgttacccccAGCTTCTGACCGTT‑‑GCT‑ATGGGCAAAAT‑CGCCGACCA‑AATGGAAAAACCG‑AGGCAACAAC‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTTACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCa                                                        <  2:82610/236‑1 (MQ=255)
                                                                                                                                                                                                     ataTTTCCGCCAGCGTCTGACTGTT‑‑GCTAATGGGAAAA‑T‑CCCCGAACACATGGCAAAACCCG‑AGG‑ACAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATcac                                                   <  1:72984/248‑1 (MQ=255)
                                                                                                                                                                                                         ttGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCt                                              >  2:14571/1‑251 (MQ=255)
                                                                                                                                                                                                               cAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAAtt                                        >  2:394793/1‑251 (MQ=255)
                                                                                                                                                                                                                       gACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAACCCCG‑AGG‑CTCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGc                                 >  1:69306/1‑250 (MQ=255)
                                                                                                                                                                                                                        aCCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCAt                               >  1:23411/1‑251 (MQ=255)
                                                                                                                                                                                                                           gTT‑‑GCTCATGGGAAAA‑T‑CCCCGAACAAAGGGCAAAAACCG‑AGA‑AAAAAC‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTgg                            <  1:438778/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggggaaaaa‑T‑CCCCGATCAAATTGCAAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACa                     <  1:514592/245‑1 (MQ=255)
                                                                                                                                                                                                                                       ggCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGACATTGGTGAAACAGcc                  >  1:542001/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aaccGCCGATCACATTGCGAAACCGG‑AGA‑AAAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTc              <  1:557181/248‑1 (MQ=255)
                                                                                                                                                                                                                                                         aaTTGCAAAACCCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGt  <  1:25364/250‑1 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
CGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGG‑CGTGCCGCCCATTTTTTT‑CACC‑GCTGACCA‑GCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTC‑ACCGATCAA‑CCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTT‑‑GCTCATGGGCAAA‑T‑CGCCGATCACATTGCGAAAACCG‑AGG‑ATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGT  >  NC_000913/4465542‑4466031

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