Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA53.6% ‑0.5 / 58.1 56L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (17/9);  new base A (0/30);  total (17/39)
Fisher's exact test for biased strand distribution p-value = 3.19e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-03

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357742‑358027
                                                                                                                                                           |                                                                                                                                            
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CCACCGCGCATTTATCCAACGCCA‑TGTGCCAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa                                                                                                                                            >  1:571782/1‑151 (MQ=255)
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa                                                                                                                                            >  2:1644016/1‑151 (MQ=255)
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa                                                                                                                                            >  1:439492/1‑151 (MQ=255)
     gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccg                                                                                                                                       >  2:517849/1‑151 (MQ=255)
     gCGTAGCGCACCACGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGGGATT‑CATTTCCGCCAA‑CAACCCCCCATGTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCGGTCTGGCTGCTATTGTT‑ACTCCGCGACAAAAccg                                                                                                                                       >  1:1690266/1‑151 (MQ=17)
        tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgccg                                                                                                                                    >  2:1279743/1‑151 (MQ=255)
                 aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCa                                                                                                                           >  2:1161824/1‑151 (MQ=255)
                  ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCag                                                                                                                          >  1:1341210/1‑151 (MQ=255)
                   gaatttCCCGCCGCACCCCTGGGGCGGCGGCTCCCGTGTTT‑ATTTTCCCCCAA‑CACCCGCGTTTTTATCCAACGCCT‑TGTGAAAAAA‑T‑GCCTCCGG‑GGCGGCTTTCTCCCAGCTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTcc                                                                                                                            <  2:1207122/144‑1 (MQ=255)
                   cattccccACCCGCAGACTTTTGTCAGCGGCTCCGGTATTT‑CTTTCCCCCCAA‑CACCCGCGCTTTTACCCACCCCAT‑TGTCAAAAAA‑T‑CCCTCCGC‑GGCGCCTTTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCagt                                                                                                                         <  1:1192672/145‑1 (MQ=255)
                      aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                      >  1:1649749/1‑151 (MQ=255)
                        cccccccGCGGCCCGGTGCCAGCGCCCCCCGTTTTT‑CTTTCCCCCCAA‑CACCCGCGTTTTTACCCACCCCCA‑TGTGAAAAAA‑T‑GCCTCCGG‑GGCGGCTTTCCCCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                    <  2:1457282/146‑1 (MQ=255)
                          cacccGCAGCCCTGTGCCAGGGGCTCCGGTTATT‑CTTTCCCCCAAA‑CACCCGCCCTTTTACCCACCCCCA‑TGGGAAAAAA‑T‑CCTTCCGC‑GGCGGCTTTCTGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                                 <  2:1475505/148‑1 (MQ=255)
                            gCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gg‑g                                                                                                              >  1:1013436/1‑151 (MQ=255)
                                  gACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCCACAAAACCGCCGTTCTCTCCAGTAGGGGC‑gg‑ggggaga                                                                                                        >  2:1085092/1‑146 (MQ=255)
                                   ccgtttgtCGGCGCCTCCCTTGTTT‑CTTTCCCCCCAA‑CAACCGCACTTTAACCCACCGCCA‑TGTAAAAAAA‑T‑GCTTTCGG‑GGGGCCTGTCTGCAATCTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑agag                                                                                                      <  2:94825/146‑1 (MQ=255)
                                   aCCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGGC‑agag                                                                                                      >  2:1382828/1‑151 (MQ=255)
                                       gtgtCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gg‑ggggcgggaaaa                                                                                                   >  1:238201/1‑141 (MQ=255)
                                       gtgtCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gg‑ggggcgaggaaa                                                                                                   >  1:898773/1‑141 (MQ=255)
                                        tgtCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gg‑ggggcgaagaaaa                                                                                                  >  2:1052313/1‑140 (MQ=255)
                                         tgtcggcgcctccctGATTCCATTCCCGCCAC‑CACCCGCGATTTTACCCACCCCCT‑TGTGAAAAAA‑G‑CCCTTCGC‑GGCGGTTGTCTCCCGGCTTTTTGC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAAt                                                                                                  <  2:916136/136‑1 (MQ=255)
                                          gcAGCGGCTCCGGTGTTT‑C‑TTTCCCCCAACCACCCGCGCTTTTACCCACCGCCA‑TGTAAAAAAT‑T‑GCCTCCGC‑GGGGGCGGTCGCCCATTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAg                                                                                                <  2:934236/149‑1 (MQ=255)
                                          ccggCGGCTACCGGTATT‑CTTTCCCCCCAA‑CACCCGCGCTTTTACCAACCGCAA‑GGTAAAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGc                                                                                               <  2:857561/148‑1 (MQ=255)
                                              cGGCTCCGGTAATT‑CTTTCCCCCAAA‑CACCCGCGCATTTACCCAACCCCA‑GGTGAAAAAA‑T‑GCCTCCGG‑GGGGGCGGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCtt                                                                                             <  2:189769/149‑1 (MQ=255)
                                                 ctcccGTGTTT‑CTTTCCCCCAAA‑CACCCGCGCTTTTACCCACCGCCA‑TGTGAAAAAA‑T‑GCCTCCGG‑GGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGc                                                                                        <  2:1003489/148‑1 (MQ=255)
                                                  tccggtgtTT‑CTTTCCCCCCAA‑C‑ACCCCCCATTTTACCCACCGCA‑AGTGGAAAAA‑TGGCCTCGGG‑GGGGGCTGTCGCCCAGTTTTTTTC‑CCCCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGcc                                                                                       <  2:595868/147‑1 (MQ=255)
                                                  caCCGTGATT‑CATTCCCCCCAA‑CACCCGC‑CTTTTTCCCCACCCCAGTGTGAAAAAA‑T‑CCCTTCGC‑GGCGGTTTTCTCCCGGTTATTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGcc                                                                                       <  1:489991/150‑1 (MQ=255)
                                                    ccGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGG‑GGCA‑AGGGAAAAAATTTTAAAcc                                                                                       >  2:17377/1‑150 (MQ=25)
                                                    ccGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGGG‑AGAGgaaaaaagtttaaacc                                                                                      >  2:879901/1‑137 (MQ=25)
                                                        taTT‑CTTTCCCGCCAA‑CAACCGCCCTTTTACCCACCGCCA‑TTTGCAAAAA‑T‑CCCTCCGC‑GGCGCTTGTCGGCCGGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCAt                                                                                 <  1:241940/150‑1 (MQ=255)
                                                          tt‑CTTTCCCCCCAA‑CACCCGCCCTTTTACCCACCGCAA‑TGTGAAAAAA‑G‑CCCTTCGG‑GGGGGTTGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCAtgt                                                                               <  2:1482349/151‑1 (MQ=255)
                                                             ctttccccccAA‑CACCCGCCATTTTACCCAACGCAA‑TGTGAAAAAA‑T‑CCCTCCGC‑GGCGGTTGTCTGCCAGTTATTTTC‑CCGCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                             <  2:1487730/143‑1 (MQ=255)
                                                             cATTCCCCCCAA‑CACCCGCGCATTTACCCAACGCCA‑TGTGAAAAAA‑T‑GCCTCCGC‑GGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                             <  2:1197807/151‑1 (MQ=255)
                                                             cATTCCCCCCAA‑CACCCCCCCTTTTACCCAACCCCA‑TGTGAAAAAT‑T‑GCCTTCGG‑GGGGGCTGTCGGCCAGTTATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                             <  2:1425841/151‑1 (MQ=255)
                                                                  cccccac‑caCCCGCGTATTAACCCAACGCCA‑‑GGGCAAAAATT‑CCCTCCGC‑GGGGGTTGTCGGCCAGTTTTTTCC‑CCGCAGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                        <  1:24692/148‑1 (MQ=255)
                                                                  ccccaaa‑caccCGCGCTTTTCCCCACCCCCG‑TGTACAAAAA‑T‑CCCTCCGC‑GGGGGCTGTCGGCCAGTTATTTTC‑CCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                        <  2:1401319/146‑1 (MQ=255)
                                                                            acCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGCGGCAGAGAGAAA‑AGTCTCAAGCCCCTCAGTTTTTGTGGCacccccg                                                                >  1:71896/1‑145 (MQ=21)
                                                                                  cATTTTTCCACCCCCT‑TGTCCAAAAT‑T‑CCCTCCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                         <  2:1217167/151‑1 (MQ=255)
                                                                                  cATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGG‑GGCA‑AAGGAAAAAACTTTAACCCCCCCCCGTTGT‑TGGGAAGCgctgcgcggg                                                         >  1:1380470/1‑146 (MQ=25)
                                                                                    tttACCCACCGCAA‑TGTGAAAAAA‑T‑GCTTCCGG‑GGGGCCTGTCGGCCGTCTATTTTC‑CCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                       <  1:505403/151‑1 (MQ=255)
                                                                                    tttACCCAACCCCA‑TGTAAAAAAA‑T‑CCCTTCGG‑GGGGGCTGTC‑GCCAGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                       <  2:1205507/151‑1 (MQ=255)
                                                                                    ttaacccACCGCCG‑GGTGAAAAAA‑T‑CCTTTCGC‑GGGGGCTGCCCCACAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                       <  2:639640/146‑1 (MQ=255)
                                                                                        ctcaaccgctgggtgAAAAAA‑T‑GCCTGCGG‑GGCGGCTTTCGGCCAGTTATTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                    <  1:626670/139‑1 (MQ=255)
                                                                                        cccACCGCCA‑TGTGAAAAAA‑T‑CCTTTCGG‑GGCGGCTGTCTCCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                   <  2:1063106/150‑1 (MQ=255)
                                                                                            aCGCCA‑GGGGAAAAAA‑G‑CCCTTCGG‑GGGGGCTTTCGCCAAGTTTTTTTC‑CCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                               <  2:1299161/151‑1 (MQ=255)
                                                                                             cccccttgtaaaaaaa‑T‑CCCTTCGG‑GGGGGCTGTCTGCCATTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                               <  2:1654513/140‑1 (MQ=255)
                                                                                             cGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCC‑AGAGAAAAAATTTTTAACCCCCCATGTTTTTTGGGCAACGCTGCTGGtgaaaaagaggga                                              >  2:887553/1‑139 (MQ=25)
                                                                                             acccccatttaaaaaa‑a‑TCCTTCGGGGGCGGTTGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                <  2:1564697/140‑1 (MQ=255)
                                                                                                   tgtgAAAAAA‑T‑CCCTCCGG‑GGGGGCTGCCGCCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                         <  2:365101/151‑1 (MQ=255)
                                                                                                   tgaaaaaaaa‑T‑CCCTCCGC‑GGCGCCTTTCTGCCAGTTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                         <  1:1161254/146‑1 (MQ=255)
                                                                                                                              ttGTCTGCC‑GCTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:684132/150‑1 (MQ=255)
                                                                                                                              cttccTGCCAGTTTTTTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:637447/147‑1 (MQ=255)
                                                                                                                              cggccggCCAGTTTTTTTC‑CCGCCCAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:517635/145‑1 (MQ=255)
                                                                                                                              cTGTCTGCCGGTTTTTTTC‑CCGCCACAAAAAGACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:1555412/151‑1 (MQ=255)
                                                                                                                              cTGTCGGCCAGCTTTTTTC‑CCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  2:1352039/151‑1 (MQ=255)
                                                                                                                              cTGTCAGCCAGTTGTTTCC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  1:246352/151‑1 (MQ=255)
                                                                                                                                           tttttc‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg    <  2:393153/150‑1 (MQ=255)
                                                                                                                                             tttc‑CCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt  <  2:439645/151‑1 (MQ=255)
                                                                                                                                              ttccccgccgAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt  <  2:1389164/147‑1 (MQ=255)
                                                                                                                                                           |                                                                                                                                            
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATT‑CATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAA‑T‑GCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GG‑GGCA‑AGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357742‑358027

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