Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8750AT41.5% 136.3 / 47.7 94N272K (AAT→AAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (38/17);  new base T (0/39);  total (38/56)
Fisher's exact test for biased strand distribution p-value = 2.73e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357726‑358023
                                                                                                                                                             |                                                                                                                                                        
cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTa                                                                                                                                                          <  1:51514/150‑1 (MQ=255)
 tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTAtt                                                                                                                                                        >  1:65431/1‑151 (MQ=255)
   aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTAtttt                                                                                                                                                      >  3:159266/1‑151 (MQ=255)
   aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTAtttt                                                                                                                                                      >  4:1847398/1‑151 (MQ=255)
    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:1301748/150‑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: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:668592/151‑1 (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‑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‑CCGCCGCAACAAAAc                                                                                                                                    <  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‑CCGCCGCAACAAAAc                                                                                                                                    <  2:168991/151‑1 (MQ=255)
                   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAAc                                                                                                                                    <  4:1476592/151‑1 (MQ=255)
                    aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAAcc                                                                                                                                   >  4:1200541/1‑151 (MQ=255)
                       gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAAccgcc                                                                                                                                >  1:1822968/1‑151 (MQ=255)
                        tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAAccgccg                                                                                                                               >  1:1155910/1‑151 (MQ=255)
                        tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAAccgccg                                                                                                                               >  1:1170167/1‑151 (MQ=255)
                           cgcATCAGGCATTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTC‑CCGCC‑AACAA‑CCGCGCATTTATCCAC‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAAccgccg                                                                                                                               <  2:428175/148‑1 (MQ=255)
                           cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTtc                                                                                                                            >  4:128819/1‑151 (MQ=255)
                                cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                       >  4:1649681/1‑151 (MQ=255)
                                cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                       >  4:1797908/1‑151 (MQ=255)
                                        tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                               >  1:516023/1‑151 (MQ=255)
                                        tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                               >  1:1387964/1‑151 (MQ=255)
                                        cccAGCCGCGCCCCTGTGTCCGCGGCTCCCGTGATTCTTTC‑CCCCC‑AACAC‑CCGCGCTTTTACCCAC‑CCCCATGAAAAAAAA‑TCCCTC‑CGCGGC‑GGCTGTCGCCCA‑TTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                                   <  1:655595/146‑1 (MQ=255)
                                          cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                             >  1:872455/1‑151 (MQ=255)
                                          cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                             >  4:1368420/1‑151 (MQ=255)
                                          cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                             >  4:215123/1‑151 (MQ=255)
                                             cggCAGACCTGTGTCGGCGCCCCCGGTAATTCTTTC‑CCACC‑AAAAC‑CCGCTCTTTTACCCAC‑CCCAAGGTAAAAAAA‑TGCCTC‑CGCGGC‑GGCTTTCCGCCA‑GTTATTTT‑C‑CCCCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCg                                                                                                             <  2:332101/146‑1 (MQ=255)
                                             ccGAAGCCCTTTGCCCGCGGCTCCCGTAATCCTTTC‑CCCCC‑AACAC‑CCGCATTTTTACCCAA‑CCCAATGTCCAAAAG‑TCCCTT‑GGGGCG‑GGCTGTCCGCCA‑GTTTTTTC‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                             <  4:298109/148‑1 (MQ=255)
                                              cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                            >  2:213819/1‑148 (MQ=255)
                                               gCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGTCgg                                                                                                            >  3:155019/1‑147 (MQ=255)
                                               gCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                         >  4:1740870/1‑149 (MQ=255)
                                               gCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                         >  1:474329/1‑149 (MQ=255)
                                                cAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                          >  2:445094/1‑148 (MQ=255)
                                                    ccTGTGTCAGCGGCTCCCGTGATTAATTC‑CCGCC‑AACAC‑CCGCCCATTTACCAAA‑CCCCATGTGAAAAAA‑TGCTTC‑CGGGGC‑GGTTTTCGCCCA‑GCTTTTTT‑C‑CCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgag‑a                                                                                                  <  2:540625/151‑1 (MQ=255)
                                                        tgtCGGGGGCTCCCGGGTTTC‑TTT‑CCCCC‑AAAAACCCGCGAATTTACCAAA‑CGCCTTGTGAAAAAG‑GGCTTG‑CGGGGG‑GGTTGTCGGCCG‑GTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aaa                                                                                                <  4:985611/149‑1 (MQ=255)
                                                        tgtCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGagagga                                                                                                  >  4:1265236/1‑146 (MQ=255)
                                                           caccgcctcccttGATTCATTC‑CCCCA‑AACAC‑CCGCACTTTTCCCAAC‑CCCCATGTAAAAAAA‑TCCCTC‑CGGGGC‑GGTTTTCGCCAA‑GTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aaa                                                                                                <  1:1841330/134‑1 (MQ=255)
                                                            aGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGG‑aaaaaaatt                                                                                          >  1:365808/1‑145 (MQ=255)
                                                               gcctcccGTAATTCTTTC‑CCCCC‑AACAC‑CCGCGCTTTTCCCCAA‑CCCCAGGTACAAAAA‑GCCTTT‑CGCGGC‑GGCTTTCTCCCA‑GTTATTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTa                                                                                        <  2:120092/145‑1 (MQ=255)
                                                               gcctcccGTAATTCATTC‑CCGCC‑AACAC‑CCGCGCATTTACCCAA‑CGCCTTGTAAAAAAT‑TCCCTT‑CGGGGG‑GTCTGTCGCCCA‑GTTTTTTT‑C‑CCGCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                       <  4:30704/131‑1 (MQ=255)
                                                                gCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑gaaaaaatttt                                                                                        >  4:63381/1‑141 (MQ=255)
                                                                 cTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAG‑AAAAAACttt                                                                                         >  4:737025/1‑147 (MQ=255)
                                                                  tACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGAAGCG‑aaaaaaaatttaaac                                                                                    >  4:34590/1‑139 (MQ=25)
                                                                   aCCGTGTTTCATTC‑CCCCC‑ACCAC‑CCGCTATTTTACCCAC‑CGCCTTGTAAAAAAA‑TGCCTC‑CGGGGC‑GGCTGTCCGCCA‑GCTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTa                                                                                        <  4:1596665/146‑1 (MQ=255)
                                                                    ccGTGTTTCTTTC‑CCGCC‑AAC‑A‑CCGGCCTTTTTCCCAACCGCCTTTTG‑AAAAA‑CCCCTTGCGGGGG‑GGTTGCCTGCCG‑GTTATTTT‑T‑CCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                    <  4:1862722/149‑1 (MQ=255)
                                                                     ccgtaatCATTT‑CCCCCAAACAC‑CCGCGTTTTTATCCAA‑CCCCTTTTAAAAAAG‑TGCCTT‑GGGGGC‑GGCGGTCGGCCG‑CCTATTTT‑T‑CCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAg                                                                                             <  4:1203910/134‑1 (MQ=255)
                                                                        gATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑aaaaaaattttaatcc                                                                                   >  1:1659499/1‑133 (MQ=25)
                                                                         aatCATTTCCCGCC‑AACAC‑CCCCGAATTTCCCCAC‑CCCCATGTAAAAAAG‑TGCCTT‑CGCGGC‑GGCTTTCTCCCA‑GTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                       <  4:376781/125‑1 (MQ=255)
                                                                          ttCTTTC‑CCGCC‑AACAC‑CCGCGCATTTACCCAC‑CGCCATGTGAAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCGGCCA‑GTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg                                                                             <  2:109690/150‑1 (MQ=255)
                                                                           tttatccccGCA‑AACAC‑CCGCGTTTTTACCAAA‑CCCCATTGGAAAAAT‑TGCCTC‑CGGGGG‑GGTTGTCTGCAG‑GTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAt                                                                              <  1:1148803/142‑1 (MQ=255)
                                                                            cATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑GAAAAAATTTAAACCCCTcc                                                                               >  2:397371/1‑145 (MQ=25)
                                                                               ttcccGCA‑AACAC‑CCGCGCATTTACCCAC‑CGCCTTGTGAAAAAA‑TCCCTT‑CGGGGC‑GGCTTTCTGCCA‑GTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgttgt                                                                         <  4:1226735/147‑1 (MQ=255)
                                                                               tccccccc‑AACCC‑CCGCGCTTTTACCAAC‑CGCCATGGCAAAAAA‑TCCCTT‑CGGGGG‑GGCTGTCTCCCA‑TCTTTTTT‑C‑CCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtg                                                                             <  1:958497/140‑1 (MQ=255)
                                                                                  cccccaaacac‑cCGCGCTTTTACCCAA‑CGCAATGTAAAAAAA‑TCCCTC‑CGCGGG‑GGCTGTCTCCCA‑GTTTTTTT‑C‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgt                                                                            <  1:997894/139‑1 (MQ=255)
                                                                                        aacac‑cCGCGCTTTTACCCAC‑CG‑CTTGGGAAAAAAGTCCCTC‑CGGGGC‑GGTTTTC‑GCCA‑GTTTTTTTCC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAta                                                                <  1:19699/151‑2 (MQ=255)
                                                                                         accacccccgcTTTTATCCAC‑CGCCATGTGAAAAAA‑TCCCTG‑CGGGGC‑GGTTGTCTGCCA‑TTTTTTTT‑T‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                   <  1:825838/140‑1 (MQ=255)
                                                                                            a‑cCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAAAATTAAACCCGCCATGTTTTGTggggatagc                                                              >  2:171346/1‑142 (MQ=25)
                                                                                              ccGCGTATTTACCCAA‑CCCCATGTAAAAAAA‑TCCCTG‑CGCGGC‑GGCTGTCGCCCA‑TTTTTTTT‑C‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                            <  1:1247173/150‑1 (MQ=255)
                                                                                              ccGCGCTTTTCTCCAC‑CCCCATGTAAAAAAA‑GCCCTC‑CGCGGC‑GGCTGTCTGCCA‑CCTTTTTC‑C‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                            <  4:326564/150‑1 (MQ=255)
                                                                                              ccGCGCATTTACCCAA‑CGCCTTGTGCAAAAA‑TGCCTC‑CGCGGC‑GGCTGTCTGCCA‑GCTTTTTT‑C‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGCGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                             <  2:213984/149‑1 (MQ=255)
                                                                                               cgcgTATTTACCCAC‑CTCCTTGTGAAAAAA‑TCCCTC‑CGGGGC‑GGCTTCTTGCCA‑GTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg‑ct                                                         <  4:727028/151‑1 (MQ=255)
                                                                                               cgcgCTTTTACCCAC‑CGCCTTGTGAAAAAT‑TCCCTG‑GGGGGG‑GGCTGTCTGCCA‑CTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg‑ct                                                         <  4:1273125/151‑1 (MQ=255)
                                                                                                  gcttttACCAAA‑CGCAATGTGAAAAAA‑TCCTTG‑CGGGGG‑GGTTGTCTGCCA‑GCTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctg‑ctg                                                        <  3:237996/146‑1 (MQ=255)
                                                                                                      ttATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑AAAAAAAATTAAAACCCACACGTTTTGTGGCGATCgcag‑c                                                          >  4:141731/1‑143 (MQ=25)
                                                                                                        acccAC‑CGCCATGTGAAAAAA‑TGCCTC‑CGCGGC‑GGCTGTCTGCCA‑GTTTTTTT‑C‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACaa                                                 <  4:65450/148‑1 (MQ=255)
                                                                                                          ccAC‑CGCCTTGCGAAAAAA‑TGCCTC‑CGCGGG‑GGTTGTCGGCCA‑TCTTTTTT‑C‑CCCCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACa                                                  <  4:1147551/147‑1 (MQ=255)
                                                                                                           ccaccgcattgtaaaaaat‑TCCCTC‑CGGGGC‑GGCTGTC‑GCCATGTTTTTTT‑C‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgt                                              <  1:1235717/137‑1 (MQ=255)
                                                                                                           caaccccaatggcaaaaaa‑TCCTTC‑CGGGGG‑GGCTGTCTGCCA‑TTTTTTTT‑C‑CCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgt                                              <  1:381981/137‑1 (MQ=255)
                                                                                                           cAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAACTTTAAACCC‑CCATTTTGTGGGGGGATCGGTGCCTGGTAaaaa                                                 >  4:1653051/1‑148 (MQ=21)
                                                                                                             accGCAATGTAAAAAAG‑GGCCTC‑CGGGGC‑GGCTGTCGCCCG‑GATTTTTC‑C‑CCCCGGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgtg                                             <  1:865763/148‑1 (MQ=255)
                                                                                                             aacgccttggaaaaaaa‑TGCCTC‑CGGGGC‑GGCTGTCTGCCA‑TTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgtg                                             <  1:1242945/138‑1 (MQ=255)
                                                                                                              cccccATGTAAAAAAA‑CGCCTC‑GGGGCC‑GTCTGCCTCCCA‑CTTATTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtg                                               <  3:478169/144‑1 (MQ=255)
                                                                                                               ccctaggaaaaaaaa‑TGCCTC‑CGGGGG‑GGCTGTCGGCCA‑TTTTTTTT‑C‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAAtgtg                                             <  1:443771/138‑1 (MQ=255)
                                                                                                                 ccttgtgAAAAAT‑TGCCTC‑CGCGGC‑GGTTGTCGCCCA‑GTTATTTT‑C‑CCCCCTCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAAc                                                   <  4:1673842/137‑1 (MQ=255)
                                                                                                                 ccagggaaaaaaa‑CGCCTC‑CGCGGC‑GGCTGTCGGCCG‑TTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGaaa                                         <  1:262835/142‑1 (MQ=255)
                                                                                                                         aaaaa‑aGGCTT‑CGGGGGGGGCGGTGTGCAA‑GTTTTTTT‑C‑CCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGa                                     <  3:241042/147‑1 (MQ=255)
                                                                                                                         aaaaa‑TCCCTG‑GGGGGG‑GGCTGT‑TGGCA‑ACTTTTTT‑TCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTg                                 <  4:1281548/150‑1 (MQ=255)
                                                                                                                         aaaaa‑TCCCTC‑CGGGGG‑GCCTGTCGCCCA‑GTTTTTTT‑C‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAg                                   <  1:1715290/148‑1 (MQ=255)
                                                                                                                          aaaa‑TCCCTC‑CGCGGC‑GGCTGTCTGCAA‑GCTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTgcg                               <  4:80577/151‑1 (MQ=255)
                                                                                                                            aaatggcctggggggc‑ggcTGCCTGCCA‑GTTTTTTT‑T‑CCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTgcg                               <  4:1765605/138‑1 (MQ=255)
                                                                                                                                 ccTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAAAATTAAAACCCCCCTTTTGTGGTGGAGCCGCTG‑CTGGGTAAAATGTGGGAAGGGGAgggggggg                            >  1:1742480/1‑141 (MQ=255)
                                                                                                                                   tt‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAGAATTTAACACCCCCTGTTGTGTGGGGATCGCCG‑CTGGGGAAAAAGTGGGAAggggggggggggg                            >  4:1548243/1‑136 (MQ=255)
                                                                                                                                       ggggc‑ggcTGTCTCCCA‑GTTTTTTT‑C‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                         <  3:298147/144‑1 (MQ=255)
                                                                                                                                         ggc‑ggcTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAAATTTAAACCC‑CCCTGTTGTTTGGGAATCGCCG‑CCGGTGACAAAGAGGGAAGGGGAGGGGGGGGGATCaacccc                  >  4:573556/1‑146 (MQ=255)
                                                                                                                                          ggggcctttCTGCCA‑TTTTTTTC‑C‑CCGCCGCAAAAAGACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                        <  1:1754900/137‑1 (MQ=255)
                                                                                                                                           gcggttttcgcccA‑GCTATTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                          <  4:1920676/131‑1 (MQ=255)
                                                                                                                                                    ggCCA‑GTTTTTTT‑T‑CCGCAGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  4:573594/145‑1 (MQ=255)
                                                                                                                                                      ccA‑GTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  4:21090/144‑1 (MQ=255)
                                                                                                                                                      ccA‑GTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa      <  4:1827703/151‑1 (MQ=255)
                                                                                                                                                             aTTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAG‑AAAAAAATTTTAACACCCCCTTTTGTGTGGGAATCGCCG‑CTGGGAAAAAAGTGGGAAAGGGAGGGGGGGGAAAAAACCCCCGAACCGCGCGAAttt  >  4:328460/1‑148 (MQ=255)
                                                                                                                                                             |                                                                                                                                                        
CTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCC‑AACAA‑CCGCGCATTTATCCAA‑CGCCATGTGCAAAAA‑TGCCTT‑CGCGGC‑GGCTGTCTGCCA‑GCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTG‑CTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357726‑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: ‑