Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC44.0% 110.5 / 38.7 91K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (36/15);  new base C (0/40);  total (36/55)
Fisher's exact test for biased strand distribution p-value = 1.25e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357730‑358023
                                                                                                                                                               |                                                                                                                                                   
aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTAtttt‑t                                                                                                                                                     <  4:1441627/151‑1 (MQ=255)
   cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑cc                                                                                                                                                  <  1:806650/150‑1 (MQ=255)
   cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑cc                                                                                                                                                  <  1:1301748/150‑1 (MQ=255)
   cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑cc                                                                                                                                                  >  1:225529/1‑150 (MQ=255)
    gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑ccgc                                                                                                                                                >  4:1503242/1‑151 (MQ=255)
         ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGcaa                                                                                                                                           <  4:1694009/151‑1 (MQ=255)
         ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGcaa                                                                                                                                           <  4:668592/151‑1 (MQ=255)
         ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGca                                                                                                                                            >  1:156500/1‑150 (MQ=255)
          cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGcaac                                                                                                                                          >  4:1247700/1‑151 (MQ=255)
              tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACaaa                                                                                                                                       >  1:1323063/1‑150 (MQ=255)
               aaGCGTGGCGCATCAGGCAATTCCAGCCGCAGCCCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑c                                                                                                                                    <  1:168558/151‑1 (MQ=255)
               aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACaaaa                                                                                                                                      <  4:1272895/150‑1 (MQ=255)
               aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑c                                                                                                                                    <  4:1476592/151‑1 (MQ=255)
               aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑c                                                                                                                                    <  2:168991/151‑1 (MQ=255)
                aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑cc                                                                                                                                   >  4:1200541/1‑151 (MQ=255)
                   gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑ccgcc                                                                                                                                >  1:1822968/1‑151 (MQ=255)
                    tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑ccgccg                                                                                                                               >  1:1155910/1‑151 (MQ=255)
                    tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑ccgccg                                                                                                                               >  1:1170167/1‑151 (MQ=255)
                       cgcATCAGGCATTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTC‑CCGCC‑AACAA‑CCGCGCATTTATCCAC‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑ccgccg                                                                                                                               <  2:428175/148‑1 (MQ=255)
                       cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTtc                                                                                                                            >  4:128819/1‑151 (MQ=255)
                            cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTcc                                                                                                                       >  4:1649681/1‑151 (MQ=255)
                            cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTcc                                                                                                                       >  4:1797908/1‑151 (MQ=255)
                                    tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTgg                                                                                                               >  1:1387964/1‑151 (MQ=255)
                                    tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTgg                                                                                                               >  1:516023/1‑151 (MQ=255)
                                    cccAGCCGCGCCCCTGTGTCCGCGGCTCCCGTGATTCTTTC‑CCCCC‑AACAC‑CCGCGCTTTTACCCAC‑CCCCATGAAAAAAAA‑TCCCTC‑CGCGGC‑GGCTGTCGCCCA‑TTTTTTTT‑C‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCagta                                                                                                                   <  1:655595/146‑1 (MQ=255)
                                      cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCg                                                                                                             >  4:1368420/1‑151 (MQ=255)
                                      cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCg                                                                                                             >  4:215123/1‑151 (MQ=255)
                                      cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCg                                                                                                             >  1:872455/1‑151 (MQ=255)
                                         cggCAGACCTGTGTCGGCGCCCCCGGTAATTCTTTC‑CCACC‑AAAAC‑CCGCTCTTTTACCCAC‑CCCAAGGTAAAAAAA‑TGCCTC‑CGCGGC‑GGCTTTCCGCCA‑GTTATTTT‑C‑CCCCCAAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCg                                                                                                             <  2:332101/146‑1 (MQ=255)
                                         ccGAAGCCCTTTGCCCGCGGCTCCCGTAATCCTTTC‑CCCCC‑AACAC‑CCGCATTTTTACCCAA‑CCCAATGTCCAAAAG‑TCCCTT‑GGGGCG‑GGCTGTCCGCCA‑GTTTTTTC‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCg                                                                                                             <  4:298109/148‑1 (MQ=255)
                                          cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                            >  2:213819/1‑148 (MQ=255)
                                           gCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGTCgg                                                                                                            >  3:155019/1‑147 (MQ=255)
                                           gCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                         >  1:474329/1‑149 (MQ=255)
                                           gCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                         >  4:1740870/1‑149 (MQ=255)
                                            cAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                          >  2:445094/1‑148 (MQ=255)
                                                ccTGTGTCAGCGGCTCCCGTGATTAATTC‑CCGCC‑AACAC‑CCGCCCATTTACCAAA‑CCCCATGTGAAAAAA‑TGCTTC‑CGGGGC‑GGTTTTCGCCCA‑GCTTTTTT‑C‑CCGCCAAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgag‑a                                                                                                  <  2:540625/151‑1 (MQ=255)
                                                    tgtCGGGGGCTCCCGGGTTTC‑TTT‑CCCCC‑AAAAACCCGCGAATTTACCAAA‑CGCCTTGTGAAAAAG‑GGCTTG‑CGGGGG‑GGTTGTCGGCCG‑GTTTTTTT‑C‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aaa                                                                                                <  4:985611/149‑1 (MQ=255)
                                                    tgtCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGagagga                                                                                                  >  4:1265236/1‑146 (MQ=255)
                                                       caccgcctcccttGATTCATTC‑CCCCA‑AACAC‑CCGCACTTTTCCCAAC‑CCCCATGTAAAAAAA‑TCCCTC‑CGGGGC‑GGTTTTCGCCAA‑GTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aaa                                                                                                <  1:1841330/134‑1 (MQ=255)
                                                        aGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGG‑aaaaaaatt                                                                                          >  1:365808/1‑145 (MQ=255)
                                                           gcctcccGTAATTCTTTC‑CCCCC‑AACAC‑CCGCGCTTTTCCCCAA‑CCCCAGGTACAAAAA‑GCCTTT‑CGCGGC‑GGCTTTCTCCCA‑GTTATTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTa                                                                                        <  2:120092/145‑1 (MQ=255)
                                                           gcctcccGTAATTCATTC‑CCGCC‑AACAC‑CCGCGCATTTACCCAA‑CGCCTTGTAAAAAAT‑TCCCTT‑CGGGGG‑GTCTGTCGCCCA‑GTTTTTTT‑C‑CCGCCCCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                       <  4:30704/131‑1 (MQ=255)
                                                            gCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑gaaaaaatttt                                                                                        >  4:63381/1‑141 (MQ=255)
                                                             cTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAG‑AAAAAACttt                                                                                         >  4:737025/1‑147 (MQ=255)
                                                              tACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGTGGAAGCG‑aaaaaaaatttaaac                                                                                    >  4:34590/1‑139 (MQ=25)
                                                               aCCGTGTTTCATTC‑CCCCC‑ACCAC‑CCGCTATTTTACCCAC‑CGCCTTGTAAAAAAA‑TGCCTC‑CGGGGC‑GGCTGTCCGCCA‑GCTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTa                                                                                        <  4:1596665/146‑1 (MQ=255)
                                                                ccGTGTTTCTTTC‑CCGCC‑AAC‑A‑CCGGCCTTTTTCCCAACCGCCTTTTG‑AAAAA‑CCCCTTGCGGGGG‑GGTTGCCTGCCG‑GTTATTTT‑T‑CCCCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                    <  4:1862722/149‑1 (MQ=255)
                                                                 ccgtaatCATTT‑CCCCCAAACAC‑CCGCGTTTTTATCCAA‑CCCCTTTTAAAAAAG‑TGCCTT‑GGGGGC‑GGCGGTCGGCCG‑CCTATTTT‑T‑CCCCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAg                                                                                             <  4:1203910/134‑1 (MQ=255)
                                                                    gATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑aaaaaaattttaatcc                                                                                   >  1:1659499/1‑133 (MQ=25)
                                                                     aatCATTTCCCGCC‑AACAC‑CCCCGAATTTCCCCAC‑CCCCATGTAAAAAAG‑TGCCTT‑CGCGGC‑GGCTTTCTCCCA‑GTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                       <  4:376781/125‑1 (MQ=255)
                                                                      ttCTTTC‑CCGCC‑AACAC‑CCGCGCATTTACCCAC‑CGCCATGTGAAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCGGCCA‑GTTTTTTT‑C‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg                                                                             <  2:109690/150‑1 (MQ=255)
                                                                       tttatccccGCA‑AACAC‑CCGCGTTTTTACCAAA‑CCCCATTGGAAAAAT‑TGCCTC‑CGGGGG‑GGTTGTCTGCAG‑GTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAt                                                                              <  1:1148803/142‑1 (MQ=255)
                                                                        cATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑GAAAAAATTTAAACCCCTcc                                                                               >  2:397371/1‑145 (MQ=25)
                                                                           ttcccGCA‑AACAC‑CCGCGCATTTACCCAC‑CGCCTTGTGAAAAAA‑TCCCTT‑CGGGGC‑GGCTTTCTGCCA‑GTTTTTTT‑C‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgttgt                                                                         <  4:1226735/147‑1 (MQ=255)
                                                                           tccccccc‑AACCC‑CCGCGCTTTTACCAAC‑CGCCATGGCAAAAAA‑TCCCTT‑CGGGGG‑GGCTGTCTCCCA‑TCTTTTTT‑C‑CCCCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg                                                                             <  1:958497/140‑1 (MQ=255)
                                                                              cccccaaacac‑cCGCGCTTTTACCCAA‑CGCAATGTAAAAAAA‑TCCCTC‑CGCGGG‑GGCTGTCTCCCA‑GTTTTTTT‑C‑CCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgt                                                                            <  1:997894/139‑1 (MQ=255)
                                                                                    aacac‑cCGCGCTTTTACCCAC‑CG‑CTTGGGAAAAAAGTCCCTC‑CGGGGC‑GGTTTTC‑GCCA‑GTTTTTTTCC‑CCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAta                                                                <  1:19699/151‑2 (MQ=255)
                                                                                     accacccccgcTTTTATCCAC‑CGCCATGTGAAAAAA‑TCCCTG‑CGGGGC‑GGTTGTCTGCCA‑TTTTTTTT‑T‑CCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                   <  1:825838/140‑1 (MQ=255)
                                                                                        a‑cCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAAAATTAAACCCGCCATGTTTTGTggggatagc                                                              >  2:171346/1‑142 (MQ=25)
                                                                                          ccGCGTATTTACCCAA‑CCCCATGTAAAAAAA‑TCCCTG‑CGCGGC‑GGCTGTCGCCCA‑TTTTTTTT‑C‑CCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                            <  1:1247173/150‑1 (MQ=255)
                                                                                          ccGCGCTTTTCTCCAC‑CCCCATGTAAAAAAA‑GCCCTC‑CGCGGC‑GGCTGTCTGCCA‑CCTTTTTC‑C‑CCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                            <  4:326564/150‑1 (MQ=255)
                                                                                          ccGCGCATTTACCCAA‑CGCCTTGTGCAAAAA‑TGCCTC‑CGCGGC‑GGCTGTCTGCCA‑GCTTTTTT‑C‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGCGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                             <  2:213984/149‑1 (MQ=255)
                                                                                           cgcgTATTTACCCAC‑CTCCTTGTGAAAAAA‑TCCCTC‑CGGGGC‑GGCTTCTTGCCA‑GTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg‑ct                                                         <  4:727028/151‑1 (MQ=255)
                                                                                           cgcgCTTTTACCCAC‑CGCCTTGTGAAAAAT‑TCCCTG‑GGGGGG‑GGCTGTCTGCCA‑CTTTTTTT‑C‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg‑ct                                                         <  4:1273125/151‑1 (MQ=255)
                                                                                              gcttttACCAAA‑CGCAATGTGAAAAAA‑TCCTTG‑CGGGGG‑GGTTGTCTGCCA‑GCTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctg‑ctg                                                        <  3:237996/146‑1 (MQ=255)
                                                                                                  ttATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑AAAAAAAATTAAAACCCACACGTTTTGTGGCGATCgcag‑c                                                          >  4:141731/1‑143 (MQ=25)
                                                                                                    acccAC‑CGCCATGTGAAAAAA‑TGCCTC‑CGCGGC‑GGCTGTCTGCCA‑GTTTTTTT‑C‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACaa                                                 <  4:65450/148‑1 (MQ=255)
                                                                                                      ccAC‑CGCCTTGCGAAAAAA‑TGCCTC‑CGCGGG‑GGTTGTCGGCCA‑TCTTTTTT‑C‑CCCCCACAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACa                                                  <  4:1147551/147‑1 (MQ=255)
                                                                                                       ccaccgcattgtaaaaaat‑TCCCTC‑CGGGGC‑GGCTGTC‑GCCATGTTTTTTT‑C‑CCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgt                                              <  1:1235717/137‑1 (MQ=255)
                                                                                                       caaccccaatggcaaaaaa‑TCCTTC‑CGGGGG‑GGCTGTCTGCCA‑TTTTTTTT‑C‑CCCCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgt                                              <  1:381981/137‑1 (MQ=255)
                                                                                                       cAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAACTTTAAACCC‑CCATTTTGTGGGGGGATCGGTGCCTGGTAaaaa                                                 >  4:1653051/1‑148 (MQ=21)
                                                                                                         accGCAATGTAAAAAAG‑GGCCTC‑CGGGGC‑GGCTGTCGCCCG‑GATTTTTC‑C‑CCCCGGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgtg                                             <  1:865763/148‑1 (MQ=255)
                                                                                                         aacgccttggaaaaaaa‑TGCCTC‑CGGGGC‑GGCTGTCTGCCA‑TTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgtg                                             <  1:1242945/138‑1 (MQ=255)
                                                                                                          cccccATGTAAAAAAA‑CGCCTC‑GGGGCC‑GTCTGCCTCCCA‑CTTATTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtg                                               <  3:478169/144‑1 (MQ=255)
                                                                                                           ccctaggaaaaaaaa‑TGCCTC‑CGGGGG‑GGCTGTCGGCCA‑TTTTTTTT‑C‑CCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgtg                                             <  1:443771/138‑1 (MQ=255)
                                                                                                             ccttgtgAAAAAT‑TGCCTC‑CGCGGC‑GGTTGTCGCCCA‑GTTATTTT‑C‑CCCCCTCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAAc                                                   <  4:1673842/137‑1 (MQ=255)
                                                                                                             ccagggaaaaaaa‑CGCCTC‑CGCGGC‑GGCTGTCGGCCG‑TTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGaaa                                         <  1:262835/142‑1 (MQ=255)
                                                                                                                     aaaaa‑aGGCTT‑CGGGGGGGGCGGTGTGCAA‑GTTTTTTT‑C‑CCGCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGa                                     <  3:241042/147‑1 (MQ=255)
                                                                                                                     aaaaa‑TCCCTG‑GGGGGG‑GGCTGT‑TGGCA‑ACTTTTTT‑TCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTg                                 <  4:1281548/150‑1 (MQ=255)
                                                                                                                     aaaaa‑TCCCTC‑CGGGGG‑GCCTGTCGCCCA‑GTTTTTTT‑C‑CCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAg                                   <  1:1715290/148‑1 (MQ=255)
                                                                                                                      aaaa‑TCCCTC‑CGCGGC‑GGCTGTCTGCAA‑GCTTTTTT‑C‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTgcg                               <  4:80577/151‑1 (MQ=255)
                                                                                                                        aaatggcctggggggc‑ggcTGCCTGCCA‑GTTTTTTT‑T‑CCCCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTgcg                               <  4:1765605/138‑1 (MQ=255)
                                                                                                                             ccTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAAAATTAAAACCCCCCTTTTGTGGTGGAGCCGCTG‑CTGGGTAAAATGTGGGAAGGGGAgggggggg                            >  1:1742480/1‑141 (MQ=255)
                                                                                                                               tt‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAGAATTTAACACCCCCTGTTGTGTGGGGATCGCCG‑CTGGGGAAAAAGTGGGAAggggggggggggg                            >  4:1548243/1‑136 (MQ=255)
                                                                                                                                   ggggc‑ggcTGTCTCCCA‑GTTTTTTT‑C‑CCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                         <  3:298147/144‑1 (MQ=255)
                                                                                                                                     ggc‑ggcTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAAATTTAAACCC‑CCCTGTTGTTTGGGAATCGCCG‑CCGGTGACAAAGAGGGAAGGGGAGGGGGGGGGATCaacccc                  >  4:573556/1‑146 (MQ=255)
                                                                                                                                      ggggcctttCTGCCA‑TTTTTTTC‑C‑CCGCCGCAAAAAGA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                        <  1:1754900/137‑1 (MQ=255)
                                                                                                                                       gcggttttcgcccA‑GCTATTTT‑C‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                          <  4:1920676/131‑1 (MQ=255)
                                                                                                                                                ggCCA‑GTTTTTTT‑T‑CCGCAGAAAAAAAT‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  4:573594/145‑1 (MQ=255)
                                                                                                                                                cccccttttttatt‑t‑CCGCCCCAAAAAAT‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc            <  2:159273/136‑1 (MQ=255)
                                                                                                                                                  ccA‑GTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  4:21090/144‑1 (MQ=255)
                                                                                                                                                  ccA‑GTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa      <  4:1827703/151‑1 (MQ=255)
                                                                                                                                                   ccatttttttt‑c‑CCCCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg         <  3:415219/141‑1 (MQ=255)
                                                                                                                                                     agttttttt‑c‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGAGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa     <  4:365935/145‑1 (MQ=255)
                                                                                                                                                         aTTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAG‑AAAAAAATTTTAACACCCCCTTTTGTGTGGGAATCGCCG‑CTGGGAAAAAAGTGGGAAAGGGAGGGGGGGGAAAAAACCCCCGAACCGCGCGAAttt  >  4:328460/1‑148 (MQ=255)
                                                                                                                                                           ttt‑t‑CCCGCGCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg  <  4:1184411/148‑1 (MQ=255)
                                                                                                                                                               |                                                                                                                                                   
ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357730‑358023

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