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

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357730‑358027
                                                                                                                                                         |                                                                                                                                                       
aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                         <  4:1066164/151‑1 (MQ=255)
  tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                       >  2:189856/1‑151 (MQ=255)
   cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                      <  4:169584/151‑1 (MQ=255)
    gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                                     <  4:141328/151‑1 (MQ=255)
     tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                    >  4:1328856/1‑151 (MQ=255)
      aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                   <  1:1292118/151‑1 (MQ=255)
          cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                               >  1:1453369/1‑151 (MQ=255)
            gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                              <  1:1055590/150‑1 (MQ=255)
               aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                          <  4:311546/151‑1 (MQ=255)
                aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAATTG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                         >  1:960258/1‑151 (MQ=255)
                  cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                        <  4:842140/150‑1 (MQ=255)
                    tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                     >  1:80790/1‑151 (MQ=255)
                     aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                                    >  1:841046/1‑151 (MQ=255)
                       cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                  >  4:1024163/1‑151 (MQ=255)
                       cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                  >  2:174553/1‑151 (MQ=255)
                                 aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                          >  4:709733/1‑149 (MQ=255)
                                   ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                      >  4:1091658/1‑151 (MQ=255)
                                     ccccccGCAGCCGTGTGTCGGGGCCTCCCGTGTTTCTTTC‑CCCCCAACACCCGCGCTTTTATCAAACGCCTTGTAAAAAATTC‑CCTC‑CGGGGCGGTTGTCGGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                        <  3:180308/143‑1 (MQ=255)
                                      cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                   >  4:643699/1‑151 (MQ=255)
                                       aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                  >  4:265081/1‑151 (MQ=255)
                                       aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                  >  1:1248575/1‑151 (MQ=255)
                                         ccGCAGCCTTGTGTCGGCGCCTCCCGTTTTTCTTTC‑CCGCAAACACCCGCGTTTTTACCCACCGCCGTGTAAAAAAATG‑CCTC‑CGCGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                 <  4:992550/150‑1 (MQ=255)
                                         ccGCAGCCTTGTGTCAGCGGCTCCCGTGTTTCTTTC‑CCCCCAACACCCGCGCTTTTCCCCACCCCCTTGTGCAAAAATC‑CCTC‑CGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                   <  1:983614/148‑1 (MQ=255)
                                         ccGCAGCCCTGTGTCAGCGGCTACCGTGTTTCTTTC‑CCGCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATC‑CCTG‑CGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                  <  1:994694/149‑1 (MQ=255)
                                         ccGCAGCCCTGTGCCGGCGCCCCCCGTGTTTCTTTC‑CCCCAAACACCCGCGTTTTTACCCACCCCCATGTAAAAAATTG‑CCTC‑CGGGGGGGCGGCCGCCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                 <  4:292756/150‑1 (MQ=255)
                                         ccGCAGACTTGTGTCGGGGCCTCCCGTGATTCTTTC‑CCCCCAACACCCGCCCTTTTACCCAACGCCTTGTG‑AAAAATT‑CCTTCCGGGGCGGTTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                  <  1:56401/149‑1 (MQ=255)
                                          cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                               >  4:496844/1‑150 (MQ=255)
                                          cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                >  1:372179/1‑150 (MQ=255)
                                            cagccgtgtgCCAGCGCCTCCCGTGATTCTTCC‑CCCCCAACACCCGCGCTTTTCCCCAACGCCATGTGAAAAAATG‑CCTC‑CGGGGCGTCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                             <  4:1232255/145‑1 (MQ=255)
                                              gcctgtgcccAGCGGCAACCGTGTTTCTTTC‑CCGCAAACACCCGCACTTTTACCAAACGCCAGGT‑AAAAAATGCCCGC‑CGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                             <  4:1202687/140‑1 (MQ=255)
                                              gcccTGTGCCAGGGCCTCCCGTGTTTCTTTC‑CCCCCAACACCCGCGCTTTTACCCACCGCCTTGTGAAAAAATC‑CTTC‑CGGGGGGGCTGTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                           <  1:156131/149‑1 (MQ=255)
                                               aCCTGTGTCAGCGGCTCCCGTGATTCATTT‑CCCCCAACACCCGCGATTTTACCCAACCCCATGTCAAAAAATG‑CCTT‑CGCGGGGGCTGCCTCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                  <  1:161173/143‑1 (MQ=255)
                                                ccttgtccggCGGCTCCCGTGATTCATTC‑CCCCCAACACCCGCGCTTTTACCCACCCCCTTGTGCAAAATTC‑CCTC‑CGCGGCGGTTGCCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                            <  1:569720/139‑1 (MQ=255)
                                                 cTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggggag                                                                                                          >  4:15904/1‑143 (MQ=255)
                                                 cTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                >  1:878203/1‑143 (MQ=255)
                                                  tgtgTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                               >  3:31396/1‑142 (MQ=255)
                                                    ttccAGCGGCTCCGTTGATTATTTC‑CCGCCAACACCCGCGCATTTACCCACCCCCTTGTGAAAAATCG‑CCTC‑GGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                              <  1:882375/139‑1 (MQ=255)
                                                    tgccAGCGCCTCCCGTAATTCTTTC‑CCCCCAACACCCGCCCTTTTACCCACCCCAATGTGAAAAAATG‑CCTC‑CGCGGCGGCTGTCTCCCACTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                      <  1:924682/147‑1 (MQ=255)
                                                     gtCGGCGGCTCCCGTTATTCTTTC‑CCGCCAACCCCCGCGCATTTACCCACCGCCATGTGCAAAAAGC‑CCTC‑CGGGGGGGCTTTCTGCCGGTTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                    <  4:221312/151‑1 (MQ=255)
                                                       cAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcagcgcaaa                                                                                                      >  4:297361/1‑137 (MQ=255)
                                                         gCGGCTCCCGTGTTTCTTTC‑CCGCCACCACCCGCGCTTTTACCCACCGCCTTGTAAAAAAATG‑CTTC‑GGGGGCGGCTGTCTCCGAGTTTTTTTCCCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                    <  1:618516/147‑1 (MQ=255)
                                                         gCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAA‑Cttt                                                                                               >  1:1373059/1‑151 (MQ=255)
                                                         gCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGgaaaaagatt                                                                                                >  4:745333/1‑144 (MQ=255)
                                                           ggCTCCCGTGATTCATTC‑CCGCCAACACCCGCGCTTTTACCCACCGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                    <  1:52958/145‑1 (MQ=25)
                                                            gCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAA‑ATTTaa                                                                                             >  1:903117/1‑150 (MQ=255)
                                                               aCCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAAGAAAAAGCCTTTAAAcc                                                                                          >  1:77570/1‑151 (MQ=255)
                                                                ccGTGTTTCTTTC‑CCCCCAACACCCGCGCTTTTACCCACCCCCATTTGAAAAATTC‑CCTC‑CGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCg                                                                                        <  1:262317/151‑1 (MQ=255)
                                                                ccGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCGGAGaaaaaaattttaacccc                                                                                         >  4:268662/1‑137 (MQ=25)
                                                                 ccgtaTTCCTTTCCCGCCACCAACCGCCATTTTACCCACCGCCATGTCAAAAAATG‑CCTC‑CGCGGGGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGccc                                                                                         <  4:313013/146‑1 (MQ=255)
                                                                    taattcttcccccccAACACCCGCGCTTTTACCCACCCCCTTGTAAAAAAAGC‑CCTC‑CGCGGGGGCTGCCTGCCAGTTTTTTTCCCGCGAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑Ct                                                                                                 <  4:723792/126‑1 (MQ=255)
                                                                      ttCTTTC‑CCCCCAACACCCGCGCTTTTACCCACCCCCTTGTGAAAAAATG‑CTTC‑CGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCAtg                                                                                   <  1:324729/150‑1 (MQ=255)
                                                                           ttccccccAACACCCGCGCTTTTACCCACCCCCTTGTGAAAAAATC‑CCTC‑CGGGGCGGTTGTCGCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCAtgt                                                                                  <  1:663391/141‑1 (MQ=255)
                                                                             ccccccaccaCCCGCGCTTTCACCCAACCCAATTTGAAAAAATG‑CCGC‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCAtgttgt                                                                               <  4:931769/144‑1 (MQ=255)
                                                                                        ccgcccttttccccACCGCCATTTCAAAAAATC‑CTTC‑CGCGGGGGCTTTCTCCCATTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCgctg                                                                 <  1:1440791/138‑1 (MQ=255)
                                                                                         cgctcttttACCCAACGCAATGTGAAAAAGTC‑CCTT‑GGCGGCGGCGTTCGGCCATTTTTTTTCCCGCCCCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCgctgct                                                               <  4:1030254/145‑1 (MQ=255)
                                                                                           cgcTTTTACCCACCCCCATGTAAAAAAATC‑CCTC‑CGGGGGGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGctgct‑‑g                                                            <  1:1295545/150‑1 (MQ=255)
                                                                                               tttATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAG‑CATTTAACCCCCCCTGTTGTG‑TGGCAATCgc                                                                   >  1:399145/1‑141 (MQ=37)
                                                                                               tttACCCACCCCCATGTAAAAAATTC‑CCCC‑CGGGGGGGCTGTCTCCCATTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGt                                                          <  4:146337/148‑1 (MQ=255)
                                                                                                        cGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAA‑CATTAAACCCCCCATCTTTTG‑TGGGAAACGCtggt‑‑gggtaaaaattg                                                 >  1:1154509/1‑138 (MQ=25)
                                                                                                             tgtgCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGGGAAAAAA‑CATTTAACCCCTCCTTTTGTG‑TGGGGAACGCTGCT‑‑GGGGACAAAATGGGAAggg                                          >  1:568071/1‑150 (MQ=255)
                                                                                                                  aaaaaTG‑CCTT‑CGGGGCGGTTGTCCGCCACTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAg                                       <  1:462100/148‑1 (MQ=255)
                                                                                                                  aaaaaTC‑CCTT‑GGCGGCGGCTGTCGCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGaa                                        <  4:1336957/147‑1 (MQ=255)
                                                                                                                   aaaaTG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAATAT‑ATTTTACACCCCCTTGTTTTG‑TGGGCATCCCCGCTGGGGTAAAAATG‑GGAAAGGGagtggg                                    >  2:19597/1‑146 (MQ=255)
                                                                                                                    aaaTG‑CCTC‑CGGGGGGGCTGTCGCCCATTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTgcgg                                  <  4:429547/151‑1 (MQ=255)
                                                                                                                    aaaTC‑CCTT‑CGCGGCGGCTGTCTCCCAGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTgcgg                                  <  1:1350146/151‑1 (MQ=255)
                                                                                                                     aaTG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAA‑ATTTTAACCCCTCCTGTTGTG‑TGGGAATCGCTGCT‑‑GGTGAAAAAGTGGAAAGGGGAGGgcggg                                 >  1:842602/1‑150 (MQ=255)
                                                                                                                      aTG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAA‑CATTAAACCCCTCCTGTTTTGTTGGGAACCG‑TGCT‑‑GGGTAAAATTTGGGAAGGGGAgggggggg                                >  3:91003/1‑144 (MQ=255)
                                                                                                                          ccTT‑GGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGa                               <  1:888597/149‑1 (MQ=255)
                                                                                                                                  ggcggcTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTgcggcg                                <  2:93248/141‑1 (MQ=255)
                                                                                                                                  ggcgcttttCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGaa                              <  4:1252503/135‑1 (MQ=255)
                                                                                                                                          tctcCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa                   <  1:645096/146‑1 (MQ=255)
                                                                                                                                            tgcccttttttttcCCGCCCAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                     <  4:807373/133‑1 (MQ=255)
                                                                                                                                             cccattttttttcCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg             <  3:103474/141‑1 (MQ=255)
                                                                                                                                              caagtttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                  <  1:1193358/136‑1 (MQ=255)
                                                                                                                                                   ttttttcCCCCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg      <  4:28885/148‑1 (MQ=255)
                                                                                                                                                        ttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGAGCCAGAGAAAAAG‑ATTTTAAAACCTCAACTTGTG‑TGGGAATCCCTGCT‑‑GGGAACAAAGTGGGAAGGGGAGGGCGGCGGAACAGCTCCCGAACCCCGCCCATTggggt  >  2:135547/1‑145 (MQ=255)
                                                                                                                                                         |                                                                                                                                                       
ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357730‑358027

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