Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC42.4% 123.5 / 26.3 85K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (38/11);  new base C (0/36);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 2.06e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA  >  NC_000913/357733‑358024
                                                                                                                                                       |                                                                                                                                                     
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTccg                                                                                                                                                    >  1:580475/1‑151 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTccg                                                                                                                                                    >  2:1938336/1‑151 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTccg                                                                                                                                                    >  2:1893546/1‑151 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTc                                                                                                                                                      <  2:1520880/149‑1 (MQ=255)
 gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAGATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTccgc                                                                                                                                                   <  2:928441/151‑1 (MQ=255)
  tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTccgcc                                                                                                                                                  >  2:1456459/1‑151 (MQ=255)
   aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTccgcc                                                                                                                                                  <  2:1894766/150‑1 (MQ=255)
    ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTccgcc                                                                                                                                                  >  1:170095/1‑149 (MQ=255)
     gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCcgccgc                                                                                                                                                >  2:1008878/1‑150 (MQ=255)
      ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGca                                                                                                                                               <  1:769636/150‑1 (MQ=255)
       cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGcaac                                                                                                                                             >  2:560089/1‑151 (MQ=255)
        tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGcaac                                                                                                                                             >  2:1489631/1‑150 (MQ=255)
          aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGcaac‑aa                                                                                                                                          >  1:1510675/1‑150 (MQ=255)
          aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑aaa                                                                                                                                         <  2:965819/151‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑aaaa                                                                                                                                        <  1:548192/150‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAAc                                                                                                                                       <  1:1302530/151‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAAc                                                                                                                                       >  1:1090089/1‑151 (MQ=255)
             aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAAcc                                                                                                                                      <  2:570887/151‑1 (MQ=255)
              gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAAccg                                                                                                                                     >  1:1742135/1‑151 (MQ=255)
                gtggcgcATCAGGCAATTCCAGCCGCAGACCTGTGCCAGCGGCTACCGTGATTCATTCCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAAccgcc                                                                                                                                   <  1:1019065/148‑1 (MQ=255)
                      cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAATTGCTTT‑CGGGCCGGTGGCCTCCA‑AGC‑AATTTTCCGGCGCAAA‑CAAAACGCCGGtc                                                                                                                               >  2:516288/1‑148 (MQ=14)
                       aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTtctct                                                                                                                            >  2:1840328/1‑151 (MQ=255)
                        tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTtctctc                                                                                                                           >  1:1741195/1‑151 (MQ=255)
                          aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTcc                                                                                                                          >  2:164634/1‑150 (MQ=255)
                          aGGCAATTCCAACCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACTACCGCACATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAATCGCCGTTCTCTcc                                                                                                                          >  1:1420711/1‑150 (MQ=255)
                                tgccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                   >  1:275211/3‑151 (MQ=255)
                                  ccccccccAGACTTGTGCCAGCGCCTCCCGTATTTCTTTCCCCCC‑AACACCCGCGCTTTTCCCCAC‑CCCCTTGTGAAAAAAGCCCTT‑CGGGGCGGCTGTCTCCC‑AGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCag                                                                                                                        <  1:980264/137‑1 (MQ=255)
                                  ccAGCCGCTGTCCTGTTCCAGCGGCTCCCTTATTTCTTTCCCCCA‑AACACCCGCGCATTTCCCCAC‑CCCATTGTAAAAAAATGCCTC‑CGCGGCGGCTGTCTCCC‑AGCTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                 <  1:167275/151‑1 (MQ=255)
                                  ccAGCCGCAGCCCTGTGCCAGGGCCCCCCGTTATTCTTTCCCCCC‑AACACCCGCGCTTTTATCCAC‑CCCCTTGTGAAAAAAGGCCTC‑CGCGGCGGCTGTCTGCA‑AGTTATTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                   <  2:392517/149‑1 (MQ=255)
                                  ccAGCCGAAGACTTGTGTCAGGGGCTCCGGGGATTCATTCCCGCA‑AACACCCGCGCTTTTACCCAA‑CGCCTTGTGCAAAAATGCCTC‑CGGGGGGGCTGTCTCCC‑AGTTATTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                   <  2:933101/149‑1 (MQ=255)
                                   cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑g                                                                                                              >  2:1054622/1‑151 (MQ=255)
                                   cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑g                                                                                                              >  1:358664/1‑151 (MQ=255)
                                   cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑g                                                                                                              >  2:289151/1‑151 (MQ=255)
                                   cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑g                                                                                                              >  2:620202/1‑151 (MQ=255)
                                      ccacagccgtgtgCCAGCGGCTCCCGGGATTCTTTCCCCCC‑AACACCCGCGCTTTTACCCAC‑CGCCATGTGCAAAATTGCCTT‑CGCGGCGCCTGTCTGCC‑AGTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggg                                                                                                           <  1:74405/142‑1 (MQ=255)
                                      ccGCAGCCCTGTGTCAGCGGCTACCGTGATTCTTTCCCGCC‑AACACCCGCGATTTTACCCAC‑CGCCATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGTTTTTTTTCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑ggg                                                                                                            <  2:1664158/150‑1 (MQ=21)
                                       cgcacccgtgttTCAGCGGCTCCCTTTTTTCTTTTCCCCC‑AACACCCCCCCATTTCCCCAA‑CGCCTTTTGAAAAACTGCCTC‑CGGGGCGGTTTTCTGCC‑AGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggg                                                                                                           <  2:490702/142‑1 (MQ=255)
                                       cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑ggg                                                                                                            >  1:1840195/1‑149 (MQ=255)
                                         cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑ggcgg                                                                                                          >  2:1385238/1‑146 (MQ=255)
                                          aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑ggggg                                                                                                          >  2:1828828/1‑147 (MQ=255)
                                           gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggg                                                                                                           >  2:276482/1‑146 (MQ=255)
                                             ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGagag                                                                                                     >  1:679009/1‑150 (MQ=255)
                                              cTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGCAgaga                                                                                                    >  2:1090180/1‑150 (MQ=255)
                                               tgtgTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggg                                                                                                           >  1:1927349/1‑142 (MQ=255)
                                                     aGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggggcgaggagaaaaa                                                                                              >  1:1050189/1‑136 (MQ=255)
                                                       cGGCTCCCGTGTTTCATCCCCCCC‑AACACCCGCGCTTTTCCCCAA‑CCCCTTGTGAAAAAATGCTTG‑CGGGGCGGCTGTCTCCC‑AGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGaa                                                                                                   <  1:1963481/142‑1 (MQ=255)
                                                        ggCTCCCGTAATTCATCCCCCCC‑AACACCCGCGCTTTTACCCAC‑CCCCATGTGAAAAAATCCCTC‑CGCGGCGGCTGTCGGCA‑AGCTTTTTTCCCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTa                                                                                         <  1:1714806/150‑1 (MQ=255)
                                                        ggCTACCGTGATTCTTTCCCGCC‑AACACCCGCGCATTTACCCAC‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTTTCTGCC‑AGTTATTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTaa                                                                                        <  2:90508/151‑1 (MQ=25)
                                                            ccctgttttCTTTCCCGCC‑AACACCCG‑GTATTTTACCAACCCCCATGTGAAAAAAGCCCTC‑CGCGGCGCCTTTCTGCA‑AGTTTTTTTCCCCCCCAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGccc                                                                                    <  1:1029627/144‑1 (MQ=255)
                                                             ccgtatttCATTCCCCCC‑AACACCCGCCCTTTTACCCAC‑CCCCTTGTGAAAAATGGCCTC‑CGCGGCGGCTGTCTGCC‑AGTTATTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGc                                                                                      <  1:885827/142‑1 (MQ=255)
                                                              ccggatttctttCCCCCAAACACCCGCGTTTTTACCCAC‑CCCCAGGTAAAAAATTCCTTT‑CGGGGCGGCTGTCGCCC‑AGTTTTTTTCCCGCCCAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGc                                                                                      <  1:1907464/139‑1 (MQ=255)
                                                                  aTTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGAGaaaaaaattttaaacc                                                                                     >  2:760795/1‑132 (MQ=25)
                                                                    tCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAAAAAC‑CTTAAAACCCTCCtg‑tt                                                                           >  2:984516/1‑151 (MQ=25)
                                                                    ccATTCCCCCA‑AACACCCCCGAATTTACCAAC‑CGCAATGT‑AAAAAAAGCCCT‑CGGGGGGGGCTTCCGGC‑AGGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCAt                                                                               <  1:1008859/147‑1 (MQ=255)
                                                                        ttCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGAGAAAAAAC‑CTTTAACCCGTCCCT‑Ttgtgt                                                                       >  1:1181364/1‑151 (MQ=25)
                                                                          ccccc‑AACACCCGCGCATTTCCCAAC‑CGCCATGTAAAAAAATCCCTT‑CGGGGGGGCTGTCGGCC‑AGTTTTTTTCCCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCAtg‑t                                                                            <  1:605189/141‑1 (MQ=255)
                                                                          ccccc‑AACACCCGCCCATTTACCCAC‑CGCCATGTGAAAAAATCCCTT‑CGCGGCGGCTTTCTGCC‑AGCTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑Ttgtgt                                                                       <  2:1397849/146‑1 (MQ=255)
                                                                             cc‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGCAAG‑GAAAAAGA‑TTTAAGCCGGTc                                                                                 >  2:1384590/1‑138 (MQ=25)
                                                                                aacaccCGCGCTTTTACCCAC‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTTTTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATc                                                                <  2:756120/151‑1 (MQ=38)
                                                                                       gcccattaccccAA‑CGCCAGGTCAAAAAAGCCCTG‑CGCGGCGGCTGCCGGCC‑AGTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCgctgct                                                          <  1:606108/140‑1 (MQ=255)
                                                                                              tATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAAAGAAAAAACATTTAAACCC‑CCATG‑TTTTGGGGGGATCGCTGCTGGtgaaa                                                   >  2:397757/1‑146 (MQ=25)
                                                                                                 ccAC‑CGCCATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCTGCC‑AGTTGTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAAtgt                                               <  2:1963078/151‑1 (MQ=255)
                                                                                                  ccacccccgtgtaaaaaaaGCCCTC‑CGCGGCGGCTGTCTGCC‑AGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAAtg                                                <  2:1085837/136‑1 (MQ=255)
                                                                                                   aa‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGCGAAAAAAA‑ATTAAACCCCCCCTG‑TTGTTGGGGCATCGCTGCTggggaca                                                   >  1:236927/1‑140 (MQ=21)
                                                                                                    acccccATGTGC‑AAAATGCCTTCCGCGGCGTCTGTCGGCC‑AGTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                             <  1:1268688/147‑1 (MQ=255)
                                                                                                      cccaaTGTGAAAAAATGCCTC‑CGCGGCGGTTGTCTGCC‑GGCTATTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                            <  2:1474553/146‑1 (MQ=255)
                                                                                                        ccATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGAAGAGAAAAAAC‑CTTAATCCCCGCCTGCTTTTGTGG‑AATCGgtgctgtggggtcaat                                               >  1:854643/1‑135 (MQ=25)
                                                                                                         ccgtgtgaaaaaaGCCTTCCGGGGCGGCTGTCGGCC‑AGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                         <  1:1673284/144‑1 (MQ=255)
                                                                                                            gtgAAAAAATCCCTT‑CGCGGCGGCTGTCTGCC‑AGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                     <  2:1717033/151‑1 (MQ=255)
                                                                                                                aaaatTCCCTC‑CGGGGCGCTTGTCCCCC‑AGTTTTTTTCCCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                  <  1:1456434/150‑1 (MQ=255)
                                                                                                                aaaaaTTCCTT‑GGCGGGGGCGGGCTGCCGACCTATTTTCCCCCGGAAAA‑AAACCCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                          <  2:815717/143‑1 (MQ=255)
                                                                                                                aaaaaTCCTTC‑CGGGGGGGCTTTCTCCC‑GGTTTTTTTCCCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                 <  1:722030/151‑1 (MQ=255)
                                                                                                                aaaaaTCCCTG‑CGCGGCGCCTTTCTGCC‑GGTTATTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                   <  1:507422/149‑1 (MQ=255)
                                                                                                                 aaaaTGCCTT‑CGCGGGGGCTGTCTCCC‑ATCTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                <  1:1443590/151‑1 (MQ=255)
                                                                                                                       ccTT‑CGCGGCGGTTGTCTGCC‑AGCTTTTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                             <  2:394529/148‑1 (MQ=255)
                                                                                                                       ccTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGAGAAAAAAA‑TTTTAAACCCCCCTT‑TTTTGTGGGGATTGCCGGTGGGAAAAAAGTGggaagggaaggggggg                              >  1:1902162/1‑141 (MQ=255)
                                                                                                                          tccggggcgccTGTCTGCC‑AGTTATTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                         <  2:1087478/145‑1 (MQ=255)
                                                                                                                          t‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGAGAAAAAAA‑TTTTAACCCCCCCTG‑TTTTGTGGGGAATGCCGCTGGGTACAAAGTGGGAAGGGGAgggggggggaac                         >  2:1811705/1‑138 (MQ=255)
                                                                                                                           ccgcgGCGCTTTTCTGCA‑GGTTTTTTTCCCGCCAAAAA‑AAATCCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                          <  1:226325/147‑1 (MQ=255)
                                                                                                                               ggcggcTTCCTGCC‑ATTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                     <  1:1389652/149‑1 (MQ=255)
                                                                                                                                  gcgtttCTGCC‑GGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                      <  1:470864/140‑1 (MQ=255)
                                                                                                                                        cggCC‑ACCTTTTTCCCCGCCGCACA‑AACACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc             <  2:1002404/146‑1 (MQ=255)
                                                                                                                                         tccc‑AGTTTTTTTCCCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc             <  2:572098/145‑1 (MQ=255)
                                                                                                                                         ctcccggtttttttcCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc              <  1:1233884/137‑1 (MQ=255)
                                                                                                                                           cc‑AGTTTTTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg            <  2:553765/146‑1 (MQ=255)
                                                                                                                                           cc‑AGTTTTTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt        <  2:772334/150‑1 (MQ=255)
                                                                                                                                            ccgtcttttttcCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt        <  1:82493/143‑1 (MQ=255)
                                                                                                                                            ccggtttttttcCCCCCGAAAA‑AAACCCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc           <  1:18721/140‑1 (MQ=255)
                                                                                                                                              agtttttttcCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt     <  1:331980/146‑1 (MQ=255)
                                                                                                                                                   ttttcCCCCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa  <  1:356676/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA  >  NC_000913/357733‑358024

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