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

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357733‑358026
                                                                                                                                                      |                                                                                                                                                    
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                   <  1:1955930/151‑1 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                   <  1:287073/151‑1 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                    <  1:987141/150‑1 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                    <  2:627353/150‑1 (MQ=255)
  tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                 >  1:361430/1‑151 (MQ=255)
  tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                                  >  2:1743184/1‑150 (MQ=255)
   aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                <  2:878369/151‑1 (MQ=255)
   aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                 >  2:1140213/1‑150 (MQ=255)
         gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                           >  2:2026022/1‑150 (MQ=255)
          aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                         <  1:1243860/151‑1 (MQ=255)
          aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                         <  1:1848082/151‑1 (MQ=255)
           tAAGCGTAGCGCTTCAGGCAATTCCAGCCGCGGACCTGTGTCAGCGGCTACCGGGATTCATTTCCGCCAACAC‑CCGCGCATTCATCCAACGCCATGTGCAAAAAT‑GCATT‑CGCGTCGGCTGTCTGCCAGCTA‑TTTTCCGCCGcaacaa                                                                                                                                          <  2:1237731/148‑1 (MQ=14)
           tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                         >  2:768295/1‑150 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                        <  2:1258916/150‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                       >  1:1792281/1‑151 (MQ=255)
             aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                      >  2:352727/1‑151 (MQ=255)
                gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                   >  1:727587/1‑151 (MQ=255)
                gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgac                                                                                                                                   >  1:1616099/1‑149 (MQ=255)
                   gcgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGtt                                                                                                                                >  1:1635819/1‑151 (MQ=255)
                     gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                              >  1:668219/1‑151 (MQ=255)
                      cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                             >  1:1114480/1‑151 (MQ=255)
                      cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                             >  1:1683744/1‑151 (MQ=255)
                         cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                          >  2:1068227/1‑151 (MQ=255)
                               aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                    >  1:281082/1‑151 (MQ=255)
                               aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                     >  2:1875817/1‑150 (MQ=255)
                                tccccgCCGCAGCCCTGTGTCGGCGGCTCCCGTGTTTCTTTCCCCCCAACAC‑CCGCGCTTTTACCCAACGCCTTGTGAAAAAAT‑CCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagta                                                                                                                      <  1:1319676/143‑1 (MQ=255)
                                   cAGCCGCAGACCTGTGTCAGAGGCTAACGTGATTCATTTCCGAAAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAACT‑GCCTT‑CGCGGACGCTGTCTGCCAGCTATTTTTCAGCAGCAACAACACCGCCGTTCACTCCAGTAGTGGc                                                                                                                 >  2:2139920/1‑150 (MQ=255)
                                              cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggg                                                                                                           >  2:1700009/1‑143 (MQ=255)
                                               tgtgCCGGCGCCTCCCGTGATTATTTCCCCCCAACAC‑CCGCGCTTTTACCCACCGCCTTGTGAAAAAAT‑GCCTC‑CGCGGCGGTGGTCGGCCGGTTTTTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCaa                                                                                                         <  2:1704582/145‑1 (MQ=255)
                                                    cggcgcCTCCCGTCTTTCATTCCCCCCAACAC‑CCCCGCATTTACCCACCCCGATGGGAAAAAGG‑GCCTC‑GGCGCCGGCTCTCTCCCGTTTATTTCCCGGCCGCACAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAg                                                                                               <  1:1490716/147‑1 (MQ=255)
                                                     aGCGGCTACCGTGAATCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGAGGCTGTCTGCCAGCTCTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggggggaagca                                                                                                   >  2:1358711/1‑136 (MQ=25)
                                                      gcgcCTCCCGTGTTTCATTCCCCCCAACAC‑CCGCCCTTTTACCCACCGCCATGTGAAAAAAG‑GCCTC‑CGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCttt                                                                                           <  1:36197/151‑1 (MQ=255)
                                                      gCGGCTCCCGTGTTTTTTTCCCCCAAACAC‑CCGCGCTTTTACCCAACGCCATGTAAAAAAAT‑CCCTT‑CGGGGGGGCTGTCTCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCt                                                                                             <  2:2131483/149‑1 (MQ=255)
                                                        ggCTCCCTTGATTCATTCCCGCCAACAC‑CCGCGTTTTTACCCACCCCATTGTGAAAAAAT‑GCCTC‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCt                                                                                             <  1:420232/147‑1 (MQ=255)
                                                        ggCTCCCGTATTTCTTTCCCGCCAACAC‑CCGCGCTTTTATCCACCCCCATGTGAAAAAAT‑GCCTC‑CGGGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTa                                                                                          <  1:640800/150‑1 (MQ=255)
                                                         gCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggtgccggggaata                                                                                                 >  2:1180879/1‑131 (MQ=255)
                                                          atACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGTAcgagggaaagatct                                                                                            >  2:1299634/2‑137 (MQ=25)
                                                             cggTGATTCATCCCCCCCAACAC‑CCGCGAATTAATCCACCGCATTTAGAAAAATT‑GCCTC‑CGCGGCGGCTGTCCCCCAGCTTTTTTCCCGCGGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGccc                                                                                     <  1:1910652/148‑1 (MQ=255)
                                                             ccGTGTTTCATTCCCGCCAACAA‑CCGCCATTTTACCCAACCCCTTGTAAAAATTT‑GCCTC‑CGGGGCGCCTGTTTCCCAGCT‑TTTTCCCCCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGccc                                                                                     <  1:2177338/149‑1 (MQ=255)
                                                             ccGTGTTTC‑TTTCCCCCAAAAACCCGCCCTTTTACCCACCGCCATGTCAAAAAAT‑CCCTT‑CGCGGCGGTTTTCTCCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAg                                                                                        <  1:972508/147‑1 (MQ=255)
                                                               gtttttCTTTCCCCCCACCAC‑CCGCGCTTTTATCCACCCCCTTGTGAAAAAAG‑GCCTT‑GGGGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGt                                                                                   <  1:482270/146‑1 (MQ=255)
                                                                tgttcctttccccccAACAC‑CCGCCCTTTCCCCCACCGCCTTGTGAAAAAGC‑GCCTG‑CGCGGCGGTTTTCGCCCAGTTTTTTTCCCGCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGt                                                                                   <  2:1097151/136‑1 (MQ=255)
                                                                 gATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAG‑AAAAAACTTTaatacc                                                                                     >  1:146954/1‑142 (MQ=25)
                                                                       tttCCCCCAACAC‑CCGCCCTTTTCACCACCCCCATGTGAAAAAAT‑CCCTT‑CGGGGCGGTTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgttg                                                                            <  2:560526/149‑1 (MQ=255)
                                                                       ttccccccAACAC‑CCCCCCTTTTACCCACCGCCATGCGAAAAAAT‑CCCTT‑GGCGGCGGCTGTCGGCCAGTTATTTTGCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGt                                                                                   <  2:1653482/136‑1 (MQ=255)
                                                                        ttCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAAAAACTTTTAACCCCCCCTGTtttgt                                                                         >  1:267943/1‑151 (MQ=25)
                                                                          cccccaccaC‑CCGCGCATTTATCCACCGCCTTGTGAAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTtgtgt                                                                         <  1:464253/146‑1 (MQ=255)
                                                                               aacac‑cCGCCCTTTTACCCACCGCCTTGTGAAAAAAG‑CCCTT‑CGCGGCGGTTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                     <  1:1032545/148‑1 (MQ=255)
                                                                                  aa‑cCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAAAAACTTTAAACCCGTCATGTTTTGTGGCaa                                                                    >  1:382756/1‑146 (MQ=255)
                                                                                      cgcgCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAGGG‑AAAAGACTTTAAACCC‑TCATGGTGTGTggggttag                                                                 >  2:333337/1‑140 (MQ=25)
                                                                                       gccaaTTTACCCAACCCCTTGT‑CAAAAATCCCCTT‑CGGGGCGGTTTTCCGCCGTCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                            <  1:1508304/146‑1 (MQ=255)
                                                                                            tttATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAAAG‑AAAAAAATTTTAAACCCCCCTGTTTTGTGGGGAACGCTGCTgg                                                          >  1:787669/1‑147 (MQ=21)
                                                                                            tttATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCGAGAG‑AAAAAAATTTTAAACCCCCATGTTTTGTGGGGAAAGCTGCTgg                                                          >  1:465384/1‑147 (MQ=25)
                                                                                            tttATCCAACCCCATGTAAAAAAAT‑CCCTT‑CGCGGCGGCTGTCTCCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                          <  1:342028/147‑1 (MQ=255)
                                                                                            tttACCCACCGCCATGTCAAAAAAT‑CCCTC‑CGGGGCGGTTGTCTGCCAGTTATTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                           <  1:904170/146‑1 (MQ=255)
                                                                                            tttACCCACCCCCATGTGAAAAAAT‑GCTTC‑CGCGGCGGCTTTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                        <  1:1969844/149‑1 (MQ=255)
                                                                                            tttACCCAACGCCATGTGAAAAATT‑CCCTC‑GGGGGCGGTTGTCTCCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                            <  1:962394/145‑1 (MQ=255)
                                                                                                 ccAACCCCATGTGCAAAAAT‑CCCTC‑CGCGGCGGTTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                 <  1:1429417/151‑1 (MQ=255)
                                                                                                  caccccaagggcaaaaaaT‑CCCTC‑CGGGGGGGCTTTCTCCCAGGTTTTTTCCCCCACCAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                     <  2:721553/133‑1 (MQ=255)
                                                                                                  cAACCCCATGTGAAAAAAG‑TCCTT‑CCGGGCGTGTTGCCCCCATGTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                     <  2:2177442/146‑1 (MQ=255)
                                                                                                    aCGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAG‑AAAAAAATTTAAACCCCGCCTGTTTTGTGGGGAACGCTGCTGGGAaaaaaaggg                                               >  1:887493/1‑145 (MQ=25)
                                                                                                     cccctggtaaaaaaaT‑GCCTC‑CGGGGGGGTTGCCTGCAAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                   <  2:468096/135‑1 (MQ=255)
                                                                                                     cGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAAAG‑AAAAAAAATTATAACCCCCCTCTTGTGTGGGGAACGCTGCTGGTGACaaaatgg                                               >  1:1242600/1‑144 (MQ=21)
                                                                                                       ccttgtaaaaaatT‑GCCTG‑CGGGCGGGCTGCCGGCAAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                <  2:626541/138‑1 (MQ=255)
                                                                                                        cATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAG‑AAAAAAATTTTATCCCCTCCTGTTGTGTGGGGAACGCTGCTGGTGACAAtgtg                                                >  1:1613105/1‑145 (MQ=37)
                                                                                                             gaaaaaaa‑GCCTTCCGGGGCGGCTGTCTCCCGTTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCGTGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                      <  2:2076678/149‑1 (MQ=255)
                                                                                                                aaaaT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAACCTTTAAGCCC‑CCCTGTTGTTTTGGGATCCCTGCTGGTGAAAAATTGGGAAGGGGAggggg                                  >  1:810743/1‑147 (MQ=255)
                                                                                                                     tcccTT‑CGCGGGGGTTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                               <  1:1498079/148‑1 (MQ=255)
                                                                                                                       ccTC‑CGCGGCGCCTTTCGCCAAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                            <  2:434068/151‑1 (MQ=255)
                                                                                                                        cTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAGTTTTAACCCCCCCA‑GTGTTGTGGGGATAGCTGCTTGGGAAAAAGTGGGGAGGGAAGAgcgg                                 >  1:2049355/1‑145 (MQ=255)
                                                                                                                            cgcgGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAG‑AAAAAAATTTTAAACCCCCCTTTTTTGTGGGGATACCTGCTGGGTACAAAGTGGGAAggggaaggggggggatc                           >  1:307131/1‑134 (MQ=255)
                                                                                                                               ggcggcTGTCCCCCGGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                            <  2:1993011/144‑1 (MQ=255)
                                                                                                                                gcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACTGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                    <  2:89198/151‑1 (MQ=255)
                                                                                                                                     tttcgccgagtttttttcCCGCCAAAACAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                    <  2:1370889/133‑1 (MQ=255)
                                                                                                                                        ctcccAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGT‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                 <  2:506444/143‑1 (MQ=255)
                                                                                                                                           ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc             <  1:1550277/147‑1 (MQ=255)
                                                                                                                                           ccAGTTTTTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt          <  1:932212/150‑1 (MQ=255)
                                                                                                                                           ccAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                    <  1:1800706/140‑1 (MQ=255)
                                                                                                                                           ccAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGGG‑AAAAAAACTTAAACCCCCCCCGTTCTGTGGCGATCGCTCCGGGGGACCAAGTGGGAAGGGGAGGGCGGCGGAACAGCTccccaaacccgc           >  2:814227/1‑139 (MQ=255)
                                                                                                                                                  tttttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAG‑AAAAAACTTTTAACCCCCCCTTTTTTGTGGGGATCGCTGCTGGGAAAAAAGTGGAAAAGGAAGGGGGGGGGAACAGCCCCCCAACCGCGCCAATTTAg   >  1:135170/1‑150 (MQ=255)
                                                                                                                                                  ttttcCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc               <  2:1308814/138‑1 (MQ=255)
                                                                                                                                                  ttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa    <  1:395976/149‑1 (MQ=255)
                                                                                                                                                   tttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt  <  1:1516533/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357733‑358026

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