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

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA  >  NC_000913/357733‑358024
                                                                                                                                                       |                                                                                                                                                     
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                    >  1:580475/1‑151 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                    >  2:1938336/1‑151 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                    >  2:1893546/1‑151 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                                      <  2:1520880/149‑1 (MQ=255)
 gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAGATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                                   <  2:928441/151‑1 (MQ=255)
  tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                  >  2:1456459/1‑151 (MQ=255)
   aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                  <  2:1894766/150‑1 (MQ=255)
    ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                  >  1:170095/1‑149 (MQ=255)
     gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                                                                                >  2:1008878/1‑150 (MQ=255)
      ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                               <  1:769636/150‑1 (MQ=255)
       cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                             >  2:560089/1‑151 (MQ=255)
        tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                             >  2:1489631/1‑150 (MQ=255)
          aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                          >  1:1510675/1‑150 (MQ=255)
          aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                         <  2:965819/151‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                        <  1:548192/150‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                       <  1:1302530/151‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                       >  1:1090089/1‑151 (MQ=255)
             aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                      <  2:570887/151‑1 (MQ=255)
              gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                     >  1:1742135/1‑151 (MQ=255)
                gtggcgcATCAGGCAATTCCAGCCGCAGACCTGTGCCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                   <  1:1019065/148‑1 (MQ=255)
                      cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAATTGCTTTCGGGCCGGTGGCCTCCAAGCAATTTTCCGGCGCAAACAAAACGCCGGtc                                                                                                                               >  2:516288/1‑148 (MQ=14)
                       aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                                                                                                                            >  2:1840328/1‑151 (MQ=255)
                        tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                                                           >  1:1741195/1‑151 (MQ=255)
                          aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                          >  2:164634/1‑150 (MQ=255)
                          aGGCAATTCCAACCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACTACCGCACATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAATCGCCGTTCTCTcc                                                                                                                          >  1:1420711/1‑150 (MQ=255)
                                tgccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                   >  1:275211/3‑151 (MQ=255)
                                  ccccccccAGACTTGTGCCAGCGCCTCCCGTATTTCTTTCCCCCCAACACCCGCGCTTTTCCCCACCCCCTTGTGAAAAAAGCCCTTCGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCag                                                                                                                        <  1:980264/137‑1 (MQ=255)
                                  ccAGCCGCTGTCCTGTTCCAGCGGCTCCCTTATTTCTTTCCCCCAAACACCCGCGCATTTCCCCACCCCATTGTAAAAAAATGCCTCCGCGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                 <  1:167275/151‑1 (MQ=255)
                                  ccAGCCGCAGCCCTGTGCCAGGGCCCCCCGTTATTCTTTCCCCCCAACACCCGCGCTTTTATCCACCCCCTTGTGAAAAAAGGCCTCCGCGGCGGCTGTCTGCAAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                   <  2:392517/149‑1 (MQ=255)
                                  ccAGCCGAAGACTTGTGTCAGGGGCTCCGGGGATTCATTCCCGCAAACACCCGCGCTTTTACCCAACGCCTTGTGCAAAAATGCCTCCGGGGGGGCTGTCTCCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                   <  2:933101/149‑1 (MQ=255)
                                   cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑g                                                                                                              >  2:1054622/1‑151 (MQ=255)
                                   cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑g                                                                                                              >  1:358664/1‑151 (MQ=255)
                                   cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑g                                                                                                              >  2:289151/1‑151 (MQ=255)
                                   cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑g                                                                                                              >  2:620202/1‑151 (MQ=255)
                                      ccacagccgtgtgCCAGCGGCTCCCGGGATTCTTTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGCAAAATTGCCTTCGCGGCGCCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggg                                                                                                           <  1:74405/142‑1 (MQ=255)
                                      ccGCAGCCCTGTGTCAGCGGCTACCGTGATTCTTTCCCGCCAACACCCGCGATTTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑ggg                                                                                                            <  2:1664158/150‑1 (MQ=21)
                                       cgcacccgtgttTCAGCGGCTCCCTTTTTTCTTTTCCCCCAACACCCCCCCATTTCCCCAACGCCTTTTGAAAAACTGCCTCCGGGGCGGTTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggg                                                                                                           <  2:490702/142‑1 (MQ=255)
                                       cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑ggg                                                                                                            >  1:1840195/1‑149 (MQ=255)
                                         cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑ggcgg                                                                                                          >  2:1385238/1‑146 (MQ=255)
                                          aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑ggggg                                                                                                          >  2:1828828/1‑147 (MQ=255)
                                           gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggg                                                                                                           >  2:276482/1‑146 (MQ=255)
                                             ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGagag                                                                                                     >  1:679009/1‑150 (MQ=255)
                                              cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGCAgaga                                                                                                    >  2:1090180/1‑150 (MQ=255)
                                               tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggg                                                                                                           >  1:1927349/1‑142 (MQ=255)
                                                     aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggggcgaggagaaaaa                                                                                              >  1:1050189/1‑136 (MQ=255)
                                                       cGGCTCCCGTGTTTCATCCCCCCCAACACCCGCGCTTTTCCCCAACCCCTTGTGAAAAAATGCTTGCGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGaa                                                                                                   <  1:1963481/142‑1 (MQ=255)
                                                        ggCTCCCGTAATTCATCCCCCCCAACACCCGCGCTTTTACCCACCCCCATGTGAAAAAATCCCTCCGCGGCGGCTGTCGGCAAGCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTa                                                                                         <  1:1714806/150‑1 (MQ=255)
                                                        ggCTACCGTGATTCTTTCCCGCCAACACCCGCGCATTTACCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTTTCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTaa                                                                                        <  2:90508/151‑1 (MQ=25)
                                                            ccctgttttCTTTCCCGCCAACACCCGGTATTTTACCAACCCCCATGTGAAAAAAGCCCTCCGCGGCGCCTTTCTGCAAGTTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGccc                                                                                    <  1:1029627/144‑1 (MQ=255)
                                                             ccgtatttCATTCCCCCCAACACCCGCCCTTTTACCCACCCCCTTGTGAAAAATGGCCTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGc                                                                                      <  1:885827/142‑1 (MQ=255)
                                                              ccggatttctttCCCCCAAACACCCGCGTTTTTACCCACCCCCAGGTAAAAAATTCCTTTCGGGGCGGCTGTCGCCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGc                                                                                      <  1:1907464/139‑1 (MQ=255)
                                                                  aTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGAGaaaaaaattttaaacc                                                                                     >  2:760795/1‑132 (MQ=25)
                                                                    tCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAAAAACCTTAAAACCCTCCtgtt                                                                           >  2:984516/1‑151 (MQ=25)
                                                                    ccATTCCCCCAAACACCCCCGAATTTACCAACCGCAATGTAAAAAAAGCCCTCGGGGGGGGCTTCCGGCAGGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAt                                                                               <  1:1008859/147‑1 (MQ=255)
                                                                        ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGAGAAAAAACCTTTAACCCGTCCCTTtgtgt                                                                       >  1:1181364/1‑151 (MQ=25)
                                                                          cccccAACACCCGCGCATTTCCCAACCGCCATGTAAAAAAATCCCTTCGGGGGGGCTGTCGGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                            <  1:605189/141‑1 (MQ=255)
                                                                          cccccAACACCCGCCCATTTACCCACCGCCATGTGAAAAAATCCCTTCGCGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgt                                                                       <  2:1397849/146‑1 (MQ=255)
                                                                             ccAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGCAAGGAAAAAGATTTAAGCCGGTc                                                                                 >  2:1384590/1‑138 (MQ=25)
                                                                                aacaccCGCGCTTTTACCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                <  2:756120/151‑1 (MQ=38)
                                                                                       gcccattaccccAACGCCAGGTCAAAAAAGCCCTGCGCGGCGGCTGCCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                          <  1:606108/140‑1 (MQ=255)
                                                                                              tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAAAGAAAAAACATTTAAACCCCCATGTTTTGGGGGGATCGCTGCTGGtgaaa                                                   >  2:397757/1‑146 (MQ=25)
                                                                                                 ccACCGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTGTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                               <  2:1963078/151‑1 (MQ=255)
                                                                                                  ccacccccgtgtaaaaaaaGCCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                <  2:1085837/136‑1 (MQ=255)
                                                                                                   aaCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGCGAAAAAAAATTAAACCCCCCCTGTTGTTGGGGCATCGCTGCTggggaca                                                   >  1:236927/1‑140 (MQ=21)
                                                                                                    acccccATGTGCAAAATGCCTTCCGCGGCGTCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                             <  1:1268688/147‑1 (MQ=255)
                                                                                                      cccaaTGTGAAAAAATGCCTCCGCGGCGGTTGTCTGCCGGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                            <  2:1474553/146‑1 (MQ=255)
                                                                                                        ccATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGAAGAGAAAAAACCTTAATCCCCGCCTGCTTTTGTGGAATCGgtgctgtggggtcaat                                               >  1:854643/1‑135 (MQ=25)
                                                                                                         ccgtgtgaaaaaaGCCTTCCGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                         <  1:1673284/144‑1 (MQ=255)
                                                                                                            gtgAAAAAATCCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                     <  2:1717033/151‑1 (MQ=255)
                                                                                                                aaaatTCCCTCCGGGGCGCTTGTCCCCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                  <  1:1456434/150‑1 (MQ=255)
                                                                                                                aaaaaTTCCTTGGCGGGGGCGGGCTGCCGACCTATTTTCCCCCGGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                          <  2:815717/143‑1 (MQ=255)
                                                                                                                aaaaaTCCTTCCGGGGGGGCTTTCTCCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                 <  1:722030/151‑1 (MQ=255)
                                                                                                                aaaaaTCCCTGCGCGGCGCCTTTCTGCCGGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                   <  1:507422/149‑1 (MQ=255)
                                                                                                                 aaaaTGCCTTCGCGGGGGCTGTCTCCCATCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                <  1:1443590/151‑1 (MQ=255)
                                                                                                                       ccTTCGCGGCGGTTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                             <  2:394529/148‑1 (MQ=255)
                                                                                                                       ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGAGAAAAAAATTTTAAACCCCCCTTTTTTGTGGGGATTGCCGGTGGGAAAAAAGTGggaagggaaggggggg                              >  1:1902162/1‑141 (MQ=255)
                                                                                                                          tccggggcgccTGTCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                         <  2:1087478/145‑1 (MQ=255)
                                                                                                                          tCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGAGAAAAAAATTTTAACCCCCCCTGTTTTGTGGGGAATGCCGCTGGGTACAAAGTGGGAAGGGGAgggggggggaac                         >  2:1811705/1‑138 (MQ=255)
                                                                                                                           ccgcgGCGCTTTTCTGCAGGTTTTTTTCCCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                          <  1:226325/147‑1 (MQ=255)
                                                                                                                               ggcggcTTCCTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                     <  1:1389652/149‑1 (MQ=255)
                                                                                                                                  gcgtttCTGCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                      <  1:470864/140‑1 (MQ=255)
                                                                                                                                        cggCCACCTTTTTCCCCGCCGCACAAACACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc             <  2:1002404/146‑1 (MQ=255)
                                                                                                                                         tcccAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc             <  2:572098/145‑1 (MQ=255)
                                                                                                                                         ctcccggtttttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc              <  1:1233884/137‑1 (MQ=255)
                                                                                                                                           ccAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg            <  2:553765/146‑1 (MQ=255)
                                                                                                                                           ccAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt        <  2:772334/150‑1 (MQ=255)
                                                                                                                                            ccgtcttttttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt        <  1:82493/143‑1 (MQ=255)
                                                                                                                                            ccggtttttttcCCCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc           <  1:18721/140‑1 (MQ=255)
                                                                                                                                              agtttttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt     <  1:331980/146‑1 (MQ=255)
                                                                                                                                                   ttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa  <  1:356676/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA  >  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: