Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA47.2% 67.8 / 68.4 91L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (33/15);  new base A (0/43);  total (33/58)
Fisher's exact test for biased strand distribution p-value = 2.07e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAAT  >  NC_000913/357740‑358036
                                                                                                                                                                |                                                                                                                                                      
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                                                                                                                                                  <  2:168991/151‑1 (MQ=255)
     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑c                                                                                                                                                  <  4:1476592/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                                                                                                                           >  1:872455/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                                                                                                                           >  4:1368420/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                                                                                                                       >  4:1740870/1‑149 (MQ=255)
                                 gCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                                       >  1:474329/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‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAga‑g‑a                                                                                                               <  2:540625/151‑1 (MQ=255)
                                          tgtCGGGGGCTCCCGGGTTTC‑TTT‑CCCCC‑AAAAACCCGCGAATTTACCAAA‑CGCCTTGTGAAAAAG‑GGCTTG‑CGGGGG‑GGTTGTCGGCCG‑GTTTTTTT‑C‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑aaa                                                                                                             <  4:985611/149‑1 (MQ=255)
                                          tgtCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGaga‑gga                                                                                                               >  4:1265236/1‑146 (MQ=255)
                                             caccgcctcccttGATTCATTC‑CCCCA‑AACAC‑CCGCACTTTTCCCAAC‑CCCCATGTAAAAAAA‑TCCCTC‑CGGGGC‑GGTTTTCGCCAA‑GTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑aaa                                                                                                             <  1:1841330/134‑1 (MQ=255)
                                              aGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGG‑G‑aaaaaaatt                                                                                                       >  1:365808/1‑145 (MQ=255)
                                                 gcctcccGTAATTCTTTC‑CCCCC‑AACAC‑CCGCGCTTTTCCCCAA‑CCCCAGGTACAAAAA‑GCCTTT‑CGCGGC‑GGCTTTCTCCCA‑GTTATTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑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‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑gaaaaaatttt                                                                                                     >  4:63381/1‑141 (MQ=255)
                                                   cTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGA‑G‑AAAAAACttt                                                                                                      >  4:737025/1‑147 (MQ=255)
                                                    tACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGTGGAAGC‑G‑aaaaaaaatttaaac                                                                                                 >  4:34590/1‑139 (MQ=25)
                                                     aCCGTGTTTCATTC‑CCCCC‑ACCAC‑CCGCTATTTTACCCAC‑CGCCTTGTAAAAAAA‑TGCCTC‑CGGGGC‑GGCTGTCCGCCA‑GCTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTa                                                                                                     <  4:1596665/146‑1 (MQ=255)
                                                      ccGTGTTTCTTTC‑CCGCC‑AAC‑A‑CCGGCCTTTTTCCCAACCGCCTTTTG‑AAAAA‑CCCCTTGCGGGGG‑GGTTGCCTGCCG‑GTTATTTT‑T‑CCCCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGcc                                                                                                 <  4:1862722/149‑1 (MQ=255)
                                                       ccgtaatCATTT‑CCCCCAAACAC‑CCGCGTTTTTATCCAA‑CCCCTTTTAAAAAAG‑TGCCTT‑GGGGGC‑GGCGGTCGGCCG‑CCTATTTT‑T‑CCCCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAg                                                                                                          <  4:1203910/134‑1 (MQ=255)
                                                          gATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑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‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtg                                                                                          <  2:109690/150‑1 (MQ=255)
                                                             tttatccccGCA‑AACAC‑CCGCGTTTTTACCAAA‑CCCCATTGGAAAAAT‑TGCCTC‑CGGGGG‑GGTTGTCTGCAG‑GTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAt                                                                                           <  1:1148803/142‑1 (MQ=255)
                                                              cATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑GAAAAAATTTAAACCCCTcc                                                                                            >  2:397371/1‑145 (MQ=25)
                                                                 ttcccGCA‑AACAC‑CCGCGCATTTACCCAC‑CGCCTTGTGAAAAAA‑TCCCTT‑CGGGGC‑GGCTTTCTGCCA‑GTTTTTTT‑C‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgttgt                                                                                      <  4:1226735/147‑1 (MQ=255)
                                                                 tccccccc‑AACCC‑CCGCGCTTTTACCAAC‑CGCCATGGCAAAAAA‑TCCCTT‑CGGGGG‑GGCTGTCTCCCA‑TCTTTTTT‑C‑CCCCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtg                                                                                          <  1:958497/140‑1 (MQ=255)
                                                                    cccccaaacac‑cCGCGCTTTTACCCAA‑CGCAATGTAAAAAAA‑TCCCTC‑CGCGGG‑GGCTGTCTCCCA‑GTTTTTTT‑C‑CCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgt                                                                                         <  1:997894/139‑1 (MQ=255)
                                                                          aacac‑cCGCGCTTTTACCCAC‑CG‑CTTGGGAAAAAAGTCCCTC‑CGGGGC‑GGTTTTC‑GCCA‑GTTTTTTTCC‑CCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAta                                                                             <  1:19699/151‑2 (MQ=255)
                                                                           accacccccgcTTTTATCCAC‑CGCCATGTGAAAAAA‑TCCCTG‑CGGGGC‑GGTTGTCTGCCA‑TTTTTTTT‑T‑CCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                                <  1:825838/140‑1 (MQ=255)
                                                                              a‑cCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑AAAAAAAATTAAACCCGCCATGTTTTGTggggatagc                                                                           >  2:171346/1‑142 (MQ=25)
                                                                                ccGCGTATTTACCCAA‑CCCCATGTAAAAAAA‑TCCCTG‑CGCGGC‑GGCTGTCGCCCA‑TTTTTTTT‑C‑CCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                                         <  1:1247173/150‑1 (MQ=255)
                                                                                ccGCGCTTTTCTCCAC‑CCCCATGTAAAAAAA‑GCCCTC‑CGCGGC‑GGCTGTCTGCCA‑CCTTTTTC‑C‑CCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                                         <  4:326564/150‑1 (MQ=255)
                                                                                ccGCGCATTTACCCAA‑CGCCTTGTGCAAAAA‑TGCCTC‑CGCGGC‑GGCTGTCTGCCA‑GCTTTTTT‑C‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGCGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                                          <  2:213984/149‑1 (MQ=255)
                                                                                 cgcgTATTTACCCAC‑CTCCTTGTGAAAAAA‑TCCCTC‑CGGGGC‑GGCTTCTTGCCA‑GTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg‑ct                                                                      <  4:727028/151‑1 (MQ=255)
                                                                                 cgcgCTTTTACCCAC‑CGCCTTGTGAAAAAT‑TCCCTG‑GGGGGG‑GGCTGTCTGCCA‑CTTTTTTT‑C‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg‑ct                                                                      <  4:1273125/151‑1 (MQ=255)
                                                                                    gcttttACCAAA‑CGCAATGTGAAAAAA‑TCCTTG‑CGGGGG‑GGTTGTCTGCCA‑GCTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctg‑ctg                                                                     <  3:237996/146‑1 (MQ=255)
                                                                                        ttATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑G‑AAAAAAAATTAAAACCCACACGTTTTGTGGCGATCgcag‑c                                                                       >  4:141731/1‑143 (MQ=25)
                                                                                          acccAC‑CGCCATGTGAAAAAA‑TGCCTC‑CGCGGC‑GGCTGTCTGCCA‑GTTTTTTT‑C‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACaa                                                              <  4:65450/148‑1 (MQ=255)
                                                                                            ccAC‑CGCCTTGCGAAAAAA‑TGCCTC‑CGCGGG‑GGTTGTCGGCCA‑TCTTTTTT‑C‑CCCCCACAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACa                                                               <  4:1147551/147‑1 (MQ=255)
                                                                                             ccaccgcattgtaaaaaat‑TCCCTC‑CGGGGC‑GGCTGTC‑GCCATGTTTTTTT‑C‑CCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgt                                                           <  1:1235717/137‑1 (MQ=255)
                                                                                             caaccccaatggcaaaaaa‑TCCTTC‑CGGGGG‑GGCTGTCTGCCA‑TTTTTTTT‑C‑CCCCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgt                                                           <  1:381981/137‑1 (MQ=255)
                                                                                             cAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAAACTTTAAACCC‑CCATTTTGTGGGGGGATCGGTGCCTGGTAaaaa                                                              >  4:1653051/1‑148 (MQ=21)
                                                                                               accGCAATGTAAAAAAG‑GGCCTC‑CGGGGC‑GGCTGTCGCCCG‑GATTTTTC‑C‑CCCCGGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgtg                                                          <  1:865763/148‑1 (MQ=255)
                                                                                               aacgccttggaaaaaaa‑TGCCTC‑CGGGGC‑GGCTGTCTGCCA‑TTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgtg                                                          <  1:1242945/138‑1 (MQ=255)
                                                                                                cccccATGTAAAAAAA‑CGCCTC‑GGGGCC‑GTCTGCCTCCCA‑CTTATTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtg                                                            <  3:478169/144‑1 (MQ=255)
                                                                                                 ccctaggaaaaaaaa‑TGCCTC‑CGGGGG‑GGCTGTCGGCCA‑TTTTTTTT‑C‑CCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgtg                                                          <  1:443771/138‑1 (MQ=255)
                                                                                                   ccttgtgAAAAAT‑TGCCTC‑CGCGGC‑GGTTGTCGCCCA‑GTTATTTT‑C‑CCCCCTCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAAc                                                                <  4:1673842/137‑1 (MQ=255)
                                                                                                   ccagggaaaaaaa‑CGCCTC‑CGCGGC‑GGCTGTCGGCCG‑TTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGaaa                                                      <  1:262835/142‑1 (MQ=255)
                                                                                                           aaaaa‑aGGCTT‑CGGGGGGGGCGGTGTGCAA‑GTTTTTTT‑C‑CCGCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGa                                                  <  3:241042/147‑1 (MQ=255)
                                                                                                           aaaaa‑TCCCTG‑GGGGGG‑GGCTGT‑TGGCA‑ACTTTTTT‑TCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTg                                              <  4:1281548/150‑1 (MQ=255)
                                                                                                           aaaaa‑TCCCTC‑CGGGGG‑GCCTGTCGCCCA‑GTTTTTTT‑C‑CCGCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAg                                                <  1:1715290/148‑1 (MQ=255)
                                                                                                            aaaa‑TCCCTC‑CGCGGC‑GGCTGTCTGCAA‑GCTTTTTT‑C‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTgcg                                            <  4:80577/151‑1 (MQ=255)
                                                                                                              aaatggcctggggggc‑ggcTGCCTGCCA‑GTTTTTTT‑T‑CCCCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTgcg                                            <  4:1765605/138‑1 (MQ=255)
                                                                                                                   ccTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAAAATTAAAACCCCCCTTTTGTGGTGGAGCCGCTG‑CTGGGTAAAATGTGGGAAGGGGAgggggggg                                         >  1:1742480/1‑141 (MQ=255)
                                                                                                                     tt‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAGAATTTAACACCCCCTGTTGTGTGGGGATCGCCG‑CTGGGGAAAAAGTGGGAAggggggggggggg                                         >  4:1548243/1‑136 (MQ=255)
                                                                                                                         ggggc‑ggcTGTCTCCCA‑GTTTTTTT‑C‑CCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                                      <  3:298147/144‑1 (MQ=255)
                                                                                                                           ggc‑ggcTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAAAATTTAAACCC‑CCCTGTTGTTTGGGAATCGCCG‑CCGGTGACAAAGAGGGAAGGGGAGGGGGGGGGATCaacccc                               >  4:573556/1‑146 (MQ=255)
                                                                                                                            ggggcctttCTGCCA‑TTTTTTTC‑C‑CCGCCGCAAAAAGA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                                     <  1:1754900/137‑1 (MQ=255)
                                                                                                                             gcggttttcgcccA‑GCTATTTT‑C‑CCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                                       <  4:1920676/131‑1 (MQ=255)
                                                                                                                                      ggCCA‑GTTTTTTT‑T‑CCGCAGAAAAAAAT‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                          <  4:573594/145‑1 (MQ=255)
                                                                                                                                      cccccttttttatt‑t‑CCGCCCCAAAAAAT‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                         <  2:159273/136‑1 (MQ=255)
                                                                                                                                        ccA‑GTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                          <  4:21090/144‑1 (MQ=255)
                                                                                                                                        ccA‑GTTTTTTT‑C‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa                   <  4:1827703/151‑1 (MQ=255)
                                                                                                                                         ccatttttttt‑c‑CCCCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg                      <  3:415219/141‑1 (MQ=255)
                                                                                                                                           agttttttt‑c‑CCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGAGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa                  <  4:365935/145‑1 (MQ=255)
                                                                                                                                               aTTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAA‑G‑AAAAAAATTTTAACACCCCCTTTTGTGTGGGAATCGCCG‑CTGGGAAAAAAGTGGGAAAGGGAGGGGGGGGAAAAAACCCCCGAACCGCGCGAAttt               >  4:328460/1‑148 (MQ=255)
                                                                                                                                                 ttt‑t‑CCCGCGCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg               <  4:1184411/148‑1 (MQ=255)
                                                                                                                                                    ttcccgccgAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt           <  4:420549/147‑1 (MQ=255)
                                                                                                                                                      cccgccaCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg         <  1:1846678/149‑1 (MQ=255)
                                                                                                                                                       ccgccgCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc          <  4:1524482/148‑1 (MQ=255)
                                                                                                                                                          ccgcAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGG‑AAAAAAATTTAACCCCG‑CCTTTTTTGTGGGGAACGCCG‑CTGGGGAAAAAGTGGGAAGGGGAGGGGGGGGGGTCAACCCCCGCACCGCGCCCATTGAGTTTGcc       >  1:39100/1‑148 (MQ=255)
                                                                                                                                                              accaaAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAAt  <  4:29813/147‑1 (MQ=255)
                                                                                                                                                                caaAA‑CCGCCGTTCTCTCCAGTAGTGGCGGAGCGAAA‑G‑AAAAAGATTTTAAACCCCCATGTTGTGTGGGGATTGGTC‑GTGGGGAAAAAGTGGGAAGGGGAGGGGGGGGGATAAGTTCCCGAAACCCGCCAATTGAGTTTGtctc     >  4:106799/1‑144 (MQ=255)
                                                                                                                                                                |                                                                                                                                                      
CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAAT  >  NC_000913/357740‑358036

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