Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8730CT42.5% 105.2 / 41.0 80S273N (AGC→AAC) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (33/13);  new base T (0/34);  total (33/47)
Fisher's exact test for biased strand distribution p-value = 4.64e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGA  >  NC_000913/357725‑358010
                                                                                                                                                       |                                                                                                                                             
aCTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTa                                                                                                                                             <  1:955844/151‑1 (MQ=255)
 cTGACCTCGTAGTCTTGATACGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTa                                                                                                                                             <  4:584492/150‑1 (MQ=255)
 cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTa                                                                                                                                             <  4:1188186/150‑1 (MQ=255)
 cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTAt                                                                                                                                            >  4:602084/1‑151 (MQ=255)
  tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTAtt                                                                                                                                           >  1:1204611/1‑151 (MQ=255)
   gAACTCGTAGGCCTGCTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGAGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGACATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTTTGCCAGCTAttt                                                                                                                                          >  1:1252609/1‑151 (MQ=255)
   gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTAttt                                                                                                                                          <  4:621991/151‑1 (MQ=255)
   gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTAt                                                                                                                                            <  1:319183/149‑1 (MQ=255)
     aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                        <  4:1066164/151‑1 (MQ=255)
     aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTAtttt                                                                                                                                         <  4:331863/150‑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)
                                                                                                                                                       |                                                                                                                                             
ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTT‑CCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAG‑CTTTAAGCCCGTCATGTTGTG‑TGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGA  >  NC_000913/357725‑358010

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: ‑